From cassar at PILOTBOOK.EU  Fri Jan  2 13:05:24 2009
From: cassar at PILOTBOOK.EU (Pierre Cassar)
Date: Fri, 2 Jan 2009 19:05:24 +0100
Subject: g2ctl fatal error
Message-ID: <FRI.2.JAN.2009.190524.0100.>

Hi All and Happy New Year.

I was trying to use gens grib2 data available at:
ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/gens/prod

geavg.t12z.pgrb2af06
geavg.t12z.pgrb2af12
geavg.t12z.pgrb2af18
...

Time ago those file were in grib1 format and I had no problem to deal with them using Grads 1.8:

e.g.
perl grib2ctl.pl geavg.t12z.pgrbaf06 geavg.t12z.pgrbaf06.idx >geavg.t12z.pgrbaf06.ctl
gribmap -i geavg.t12z.pgrbaf06.ctl
etc...

Now if I issue (Windows XP, Vista, Freebsd 6.3, Freebsd 7.0):
perl g2ctl.pl geavg.t12z.pgrbaf06 geavg.t12z.pgrbaf06.idx >geavg.t12z.pgrbaf06.ctl

I only obtain:

*** FATAL ERROR: unknown option -ctl_ens ***

Big problem:
  either geavg.t12z.pgrbaf06 is missing or not a grib file
  or wgrib2 is not on your path or wgrib2 is too old
  or can not write to C:\Users\Cassar\AppData\Local\Temp\g7224.tmp (full disk or permissions).

(On Freebsd only changes the location of tmp dir)

I have these versions:

g2ctl.pl 0.0.4k
wgrib2 v0.1.7.2.oga.1 (on windows)
wgrib2 v0.1.7.1.oga.3 (on freebsd)

and wgrib2 is in my path.

Any ideas, suggestions?

Thank you.

Pierre Cassar.


From matt.alonso at GMAIL.COM  Fri Jan  2 13:35:24 2009
From: matt.alonso at GMAIL.COM (Matt Alonso)
Date: Fri, 2 Jan 2009 13:35:24 -0500
Subject: projection information
Message-ID: <FRI.2.JAN.2009.133524.0500.>

Hi All,

Does anyone have - or know where I could find - more detailed information
about the projections GrADS uses?  I was attempting to use another GIS
program (ArcGIS) to come up w/some graphics to use as basemaps in GrADS but
I have not been able to find a projection in ArcGIS that seems to line up
with GrADS' nps or lambert projections.  Has anyone been able to do
something similar to this?  Or does anyone have any more detailed
information that might point me in the right direction?

Thanks in advance.

Cheers,
Matt
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090102/d8a7427c/attachment.html 

From arlindo.dasilva at GMAIL.COM  Mon Jan  5 00:30:28 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 5 Jan 2009 00:30:28 -0500
Subject: New Release GrADS 2.0.a4
In-Reply-To: <495BB2F5.6020404@scotiaweather.com>
Message-ID: <MON.5.JAN.2009.003028.0500.>

On Wed, Dec 31, 2008 at 12:59 PM, mmacleod <mmacleod at scotiaweather.com>wrote:

> Hi Jennifer.
>
> Happy New Year!
>
> I have checked the COLA web site and see that there is a Ubuntu compile
> for  version 2.0.a3 but not for 2.0.a4.  We use the Ubuntu O/S on our
> system.  Do you have a time-line for a compile of 2.0.a4 under Ubuntu?
>
> I also use a Windows platform for some testing of model runs.  Will
> there be a Windows XP  compile of a4 in the near future?
>

I have uploaded a Win32 build  to the OpenGrADS download area at
Sourceforge:

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=290254

Notice that this build has been labeled as v2.0.a4.oga.1 (meaning v2.0.a4 +
OpenGrADS patch 1.)  In particular, the sources have been slightly patched
for Windows; consult the ChangeLog for additional details. For any practical
purposes it should have the same functionality as the official COLA v2.0.a4
release.

The OpenGrADS distribution of the Windows version of GrADS comes in two
flavors. Most users will be interested in the Win32 Superpack, which bundles
the GrADS binaries, and all else that is needed to run GrADS on Windows. (A
relocatable zip version is also available for those users who do not have
the necessary privileges to install software.) For those experienced users
working on a Cygwin environment, a lighter weight version is also available,
although it requires the same level of setup as on a Linux platform.  For
 additional information, consult the wiki:


http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows

Enjoy,

       Arlindo


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/696e2822/attachment.html 

From arlindo.dasilva at GMAIL.COM  Mon Jan  5 00:54:34 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 5 Jan 2009 00:54:34 -0500
Subject: New Release GrADS 2.0.a4
In-Reply-To: <495BB2F5.6020404@scotiaweather.com>
Message-ID: <MON.5.JAN.2009.005434.0500.>

On Wed, Dec 31, 2008 at 12:59 PM, mmacleod <mmacleod at scotiaweather.com>wrote:

> Hi Jennifer.
>
> Happy New Year!
>
> I have checked the COLA web site and see that there is a Ubuntu compile
> for  version 2.0.a3 but not for 2.0.a4.  We use the Ubuntu O/S on our
> system.  Do you have a time-line for a compile of 2.0.a4 under Ubuntu?
>

I am not sure why you need a Ubuntu specific build but the x86_64 build that
I posted on Sourgeforge

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032

happened to be made on Ubuntu 8.04. However, the i686 build which was made
on a RedHat machine works just fine on Ubuntu (if running i686 on 64-bit
Ubuntu be sure to install the ia32 shared libraries --- these are not
installed by default on Ubuntu.)

   Arlindo


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/bd8fa19d/attachment.html 

From arlindo.dasilva at GMAIL.COM  Mon Jan  5 01:08:31 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 5 Jan 2009 01:08:31 -0500
Subject: New Release GrADS 2.0.a4
In-Reply-To: <CA13161D-2E09-45BD-A5A1-F5DD4EAF8EA0@cola.iges.org>
Message-ID: <MON.5.JAN.2009.010831.0500.>

All,
   I have posted builds for FreeBSD and universal binaries for Mac OS X on
the OpenGrADS download area at Sourceforge:

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032

  I also included my Linux x86_64/i686 builds along with patched sources
(which build out of the box with the OpenGrADS supplibs). The FreeBSD build
was made of FreeBSD v6.3; they should work fine on FreeBSD 7.0, let me know
otherwise.

Notice that this build has been labeled as v2.0.a4.oga.1 (meaning v2.0.a4 +
OpenGrADS patch 1.)  In particular, the sources have been slightly patched
for building with the OpenGrADS supplibs; consult the ChangeLog for
additional details. For any practical purposes it should have the same
functionality as the official COLA v2.0.a4 release.

Question: Does anybody care for IRIX64 and IA64 (Altix) builds?

As I mentioned on a previous posting, the Windows build can be found here:

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=290254

  Enjoy!

     Arlindo


On Wed, Dec 31, 2008 at 12:03 PM, Jennifer Adams <jma at cola.iges.org> wrote:

> Dear All,
> I have posted a new version of GrADS: 2.0.a4. This is primarily a bug-fix
> release. The GrADS downloads web page (*
> http://iges.org/grads/downloads.html*<http://iges.org/grads/downloads.html>)
> has been updated with links to the source code and pre-compiled binaries for
> a few OS flavors. Here are highlights of what is new in 2.0.a4:
>
>     Bugs:
>     * Fixed function stnave(): error occurred when setting undef mask for
> result
>       if some reports were missing. Led to incorrect averages!!
>     * Fixed function asum(): error occurred when setting undef mask for
> final result.
>       Led to answer being always undefined!!
>     * Fixed seg fault for non-templated netcdf files when X dim wraps the
> globe
>       and other dims are fixed.
>     * No precision errors for some hourly GDS data sets.
>
>     Misc:
>     * 'set frame circle' now works for orthograhic projections
>     * Support for grib2 grid 204 (curvilinear orthogonal)
>     * The color of a constant grid is whatever 'set ccolor' is, instead of
> being fixed red
>     * The 'q ctlinfo' and 'q file' output is enhanced and more accurate.
>       (This one helps GDS 2.0 know when a data set is 4D or 5D.)
>     * Create postscript output from inside GrADS using 'print' command
> without 'enable print'
>
>
>
> --
> Jennifer M. Adams
> IGES/COLA
> 4041 Powder Mill Road, Suite 302
> Calverton, MD 20705
> jma at cola.iges.org
>
>
>
>


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/d1c4f1b0/attachment.html 

From meteo at SPORT123.IT  Mon Jan  5 04:43:23 2009
From: meteo at SPORT123.IT (Valerio Guerra)
Date: Mon, 5 Jan 2009 10:43:23 +0100
Subject: New Release GrADS 2.0.a4
Message-ID: <MON.5.JAN.2009.104323.0100.>

Arlindo da Silva wrote:

> The FreeBSD build was made of FreeBSD v6.3; they should work fine on FreeBSD 7.0, let me know otherwise.

Thank you Arlindo for the new Freebsd build.

I've tested the batch version of gradsdap on Freebsd 6.1 and it's ok.
I also tried a test on Freebsd 7.0 and it was unsuccessful:

$ gradsdap -bl
/libexec/ld-elf.so.1: /usr/X11R6/lib/libXext.so.6: unsupported file layout

It seems to me strange: I do not use X ???

Any idea?

Thank you again and Happy New Year.

Valerio Guerra


From cassar at PILOTBOOK.EU  Mon Jan  5 05:47:53 2009
From: cassar at PILOTBOOK.EU (Pierre Cassar)
Date: Mon, 5 Jan 2009 11:47:53 +0100
Subject: g2ctl fatal error
Message-ID: <MON.5.JAN.2009.114753.0100.>

Wesley Ebisuzaki wrote:

> Hi Pierre,
>
> You need a new version of wgrib2. The -ctl_ens option was added with v0.1.7.7g 10/2008.
>
> Happy new year,
> Wesley

Thank you Wesley for your note.

Unfortunately I tested wgrib2.v1.7.7d from
ftp://ftp.cpc.ncep.noaa.gov/wd51we/wgrib2 on Windows Vista but it
always crashes (UAC on/off, Administrator on/off, etc).
The latest windows build running on Vista is the 0.1.7.2.oga.1.
Where can I find wgrib2 version 1.7.7d (or later) running on Windows
Vista?
Moreover I was not able to find a Freebsd 6.3 (and possibly 7.0) build
for version 1.7.7d or later. My latest working wgrib2 build for
Freebsd is the 0.1.7.1.oga.3.

Thank you.

Pierre Cassar.


From meteo at SPORT123.IT  Mon Jan  5 06:05:19 2009
From: meteo at SPORT123.IT (Valerio Guerra)
Date: Mon, 5 Jan 2009 12:05:19 +0100
Subject: wgrib2 on Freebsd
Message-ID: <MON.5.JAN.2009.120519.0100.>

Pierre Cassar wrote:

> Moreover I was not able to find a Freebsd 6.3 (and possibly 7.0) build for version 1.7.7d or later. My latest working wgrib2 build for Freebsd is the 0.1.7.1.oga.3

I'm also interested in an updated wgrib2 build for FreeBSD 6.x and in the near future for FreeBSD 7.0. I tried to compile it by myself but the effort was too big for my sysadmin skills :-(

Best regards.

Valerio Guerra


From arlindo.dasilva at GMAIL.COM  Mon Jan  5 08:14:27 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 5 Jan 2009 08:14:27 -0500
Subject: New Release GrADS 2.0.a4
In-Reply-To: <PC187020090105104323074682832701@Poseidon>
Message-ID: <MON.5.JAN.2009.081427.0500.>

On Mon, Jan 5, 2009 at 4:43 AM, Valerio Guerra <meteo at sport123.it> wrote:

> Arlindo da Silva wrote:
>
> > The FreeBSD build was made of FreeBSD v6.3; they should work fine on
> FreeBSD 7.0, let me know otherwise.
>
> Thank you Arlindo for the new Freebsd build.
>
> I've tested the batch version of gradsdap on Freebsd 6.1 and it's ok.
> I also tried a test on Freebsd 7.0 and it was unsuccessful:
>
> $ gradsdap -bl
> /libexec/ld-elf.so.1: /usr/X11R6/lib/libXext.so.6: unsupported file layout
>
> It seems to me strange: I do not use X ???
>

It is strange. Even if you do not use X you need the X libraries to start
GrADS.


>
> Any idea?
>

I just uploaded a FreeBSD 7.0 build; see if that works.



>
> Thank you again and Happy New Year.
>

Happy New Year!

   Arlindo


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/fcfadfd9/attachment.html 

From meteo at SPORT123.IT  Mon Jan  5 10:24:00 2009
From: meteo at SPORT123.IT (Valerio Guerra)
Date: Mon, 5 Jan 2009 16:24:00 +0100
Subject: New Release GrADS 2.0.a4
Message-ID: <MON.5.JAN.2009.162400.0100.>

Arlindo da Silva wrote:

> I just uploaded a FreeBSD 7.0 build; see if that works.

thank you Arlindo but unfortunately I get the same error of your Freebsd 6.3 build:

$ gradsdap -bl
/libexec/ld-elf.so.1: /usr/X11R6/lib/libXext.so.6: unsupported file layout

It might be a problem tied to the 64 bit architecture of the Freebsd 7.0 test machine?
Any ideas, suggestion?

Currently I use Freebsd 6.1 and Grads 2.0.a4 with no problems but my sysadmin told me that in the near future we have to upgrade to Freebsd 7 (64 bit) so I'm quite worried...

Thank you very much.

Valerio Guerra


From brian.walsh at AMEC.COM  Mon Jan  5 12:58:55 2009
From: brian.walsh at AMEC.COM (Walsh, Brian J)
Date: Mon, 5 Jan 2009 12:58:55 -0500
Subject: Filled Contours
Message-ID: <MON.5.JAN.2009.125855.0500.>


Hello,
 
I am trying to do a plot of 3 variables; one is shaded (dew point
depression), the second is dashed contours (temperature) and the third
is area of freezing rain (solid contour).  Is there any function in
grads where I can fill the solid contour with something such as dots,
stripes, lines, etc. or is this not possible since I already have one
variable being shaded?
 
Thanks,
 
Brian Walsh


The information contained in this e-mail is intended only for the individual or entity to whom it is addressed. 
Its contents (including any attachments) may contain confidential and/or privileged information. 
If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents. 
If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/920c29c8/attachment.html 

From brian.walsh at AMEC.COM  Mon Jan  5 13:06:26 2009
From: brian.walsh at AMEC.COM (Walsh, Brian J)
Date: Mon, 5 Jan 2009 13:06:26 -0500
Subject: Multiple Parameters on one Plot
Message-ID: <MON.5.JAN.2009.130626.0500.>


I am trying to do a plot of 3 variables; one is shaded (dew point
depression), the second is dashed contours (temperature) and the third
is area of freezing rain (solid contour).  Is there any function in
grads where I can fill the solid contour with something such as dots,
stripes, lines, etc. or is this not possible since I already have one
variable being shaded?
 
Thanks,
 


The information contained in this e-mail is intended only for the individual or entity to whom it is addressed. 
Its contents (including any attachments) may contain confidential and/or privileged information. 
If you are not an intended recipient you must not use, disclose, disseminate, copy or print its contents. 
If you receive this e-mail in error, please notify the sender by reply e-mail and delete and destroy the message.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/8afb0f6a/attachment.html 

From matias at TROPOSFERA.ATMOSFCU.UNAM.MX  Mon Jan  5 07:24:31 2009
From: matias at TROPOSFERA.ATMOSFCU.UNAM.MX (Matias Mendez)
Date: Mon, 5 Jan 2009 07:24:31 -0500
Subject: Filled Contours
In-Reply-To: <D688BDE22B02984496332260B9E36A8F0A7BE22B@ATT1-EX1.am.int.amec.com>
Message-ID: <MON.5.JAN.2009.072431.0500.>

Hello Brian,

Check the following link, where you can get some useful GrADS scripts:

http://wind.geophys.tohoku.ac.jp/~kodamail/wiki/wiki.cgi/GrADS/script?lang=en

You must use hatch.gs to draw hatch, you can specify value range to be
drawn, angle, density and hatch interval.

Good luck,

Matias



El lun, 05-01-2009 a las 12:58 -0500, Walsh, Brian J escribi?:
> Hello,
>
> I am trying to do a plot of 3 variables; one is shaded (dew point
> depression), the second is dashed contours (temperature) and the third
> is area of freezing rain (solid contour).  Is there any function in
> grads where I can fill the solid contour with something such as dots,
> stripes, lines, etc. or is this not possible since I already have one
> variable being shaded?
>
> Thanks,
>
> Brian Walsh
> The information contained in this e-mail is intended only for the
> individual or entity to whom it is addressed.
> Its contents (including any attachments) may contain confidential
> and/or privileged information.
> If you are not an intended recipient you must not use, disclose,
> disseminate, copy or print its contents.
> If you receive this e-mail in error, please notify the sender by reply
> e-mail and delete and destroy the message.
>


From msponsler at COMCAST.NET  Mon Jan  5 14:12:50 2009
From: msponsler at COMCAST.NET (Mark Sponsler)
Date: Mon, 5 Jan 2009 19:12:50 +0000
Subject: New Release GrADS 2.0.a4
Message-ID: <MON.5.JAN.2009.191250.0000.>

Many thanks Arlindo - Just our of curiosity - does this version include the extensions that are in 1.9?
--
Thanks,
Mark

-------------- Original message --------------
From: Arlindo da Silva <arlindo.dasilva at GMAIL.COM>

On Wed, Dec 31, 2008 at 12:59 PM, mmacleod <mmacleod at scotiaweather.com> wrote:

Hi Jennifer.

Happy New Year!

I have checked the COLA web site and see that there is a Ubuntu compile
for  version 2.0.a3 but not for 2.0.a4.  We use the Ubuntu O/S on our
system.  Do you have a time-line for a compile of 2.0.a4 under Ubuntu?

I also use a Windows platform for some testing of model runs.  Will
there be a Windows XP  compile of a4 in the near future?



I have uploaded a Win32 build  to the OpenGrADS download area at Sourceforge:


https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=290254



Notice that this build has been labeled as v2.0.a4.oga.1 (meaning v2.0.a4 + OpenGrADS patch 1.)  In particular, the sources have been slightly patched for Windows; consult the ChangeLog for additional details. For any practical purposes it should have the same functionality as the official COLA v2.0.a4 release.


The OpenGrADS distribution of the Windows version of GrADS comes in two flavors. Most users will be interested in the Win32 Superpack, which bundles the GrADS binaries, and all else that is needed to run GrADS on Windows. (A relocatable zip version is also available for those users who do not have the necessary privileges to install software.) For those experienced users working on a Cygwin environment, a lighter weight version is also available, although it requires the same level of setup as on a Linux platform.  For  additional information, consult the wiki:


   http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows



Enjoy,


       Arlindo



--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/7c9316da/attachment.html 

From hrc at EAS.GATECH.EDU  Mon Jan  5 15:06:57 2009
From: hrc at EAS.GATECH.EDU (Hai-Ru Chang)
Date: Mon, 5 Jan 2009 15:06:57 -0500
Subject: netcdf in grads-2.0a4
Message-ID: <MON.5.JAN.2009.150657.0500.>

Hello,
            I tried to compile grads-2.0a4 on my opteron server
running Suse 10.3.  But the configure refused to include netcdf
support with following messages:

Checking in system locations for libraries to support netcdf ...
checking for udunits.h... (cached) yes
checking for utInit in -ludunits... (cached) yes
checking for nc_open in -lnetcdf... no
checking for netcdf.h with compiler... yes
checking for netcdf.h with preprocessor... yes
checking for netcdf 3 interface... yes
- netcdf disabled


As you can see the program can find netcdf lib.  But it can't find
'nc_open'.   I use gcc to compile grads but my netcdf was compile
with icc.  Will this cause problem?  How can I fix it?  I need netcdf
compiled by icc for other applications.





Dr. Hai-Ru Chang
School of Earth and Atmospheric Sciences
Georgia Institute of Technology
Atlanta, GA 30332-0340


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/df3ef2d6/attachment.html 

From arlindo.dasilva at GMAIL.COM  Mon Jan  5 15:46:51 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 5 Jan 2009 15:46:51 -0500
Subject: netcdf in grads-2.0a4
In-Reply-To: <4734CFFD-43B1-4C75-A799-5B7992FBB344@eas.gatech.edu>
Message-ID: <MON.5.JAN.2009.154651.0500.>

On Mon, Jan 5, 2009 at 3:06 PM, Hai-Ru Chang <hrc at eas.gatech.edu> wrote:

> Hello,           I tried to compile grads-2.0a4 on my opteron server
> running Suse 10.3.  But the configure refused to include netcdf support with
> following messages:
>
> Checking in system locations for libraries to support netcdf ...
> checking for udunits.h... (cached) yes
> checking for utInit in -ludunits... (cached) yes
> checking for nc_open in -lnetcdf... no
> checking for netcdf.h with compiler... yes
> checking for netcdf.h with preprocessor... yes
> checking for netcdf 3 interface... yes
> - netcdf disabled
>
>
> As you can see the program can find netcdf lib.  But it can't find
> 'nc_open'.
>

Most likely it is not finding a required icc library....


> I use gcc to compile grads but my netcdf was compile with icc.  Will this
> cause problem?  How can I fix it?  I need netcdf compiled by icc for other
> applications.
>

Try compiling grads with icc/icpc

export CC=icc
export CXX=icpc
./configure

  Arlindo


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090105/52996754/attachment.html 

From ela at COLA.IGES.ORG  Mon Jan  5 17:22:39 2009
From: ela at COLA.IGES.ORG (Eric Altshuler)
Date: Mon, 5 Jan 2009 17:22:39 -0500
Subject: netcdf in grads-2.0a4
In-Reply-To: <1153612368.15941231194144587.JavaMail.root@mail.iges.org>
Message-ID: <MON.5.JAN.2009.172239.0500.>

Dear Dr. Chang,

You can build multiple netcdf installations using different compilers, just put them in different directories (--prefix argument to configure script for netcdf). I have one netcdf build using ifort/icc and another using gfortran/gcc. Then, when building other codes that use netcdf, you can specify the path to the netcdf build you want to use (usually through a NETCDF environment variable or by explicitly specifying the lib and include paths).

Best regards,

Eric L. Altshuler
Assistant Research Scientist
Center for Ocean-Land-Atmosphere Studies
4041 Powder Mill Road, Suite 302
Calverton, MD 20705-3106
USA

E-mail: ela at cola.iges.org
Phone: (301) 902-1257
Fax: (301) 595-9793

From: "Hai-Ru Chang" <hrc at EAS.GATECH.EDU>
To: GRADSUSR at LIST.CINECA.IT
Sent: Monday, January 5, 2009 3:06:57 PM GMT -05:00 US/Canada Eastern
Subject: netcdf in grads-2.0a4

Hello,
            I tried to compile grads-2.0a4 on my opteron server
running Suse 10.3.  But the configure refused to include netcdf
support with following messages:

Checking in system locations for libraries to support netcdf ...
checking for udunits.h... (cached) yes
checking for utInit in -ludunits... (cached) yes
checking for nc_open in -lnetcdf... no
checking for netcdf.h with compiler... yes
checking for netcdf.h with preprocessor... yes
checking for netcdf 3 interface... yes
- netcdf disabled


As you can see the program can find netcdf lib.  But it can't find
'nc_open'.   I use gcc to compile grads but my netcdf was compile
with icc.  Will this cause problem?  How can I fix it?  I need netcdf
compiled by icc for other applications.





Dr. Hai-Ru Chang
School of Earth and Atmospheric Sciences
Georgia Institute of Technology
Atlanta, GA 30332-0340


From preuter at LABRI.FR  Tue Jan  6 06:57:53 2009
From: preuter at LABRI.FR (Patrick Reuter)
Date: Tue, 6 Jan 2009 12:57:53 +0100
Subject: Nearest defined neighbor function
Message-ID: <TUE.6.JAN.2009.125753.0100.>

Dear all,

I have a variable which is only defined over the sea, elsewhere there
are undefined values. Is there a function in Grads that can draw me a
map where at locations of undefined values, instead of leaving it
blank, draws me the value of the nearest defined value of the
8-neighborhood ? I tried smth9, but it does not exactly what I want.

Can anybody help?

Thanks in advance

  Patrick


From preuter at LABRI.FR  Tue Jan  6 08:39:09 2009
From: preuter at LABRI.FR (Patrick Reuter)
Date: Tue, 6 Jan 2009 14:39:09 +0100
Subject: 2nd try : Nearest defined neighbor function
Message-ID: <TUE.6.JAN.2009.143909.0100.>

Sorry for multiples, message was somehow rejected ..

Dear all,

I have a variable which is only defined over the sea, elsewhere there
are undefined values. Is there a function in Grads that can draw me a
map where at locations of undefined values, instead of leaving it
blank, draws me the value of the nearest defined value of the
8-neighborhood ? I know about the existence of 'smth9', but it does
not exactly what I want, since it smoothes out, and 'const', since you
have to specify a default value which I don't know since it should to
be determined from the nearest neighbor !

Can anybody help?

Thanks in advance

   Patrick


From romero619 at HOTMAIL.COM  Tue Jan  6 11:47:16 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Tue, 6 Jan 2009 08:47:16 -0800
Subject: 2nd try : Nearest defined neighbor function
In-Reply-To: <20090106143909.10134n0imhcbuebh@webmail.labri.fr>
Message-ID: <TUE.6.JAN.2009.084716.0800.>


check out the const function in the grads documentation.
it allows you to set a specific value/set of values to a constant value, i.e., in your case, set all undefined values to zero.
'd const(htsgwsfc,0,-u)'
 
'htsgwsfc' is ocean wave height.
this would plot the variable 'htsgwsfc', and set all areas where 'htsgwsfc' to zero, near the coast or overland.Note that this fades the value to zero near the coast, so its not extremely accurate, since the values near the coast might actually be higher than shown when using the const function. Example: an area near the coast where the closest-defined value is '3'; the "nearest neighbor" value should also be '3', but the const function fades to zero near the coast, so the value might be plotted as '1', '2', or near '0' instead
 
hope this helps...
> Date: Tue, 6 Jan 2009 14:39:09 +0100> From: preuter at LABRI.FR> Subject: 2nd try : Nearest defined neighbor function> To: GRADSUSR at LIST.CINECA.IT> > Sorry for multiples, message was somehow rejected ..> > Dear all,> > I have a variable which is only defined over the sea, elsewhere there> are undefined values. Is there a function in Grads that can draw me a> map where at locations of undefined values, instead of leaving it> blank, draws me the value of the nearest defined value of the> 8-neighborhood ? I know about the existence of 'smth9', but it does> not exactly what I want, since it smoothes out, and 'const', since you> have to specify a default value which I don't know since it should to> be determined from the nearest neighbor !> > Can anybody help?> > Thanks in advance> > Patrick
_________________________________________________________________
It?s the same Hotmail?. If by ?same? you mean up to 70% faster.
http://windowslive.com/online/hotmail?ocid=TXT_TAGLM_WL_hotmail_acq_broad1_122008
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090106/15430648/attachment.html 

From romero619 at HOTMAIL.COM  Tue Jan  6 11:52:59 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Tue, 6 Jan 2009 08:52:59 -0800
Subject: trouble plotting latest gshhs_land.shp ver.1.6 using shape
 extension
Message-ID: <TUE.6.JAN.2009.085259.0800.>


Hello,
I just downloaded the latest version of the gshhs_* shapefiles, version 1.6, since this latest version uses polygons (the previous versions used 'lines') and can be used to do filled, hi-res basemaps.
 
here's the url to the latest versions:
ftp://ftp.ihg.uni-duisburg.de/GIS/GISData/GSHHS/
 
I've tried plotting the gshhs_land.shp file in opengrads 1.9 on both windows&linux, and this shapefile is crashing grads in both cases.
the gshhs_lakes.shp file works fine, but the land file crashes.
can anyone else please test&confirm this?
 
Thanks,
P.Romero
_________________________________________________________________
Send e-mail anywhere. No map, no compass.
http://windowslive.com/oneline/hotmail?ocid=TXT_TAGLM_WL_hotmail_acq_anywhere_122008
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090106/21670716/attachment.html 

From smcmillan at PLANALYTICS.COM  Tue Jan  6 11:45:28 2009
From: smcmillan at PLANALYTICS.COM (Stephen R McMillan)
Date: Tue, 6 Jan 2009 10:45:28 -0600
Subject: 2nd try : Nearest defined neighbor function
In-Reply-To: <20090106143909.10134n0imhcbuebh@webmail.labri.fr>
Message-ID: <TUE.6.JAN.2009.104528.0600.>

Patrick,
What if an undefined point has more than one "nearest neighbor" and are of 
different values?  Do you want to assign an average of these, the highest, 
lowest, etc.?  I'm assuming each gridpoint has up to eight nearest 
neighbors in complex coastal areas.  For example, you could get the max 
nearest neighbor using something like this:

'define newvar= max(max(var,x=-152,x=-150),y=262,y=264)'

Then assign result to point (in this case, x=-151, y=263) using the 'set 
defval var...' command.

Similarly, you could use 'define var= min(min...' to get min value of 
nearest 8-neighbor or the "aave" command to get average of nearest 
neighbor.

Hope this helps....
Stephen McMillan




Patrick Reuter <preuter at LABRI.FR> 
Sent by: GRADSUSR at LIST.CINECA.IT
01/06/2009 07:39 AM
Please respond to
GRADSUSR at LIST.CINECA.IT


To
GRADSUSR at LIST.CINECA.IT
cc

Subject
2nd try : Nearest defined neighbor function





Sorry for multiples, message was somehow rejected ..

Dear all,

I have a variable which is only defined over the sea, elsewhere there
are undefined values. Is there a function in Grads that can draw me a
map where at locations of undefined values, instead of leaving it
blank, draws me the value of the nearest defined value of the
8-neighborhood ? I know about the existence of 'smth9', but it does
not exactly what I want, since it smoothes out, and 'const', since you
have to specify a default value which I don't know since it should to
be determined from the nearest neighbor !

Can anybody help?

Thanks in advance

Patrick




***************************************************
The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary.  If you are not the intended recipient, you may not review, copy or distribute this message.  If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090106/ecdde99f/attachment.html 

From romero619 at HOTMAIL.COM  Tue Jan  6 12:05:56 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Tue, 6 Jan 2009 09:05:56 -0800
Subject: 2nd try : Nearest defined neighbor function
In-Reply-To: <OF0619735A.5BE0BD5C-ON85257536.00558A3C-86257536.005C267F@planalytics.com>
Message-ID: <TUE.6.JAN.2009.090556.0800.>


it would be very cool if grads included "nearest neighbor" functions.
maybe a UDF? how hard/complicated would it be to setup the math for each case: max, min, average? would a cause plotting to be somewhat slower/processor intensive?



Date: Tue, 6 Jan 2009 10:45:28 -0600From: smcmillan at PLANALYTICS.COMSubject: Re: 2nd try : Nearest defined neighbor functionTo: GRADSUSR at LIST.CINECA.ITPatrick,What if an undefined point has more than one "nearest neighbor" and are of different values?  Do you want to assign an average of these, the highest, lowest, etc.?  I'm assuming each gridpoint has up to eight nearest neighbors in complex coastal areas.  For example, you could get the max nearest neighbor using something like this:'define newvar= max(max(var,x=-152,x=-150),y=262,y=264)'Then assign result to point (in this case, x=-151, y=263) using the 'set defval var...' command.Similarly, you could use 'define var= min(min...' to get min value of nearest 8-neighbor or the "aave" command to get average of nearest neighbor.Hope this helps....Stephen McMillan



Patrick Reuter <preuter at LABRI.FR>Sent by: GRADSUSR at LIST.CINECA.IT
01/06/2009 07:39 AM




Please respond toGRADSUSR at LIST.CINECA.IT





To
GRADSUSR at LIST.CINECA.IT


cc



Subject
2nd try : Nearest defined neighbor function




Sorry for multiples, message was somehow rejected ..Dear all,I have a variable which is only defined over the sea, elsewhere thereare undefined values. Is there a function in Grads that can draw me amap where at locations of undefined values, instead of leaving itblank, draws me the value of the nearest defined value of the8-neighborhood ? I know about the existence of 'smth9', but it doesnot exactly what I want, since it smoothes out, and 'const', since youhave to specify a default value which I don't know since it should tobe determined from the nearest neighbor !Can anybody help?Thanks in advancePatrick***************************************************The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary. If you are not the intended recipient, you may not review, copy or distribute this message. If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.***************************************************
_________________________________________________________________
It?s the same Hotmail?. If by ?same? you mean up to 70% faster.
http://windowslive.com/online/hotmail?ocid=TXT_TAGLM_WL_hotmail_acq_broad1_122008
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090106/29d5a093/attachment.html 

From arlindo.dasilva at GMAIL.COM  Tue Jan  6 13:30:17 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Tue, 6 Jan 2009 13:30:17 -0500
Subject: trouble plotting latest gshhs_land.shp ver.1.6 using shape
 extension
In-Reply-To: <BAY117-W3AD1358FB638CA7C1F10FF9DE0@phx.gbl>
Message-ID: <TUE.6.JAN.2009.133017.0500.>

On Tue, Jan 6, 2009 at 11:52 AM, Pablo Romero <romero619 at hotmail.com> wrote:

>  Hello,
> I just downloaded the latest version of the gshhs_* shapefiles, version
> 1.6, since this latest version uses polygons (the previous versions used
> 'lines') and can be used to do filled, hi-res basemaps.
>
> here's the url to the latest versions:
> *ftp://ftp.ihg.uni-duisburg.de/GIS/GISData/GSHHS/*<ftp://ftp.ihg.uni-duisburg.de/GIS/GISData/GSHHS/>
>
> I've tried plotting the gshhs_land.shp file in opengrads 1.9 on both
> windows&linux, and this shapefile is crashing grads in both cases.
> the gshhs_lakes.shp file works fine, but the land file crashes.
> can anyone else please test&confirm this?
>
>

Can you send a sample script reproducing the problem?



>
> Thanks,
> P.Romero
>
>
> ------------------------------
> Send e-mail anywhere. No map, no compass. Get your Hotmail(R) account now.<http://windowslive.com/oneline/hotmail?ocid=TXT_TAGLM_WL_hotmail_acq_anywhere_122008>
>



--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090106/0b320aa9/attachment.html 

From teddyallen at YAHOO.COM  Tue Jan  6 18:36:31 2009
From: teddyallen at YAHOO.COM (Teddy Allen)
Date: Tue, 6 Jan 2009 15:36:31 -0800
Subject: GrADS output to textfile?
Message-ID: <TUE.6.JAN.2009.153631.0800.>

Dear GrADS users,
I have seen traces of information leading to the creation of .txt files from GrADS output generated from NetCDF files. Unfortunately, I am having a bit of trouble trying to accomplish this. One method says to:
?
'set gxout fwrite'
'set fwrite file.txt'
'd [variable]'
'disable fwrite'
?
While another points to a script from the cookbook: http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file
?
?
I also try :
set gxout fwrite
d [variable]
disable fwrite
?
GrADS? then sends an output to my hard-drive, but it is not in an ascii .txt format.
Does anybody have a simple method to produce text output files from GrADS? I am not having any success in the above options. I am using GrADS 2.0 on a windows XP system. 
?
Thank you very much.
Sincerely,
teddy 


?

web: http://personal.ecu.edu/tla0923/HOME.html
"To move quickly, go alone. To move far, go together."
African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? *
?


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090106/4b49769a/attachment.html 

From smcmillan at PLANALYTICS.COM  Wed Jan  7 01:35:23 2009
From: smcmillan at PLANALYTICS.COM (Stephen R McMillan)
Date: Wed, 7 Jan 2009 00:35:23 -0600
Subject: GrADS output to textfile?
In-Reply-To: <456055.26544.qm@web51509.mail.re2.yahoo.com>
Message-ID: <WED.7.JAN.2009.003523.0600.>

Teddy,
The fwrite command produces binary, not text output.  To produce ascii 
text output of what you're currently displaying, use 'set gxout print'. 
You control the output with 'set prnopts....'  Display your data, then 
write to file using the write command.  Here's a simple example (assuming 
file is open and appropriately dimensioned):

'set gxout print'
'set prnopts %0.1f 10 1'
'd air-273.16'
dummy=write('c:/temp/output.txt',result)

Your output might look something like this:
Printing Grid -- 77 Values -- Undef = 32766
11.1 10.7 8.4 9.0 12.9 15.4 16.3 17.1 17.5 17.5 
17.6 8.3 8.6 7.0 6.3 8.2 10.3 11.7 12.4 13.0 
13.8 13.9 7.0 7.8 6.8 4.2 3.1 4.7 7.0 7.9
-2.2 0.1 2.4 3.2 3.5 0.8 0.2 -2.3 -4.3 -4.7

You could also output one grid value at a time without using 'set gxout 
print' and 'set prnopts.' The opengrads cookbook listing you referenced 
goes into more detail on the above and alternate methods. Additional 
references:
set gxout print:  http://grads.iges.org/grads/gadoc/gradcomdsetgxout.html
set prnopts:  http://grads.iges.org/grads/gadoc/gradcomdsetprnopts.html
write command:  http://grads.iges.org/grads/gadoc/script.html#intrinsic

Stephen McMillan





Teddy Allen <teddyallen at YAHOO.COM> 
Sent by: GRADSUSR at LIST.CINECA.IT
01/06/2009 05:36 PM
Please respond to
teddyallen at yahoo.com


To
GRADSUSR at LIST.CINECA.IT
cc

Subject
GrADS output to textfile?






Dear GrADS users,
I have seen traces of information leading to the creation of .txt files 
from GrADS output generated from NetCDF files. Unfortunately, I am having 
a bit of trouble trying to accomplish this. One method says to:
 
'set gxout fwrite'
'set fwrite file.txt'
'd [variable]'
'disable fwrite'
 
While another points to a script from the cookbook: 
http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file
 
 
I also try :
set gxout fwrite
d [variable]
disable fwrite
 
GrADS  then sends an output to my hard-drive, but it is not in an ascii 
.txt format.
Does anybody have a simple method to produce text output files from GrADS? 
I am not having any success in the above options. I am using GrADS 2.0 on 
a windows XP system. 
 
Thank you very much.
Sincerely,
teddy 


 
web: http://personal.ecu.edu/tla0923/HOME.html
"To move quickly, go alone. To move far, go together."
African 
Proverb                                                                                                        
*
 




***************************************************
The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary.  If you are not the intended recipient, you may not review, copy or distribute this message.  If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/db67d703/attachment.html 

From preuter at LABRI.FR  Wed Jan  7 03:55:51 2009
From: preuter at LABRI.FR (Patrick Reuter)
Date: Wed, 7 Jan 2009 09:55:51 +0100
Subject: 2nd try : Nearest defined neighbor function
In-Reply-To: <OF0619735A.5BE0BD5C-ON85257536.00558A3C-86257536.005C267F@planalytics.com>
Message-ID: <WED.7.JAN.2009.095551.0100.>

Dear Stephen and Pablo,

that is exactly what I was looking for, thank you. I discovered defval
and it does the job. I assume that now I have to loop through the XY
grid values, and "defval" the undefined values manually for each grid
point -- until somebody knows how to make a UDF out of it.

Happy new year
    Patrick


"Stephen R McMillan" <smcmillan at PLANALYTICS.COM> a ?crit :

> Patrick,
> What if an undefined point has more than one "nearest neighbor" and are of
> different values?  Do you want to assign an average of these, the highest,
> lowest, etc.?  I'm assuming each gridpoint has up to eight nearest
> neighbors in complex coastal areas.  For example, you could get the max
> nearest neighbor using something like this:
>
> 'define newvar= max(max(var,x=-152,x=-150),y=262,y=264)'
>
> Then assign result to point (in this case, x=-151, y=263) using the 'set
> defval var...' command.
>
> Similarly, you could use 'define var= min(min...' to get min value of
> nearest 8-neighbor or the "aave" command to get average of nearest
> neighbor.
>
> Hope this helps....
> Stephen McMillan
>
>
>
>
> Patrick Reuter <preuter at LABRI.FR>
> Sent by: GRADSUSR at LIST.CINECA.IT
> 01/06/2009 07:39 AM
> Please respond to
> GRADSUSR at LIST.CINECA.IT
>
>
> To
> GRADSUSR at LIST.CINECA.IT
> cc
>
> Subject
> 2nd try : Nearest defined neighbor function
>
>
>
>
>
> Sorry for multiples, message was somehow rejected ..
>
> Dear all,
>
> I have a variable which is only defined over the sea, elsewhere there
> are undefined values. Is there a function in Grads that can draw me a
> map where at locations of undefined values, instead of leaving it
> blank, draws me the value of the nearest defined value of the
> 8-neighborhood ? I know about the existence of 'smth9', but it does
> not exactly what I want, since it smoothes out, and 'const', since you
> have to specify a default value which I don't know since it should to
> be determined from the nearest neighbor !
>
> Can anybody help?
>
> Thanks in advance
>
> Patrick
>
>
>
>
> ***************************************************
> The information contained in this e-mail message is intended only
> for the use of the recipient(s) named above and may contain
> information that is privileged, confidential, and/or proprietary.
> If you are not the intended recipient, you may not review, copy or
> distribute this message.  If you have received this communication in
> error, please notify the sender immediately by e-mail, and delete
> the original message.
> ***************************************************


From teddyallen at YAHOO.COM  Wed Jan  7 07:51:24 2009
From: teddyallen at YAHOO.COM (Teddy Allen)
Date: Wed, 7 Jan 2009 04:51:24 -0800
Subject: GrADS output to textfile?
In-Reply-To: <OF07DEB0B6.343C4622-ON85257537.001F35DD-86257537.0024338D@planalytics.com>
Message-ID: <WED.7.JAN.2009.045124.0800.>

Thank you very much Stephen, but my version of GrADS (v2 windows xp) does not recognize the WRITE command. Are you using a newer or older version than the version 2?
Thank you again for your support. I appreciate it!
teddy 



web: http://personal.ecu.edu/tla0923/HOME.html
"To move quickly, go alone. To move far, go together."
African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? *
?

--- On Wed, 1/7/09, Stephen R McMillan <smcmillan at PLANALYTICS.COM> wrote:

From: Stephen R McMillan <smcmillan at PLANALYTICS.COM>
Subject: Re: GrADS output to textfile?
To: GRADSUSR at LIST.CINECA.IT
Date: Wednesday, January 7, 2009, 1:35 AM



Teddy,
The fwrite command produces binary, not text output. ?To produce ascii text output of what you're currently displaying, use 'set gxout print'. ?You control the output with 'set prnopts....' ?Display your data, then write to file using the write command. ?Here's a simple example (assuming file is open and appropriately dimensioned):

'set gxout print'
'set prnopts %0.1f 10 1'
'd air-273.16'
dummy=write('c:/temp/output.txt',result)

Your output might look something like this:
Printing Grid -- 77 Values -- Undef = 32766
11.1 10.7 8.4 9.0 12.9 15.4 16.3 17.1 17.5 17.5 
17.6 8.3 8.6 7.0 6.3 8.2 10.3 11.7 12.4 13.0 
13.8 13.9 7.0 7.8 6.8 4.2 3.1 4.7 7.0 7.9
-2.2 0.1 2.4 3.2 3.5 0.8 0.2 -2.3 -4.3 -4.7

You could also output one grid value at a time without using 'set gxout print' and 'set prnopts.' The opengrads cookbook listing you referenced goes into more detail on the above and alternate methods. Additional references:
set gxout print: ?http://grads.iges.org/grads/gadoc/gradcomdsetgxout.html
set prnopts: ?http://grads.iges.org/grads/gadoc/gradcomdsetprnopts.html
write command: ?http://grads.iges.org/grads/gadoc/script.html#intrinsic

Stephen McMillan








Teddy Allen <teddyallen at YAHOO.COM>
Sent by: GRADSUSR at LIST.CINECA.IT
01/06/2009 05:36 PM




Please respond to
teddyallen at yahoo.com






To
GRADSUSR at LIST.CINECA.IT


cc



Subject
GrADS output to textfile?












Dear GrADS users,
I have seen traces of information leading to the creation of .txt files from GrADS output generated from NetCDF files. Unfortunately, I am having a bit of trouble trying to accomplish this. One method says to:
?
'set gxout fwrite'
'set fwrite file.txt'
'd [variable]'
'disable fwrite'
?
While another points to a script from the cookbook: http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file
?
?
I also try :
set gxout fwrite
d [variable]
disable fwrite
?
GrADS? then sends an output to my hard-drive, but it is not in an ascii .txt format.
Does anybody have a simple method to produce text output files from GrADS? I am not having any success in the above options. I am using GrADS 2.0 on a windows XP system. 
?
Thank you very much.
Sincerely,
teddy 


?
web: http://personal.ecu.edu/tla0923/HOME.html
"To move quickly, go alone. To move far, gotogether."
African Proverb??????????????????????????????????????????????????????????????????????????????????????????????????????? *
?



***************************************************
The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary. If you are not the intended recipient, you may not review, copy or distribute this message. If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.
***************************************************



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/7b291fd1/attachment.html 

From littaaj at REDIFFMAIL.COM  Wed Jan  7 10:19:33 2009
From: littaaj at REDIFFMAIL.COM (Litta AJ)
Date: Wed, 7 Jan 2009 15:19:33 -0000
Subject: help
Message-ID: <WED.7.JAN.2009.151933.0000.LITTAAJ@REDIFFMAIL.COM>

Dear members,
 ?
I am writing this mail to get an information. How can we plot black and white figures in Grads. I want to plot some of my model results in black and white. 

regards,
Litta
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/ff84983c/attachment.html 

From feudale at ICTP.IT  Wed Jan  7 10:37:11 2009
From: feudale at ICTP.IT (Laura Feudale)
Date: Wed, 7 Jan 2009 16:37:11 +0100
Subject: help
In-Reply-To: <20090107151933.55615.qmail@f4mail213.rediffmail.com>
Message-ID: <WED.7.JAN.2009.163711.0100.>

Dear Litta,
at the beginning of your script write

'set display greyscale white'
'c'

and then display your variable.
Regards,
-- Laura

-----------------------------------------------------------------
 Laura Feudale, Ph.D. (Climate Dynamics, George Mason University)
*Current Affiliation:            Earth System Physics Section
 The Abdus Salam International Centre for Theoretical Physics
 Strada Costiera 11, 34014 Trieste, ITALY
 http://www.ictp.it
 mailto: feudale at ictp.it
-----------------------------------------------------------------



Litta AJ wrote:
>
> Dear members,
>
> I am writing this mail to get an information. How can we plot black
> and white figures in Grads. I want to plot some of my model results in
> black and white.
>
> regards,
> Litta
>
>
>


From smcmillan at PLANALYTICS.COM  Wed Jan  7 11:00:06 2009
From: smcmillan at PLANALYTICS.COM (Stephen R McMillan)
Date: Wed, 7 Jan 2009 10:00:06 -0600
Subject: GrADS output to textfile?
In-Reply-To: <509415.98785.qm@web51506.mail.re2.yahoo.com>
Message-ID: <WED.7.JAN.2009.100006.0600.>

Teddy....I'm using v1.9.0-rc1 on windows xp.  I'm using 'write' from a 
script.  Doesn't v2 allow it from a script?  Stephen Mc



Teddy Allen <teddyallen at YAHOO.COM> 
Sent by: GRADSUSR at LIST.CINECA.IT
01/07/2009 06:51 AM
Please respond to
teddyallen at yahoo.com


To
GRADSUSR at LIST.CINECA.IT
cc

Subject
Re: GrADS output to textfile?






Thank you very much Stephen, but my version of GrADS (v2 windows xp) does 
not recognize the WRITE command. Are you using a newer or older version 
than the version 2?
Thank you again for your support. I appreciate it!
teddy 

web: http://personal.ecu.edu/tla0923/HOME.html
"To move quickly, go alone. To move far, go together."
African 
Proverb                                                                                                        
*
 


--- On Wed, 1/7/09, Stephen R McMillan <smcmillan at PLANALYTICS.COM> wrote:
From: Stephen R McMillan <smcmillan at PLANALYTICS.COM>
Subject: Re: GrADS output to textfile?
To: GRADSUSR at LIST.CINECA.IT
Date: Wednesday, January 7, 2009, 1:35 AM


Teddy,
The fwrite command produces binary, not text output.  To produce ascii 
text output of what you're currently displaying, use 'set gxout print'. 
 You control the output with 'set prnopts....'  Display your data, then 
write to file using the write command.  Here's a simple example (assuming 
file is open and appropriately dimensioned):

'set gxout print'
'set prnopts %0.1f 10 1'
'd air-273.16'
dummy=write('c:/temp/output.txt',result)

Your output might look something like this:
Printing Grid -- 77 Values -- Undef = 32766
11.1 10.7 8.4 9.0 12.9 15.4 16.3 17.1 17.5 17.5 
17.6 8.3 8.6 7.0 6.3 8.2 10.3 11.7 12.4 13.0 
13.8 13.9 7.0 7.8 6.8 4.2 3.1 4.7 7.0 7.9
-2.2 0.1 2.4 3.2 3.5 0.8 0.2 -2.3 -4.3 -4.7

You could also output one grid value at a time without using 'set gxout 
print' and 'set prnopts.' The opengrads cookbook listing you referenced 
goes into more detail on the above and alternate methods. Additional 
references:
set gxout print:  http://grads.iges.org/grads/gadoc/gradcomdsetgxout.html
set prnopts:  http://grads.iges.org/grads/gadoc/gradcomdsetprnopts.html
write command:  http://grads.iges.org/grads/gadoc/script.html#intrinsic

Stephen McMillan





Teddy Allen <teddyallen at YAHOO.COM>
Sent by: GRADSUSR at LIST.CINECA.IT 
01/06/2009 05:36 PM 
Please respond to
teddyallen at yahoo.com


To
GRADSUSR at LIST.CINECA.IT 
cc

Subject
GrADS output to textfile?







Dear GrADS users,
I have seen traces of information leading to the creation of .txt files 
from GrADS output generated from NetCDF files. Unfortunately, I am having 
a bit of trouble trying to accomplish this. One method says to:
 
'set gxout fwrite'
'set fwrite file.txt'
'd [variable]'
'disable fwrite'
 
While another points to a script from the cookbook: 
http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file
 
 
I also try :
set gxout fwrite
d [variable]
disable fwrite
 
GrADS  then sends an output to my hard-drive, but it is not in an ascii 
.txt format.
Does anybody have a simple method to produce text output files from GrADS? 
I am not having any success in the above options. I am using GrADS 2.0 on 
a windows XP system. 
 
Thank you very much.
Sincerely,
teddy 


 
web: http://personal.ecu.edu/tla0923/HOME.html
"To move quickly, go alone. To move far, gotogether."
African 
Proverb                                                                                                        
*
 



***************************************************
The information contained in this e-mail message is intended only for the 
use of the recipient(s) named above and may contain information that is 
privileged, confidential, and/or proprietary. If you are not the intended 
recipient, you may not review, copy or distribute this message. If you 
have received this communication in error, please notify the sender 
immediately by e-mail, and delete the original message.
***************************************************




***************************************************
The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary.  If you are not the intended recipient, you may not review, copy or distribute this message.  If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/4ce36324/attachment.html 

From romero619 at HOTMAIL.COM  Wed Jan  7 13:07:22 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Wed, 7 Jan 2009 10:07:22 -0800
Subject: aave in pcgrads19, error: function not found?
Message-ID: <WED.7.JAN.2009.100722.0800.>


its the first time Ive tried using the aave function, but I received an error in grads saying "function/comand not found". has this function been removed?
 
P.Romero
_________________________________________________________________
Windows LiveTM: Keep your life in sync. 
http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/d71bdab1/attachment.html 

From sajjad_met at HOTMAIL.COM  Wed Jan  7 14:03:55 2009
From: sajjad_met at HOTMAIL.COM (sajjad saeed)
Date: Thu, 8 Jan 2009 00:03:55 +0500
Subject: eof
In-Reply-To: <4964CC27.4090609@ictp.it>
Message-ID: <THU.8.JAN.2009.000355.0500.>


Hi all,

Is it possible to calculate EOFs using Grads? 
If yes! I will appreciate if any body have any script to calculate the eofs and share it?

Thankx

Sajjad


_________________________________________________________________
More than messages?check out the rest of the Windows Live?.
http://www.microsoft.com/windows/windowslive/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/08290e44/attachment.html 

From matias at TROPOSFERA.ATMOSFCU.UNAM.MX  Wed Jan  7 14:28:13 2009
From: matias at TROPOSFERA.ATMOSFCU.UNAM.MX (J. Matias Mendez Perez)
Date: Wed, 7 Jan 2009 13:28:13 -0600
Subject: eof
In-Reply-To: <BAY128-W30F198B64A2207644A510EFBDF0@phx.gbl>
Message-ID: <WED.7.JAN.2009.132813.0600.>

Hi Sajjad,

Matthias M?nnich has developed a tool to compute EOFs:

http://www.atmos.ucla.edu/~munnich/Grads/EOF/

Good luck,

Matias

>
> Hi all,
>
> Is it possible to calculate EOFs using Grads?
> If yes! I will appreciate if any body have any script to calculate the
> eofs and share it?
>
> Thankx
>
> Sajjad
>
>
> _________________________________________________________________
> More than messages?check out the rest of the Windows Live?.
> http://www.microsoft.com/windows/windowslive/


From sajjad_met at HOTMAIL.COM  Wed Jan  7 15:23:26 2009
From: sajjad_met at HOTMAIL.COM (sajjad saeed)
Date: Thu, 8 Jan 2009 01:23:26 +0500
Subject: eof
In-Reply-To: <fe968ef5fbdfbf37b7ce2ee1dfbd0f9f.squirrel@atmosfera.unam.mx>
Message-ID: <THU.8.JAN.2009.012326.0500.>


Dear Matias,
 
Thankyou very much.
 
Sajjad> Date: Wed, 7 Jan 2009 13:28:13 -0600> From: matias at TROPOSFERA.ATMOSFCU.UNAM.MX> Subject: Re: eof> To: GRADSUSR at LIST.CINECA.IT> > Hi Sajjad,> > Matthias M?nnich has developed a tool to compute EOFs:> > http://www.atmos.ucla.edu/~munnich/Grads/EOF/> > Good luck,> > Matias> > >> > Hi all,> >> > Is it possible to calculate EOFs using Grads?> > If yes! I will appreciate if any body have any script to calculate the> > eofs and share it?> >> > Thankx> >> > Sajjad> >> >> > _________________________________________________________________> > More than messages?check out the rest of the Windows Live?.> > http://www.microsoft.com/windows/windowslive/
_________________________________________________________________
More than messages?check out the rest of the Windows Live?.
http://www.microsoft.com/windows/windowslive/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/73420194/attachment.html 

From scchan at ATMOS.UMD.EDU  Wed Jan  7 15:05:47 2009
From: scchan at ATMOS.UMD.EDU (Steven C Chan)
Date: Wed, 7 Jan 2009 15:05:47 -0500
Subject: eof
In-Reply-To: <fe968ef5fbdfbf37b7ce2ee1dfbd0f9f.squirrel@atmosfera.unam.mx>
Message-ID: <WED.7.JAN.2009.150547.0500.>

If I have not mistaken, Grads 2.* does not support UDF. So that may not
work with your version of Grads (of course you can always install an
older Grads version in your own directory).

Steven Chan

J. Matias Mendez Perez wrote:
> Hi Sajjad,
>
> Matthias M?nnich has developed a tool to compute EOFs:
>
> http://www.atmos.ucla.edu/~munnich/Grads/EOF/
>
> Good luck,
>
> Matias
>
>
>> Hi all,
>>
>> Is it possible to calculate EOFs using Grads?
>> If yes! I will appreciate if any body have any script to calculate the
>> eofs and share it?
>>
>> Thankx
>>
>> Sajjad
>>
>>
>> _________________________________________________________________
>> More than messages?check out the rest of the Windows Live?.
>> http://www.microsoft.com/windows/windowslive/
>>


From matias at TROPOSFERA.ATMOSFCU.UNAM.MX  Wed Jan  7 16:25:13 2009
From: matias at TROPOSFERA.ATMOSFCU.UNAM.MX (J. Matias Mendez Perez)
Date: Wed, 7 Jan 2009 15:25:13 -0600
Subject: eof - fortran code
In-Reply-To: <49650B1B.5080303@atmos.umd.edu>
Message-ID: <WED.7.JAN.2009.152513.0600.>

Hello everybody,

Another option is to use a fortran code to compute EOFs and then to
display using GrADS. I have used the following code for EOFs and rotated
EOFs:

http://www.lasg.ac.cn/staff/ljp/subroutine/REOF.f

Good luck,

Matias

> If I have not mistaken, Grads 2.* does not support UDF. So that may not
> work with your version of Grads (of course you can always install an
> older Grads version in your own directory).
>
> Steven Chan
>
> J. Matias Mendez Perez wrote:
>> Hi Sajjad,
>>
>> Matthias M?nnich has developed a tool to compute EOFs:
>>
>> http://www.atmos.ucla.edu/~munnich/Grads/EOF/
>>
>> Good luck,
>>
>> Matias
>>
>>
>>> Hi all,
>>>
>>> Is it possible to calculate EOFs using Grads?
>>> If yes! I will appreciate if any body have any script to calculate the
>>> eofs and share it?
>>>
>>> Thankx
>>>
>>> Sajjad
>>>
>>>
>>> _________________________________________________________________
>>> More than messages?check out the rest of the Windows Live?.
>>> http://www.microsoft.com/windows/windowslive/
>>>
>


From arlindo.dasilva at GMAIL.COM  Wed Jan  7 17:59:08 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Wed, 7 Jan 2009 17:59:08 -0500
Subject: eof
In-Reply-To: <BAY128-W30F198B64A2207644A510EFBDF0@phx.gbl>
Message-ID: <WED.7.JAN.2009.175908.0500.>

On Wed, Jan 7, 2009 at 2:03 PM, sajjad saeed <sajjad_met at hotmail.com> wrote:

>  Hi all,
>
> Is it possible to calculate EOFs using Grads?
> If yes! I will appreciate if any body have any script to calculate the eofs
> and share it?
>

If you are willing to use Python, PyGrADS provides a function to compute
EOFs. You can find an example at the wiki:

http://opengrads.org/wiki/index.php?title=PyGrADS_Interactive_Shell#Computing_EOFs

PyGrADS works with both GrADS v1.9-rc1 and partially with GrADS v2.0.a4.
 The main restriction of PyGrADS with GrADS v2 is that you cannot import
data directly back into GrADS, but you can export data and compute EOFs, and
either do your plots right in python or write the data to a file and read it
back into GrADS.

     Arlindo

Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/53df4834/attachment.html 

From romero619 at HOTMAIL.COM  Wed Jan  7 18:28:06 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Wed, 7 Jan 2009 15:28:06 -0800
Subject: listserver rejecting messages
Message-ID: <WED.7.JAN.2009.152806.0800.>


is anyone else experiencing this problem?
_________________________________________________________________
Windows LiveTM: Keep your life in sync. 
http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t1_allup_howitworks_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/49781772/attachment.html 

From romero619 at HOTMAIL.COM  Wed Jan  7 18:30:09 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Wed, 7 Jan 2009 15:30:09 -0800
Subject: scripting question about 'define' and 'set defval'
Message-ID: <WED.7.JAN.2009.153009.0800.>


I think Im missing something obvious, but I cant feed a variable created using 'define' into the 'set defval' function...right now, Im using 'query define' in order to retrieve the value, but Im hoping there's a better, easier, more reliable way... 'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')''q define' *assuming newvar is the ONLY defined variable, I can do this...'theVal=sublin(result,1)'theVal=subwrd(theVal,2)'set defval ht 'i' 'j' 'theVal is there a simpler way to access my 'newvar' value, without having to 'query define' in order to get it?Please help, thanks...P.Romero
_________________________________________________________________
Windows LiveTM: Keep your life in sync. 
http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090107/7f5b3e80/attachment.html 

From smcmillan at PLANALYTICS.COM  Thu Jan  8 01:13:28 2009
From: smcmillan at PLANALYTICS.COM (Stephen R McMillan)
Date: Thu, 8 Jan 2009 00:13:28 -0600
Subject: scripting question about 'define' and 'set defval'
In-Reply-To: <BAY117-W135551BD13E6131B3435D6F9DF0@phx.gbl>
Message-ID: <THU.8.JAN.2009.001328.0600.>

Pablo,
I believe you should be able to do your 'set defval ht 'i' 'j' 'newvar 
immediately following the 'define newvar...' line.
Stephen



Pablo Romero <romero619 at HOTMAIL.COM> 
Sent by: GRADSUSR at LIST.CINECA.IT
01/07/2009 05:30 PM
Please respond to
GRADSUSR at LIST.CINECA.IT


To
GRADSUSR at LIST.CINECA.IT
cc

Subject
scripting question about 'define' and 'set defval'





I think Im missing something obvious, but I cant feed a variable created 
using 'define' into the 'set defval' function...
right now, Im using 'query define' in order to retrieve the value, but Im 
hoping there's a better, easier, more reliable way...
 
'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')'
'q define'
 
*assuming newvar is the ONLY defined variable, I can do this...
'theVal=sublin(result,1)
'theVal=subwrd(theVal,2)
'set defval ht 'i' 'j' 'theVal
 
is there a simpler way to access my 'newvar' value, without having to 
'query define' in order to get it?
Please help, thanks...
P.Romero

Windows LiveTM: Keep your life in sync.  Check it out.



***************************************************
The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary.  If you are not the intended recipient, you may not review, copy or distribute this message.  If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/b58412b4/attachment.html 

From romero619 at HOTMAIL.COM  Thu Jan  8 10:33:20 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Thu, 8 Jan 2009 07:33:20 -0800
Subject: scripting question about 'define' and 'set defval'
In-Reply-To: <OF2F672D56.820CF700-ON85257538.00177A60-86257538.002231D1@planalytics.com>
Message-ID: <THU.8.JAN.2009.073320.0800.>


Stephen,
I tried that initially, but it didnt work.
I even tried say newvar and to see if it recognized newvar but it just echoed "newvar" at the prompt.
can you/anyone else test and confirm this?
here are the exact lines from the script Im using...
 
'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')'
'set defval ht 'i' 'j' 'newvar*  'q define'* tmpvar=sublin(result,2)* tmpvar=subwrd(tmpvar,2)* 'set defval myvar 'i' 'j' 'tmpvar again, the above didnt work, 'set defval' throws an error because 'newvar' is empty/undefined.
P.Romero



Date: Thu, 8 Jan 2009 00:13:28 -0600From: smcmillan at PLANALYTICS.COMSubject: Re: scripting question about 'define' and 'set defval'To: GRADSUSR at LIST.CINECA.ITPablo,I believe you should be able to do your 'set defval ht 'i' 'j' 'newvar  immediately following the 'define newvar...' line.Stephen



Pablo Romero <romero619 at HOTMAIL.COM>Sent by: GRADSUSR at LIST.CINECA.IT
01/07/2009 05:30 PM




Please respond toGRADSUSR at LIST.CINECA.IT





To
GRADSUSR at LIST.CINECA.IT


cc



Subject
scripting question about 'define' and 'set defval'




I think Im missing something obvious, but I cant feed a variable created using 'define' into the 'set defval' function...right now, Im using 'query define' in order to retrieve the value, but Im hoping there's a better, easier, more reliable way... 'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')''q define' *assuming newvar is the ONLY defined variable, I can do this...'theVal=sublin(result,1)'theVal=subwrd(theVal,2)'set defval ht 'i' 'j' 'theVal is there a simpler way to access my 'newvar' value, without having to 'query define' in order to get it?Please help, thanks...P.Romero

Windows LiveTM: Keep your life in sync.  Check it out.***************************************************The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary. If you are not the intended recipient, you may not review, copy or distribute this message. If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.***************************************************
_________________________________________________________________
Windows Live?: Keep your life in sync. 
http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/ee5bc964/attachment.html 

From smcmillan at PLANALYTICS.COM  Thu Jan  8 10:37:23 2009
From: smcmillan at PLANALYTICS.COM (Stephen R McMillan)
Date: Thu, 8 Jan 2009 09:37:23 -0600
Subject: scripting question about 'define' and 'set defval'
In-Reply-To: <BAY117-W2960C24CF7FA5CBBC285A3F9DC0@phx.gbl>
Message-ID: <THU.8.JAN.2009.093723.0600.>

Where did the variable 'ht' come from?  Did you previously define it as 
your field variable?  Stephen Mc



Pablo Romero <romero619 at HOTMAIL.COM> 
Sent by: GRADSUSR at LIST.CINECA.IT
01/08/2009 09:33 AM
Please respond to
GRADSUSR at LIST.CINECA.IT


To
GRADSUSR at LIST.CINECA.IT
cc

Subject
Re: scripting question about 'define' and 'set defval'





Stephen,
I tried that initially, but it didnt work.
I even tried say newvar and to see if it recognized newvar but it just 
echoed "newvar" at the prompt.
can you/anyone else test and confirm this?
here are the exact lines from the script Im using...
 
'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')'
'set defval ht 'i' 'j' 'newvar
*  'q define'
* tmpvar=sublin(result,2)
* tmpvar=subwrd(tmpvar,2)
* 'set defval myvar 'i' 'j' 'tmpvar
 
again, the above didnt work, 'set defval' throws an error because 'newvar' 
is empty/undefined.
P.Romero


Date: Thu, 8 Jan 2009 00:13:28 -0600
From: smcmillan at PLANALYTICS.COM
Subject: Re: scripting question about 'define' and 'set defval'
To: GRADSUSR at LIST.CINECA.IT


Pablo,
I believe you should be able to do your 'set defval ht 'i' 'j' 'newvar 
 immediately following the 'define newvar...' line.
Stephen




Pablo Romero <romero619 at HOTMAIL.COM>
Sent by: GRADSUSR at LIST.CINECA.IT 01/07/2009 05:30 PM 
Please respond to
GRADSUSR at LIST.CINECA.IT




To
GRADSUSR at LIST.CINECA.IT 
cc

Subject
scripting question about 'define' and 'set defval'






I think Im missing something obvious, but I cant feed a variable created 
using 'define' into the 'set defval' function...
right now, Im using 'query define' in order to retrieve the value, but Im 
hoping there's a better, easier, more reliable way...
 
'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')'
'q define'
 
*assuming newvar is the ONLY defined variable, I can do this...
'theVal=sublin(result,1)
'theVal=subwrd(theVal,2)
'set defval ht 'i' 'j' 'theVal
 
is there a simpler way to access my 'newvar' value, without having to 
'query define' in order to get it?
Please help, thanks...
P.Romero


Windows LiveTM: Keep your life in sync.  Check it out.


***************************************************
The information contained in this e-mail message is intended only for the 
use of the recipient(s) named above and may contain information that is 
privileged, confidential, and/or proprietary. If you are not the intended 
recipient, you may not review, copy or distribute this message. If you 
have received this communication in error, please notify the sender 
immediately by e-mail, and delete the original message.
***************************************************


Windows Live?: Keep your life in sync.  Check it out.



***************************************************
The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary.  If you are not the intended recipient, you may not review, copy or distribute this message.  If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.
***************************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/0d7f4a79/attachment.html 

From jhuddleston at HUGHES.NET  Thu Jan  8 10:20:54 2009
From: jhuddleston at HUGHES.NET (John Huddleston)
Date: Thu, 8 Jan 2009 08:20:54 -0700
Subject: aura he5 data exceeds dimension
Message-ID: <THU.8.JAN.2009.082054.0700.>

Hi



I am accessing
http://acdisc.sci.gsfc.nasa.gov/opendap/Aura_OMI_Level3/OMAEROe.003/2006/OMI
-Aura_L3-OMAEROe_2006m0101_v003-2008m0815t092235.he5 with the latest v2.0a4
gradsdap.



It can open the file and the 'q file' works fine; however, when I display
some of the aerosol optical variables, there is an error:



gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds dimension
bound



Data Request Error:  Error for variable 'aerosolopticalt'

  Error ocurred at column 1

DISPLAY error:  Invalid expression

  Expression = aerosolopticalt



Is this a problem with the data or bounds limitation in the gradsdap
application?



Here is the information on the gradsdap startup



Starting X server under c:\GrADS20/Xming

Starting gradsdap under c:\GrADS20\win32 ...



Grid Analysis and Display System (GrADS) Version 2.0.a4.oga.1

Copyright (c) 1988-2008 by Brian Doty and the

Institute for Global Environment and Society (IGES)

GrADS comes with ABSOLUTELY NO WARRANTY

See file COPYRIGHT for more information



Config: v2.0.a3.oga.1 little-endian readline printim grib2 netcdf hdf4-sds
opendap-grids,stn Athena



John Huddleston















-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/27df063b/attachment.html 

From romero619 at HOTMAIL.COM  Thu Jan  8 10:46:21 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Thu, 8 Jan 2009 07:46:21 -0800
Subject: scripting question about 'define' and 'set defval'
In-Reply-To: <OF9A4CF1FE.29889469-ON85257538.0055A992-86257538.0055EB89@planalytics.com>
Message-ID: <THU.8.JAN.2009.074621.0800.>


sorry, typo...I changed the names in the script during some testing and missed in within the email.
but the problem still exists, I can't find a simple way to just access the defined variable.
here it is again...
 
  'define tmp2=aave(tmp1,x='x1',x='x2',y='y1',y='y2')'  'set defval tmp1 'i' 'j' 'tmp2  say   'set defval tmp1 'i' 'j' 'tmp2
 
the last line prints something like "set defval tmp1 -65 100 tmp2", instead of substituting the value of tmp2.
 
the variable names are tmp1 and tmp2, though its irrelevant as long as theyre consistent (I already checked that they are, in fact, correct, since I can see & access them through 'q define')



Date: Thu, 8 Jan 2009 09:37:23 -0600From: smcmillan at PLANALYTICS.COMSubject: Re: scripting question about 'define' and 'set defval'To: GRADSUSR at LIST.CINECA.ITWhere did the variable 'ht' come from?  Did you previously define it as your field variable?  Stephen Mc



Pablo Romero <romero619 at HOTMAIL.COM>Sent by: GRADSUSR at LIST.CINECA.IT
01/08/2009 09:33 AM




Please respond toGRADSUSR at LIST.CINECA.IT





To
GRADSUSR at LIST.CINECA.IT


cc



Subject
Re: scripting question about 'define' and 'set defval'




Stephen,I tried that initially, but it didnt work.I even tried say newvar and to see if it recognized newvar but it just echoed "newvar" at the prompt.can you/anyone else test and confirm this?here are the exact lines from the script Im using... 'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')''set defval ht 'i' 'j' 'newvar*  'q define'* tmpvar=sublin(result,2)* tmpvar=subwrd(tmpvar,2)* 'set defval myvar 'i' 'j' 'tmpvar again, the above didnt work, 'set defval' throws an error because 'newvar' is empty/undefined.P.Romero

Date: Thu, 8 Jan 2009 00:13:28 -0600From: smcmillan at PLANALYTICS.COMSubject: Re: scripting question about 'define' and 'set defval'To: GRADSUSR at LIST.CINECA.ITPablo,I believe you should be able to do your 'set defval ht 'i' 'j' 'newvar  immediately following the 'define newvar...' line.Stephen



Pablo Romero <romero619 at HOTMAIL.COM>Sent by: GRADSUSR at LIST.CINECA.IT01/07/2009 05:30 PM




Please respond toGRADSUSR at LIST.CINECA.IT





To
GRADSUSR at LIST.CINECA.IT


cc



Subject
scripting question about 'define' and 'set defval'




I think Im missing something obvious, but I cant feed a variable created using 'define' into the 'set defval' function...right now, Im using 'query define' in order to retrieve the value, but Im hoping there's a better, easier, more reliable way... 'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')''q define' *assuming newvar is the ONLY defined variable, I can do this...'theVal=sublin(result,1)'theVal=subwrd(theVal,2)'set defval ht 'i' 'j' 'theVal is there a simpler way to access my 'newvar' value, without having to 'query define' in order to get it?Please help, thanks...P.Romero

Windows LiveTM: Keep your life in sync.  Check it out.***************************************************The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary. If you are not the intended recipient, you may not review, copy or distribute this message. If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.***************************************************

Windows Live?: Keep your life in sync.  Check it out.***************************************************The information contained in this e-mail message is intended only for the use of the recipient(s) named above and may contain information that is privileged, confidential, and/or proprietary. If you are not the intended recipient, you may not review, copy or distribute this message. If you have received this communication in error, please notify the sender immediately by e-mail, and delete the original message.***************************************************
_________________________________________________________________
Windows Live? Hotmail?: Chat. Store. Share. Do more with mail. 
http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t1_hm_justgotbetter_howitworks_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/9c94a740/attachment.html 

From scpiper at UCSD.EDU  Thu Jan  8 10:32:25 2009
From: scpiper at UCSD.EDU (Stephen Piper)
Date: Thu, 8 Jan 2009 07:32:25 -0800
Subject: batch output for list of files and lat-lons
Message-ID: <THU.8.JAN.2009.073225.0800.>

Hi,

I want to extract annual average vertical profiles from a monthly
binary file for a list of lat-lon locations.
I have a working .ctl file that works fine if I invoke it
interactively in Mac Unix:

~piper/bin/grads
g>open my.ctl
g>set lat 45
g>set lon
g>set z 1 9
g> set gxout print
g>d ave(bioa,t=1,t=12)

Now how can I do this for a long list of lat-lons for several
different binary files instead of
typing each one in interactively?  ie., how can I invoke grads so
that I can give it a prepared
list, and don't have
to interactively enter each file and lat-lon pair?

Thanks,
Steve


From michael.l.sestak at NAVY.MIL  Thu Jan  8 10:59:39 2009
From: michael.l.sestak at NAVY.MIL (Sestak, Michael L CIV 63134)
Date: Thu, 8 Jan 2009 07:59:39 -0800
Subject: batch output for list of files and lat-lons
In-Reply-To: <7F80D088-CA93-4DD8-A9D5-39E7D7A69690@ucsd.edu>
Message-ID: <THU.8.JAN.2009.075939.0800.>

Look at grads scripting in the online manual.  You can use a while loop to read your lat lon file using thr read function and then use ordinary grads commands to make your displays.



From: Stephen Piper
Sent: Thu 1/8/2009 7:32
To: GRADSUSR at LIST.CINECA.IT
Subject: batch output for list of files and lat-lons


Hi,

I want to extract annual average vertical profiles from a monthly
binary file for a list of lat-lon locations.
I have a working .ctl file that works fine if I invoke it
interactively in Mac Unix:

~piper/bin/grads
g>open my.ctl
g>set lat 45
g>set lon
g>set z 1 9
g> set gxout print
g>d ave(bioa,t=1,t=12)

Now how can I do this for a long list of lat-lons for several
different binary files instead of
typing each one in interactively?  ie., how can I invoke grads so
that I can give it a prepared
list, and don't have
to interactively enter each file and lat-lon pair?

Thanks,
Steve
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/9877ccbb/attachment.html 

From AHassankunju at FLAGSTONERE.BM  Thu Jan  8 10:55:10 2009
From: AHassankunju at FLAGSTONERE.BM (Annes Hassankunju)
Date: Thu, 8 Jan 2009 11:55:10 -0400
Subject: scripting question about 'define' and 'set defval'
In-Reply-To: <BAY117-W2960C24CF7FA5CBBC285A3F9DC0@phx.gbl>
Message-ID: <THU.8.JAN.2009.115510.0400.>

signoff gradsusr

 

 

 


Where did the variable 'ht' come from?  Did you previously define it as
your field variable?  Stephen Mc



Pablo Romero <romero619 at HOTMAIL.COM>
Sent by: GRADSUSR at LIST.CINECA.IT

01/08/2009 09:33 AM

Please respond to
GRADSUSR at LIST.CINECA.IT

To

GRADSUSR at LIST.CINECA.IT

cc

	
Subject

Re: scripting question about 'define' and 'set defval'

 

        	



Stephen,
I tried that initially, but it didnt work.
I even tried say newvar and to see if it recognized newvar but it just
echoed "newvar" at the prompt.
can you/anyone else test and confirm this?
here are the exact lines from the script Im using...
 
'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')'
'set defval ht 'i' 'j' 'newvar
*  'q define'
* tmpvar=sublin(result,2)
* tmpvar=subwrd(tmpvar,2)
* 'set defval myvar 'i' 'j' 'tmpvar
 
again, the above didnt work, 'set defval' throws an error because
'newvar' is empty/undefined.
P.Romero

________________________________



Date: Thu, 8 Jan 2009 00:13:28 -0600
From: smcmillan at PLANALYTICS.COM
Subject: Re: scripting question about 'define' and 'set defval'
To: GRADSUSR at LIST.CINECA.IT


Pablo,
I believe you should be able to do your 'set defval ht 'i' 'j' 'newvar
immediately following the 'define newvar...' line.
Stephen




Pablo Romero <romero619 at HOTMAIL.COM>
Sent by: GRADSUSR at LIST.CINECA.IT01/07/2009 05:30 PM

Please respond to
GRADSUSR at LIST.CINECA.IT





To

GRADSUSR at LIST.CINECA.IT

cc

	
Subject

scripting question about 'define' and 'set defval'

 

        	



I think Im missing something obvious, but I cant feed a variable created
using 'define' into the 'set defval' function...
right now, Im using 'query define' in order to retrieve the value, but
Im hoping there's a better, easier, more reliable way...
 
'define newvar=aave(myvar,x='x1',x='x2',y='y1',y='y2')'
'q define'
 
*assuming newvar is the ONLY defined variable, I can do this...
'theVal=sublin(result,1)
'theVal=subwrd(theVal,2)
'set defval ht 'i' 'j' 'theVal
 
is there a simpler way to access my 'newvar' value, without having to
'query define' in order to get it?
Please help, thanks...
P.Romero

________________________________



Windows LiveTM: Keep your life in sync.  Check it out.
<http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_01200
9> 


***************************************************
The information contained in this e-mail message is intended only for
the use of the recipient(s) named above and may contain information that
is privileged, confidential, and/or proprietary. If you are not the
intended recipient, you may not review, copy or distribute this message.
If you have received this communication in error, please notify the
sender immediately by e-mail, and delete the original message.
***************************************************

________________________________



Windows Live(tm): Keep your life in sync.  Check it out.
<http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_allup_explore_01200
9> 


***************************************************
The information contained in this e-mail message is intended only for
the use of the recipient(s) named above and may contain information that
is privileged, confidential, and/or proprietary. If you are not the
intended recipient, you may not review, copy or distribute this message.
If you have received this communication in error, please notify the
sender immediately by e-mail, and delete the original message.
***************************************************


**********************************************************************
This communication contains information which is confidential and may also be legally privileged. It is for the exclusive use of the intended recipient(s). If you are not the intended recipient(s), disclosure, copying, distribution, or other use of, or action taken or omitted to be taken in reliance upon, this communication or the information in it is prohibited and maybe unlawful. If you have received this communication in error please notify the sender by return email, delete it from your system and destroy any copies.
**********************************************************************

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/081067ef/attachment.html 

From jma at COLA.IGES.ORG  Thu Jan  8 12:13:36 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Thu, 8 Jan 2009 12:13:36 -0500
Subject: aura he5 data exceeds dimension
In-Reply-To: <BE54F6D357774DA4BAADC20597A0BE8C@JOHNIYWTKBFEG5>
Message-ID: <THU.8.JAN.2009.121336.0500.>

Hi, John --
There is an extra dimension for those variables that the sdfopen code  
doesn't find because it's not 'coards compliant'. Check the ncdump  
output and you'll find:

dimensions:
         nWavelDiagnostic = 5 ;
         lat = 720 ;
         lon = 1440 ;
data:
  nWavelDiagnostic = 0, 1, 2, 3, 4 ;

So you have to tell GrADS how to map that "nWavelDiagnostic"  
dimension. I'd use EDEF for that. Plus, the variable names are very  
long, so you have to alias them to something more meaningful to you.  
Once you open the file, check the 'query attr' output for the titles  
of these variables -- they are also very long. Anyway, here's a bare- 
bones descriptor file that works (note the dset entry shouldn't wrap):

dset http://acdisc.sci.gsfc.nasa.gov/opendap/Aura_OMI_Level3/OMAEROe.003/2006/OMI-Aura_L3-OMAEROe_2006m0101_v003-2008m0815t092235.he5
dtype netcdf
undef -88888 missing_value
unpack scale_factor add_offset
xdef 1440 linear -179.875 0.25
ydef  720 linear -89.875 0.25
zdef 1 linear 0 1
tdef 1 linear 00Z01JAN0001 1mn
edef 5 names  0 1 2 3 4
vars 13
AbsorbingAerosolOpticalThicknessMW=>var1          0  e,y,x  var1
AerosolModelMW=>var2                              0  y,x    var2
AerosolOpticalThicknessMW=>var3                   0  e,y,x  var3
AerosolOpticalThicknessPassedThresholdMean=>var4  0  y,x    var4
AerosolOpticalThicknessPassedThresholdStd=>var5   0  y,x    var5
SingleScatteringAlbedoMW=>var6                    0  e,y,x  var6
SingleScatteringAlbedoPassedThresholdMean=>var7   0  y,x    var7
SingleScatteringAlbedoPassedThresholdStd=>var8    0  y,x    var8
SolarZenithAngle=>var9      0  y,x  solarzenithangle
TerrainReflectivity=>var10  0  y,x  TerrainReflectivity
UVAerosolIndex=>var11       0  y,x  UVAerosolIndex
VISAerosolIndex=>var12      0  y,x  VISAerosolIndex
ViewingZenithAngle=>var13   0  y,x  ViewingZenithAngle
endvars


Jennifer

On Jan 8, 2009, at 10:20 AM, John Huddleston wrote:

> Hi
>
> I am accessing http://acdisc.sci.gsfc.nasa.gov/opendap/Aura_OMI_Level3/OMAEROe.003/2006/OMI-Aura_L3-OMAEROe_2006m0101_v003-2008m0815t092235.he5 
>  with the latest v2.0a4 gradsdap.
>
> It can open the file and the ?q file? works fine; however, when I  
> display some of the aerosol optical variables, there is an error:
>
> gancgrid error: nc_get_vara_double failed; NetCDF: Index exceeds  
> dimension bound
>
> Data Request Error:  Error for variable 'aerosolopticalt'
>   Error ocurred at column 1
> DISPLAY error:  Invalid expression
>   Expression = aerosolopticalt
>
> Is this a problem with the data or bounds limitation in the gradsdap  
> application?
>
> Here is the information on the gradsdap startup
>
> Starting X server under c:\GrADS20/Xming
> Starting gradsdap under c:\GrADS20\win32 ...
>
> Grid Analysis and Display System (GrADS) Version 2.0.a4.oga.1
> Copyright (c) 1988-2008 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> GrADS comes with ABSOLUTELY NO WARRANTY
> See file COPYRIGHT for more information
>
> Config: v2.0.a3.oga.1 little-endian readline printim grib2 netcdf  
> hdf4-sds opendap-grids,stn Athena
>
> John Huddleston
>
>
>
>
>
>
>

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/45c7145b/attachment.html 

From timlin at ATMOS.UIUC.EDU  Thu Jan  8 14:17:13 2009
From: timlin at ATMOS.UIUC.EDU (Michael Timlin)
Date: Thu, 8 Jan 2009 13:17:13 -0600
Subject: scripting question about 'define' and 'set defval'
In-Reply-To: <BAY117-W17E98A39F096C528E896A7F9DC0@phx.gbl>
Message-ID: <THU.8.JAN.2009.131713.0600.>

An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090108/662dbea4/attachment.html 

From meteo at SPORT123.IT  Fri Jan  9 02:43:26 2009
From: meteo at SPORT123.IT (Valerio Guerra)
Date: Fri, 9 Jan 2009 08:43:26 +0100
Subject: listserver rejecting messages
Message-ID: <FRI.9.JAN.2009.084326.0100.>

Pablo Romero wrote:
> is anyone else experiencing this problem?

Well, in the last few day I received 3 message starting like this one:

> From: L-Soft list server at CINECA <LISTSERV at LIST.CINECA.IT>
> Subject: Rejected posting to GRADSUSR at LIST.CINECA.IT
>
> Your message is being returned to you unprocessed because it appears to have already  been  distributed to  the  GRADSUSR  list...

Anyway my posting to the list seemed successful.

Valerio Guerra.


From matt.alonso at GMAIL.COM  Fri Jan  9 08:17:26 2009
From: matt.alonso at GMAIL.COM (Matt Alonso)
Date: Fri, 9 Jan 2009 08:17:26 -0500
Subject: listserver rejecting messages
In-Reply-To: <BAY117-W124B4892A3EC2A21BC8A9AF9DF0@phx.gbl>
Message-ID: <FRI.9.JAN.2009.081726.0500.>

I had the same issue see the email I received below:

"Your message  is being returned to  you unprocessed because it  appears to
have
already  been  distributed to  the  GRADSUSR  list.  That  is, a  message
 with
identical text  (but possibly with different  mail headers) has been  posted
to
the list recently, either by you or by  someone else. If you have a good
reason
to resend this message to the list (for instance because you have been
notified
of a hardware failure with loss of  data), please alter the text of the
message
in some way and  resend it to the list. Note that  altering the "Subject:"
line
or adding blank  lines at the top  or bottom of the message  is not
sufficient;
you should  instead add a  sentence or  two at the  top explaining why  you
are
resending the  message, so that the  other subscribers understand why  they
are
getting two copies of the same message."



On Wed, Jan 7, 2009 at 6:28 PM, Pablo Romero <romero619 at hotmail.com> wrote:

>  is anyone else experiencing this problem?
> ------------------------------
> Windows LiveTM: Keep your life in sync. See how it works.<http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t1_allup_howitworks_012009>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090109/719739d4/attachment.html 

From arlindo.dasilva at GMAIL.COM  Fri Jan  9 09:47:39 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Fri, 9 Jan 2009 09:47:39 -0500
Subject: listserver rejecting messages
In-Reply-To: <6cc1fd310901090517y5805d06fr8a204006965a41c0@mail.gmail.com>
Message-ID: <FRI.9.JAN.2009.094739.0500.>

On Fri, Jan 9, 2009 at 8:17 AM, Matt Alonso <matt.alonso at gmail.com> wrote:

> I had the same issue see the email I received below:
>

Mii too.


>
>
> "Your message  is being returned to  you unprocessed because it  appears to
> have
> already  been  distributed to  the  GRADSUSR  list.  That  is, a  message
>  with
> identical text  (but possibly with different  mail headers) has been
>  posted to
> the list recently, either by you or by  someone else. If you have a good
> reason
> to resend this message to the list (for instance because you have been
> notified
> of a hardware failure with loss of  data), please alter the text of the
> message
> in some way and  resend it to the list. Note that  altering the "Subject:"
> line
> or adding blank  lines at the top  or bottom of the message  is not
> sufficient;
> you should  instead add a  sentence or  two at the  top explaining why  you
> are
> resending the  message, so that the  other subscribers understand why  they
> are
> getting two copies of the same message."
>
>
>
>
> On Wed, Jan 7, 2009 at 6:28 PM, Pablo Romero <romero619 at hotmail.com>wrote:
>
>>  is anyone else experiencing this problem?
>> ------------------------------
>> Windows LiveTM: Keep your life in sync. See how it works.<http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t1_allup_howitworks_012009>
>>
>
>


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090109/72d741e6/attachment.html 

From graeme at SURFLINE.COM  Fri Jan  9 19:37:38 2009
From: graeme at SURFLINE.COM (Graeme Rae)
Date: Fri, 9 Jan 2009 16:37:38 -0800
Subject: Killing the splash message
In-Reply-To: <20081201141738.90EED20A2C@mx2.cineca.it>
Message-ID: <FRI.9.JAN.2009.163738.0800.>

An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090109/d25895e5/attachment.html 

From arlindo.dasilva at GMAIL.COM  Sat Jan 10 08:34:31 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Sat, 10 Jan 2009 08:34:31 -0500
Subject: Killing the splash message
In-Reply-To: <C58D2DD2.146271%graeme@surfline.com>
Message-ID: <SAT.10.JAN.2009.083431.0500.>

On Fri, Jan 9, 2009 at 7:37 PM, Graeme Rae <graeme at surfline.com> wrote:

>
>
> Is there any way to turn off this message:
>
>
> > Grid Analysis and Display System (GrADS) Version 1.8SL11
> > Copyright (c) 1988-2001 by Brian Doty
> > Center for Ocean-Land-Atmosphere Studies
> > Institute for Global Environment and Society
> > All Rights Reserved
> >
> > Config: v1.8SL11 32-bit little-endian readline printim
> >
> > Issue 'q config' command for more information.
> >
> > GX Package Initialization: Size = 11 8.5
> > Running in Batch mode
> > No hardcopy metafile open
> > GX package terminated
>
>
>
> I'm outputting to a log file and 95% of the log is the above, the
> occasional
> error message or useful information is a little hard to see!
>

On Linux/Unix (and on Windows under cygwin) just pipe it through tail, e.g.,

grads -blc 'run somescript.gs' | tail -c +16

This will eliminate the first 16 lines of output.

   Arlindo


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090110/ebd1b9dd/attachment.html 

From j.burn at GEOS.COM  Sun Jan 11 05:59:36 2009
From: j.burn at GEOS.COM (James Burn)
Date: Sun, 11 Jan 2009 11:59:36 +0100
Subject: Access to Thredds and Opendap data
Message-ID: <SUN.11.JAN.2009.115936.0100.>

Hi

I'm wondering if it is possible to use Grads 2.0 to access data via an
Opendap server such as
http://nomads.ncdc.noaa.gov/thredds/catalog/gfs4/catalog.html

I've not seen any mention of Thredds here on the Grads list so am guessing
its not the done thing. I've Windows XP here and Opensuse Linux at home.

When I try to sdfopen one of the Grib2 files eg
http://nomads.ncdc.noaa.gov/thredds/dodsC/gfs4/200901/20090108/gfs_4_2009010
8_1800_171.grb2

I get gaopnc error and "Can't ingest SDF Metadata".

I've also tried making the appropriate .inv and .idx files locally and
opening a local .ctl file which seems more promising, but seems to fall
over when accessing the binary file via http.

Is what I am trying to do possible? Or is data access limited to dods and
local servers?  My feeling is that I'm missing a trick or two.

This is my first post, although I've been hacking together scripts for a
year now and would like to thank all those who contribute to this list for
the tips I've picked up so far!

Thanks
James Burn


From madan at EES.HOKUDAI.AC.JP  Mon Jan 12 02:54:15 2009
From: madan at EES.HOKUDAI.AC.JP (madan sigdel)
Date: Mon, 12 Jan 2009 16:54:15 +0900
Subject: unit for vertical integartion
In-Reply-To: <20090111110034.3910D1FD32@mx2.cineca.it>
Message-ID: <MON.12.JAN.2009.165415.0900.>

Dear Users

Will any expert can tell me the standard unit for moisture
convergence/divergence calculated using grads.

 First I have calculated vertical integrated moisture flux as kg/m s. I
think this unit is right as sp.humidity is kg/kg, pressure is hPa and wind
in m/s. Now using same moisture flux I want to calculate the
convergence/divergence. Waiting for your answer.

Thank you,
madan


From jma at COLA.IGES.ORG  Mon Jan 12 07:32:34 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Mon, 12 Jan 2009 07:32:34 -0500
Subject: Access to Thredds and Opendap data
In-Reply-To: <20090111110034.3910D1FD32@mx2.cineca.it>
Message-ID: <MON.12.JAN.2009.073234.0500.>

You should be using the 'gradsdap' executable to access OPeNDAP
servers. I am able to open one of these data sets with sdfopen. Beware
that there are many more coordinate variables than the usual set, and
to see all the variables you will need to use xdfopen or open with a
descriptor file to map the following to X, Y, Z, T, and E:

netcdf gfs_4_20090109_1800_171 { // format variant: 64bit
dimensions:
         time1 = 1 ;
         pressure1 = 1 ;
         lat = 361 ;
         lon = 720 ;
         pressure = 26 ;
         time = 1 ;
         pressure3 = 21 ;
         pressure_difference_layer1 = 1 ;
         pressure4 = 2 ;
         potential_vorticity_surface = 2 ;
         height_above_ground1 = 1 ;
         pressure2 = 6 ;
         pressure_difference_layer = 1 ;
         depth_below_surface_layer = 4 ;
         altitude_above_msl = 3 ;
         height_above_ground = 1 ;
         bounds_dim = 2 ;


On Jan 11, 2009, at 5:59 AM, James Burn wrote:

> Hi
>
> I'm wondering if it is possible to use Grads 2.0 to access data via an
> Opendap server such as
> http://nomads.ncdc.noaa.gov/thredds/catalog/gfs4/catalog.html
>
> I've not seen any mention of Thredds here on the Grads list so am
> guessing
> its not the done thing. I've Windows XP here and Opensuse Linux at
> home.
>
> When I try to sdfopen one of the Grib2 files eg
> http://nomads.ncdc.noaa.gov/thredds/dodsC/gfs4/200901/20090108/gfs_4_2009010
> 8_1800_171.grb2
>
> I get gaopnc error and "Can't ingest SDF Metadata".
>
> I've also tried making the appropriate .inv and .idx files locally and
> opening a local .ctl file which seems more promising, but seems to
> fall
> over when accessing the binary file via http.
>
> Is what I am trying to do possible? Or is data access limited to
> dods and
> local servers?  My feeling is that I'm missing a trick or two.
>
> This is my first post, although I've been hacking together scripts
> for a
> year now and would like to thank all those who contribute to this
> list for
> the tips I've picked up so far!
>
> Thanks
> James Burn

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090112/fec82ec5/attachment.html 

From hrc at EAS.GATECH.EDU  Mon Jan 12 12:39:27 2009
From: hrc at EAS.GATECH.EDU (Hai-Ru Chang)
Date: Mon, 12 Jan 2009 12:39:27 -0500
Subject: gd erroe in grads-2.0a4
Message-ID: <MON.12.JAN.2009.123927.0500.>

Hello,
            When I tried to compile grads-2.0a4, I met the following:

gxmeta.o: In function `gxhpng':
gxmeta.c:(.text+0x1d40): undefined reference to `gdImageCreateFromPng'
gxmeta.c:(.text+0x296c): undefined reference to `gdImageCreateFromPng'
gxmeta.c:(.text+0x2aa0): undefined reference to `gdImageCreateFromPng'
gxmeta.c:(.text+0x2ba1): undefined reference to `gdImagePng'


Cam anyone tell me how to solve it?  Thanks!





Dr. Hai-Ru Chang
School of Earth and Atmospheric Sciences
Georgia Institute of Technology
Atlanta, GA 30332-0340


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090112/0c7698fb/attachment.html 

From arlindo.dasilva at GMAIL.COM  Mon Jan 12 14:47:40 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 12 Jan 2009 14:47:40 -0500
Subject: gd erroe in grads-2.0a4
In-Reply-To: <FB553FFF-C991-4DC8-9E5F-475EFD0CD317@eas.gatech.edu>
Message-ID: <MON.12.JAN.2009.144740.0500.>

On Mon, Jan 12, 2009 at 12:39 PM, Hai-Ru Chang <hrc at eas.gatech.edu> wrote:

> Hello,           When I tried to compile grads-2.0a4, I met the following:
>
> gxmeta.o: In function `gxhpng':
> gxmeta.c:(.text+0x1d40): undefined reference to `gdImageCreateFromPng'
> gxmeta.c:(.text+0x296c): undefined reference to `gdImageCreateFromPng'
> gxmeta.c:(.text+0x2aa0): undefined reference to `gdImageCreateFromPng'
> gxmeta.c:(.text+0x2ba1): undefined reference to `gdImagePng'
>
>
> Cam anyone tell me how to solve it?  Thanks!
>

Which gd version do you have? Look at this for the required version:

    http://grads.iges.org/grads/gadoc/supplibs.html

          Arlindo


>
>
>
>
>
> Dr. Hai-Ru Chang
>
> School of Earth and Atmospheric Sciences
>
> Georgia Institute of Technology
>
> Atlanta, GA 30332-0340
>
>
>
>


--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090112/f9ae61da/attachment.html 

From maldmm5 at GMAIL.COM  Wed Jan 14 16:03:49 2009
From: maldmm5 at GMAIL.COM (Mald MM5)
Date: Thu, 15 Jan 2009 10:03:49 +1300
Subject: How to remove the boxes in the contour labels????
Message-ID: <THU.15.JAN.2009.100349.1300.>

Hi there anybody,
I have a plot where contour plot is superimposed on a shaded plot. There is
a box coming with the contour label. However, I want the box to be removed
and keep the label only. How do I achieve this?...

Thanks
Mald
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090115/1eec8315/attachment.html 

From robock at ENVSCI.RUTGERS.EDU  Wed Jan 14 16:14:18 2009
From: robock at ENVSCI.RUTGERS.EDU (Alan Robock)
Date: Wed, 14 Jan 2009 16:14:18 -0500
Subject: help with GDS
Message-ID: <WED.14.JAN.2009.161418.0500.>

Dear friends,

My student, Anthony DeAngelis, and I are trying to use GDS for the first
time and have been unsuccessful.  Could you please help us?  The error
messages are not particularly useful.  Thanks very much.

We installed GrADS and Java on our machine and tried to read cumulus
precip from a GSFC archive:

[toine379 at control5 ~]$ java -version
java version "1.4.2"
gcj (GCC) 3.4.4 20050721 (Red Hat 3.4.4-2)
Copyright (C) 2004 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is
NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE.

[toine379 at control5 ~]$ grads

Grid Analysis and Display System
(GrADS) Version 2.0.a3
Copyright (c) 1988-2008 by Brian Doty and the
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information

Config: v2.0.a3 little-endian readline printim grib2 netcdf hdf4-sds
Issue 'q config' command for more information.
Landscape mode? ('n' for portrait):
GX Package Initialization: Size = 11 8.5
ga->  sdfopen
http://nsipp.gsfc.nasa.gov:9090/dods/_expr_{/CO2c20c/ens01/CO2_prs}{tloop(aave(precu,global))}{0:360,-87:90,1000:1000,00Z01JAN2000:00Z01FEB2000}
Scanning self-describing file:
http://nsipp.gsfc.nasa.gov:9090/dods/_expr_{/CO2c20c/ens01/CO2_prs}{tloop(aave(precu,global))}{0:360,-87:90,1000:1000,00Z01JAN2000:00Z01FEB2000}
gaopnc error: nc_open failed to open file
http://nsipp.gsfc.nasa.gov:9090/dods/_expr_{/CO2c20c/ens01/CO2_prs}{tloop(aave(precu,global))}{0:360,-87:90,1000:1000,00Z01JAN2000:00Z01FEB2000}
read_metadata: gaopnc failed
gadsdf: Couldn't ingest SDF metadata.


Alan

Alan Robock, Professor II
   Director, Meteorology Undergraduate Program
   Associate Director, Center for Environmental Prediction
Department of Environmental Sciences        Phone: +1-732-932-9800 x6222
Rutgers University                                  Fax: +1-732-932-8644
14 College Farm Road                   E-mail: robock at envsci.rutgers.edu
New Brunswick, NJ 08901-8551  USA      http://envsci.rutgers.edu/~robock


From ela at COLA.IGES.ORG  Wed Jan 14 16:33:54 2009
From: ela at COLA.IGES.ORG (Eric Altshuler)
Date: Wed, 14 Jan 2009 16:33:54 -0500
Subject: help with GDS
In-Reply-To: <Pine.LNX.4.64.0901141607110.10115@envsci.rutgers.edu>
Message-ID: <WED.14.JAN.2009.163354.0500.>

Dear Dr. Robock,

To access data on a GDS, you need to use the 'gradsdap' executable, not 'grads'.

Best regards,

Eric L. Altshuler
Assistant Research Scientist
Center for Ocean-Land-Atmosphere Studies
4041 Powder Mill Road, Suite 302
Calverton, MD 20705-3106
USA

E-mail: ela at cola.iges.org
Phone: (301) 902-1257
Fax: (301) 595-9793

----- Original Message -----
From: "Alan Robock" <robock at ENVSCI.RUTGERS.EDU>
To: GRADSUSR at LIST.CINECA.IT
Sent: Wednesday, January 14, 2009 4:14:18 PM GMT -05:00 US/Canada Eastern
Subject: help with GDS

Dear friends,

My student, Anthony DeAngelis, and I are trying to use GDS for the first
time and have been unsuccessful.  Could you please help us?  The error
messages are not particularly useful.  Thanks very much.

We installed GrADS and Java on our machine and tried to read cumulus
precip from a GSFC archive:

[toine379 at control5 ~]$ java -version
java version "1.4.2"
gcj (GCC) 3.4.4 20050721 (Red Hat 3.4.4-2)
Copyright (C) 2004 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is
NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR
PURPOSE.

[toine379 at control5 ~]$ grads

Grid Analysis and Display System
(GrADS) Version 2.0.a3
Copyright (c) 1988-2008 by Brian Doty and the
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information

Config: v2.0.a3 little-endian readline printim grib2 netcdf hdf4-sds
Issue 'q config' command for more information.
Landscape mode? ('n' for portrait):
GX Package Initialization: Size = 11 8.5
ga->  sdfopen
http://nsipp.gsfc.nasa.gov:9090/dods/_expr_{/CO2c20c/ens01/CO2_prs}{tloop(aave(precu,global))}{0:360,-87:90,1000:1000,00Z01JAN2000:00Z01FEB2000}
Scanning self-describing file:
http://nsipp.gsfc.nasa.gov:9090/dods/_expr_{/CO2c20c/ens01/CO2_prs}{tloop(aave(precu,global))}{0:360,-87:90,1000:1000,00Z01JAN2000:00Z01FEB2000}
gaopnc error: nc_open failed to open file
http://nsipp.gsfc.nasa.gov:9090/dods/_expr_{/CO2c20c/ens01/CO2_prs}{tloop(aave(precu,global))}{0:360,-87:90,1000:1000,00Z01JAN2000:00Z01FEB2000}
read_metadata: gaopnc failed
gadsdf: Couldn't ingest SDF metadata.


Alan

Alan Robock, Professor II
   Director, Meteorology Undergraduate Program
   Associate Director, Center for Environmental Prediction
Department of Environmental Sciences        Phone: +1-732-932-9800 x6222
Rutgers University                                  Fax: +1-732-932-8644
14 College Farm Road                   E-mail: robock at envsci.rutgers.edu
New Brunswick, NJ 08901-8551  USA      http://envsci.rutgers.edu/~robock


From jma at COLA.IGES.ORG  Thu Jan 15 08:57:03 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Thu, 15 Jan 2009 08:57:03 -0500
Subject: New Documentation for the Ensemble Dimension
Message-ID: <THU.15.JAN.2009.085703.0500.>

Hi, Everyone --
I have added a new section to the User's Guide that explains how the
ensemble dimension has been implemented in GrADS 2.0. Additional
topics include how to create a descriptor file for an ensemble data
set, how to organize your data using templating, special
considerations for certain data formats, and some examples (with
illustrations!) using real data.

http://iges.org/grads/gadoc/ensembles.html

Real time NCEP GFS Ensemble forecasts are behind our GDS at http://monsoondata.org:9090/dods/gfsens/
Why not take a moment to fire up 'gradsdap', use 'sdfopen' to access
the latest forecast, and start playing with ensemble data. It's cool!

--Jennifer


From hmjbarbosa at GMAIL.COM  Thu Jan 15 10:17:28 2009
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Thu, 15 Jan 2009 13:17:28 -0200
Subject: Multiple Parameters on one Plot
In-Reply-To: <D688BDE22B02984496332260B9E36A8F0A7BE265@ATT1-EX1.am.int.amec.com>
Message-ID: <THU.15.JAN.2009.131728.0200.>

Dear Brian,

What you are searching for is called pattern fill, but that is not
available in grads.
The only option available for filling a plot is color shading. Not dots, lines,
stripes, etc...

I have done that however, but using gimp.


From hmjbarbosa at GMAIL.COM  Thu Jan 15 13:21:07 2009
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Thu, 15 Jan 2009 16:21:07 -0200
Subject: turn messages off
Message-ID: <THU.15.JAN.2009.162107.0200.>

Hi All,

Simple question, is it possible to turn off grads messages such as the
ones shown below?

ga->reinit
No hardcopy metafile open
All files closed; all defined objects released
All GrADS attributes have been reinitialized

Or like:

ga->sfopen some_double_precision_file
SDF:  The double precision values of the variable lon_bnds
will be automatically converted to single precision.
SDF:  The double precision values of the variable lat_bnds
will be automatically converted to single precision.

I've a script which does a lot of 'reinit's and 'sdfopen's
and I would like to have no messages shown after these
commands.

Thanks,
Henrique


From aditsuda at YAHOO.COM  Fri Jan 16 04:02:54 2009
From: aditsuda at YAHOO.COM (aditsuda jamroensan)
Date: Fri, 16 Jan 2009 01:02:54 -0800
Subject: plot pollution concentration above ground levels
Message-ID: <FRI.16.JAN.2009.010254.0800.>

Dear Grads users

     My wrfchem model give pollution concentrations
above sea levels. Is there any way I can plot the
concentrations above ground levels?

Thanks,
Ada


From kieucq at YAHOO.CO.UK  Fri Jan 16 11:48:34 2009
From: kieucq at YAHOO.CO.UK (Lemmy)
Date: Fri, 16 Jan 2009 16:48:34 +0000
Subject: New Documentation for the Ensemble Dimension
In-Reply-To: <4AF851A9-CC2E-4E05-A199-0D7C74E36B8D@cola.iges.org>
Message-ID: <FRI.16.JAN.2009.164834.0000.>

Hi Jennifer,

I have gone through the instructions of using ensemble dimension but somehow my data files are not opened properly.

Here are the ensemble of 11 data files I have:

mmm00.dat
mmm01.dat
.....
mmm11.dat

Here is my ctl file:

DSET ^mmm%e2.dat
OPTIONS template
TITLE etkf
UNDEF -99999
XDEF  51 linear 0  2.5
YDEF  21 linear 0  2.5
ZDEF  1   linear 0 1
TDEF  500  linear 00z01apr99 3hr
EDEF  11 names 1 2 3 3 4 5 6 7 8 9 10 11
vars 1
s 1 99 height
endvars

and here is the result:

Scanning description file:  baro.ctl
Data file mmm%e2.dat is open as file 1

I would expect it to open mmm01.dat,...., but it plainly opened just mmm%e2.dat. What surprised me further is that although the file mmm%e2.dat does not exist, the X-window is still popped up as if it works. I appreciate your any help very much.

Lemmy






--- On Thu, 15/1/09, Jennifer Adams <jma at COLA.IGES.ORG> wrote:
From: Jennifer Adams <jma at COLA.IGES.ORG>
Subject: New Documentation for the Ensemble Dimension
To: GRADSUSR at LIST.CINECA.IT
Date: Thursday, 15 January, 2009, 7:57 AM

Hi, Everyone --
I have added a new section to the User's Guide that explains how the
ensemble dimension has been implemented in GrADS 2.0. Additional
topics include how to create a descriptor file for an ensemble data
set, how to organize your data using templating, special
considerations for certain data formats, and some examples (with
illustrations!) using real data.

http://iges.org/grads/gadoc/ensembles.html

Real time NCEP GFS Ensemble forecasts are behind our GDS at
http://monsoondata.org:9090/dods/gfsens/
Why not take a moment to fire up 'gradsdap', use 'sdfopen' to
access
the latest forecast, and start playing with ensemble data. It's cool!

--Jennifer


From kieucq at YAHOO.CO.UK  Fri Jan 16 11:53:33 2009
From: kieucq at YAHOO.CO.UK (Lemmy)
Date: Fri, 16 Jan 2009 16:53:33 +0000
Subject: the Ensemble Dimension
In-Reply-To: <4AF851A9-CC2E-4E05-A199-0D7C74E36B8D@cola.iges.org>
Message-ID: <FRI.16.JAN.2009.165333.0000.>

I have gone through the instructions of using ensemble dimension but somehow my data files are not opened properly.

Here are the ensemble of 11 data files I have:

mmm00.dat
mmm01.dat
.....
mmm11.dat

Here is my ctl file:

DSET ^mmm%e2.dat
OPTIONS template
TITLE etkf
UNDEF -99999
XDEF  51 linear 0  2.5
YDEF  21 linear 0  2.5
ZDEF  1   linear 0 1
TDEF  500  linear 00z01apr99 3hr
EDEF  11 names 1 2 3 3 4 5 6 7 8 9 10 11
vars 1
s 1 99 height
endvars

and here is the result:

Scanning description file:  baro.ctl
Data file mmm%e2.dat is open as file 1

I would expect it to open mmm01.dat,...., but it plainly opened just mmm%e2.dat. What surprised me further is that although the file mmm%e2.dat does not exist, the X-window is still popped up as if it works. I appreciate your any help very much.

Lemmy


From jma at COLA.IGES.ORG  Fri Jan 16 12:01:49 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Fri, 16 Jan 2009 12:01:49 -0500
Subject: New Documentation for the Ensemble Dimension
In-Reply-To: <817831.68897.qm@web25218.mail.ukl.yahoo.com>
Message-ID: <FRI.16.JAN.2009.120149.0500.>

The template substitution string should be just %e, and your edef
should have the 2-digit names:
   dset ^mmm%e.dat
   edef 11 names 01 02 03 04 05 06 07 08 09 10 11
There is no formatting of numbers with %e  -- it is a straight
substitution of the ensemble name given in EDEF.

Also, since your data do not vary in Z, then the variable declaration
should be:
s 0 99 height

--Jennifer

On Jan 16, 2009, at 11:48 AM, Lemmy wrote:

> Hi Jennifer,
>
> I have gone through the instructions of using ensemble dimension but
> somehow my data files are not opened properly.
>
> Here are the ensemble of 11 data files I have:
>
> mmm00.dat
> mmm01.dat
> .....
> mmm11.dat
>
> Here is my ctl file:
>
> DSET ^mmm%e2.dat
> OPTIONS template
> TITLE etkf
> UNDEF -99999
> XDEF  51 linear 0  2.5
> YDEF  21 linear 0  2.5
> ZDEF  1   linear 0 1
> TDEF  500  linear 00z01apr99 3hr
> EDEF  11 names 1 2 3 3 4 5 6 7 8 9 10 11
> vars 1
> s 1 99 height
> endvars
>
> and here is the result:
>
> Scanning description file:  baro.ctl
> Data file mmm%e2.dat is open as file 1
>
> I would expect it to open mmm01.dat,...., but it plainly opened just
> mmm%e2.dat. What surprised me further is that although the file mmm
> %e2.dat does not exist, the X-window is still popped up as if it
> works. I appreciate your any help very much.
>
> Lemmy
>
>
>
>
>
>
> --- On Thu, 15/1/09, Jennifer Adams <jma at COLA.IGES.ORG> wrote:
> From: Jennifer Adams <jma at COLA.IGES.ORG>
> Subject: New Documentation for the Ensemble Dimension
> To: GRADSUSR at LIST.CINECA.IT
> Date: Thursday, 15 January, 2009, 7:57 AM
>
> Hi, Everyone --
> I have added a new section to the User's Guide that explains how the
> ensemble dimension has been implemented in GrADS 2.0. Additional
> topics include how to create a descriptor file for an ensemble data
> set, how to organize your data using templating, special
> considerations for certain data formats, and some examples (with
> illustrations!) using real data.
>
> http://iges.org/grads/gadoc/ensembles.html
>
> Real time NCEP GFS Ensemble forecasts are behind our GDS at
> http://monsoondata.org:9090/dods/gfsens/
> Why not take a moment to fire up 'gradsdap', use 'sdfopen' to
> access
> the latest forecast, and start playing with ensemble data. It's cool!
>
> --Jennifer

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/1af757c0/attachment.html 

From arlindo.dasilva at GMAIL.COM  Fri Jan 16 12:25:41 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Fri, 16 Jan 2009 12:25:41 -0500
Subject: How to remove the boxes in the contour labels????
In-Reply-To: <54c90bba0901141303o265f5947lc1ce5d139fd71df8@mail.gmail.com>
Message-ID: <FRI.16.JAN.2009.122541.0500.>

On Wed, Jan 14, 2009 at 4:03 PM, Mald MM5 <maldmm5 at gmail.com> wrote:

> Hi there anybody,
> I have a plot where contour plot is superimposed on a shaded plot. There is
> a box coming with the contour label. However, I want the box to be removed
> and keep the label only. How do I achieve this?...
>

Short of changing the source code, I don't think there is another way of
accomplishing it.  Why would you want to do this?  When writing gxyat I had
a bug at one point that caused those rectangles containing the labels not to
be filled. I can tell you that Brian had a good reason for doing this: it
does not look pretty, and worse, it can make the labels very difficult to
read when it intersects a continent boundary.

  Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/cc8bd728/attachment.html 

From romero619 at HOTMAIL.COM  Fri Jan 16 12:35:45 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Fri, 16 Jan 2009 09:35:45 -0800
Subject: Need Help developing WGS84 ("Web") Mercator Projection for GrADS
Message-ID: <FRI.16.JAN.2009.093545.0800.>


Hello,
Im currently reviewing, testing, & attempting to change the grads2.0 source code in order to add a Mercator projection.
For those who are familiar in GrADS code/development, I could please use some assistance.
I've narrowed the projection routines to the gwxmap.c source file.
 
I need help with:
2) how to create the "scaling" functions to correctly adjust the x,y values for the grads (internal) grid
1)determining in how many different places the scaling function needs to be applied, 
  (*since it seems there are a few different routines throughout all of the source code files where the
   scaling is performed)
3)question: once the projection routine is implemented, it *should* take care of correctly scaling&drawing everything, including maps, grids, shapefiles (using shapefile extension), symbols, etc., is this correct?
 
The functions for transforming from lat/lon to x,y are fairly simple: 
http://en.wikipedia.org/wiki/Mercator_projection 
x=lon; 
longitude remains unchanged
 
y=ln(tan((pi/4)+(lat/2))); 
latitude is variable
 
so, I believe the key is just setting up the transformation functions and correctly referring to them throughout the source code.
 
Most other graphics packages already include a Mercator projection, and it gets used alot for web mapipng/tiling, such as google earth & MS Virtual Earth.
Grads needs this...
 
Those interested, please email me.
Thanks,,
P.Romero
 
 
 
_________________________________________________________________
Windows Live? Hotmail?: Chat. Store. Share. Do more with mail. 
http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t1_hm_justgotbetter_explore_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/2328d812/attachment.html 

From arlindo.dasilva at GMAIL.COM  Fri Jan 16 12:42:43 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Fri, 16 Jan 2009 12:42:43 -0500
Subject: Multiple Parameters on one Plot
In-Reply-To: <e71668fa0901150717k2719d12aw6b841cf89678bbc@mail.gmail.com>
Message-ID: <FRI.16.JAN.2009.124243.0500.>

On Thu, Jan 15, 2009 at 10:17 AM, Henrique Barbosa <hmjbarbosa at gmail.com>wrote:

> Dear Brian,
>
> What you are searching for is called pattern fill, but that is not
> available in grads.
> The only option available for filling a plot is color shading. Not dots,
> lines,
> stripes, etc...
>
> I have done that however, but using gimp.
>
> From grads I exported two images, one with the color shaded variable
> plotted,
> and the other with a second variable plotted using a solid contour.
> From Gimp I imported each image at a different layer, keeping the contoured
> one
> on top. Erased its background to let the bottom one come visible.
> Then filled it with a pattern, and then had grads merge both layers
> and save the composed image.
>
> Quite time consuming, but it worked for me.
>
> Another option could be using PyGrads to rely on python
> graphical capabilities.... Never done it, but maybe it has
> pattern filling.
>

And yet one more option is to use gxyat to create SVG instead of PNG and pop
that into a vector graphics application such as Adobe Illustrator or
Inkscape. The advantage of vector graphics is that you do not loose any
resolution.

If you know C, and are willing to learn the cairo library, you could extend
gxyat and allow the user to substitude specific colors with specified fill
patterns. Drop a note at opengrads-devel at lists.sf.net if are interested in
working on gxyat.

       Arlindo




> Good luck,
> Henrique
>
>
> On Mon, Jan 5, 2009 at 4:06 PM, Walsh, Brian J <brian.walsh at amec.com>
> wrote:
> > I am trying to do a plot of 3 variables; one is shaded (dew point
> > depression), the second is dashed contours (temperature) and the third is
> > area of freezing rain (solid contour).  Is there any function in grads
> where
> > I can fill the solid contour with something such as dots, stripes, lines,
> > etc. or is this not possible since I already have one variable being
> shaded?
> >
> > Thanks,
> >
> > The information contained in this e-mail is intended only for the
> individual
> > or entity to whom it is addressed.
> > Its contents (including any attachments) may contain confidential and/or
> > privileged information.
> > If you are not an intended recipient you must not use, disclose,
> > disseminate, copy or print its contents.
> > If you receive this e-mail in error, please notify the sender by reply
> > e-mail and delete and destroy the message.
> >
>



--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/6defb028/attachment.html 

From aditsuda at YAHOO.COM  Fri Jan 16 12:55:24 2009
From: aditsuda at YAHOO.COM (Aditsuda Jamroensan)
Date: Fri, 16 Jan 2009 18:55:24 +0100
Subject: plot concentration above ground levels
Message-ID: <FRI.16.JAN.2009.185524.0100.>

Dear Grads users


     My wrfchem model give pollution concentrations
above sea levels. Is there any way I can plot the
concentrations above ground levels?


Thanks,
Ada


From jma at COLA.IGES.ORG  Fri Jan 16 13:25:58 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Fri, 16 Jan 2009 13:25:58 -0500
Subject: How to remove the boxes in the contour labels????
In-Reply-To: <54c90bba0901141303o265f5947lc1ce5d139fd71df8@mail.gmail.com>
Message-ID: <FRI.16.JAN.2009.132558.0500.>

You can't do what you are asking for with the current version of
GrADS, but Brian is working on an overhaul of the contour labeling for
a future release (I don't know when it will be available). The code is
nearly ready but the user interface for controlling the options hasn't
been finalized. The filled rectangles behind the contour labels will
be replaced by a gap in the contour line.

In the meanwhile, you can use 'set background' to control the color of
the filled rectangle behind the contour labels. Although the boxes
often straddle both sides of a contour line, you may be able to
diminish their visual impact by drawing each contour separately,
specifying the background colors that you used in your shaded plot.

Jennifer



On Jan 14, 2009, at 4:03 PM, Mald MM5 wrote:

> Hi there anybody,
>
> I have a plot where contour plot is superimposed on a shaded plot.
> There is a box coming with the contour label. However, I want the
> box to be removed and keep the label only. How do I achieve this?...
>
> Thanks
> Mald
>
>

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/5ad9c97f/attachment.html 

From taranaphy at YAHOO.COM  Fri Jan 16 13:39:46 2009
From: taranaphy at YAHOO.COM (Tarana Mahzabin)
Date: Fri, 16 Jan 2009 19:39:46 +0100
Subject: about cyclone track
Message-ID: <FRI.16.JAN.2009.193946.0100.>

Dear users,

I have some Lat lon value of Cyclone Eye position.
Is there any way to plot this data through GrADS...
so that I can get the track?
Any body can help me,Please.

Tarana


From jma at COLA.IGES.ORG  Fri Jan 16 13:49:52 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Fri, 16 Jan 2009 13:49:52 -0500
Subject: Multiple Parameters on one Plot
In-Reply-To: <D688BDE22B02984496332260B9E36A8F0A7BE265@ATT1-EX1.am.int.amec.com>
Message-ID: <FRI.16.JAN.2009.134952.0500.>

You can create a stippling effect for your top layer (area of freezing
rain) like this:

set gxout grid
set grid off
set gridln off
set digsiz 0.02
d const(maskout(cfrzr,cfrzr-0.5),0)

This will draw little tiny 0's at the center of every grid point that
has a cfrzr value > 0.5. You may have to play with the digsize to get
it to look just right. And you can use 'set ccolor' to control the
color of the stippling.

Jennifer





On Jan 5, 2009, at 1:06 PM, Walsh, Brian J wrote:

> I am trying to do a plot of 3 variables; one is shaded (dew point
> depression), the second is dashed contours (temperature) and the
> third is area of freezing rain (solid contour).  Is there any
> function in grads where I can fill the solid contour with something
> such as dots, stripes, lines, etc. or is this not possible since I
> already have one variable being shaded?
>
> Thanks,
>
> The information contained in this e-mail is intended only for the
> individual or entity to whom it is addressed.
> Its contents (including any attachments) may contain confidential
> and/or privileged information.
> If you are not an intended recipient you must not use, disclose,
> disseminate, copy or print its contents.
> If you receive this e-mail in error, please notify the sender by
> reply e-mail and delete and destroy the message.

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/59941e16/attachment.html 

From Charles.Seman at NOAA.GOV  Fri Jan 16 14:43:52 2009
From: Charles.Seman at NOAA.GOV (Charles Seman)
Date: Fri, 16 Jan 2009 14:43:52 -0500
Subject: about cyclone track
In-Reply-To: <20090116184045.44EBC21139@mx2.cineca.it>
Message-ID: <FRI.16.JAN.2009.144352.0500.>

Tarana,

Please find attached two GrADS scripts:
plot_hurricane_center3.gs
  uses minimum sea-level or surface pressure to identify "center" (user
defines variable name)
  can specify more than one "model"
plot_cyclone_track-interactive.gs
  uses up to 3 different parameters to find "center": surface pressure,
relative vorticity, max wind (user defines an area for finding "center"
for each time level)
  specify one "model" (can specify separate input datasets for
vorticity, max wind)

I hope this helps.
If you decide to use either of the scripts, and have questions, please
let me know.

Thanks,
Chuck

Tarana Mahzabin wrote:
> Dear users,
>
> I have some Lat lon value of Cyclone Eye position.
> Is there any way to plot this data through GrADS...
> so that I can get the track?
> Any body can help me,Please.
>
> Tarana
>

--

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
 Charles Seman                                Charles.Seman at noaa.gov
 U.S. Department of Commerce / NOAA / OAR
 Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
 201 Forrestal Road                              fax: (609) 987-5063
 Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not necessarily
reflect any position of the Government or NOAA."

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: plot_hurricane_center3.gs
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/5c2af31d/attachment.pl 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: plot_cyclone_track-interactive.gs
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/5c2af31d/attachment-0001.pl 

From Charles.Seman at NOAA.GOV  Fri Jan 16 15:07:17 2009
From: Charles.Seman at NOAA.GOV (Charles Seman)
Date: Fri, 16 Jan 2009 15:07:17 -0500
Subject: How to remove the boxes in the contour labels????
In-Reply-To: <54c90bba0901141303o265f5947lc1ce5d139fd71df8@mail.gmail.com>
Message-ID: <FRI.16.JAN.2009.150717.0500.>

Mald,

I'm not sure if this might be of some help?  Please find attached a
GrADS script function clabel.gsf which is a script to plot contour
labels interactively. Also attached is a GrADS script fig1ab_sd.gs which
uses clabel.gsf...  If you find this useful, and have questions, please
let me know...
Thanks,
Chuck

some documentation:
...in fig1ab_sd.gs, input parameters for "clabel" (need 3 parameters:
csiz = contour label size, cformat = arithmetic format, and clist =
contour values to plot (here, 4 values in "sd_level"; sd_style defines
contour line style for each contour value in sd_level)):
*
*  Define the contour levels for the sample standard deviation...
*
sd_level = '0.2 0.5 0.8 1.1'
*sd_style = ' 5   3   2   1'
sd_style = ' 1   1   1   1'
*sd_label = 'off'
*sd_label = 'on'
sd_label = 'manual' ; csiz = 0.08 ; cformat = '%3.1f'

...further down in fig1ab_sd.gs, a loop to plot standard deviation
contour line and prompts user whether to label it using "clabel()" (here
contour lines and labels are plotted one-by-one), and after each
contour, prompts user for a replot to retry contour label; goes on to
next contour level if user does not type "yes":

  n=1 ; sd_clev = subwrd(sd_level,n) ; sd_cstyle = subwrd(sd_style,n)
  while ( sd_clev != '' )
    say 'n, sd_clev = 'n', 'sd_clev
    'set gxout contour'
    'set clevs 'sd_clev
    'set cstyle 'sd_cstyle
    'set ccolor 1'
    'set xlab off'
    'set ylab off'
    'set x 1'
    'set y 1 'ny
    'set lev 'zlevs
    'set t 1'
    'd sd'
    if( sd_label = 'manual' )
      say
      say 'Contour level = 'sd_clev
      say 'Do you want to plot a contour label?'
      say 'Type y(es) or Y(es), or press "ENTER" to continue...'
      pull plot
      if( plot='y' | plot='yes' | plot='Y' | plot='Yes' )
        clabel(csiz,cformat,sd_clev)
      endif
      say
      say 'Do you want to replot panel?'
      say 'Type y(es) or Y(es), or press "ENTER" to continue...'
      pull replot
      if( replot='y' | replot='yes' | replot='Y' | replot='Yes' )
        say
        say 'plot zonal and time averaged "'var_name'"...'
        say
        'set gxout shaded'
        'set clevs '_vcl.m1
        'set ccols '_vcc.m1
        'set xlab off'
        'set ylab off'
        'set x 1'
        'set y 1 'ny
        'set lev 'zlevs
        'set t 1'
        'd ztav'
        'set gxout stat'
        'd ztav'
        min_max = sublin(result,8)
        say '   'min_max
        if( plot_min_max = 'yes' )
          xmm = xl_pvar_parea + 0.1
          ymm = yt_pvar_parea - 0.1
          'set string 1 l 4'
          'set strsiz 0.08'
          'draw string 'xmm' 'ymm' Data: 'min_max' '
        endif
        if( plot_min_temp = 'yes' )
          say
          say '   plot level of minimum temperature...'
          say
          j=1
          while ( j <= ny )
            'set y 'j
*           say '      xloc, yloc = 'xloc.j', 'yloc.j
            'set line 1'
            'draw mark 3 'xloc.j' 'yloc.j' 0.03'
            j=j+1
          endwhile
        endif
        n=0
      endif
    endif
    n=n+1
    sd_clev = subwrd(sd_level,n) ; sd_cstyle = subwrd(sd_style,n)
  endwhile



Mald MM5 wrote:
> Hi there anybody,
>
> I have a plot where contour plot is superimposed on a shaded plot.
> There is a box coming with the contour label. However, I want the box
> to be removed and keep the label only. How do I achieve this?...
>
> Thanks
> Mald
>
>

--

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
 Charles Seman                                Charles.Seman at noaa.gov
 U.S. Department of Commerce / NOAA / OAR
 Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
 201 Forrestal Road                              fax: (609) 987-5063
 Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not necessarily
reflect any position of the Government or NOAA."

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: clabel.gsf
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/0ccad345/attachment.pl 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: fig1ab_sd.gs
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090116/0ccad345/attachment-0001.pl 

From furue at HAWAII.EDU  Fri Jan 16 19:33:51 2009
From: furue at HAWAII.EDU (Ryo Furue)
Date: Fri, 16 Jan 2009 14:33:51 -1000
Subject: How to remove the boxes in the contour labels????
In-Reply-To: <54c90bba0901141303o265f5947lc1ce5d139fd71df8@mail.gmail.com>
Message-ID: <FRI.16.JAN.2009.143351.1000.>

Hi Mald,

| I have a plot where contour plot is superimposed on a shaded
| plot. There is a box coming with the contour label. However,
| I want the box to be removed and keep the label only. How do I
| achieve this?...

I'm reasonably sure that this is not what you want :-)
but here we go:

 $ egrep -v '^c0 g0 [0-9]+ [0-9]+ [0-9]+ [0-9]+ B' plot.eps > out.eps

where "$" is the shell prompt and "plot.eps" is a PostScript
file generated by gxeps from a GrADS output.

If you do this, you can see that contour lines and labels overlap
(because the boxes are gone) and the result isn't pretty.

I suppose that a more proper solution is to wait for a future
version of GrADS (see Jennifer's message) where the contour labels
will be replaced by a gap in the contour line.

Cheers,
Ryo


From madan at EES.HOKUDAI.AC.JP  Sat Jan 17 02:13:22 2009
From: madan at EES.HOKUDAI.AC.JP (madan sigdel)
Date: Sat, 17 Jan 2009 16:13:22 +0900
Subject: unit for veritical integrated moisture convergence ( Reposting)
In-Reply-To: <4970E8F5.2080306@noaa.gov>
Message-ID: <SAT.17.JAN.2009.161322.0900.>

Dear Users

Please tell me the standard unit for moisture convergence/divergence
[-hdivg(qu, qv)]calculated using grads.

 First I have calculated vertical integrated moisture flux (qu, qv) and unit
as kg/m s. I think this unit is right as specific humidity is kg/kg,
pressure is hPa and wind in m/s. After this using same moisture flux I want
to calculate the convergence/divergence. In some papers the unit is
mentioned as mm/day. How can I get it.
Please suggest me.

Thank you,
madan


From jbjamali at GMAIL.COM  Sat Jan 17 05:36:41 2009
From: jbjamali at GMAIL.COM (javad bodaghjamali)
Date: Sat, 17 Jan 2009 14:06:41 +0330
Subject: monthly and seasonal rainfall analysis
Message-ID: <SAT.17.JAN.2009.140641.0330.>

Dear Grads users,

Hi, I have some problems in do analysis monthly and seasonal from daily
data.

Would you please infrom me any comment or help in reagrd of the following
matters;

1) Is it possible to change the time increment in contorol file of grads
which has been build for daily rain data?
    fro example from mm/day to mm/month . Is it needed to use some *.gs
files or use fortan program and combine to grads.

2) I want to make a time a series of total No. of stations which has been
used in each day for our analysis. In our contorol file the No. of stations
as an variable will be produced. Is there any function that could summing
all of  the no. of station  over total grids  , could SUMG do it or not?

3) I want to make a point to point correlation between two variables which
both variables are X, Y, and Time dependent and make a spatial distribution
for correlation cefficient for two variables . Is there any *.gs files to do
it?

Thanks in advance and looking forward to your reply.

Javad BodaghJamali
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090117/5898e9d0/attachment.html 

From j.burn at GEOS.COM  Sat Jan 17 13:37:40 2009
From: j.burn at GEOS.COM (Burn James)
Date: Sat, 17 Jan 2009 19:37:40 +0100
Subject: Access to Thredds and Opendap data
Message-ID: <SAT.17.JAN.2009.193740.0100.>

Thanks, I'm embarrassed to say gradsdap was my problem! The icon setup on
the desktop here opened "grads.exe" rather than gradsdap.  Not much of an
excuse and definitely enough for me to feel pretty daft.

So now I can see individual timesteps, I'll have a play with template
option in descriptor file to access the different timesteps.

Thanks


On Mon, 12 Jan 2009 07:32:34 -0500, Jennifer Adams <jma at COLA.IGES.ORG>
wrote:

>You should be using the 'gradsdap' executable to access OPeNDAP
>servers. I am able to open one of these data sets with sdfopen. Beware
>that there are many more coordinate variables than the usual set, and
>to see all the variables you will need to use xdfopen or open with a
>descriptor file to map the following to X, Y, Z, T, and E:
>
>netcdf gfs_4_20090109_1800_171 { // format variant: 64bit
>dimensions:
>         time1 = 1 ;
>         pressure1 = 1 ;
>         lat = 361 ;
>         lon = 720 ;
>         pressure = 26 ;
>         time = 1 ;
>         pressure3 = 21 ;
>         pressure_difference_layer1 = 1 ;
>         pressure4 = 2 ;
>         potential_vorticity_surface = 2 ;
>         height_above_ground1 = 1 ;
>         pressure2 = 6 ;
>         pressure_difference_layer = 1 ;
>         depth_below_surface_layer = 4 ;
>         altitude_above_msl = 3 ;
>         height_above_ground = 1 ;
>         bounds_dim = 2 ;
>
>
>On Jan 11, 2009, at 5:59 AM, James Burn wrote:
>
>> Hi
>>
>> I'm wondering if it is possible to use Grads 2.0 to access data via an
>> Opendap server such as
>> http://nomads.ncdc.noaa.gov/thredds/catalog/gfs4/catalog.html
>>
>> I've not seen any mention of Thredds here on the Grads list so am
>> guessing
>> its not the done thing. I've Windows XP here and Opensuse Linux at
>> home.
>>
>> When I try to sdfopen one of the Grib2 files eg
>>
http://nomads.ncdc.noaa.gov/thredds/dodsC/gfs4/200901/20090108/gfs_4_2009010
>> 8_1800_171.grb2
>>
>> I get gaopnc error and "Can't ingest SDF Metadata".
>>
>> I've also tried making the appropriate .inv and .idx files locally and
>> opening a local .ctl file which seems more promising, but seems to
>> fall
>> over when accessing the binary file via http.
>>
>> Is what I am trying to do possible? Or is data access limited to
>> dods and
>> local servers?  My feeling is that I'm missing a trick or two.
>>
>> This is my first post, although I've been hacking together scripts
>> for a
>> year now and would like to thank all those who contribute to this
>> list for
>> the tips I've picked up so far!
>>
>> Thanks
>> James Burn
>
>--
>Jennifer M. Adams
>IGES/COLA
>4041 Powder Mill Road, Suite 302
>Calverton, MD 20705
>jma at cola.iges.org
>
>
>
>


From appopson at YAHOO.FR  Mon Jan 19 07:43:59 2009
From: appopson at YAHOO.FR (louis Vonder)
Date: Mon, 19 Jan 2009 12:43:59 +0000
Subject: Vertically integrated moisture flux
Message-ID: <MON.19.JAN.2009.124359.0000.>

Hi grads users,

Please? I want to know if one of you? can have a script to compute 
"Vertically integrated moisture flux" using reanalysis data?

If is yes, please can you share your script with me?

Thanks



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090119/16c5c26a/attachment.html 

From appopson at YAHOO.FR  Mon Jan 19 07:57:18 2009
From: appopson at YAHOO.FR (louis Vonder)
Date: Mon, 19 Jan 2009 12:57:18 +0000
Subject: A script: Vertically integrated moisture flux
Message-ID: <MON.19.JAN.2009.125718.0000.>

Dear Grads users,

I wannot if anybody have a grads script to compute "Vertically integrated moisture flux"?? using reanalysis data?

In the positive case, can you share your script with me?

Thanks in advance.



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090119/86ab98c1/attachment.html 

From shimon at CYCLONE.TAU.AC.IL  Tue Jan 20 00:01:34 2009
From: shimon at CYCLONE.TAU.AC.IL (Simon Krichak)
Date: Tue, 20 Jan 2009 07:01:34 +0200
Subject: Calculation of  Madden-Julian Oscillation and other
 teleconnection indices
Message-ID: <TUE.20.JAN.2009.070134.0200.>

Dear Grads users,

I want to calculate index of Madden-Julian Oscillation and several other teleconnection indices 
based on a global data archive using GrADS.
Does anybody have (and hopefully ready to share) a ready script for performing it ?
Will appreciate any help or recommendation.

Simon

  Simon Krichak
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090120/fcccd41b/attachment.html 

From mitchell at ATMOS.WASHINGTON.EDU  Tue Jan 20 00:21:08 2009
From: mitchell at ATMOS.WASHINGTON.EDU (Todd MItchell)
Date: Mon, 19 Jan 2009 21:21:08 -0800
Subject: Calculation of  Madden-Julian Oscillation and other
 teleconnection indices
In-Reply-To: <562046C5D8194C7A8B5937B640E48EE3@shimonlaptop>
Message-ID: <MON.19.JAN.2009.212108.0800.>

Hi Simon,                           Monday 19 January 2009

Wheeler and Hendon provide a MJO index at
http://www.bom.gov.au/bmrc/clfor/cfstaff/matw/maproom/RMM/

Todd Mitchell


Simon Krichak wrote:
> Dear Grads users,
>
> I want to calculate index of Madden-Julian Oscillation and several
> other teleconnection indices
> based on a global data archive using GrADS.
> Does anybody have (and hopefully ready to share) a ready script for
> performing it ?
> Will appreciate any help or recommendation.
>
> Simon
>
>
>     Simon Krichak
>


--
Joint Institute for the Study of the Atmosphere and Ocean (JISAO)
University of Washington; Box 355672, Seattle WA 98195-5672
Phone: 206 685 3786; Fax: 206 685 3397; jisao.washington.edu
Street Address: 3737 Brooklyn Ave. NE, Seattle WA 98105


From shimon at CYCLONE.TAU.AC.IL  Tue Jan 20 00:36:31 2009
From: shimon at CYCLONE.TAU.AC.IL (Simon Krichak)
Date: Tue, 20 Jan 2009 07:36:31 +0200
Subject: Calculation of  Madden-Julian Oscillation and other
 teleconnection indices
Message-ID: <TUE.20.JAN.2009.073631.0200.>

Thank you Todd,

I will look at it carefully. Still, I would like to know how to calculate
the index based on the model data I have.

Thanks again,
   Simon

----- Original Message -----
From: "Todd MItchell" <mitchell at ATMOS.WASHINGTON.EDU>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Tuesday, January 20, 2009 7:21 AM
Subject: Re: Calculation of Madden-Julian Oscillation and other
teleconnection indices


> Hi Simon,                           Monday 19 January 2009
>
> Wheeler and Hendon provide a MJO index at
> http://www.bom.gov.au/bmrc/clfor/cfstaff/matw/maproom/RMM/
>
> Todd Mitchell
>
>
> Simon Krichak wrote:
>> Dear Grads users,
>>
>> I want to calculate index of Madden-Julian Oscillation and several
>> other teleconnection indices
>> based on a global data archive using GrADS.
>> Does anybody have (and hopefully ready to share) a ready script for
>> performing it ?
>> Will appreciate any help or recommendation.
>>
>> Simon
>>
>>
>>     Simon Krichak
>>
>
>
> --
> Joint Institute for the Study of the Atmosphere and Ocean (JISAO)
> University of Washington; Box 355672, Seattle WA 98195-5672
> Phone: 206 685 3786; Fax: 206 685 3397; jisao.washington.edu
> Street Address: 3737 Brooklyn Ave. NE, Seattle WA 98105


--------------------------------------------------------------------------------



No virus found in this incoming message.
Checked by AVG - http://www.avg.com
Version: 8.0.176 / Virus Database: 270.10.9/1900 - Release Date: 18/01/2009
12:11


From romero619 at HOTMAIL.COM  Tue Jan 20 01:08:45 2009
From: romero619 at HOTMAIL.COM (P.Romero)
Date: Tue, 20 Jan 2009 00:08:45 -0600
Subject: printim directly to stdout
Message-ID: <TUE.20.JAN.2009.000845.0600.>

Would it be possible (on linux) to redirect printim's binary output directly
to stdout?

Just curious.

P.Romero

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090120/59b5257a/attachment.html 

From romero619 at HOTMAIL.COM  Tue Jan 20 12:44:00 2009
From: romero619 at HOTMAIL.COM (P.Romero)
Date: Tue, 20 Jan 2009 11:44:00 -0600
Subject: source code for solid/filled arrowhead
Message-ID: <TUE.20.JAN.2009.114400.0600.>

In case anyone is interested in having solid/filled arrowheads for 'gxout
vector', here's the source code I came up with:



//modify gagx.c source file

//@ around line 3265

//replace entire gaarrw function, and add/modify the static 'a150'/'a165'
variables beforehand



static gadouble a150 = 150.0*3.141592654/180.0;

static gadouble a165 = 165.0*3.141592654/180.0;

void gaarrw (gadouble x, gadouble y, gadouble ang, gadouble siz, gadouble
asiz) {

gadouble xx,yy;

gadouble xy[8];



  if (siz<0.0001) {

    gxmark (2, x, y, 0.01);

    return;

  }

  gxplot (x,y,3);

  xx = x+siz*cos(ang);

  yy = y+siz*sin(ang);

  gxplot (xx,yy,2);

  if (asiz==0.0) return;

  if (asiz<0.0) asiz = -1.0*asiz*siz;



  //new filled arrowhead

  xy[0]=xx;

  xy[1]=yy;

  xy[2]=xx+asiz*cos(ang+a165);

  xy[3]=yy+asiz*sin(ang+a165);

  xy[4]=xx+asiz*cos(ang-a165);

  xy[5]=yy+asiz*sin(ang-a165);

  xy[6]=xy[0];

  xy[7]=xy[1];

  gxfill(xy,4);



   //original non-filled arrowhead

  //plot arrow head lines

  //first plot one line, then move back to tip, then plot second line

  /*

  gxplot (xx+asiz*cos(ang+a150),yy+asiz*sin(ang+a150),2);

  gxplot (xx,yy,3);

  gxplot (xx+asiz*cos(ang-a150),yy+asiz*sin(ang-a150),2);

  */

}



.I tested this and it seems to work fine.

I changed the arrowhead angle to 165 degrees because the arrowheads seemed
too wide at 150, just my opinion.

I guess ideally this function should have arguments for "isFilled" and
"arrowHeadAngle", to allow the user to select standard/filled arrowhead, and
the angle used to determine the arrowhead's spread.

P.Romero

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090120/95126569/attachment.html 

From paul.earnshaw at METOFFICE.GOV.UK  Wed Jan 21 10:14:58 2009
From: paul.earnshaw at METOFFICE.GOV.UK (Paul Earnshaw)
Date: Wed, 21 Jan 2009 15:14:58 +0000
Subject: New Documentation for the Ensemble Dimension
In-Reply-To: <4AF851A9-CC2E-4E05-A199-0D7C74E36B8D@cola.iges.org>
Message-ID: <WED.21.JAN.2009.151458.0000.>

Jennifer,

Thank you for the extra documentation, it has been a great help.

There was one thing I think may need to be added, and I am sorry if this
is in the documentation and I have missed it somewhere. Could you make
it clear how GrADS identifies the ensemble dimension in SDF files?

By this I mean having either the attribute "axis" or "grads_dim" set to
"E" or "e" for the ensemble dimension.

I did eventually work it out by digging in the source code, and I think
it would be useful in the main documentation.

Cheers,
Paul

On Thu, 2009-01-15 at 08:57 -0500, Jennifer Adams wrote:
> Hi, Everyone --
> I have added a new section to the User's Guide that explains how the
> ensemble dimension has been implemented in GrADS 2.0. Additional
> topics include how to create a descriptor file for an ensemble data
> set, how to organize your data using templating, special
> considerations for certain data formats, and some examples (with
> illustrations!) using real data.
>
> http://iges.org/grads/gadoc/ensembles.html
>
> Real time NCEP GFS Ensemble forecasts are behind our GDS at http://monsoondata.org:9090/dods/gfsens/
> Why not take a moment to fire up 'gradsdap', use 'sdfopen' to access
> the latest forecast, and start playing with ensemble data. It's cool!
>
> --Jennifer
--
Dr. Paul Earnshaw   Global Model Development and Diagnostics
Met Office  FitzRoy Road  Exeter  Devon  EX1 3PB  United Kingdom
Tel: +44 (0)1392 884525  Fax: +44 (0)1392 885681
email: paul.earnshaw at metoffice.gov.uk website: www.metoffice.gov.uk

Met Office climate change predictions can now be viewed on Google Earth
http://www.metoffice.gov.uk/research/hadleycentre/google/


From paul.earnshaw at METOFFICE.GOV.UK  Wed Jan 21 10:20:06 2009
From: paul.earnshaw at METOFFICE.GOV.UK (Paul Earnshaw)
Date: Wed, 21 Jan 2009 15:20:06 +0000
Subject: New Documentation for the Ensemble Dimension
In-Reply-To: <4AF851A9-CC2E-4E05-A199-0D7C74E36B8D@cola.iges.org>
Message-ID: <WED.21.JAN.2009.152006.0000.>

Jennifer,

Thank you for the extra documentation, it has been a great help.

There was one thing I think may need to be added, and I am sorry if this
is in the documentation and I have missed it somewhere. Could you make
it clear how GrADS identifies the ensemble dimension in SDF files?

By this I mean having either the attribute "axis" or "grads_dim" set to
"E" or "e" for the ensemble dimension.

I did eventually work it out by digging in the source code, and I think
it would be useful in the main documentation.

Cheers,
Paul

--
Dr. Paul Earnshaw   Global Model Development and Diagnostics
Met Office  FitzRoy Road  Exeter  Devon  EX1 3PB  United Kingdom
Tel: +44 (0)1392 884525  Fax: +44 (0)1392 885681
email: paul.earnshaw at metoffice.gov.uk website: www.metoffice.gov.uk

Met Office climate change predictions can now be viewed on Google Earth
http://www.metoffice.gov.uk/research/hadleycentre/google/


From jma at COLA.IGES.ORG  Wed Jan 21 10:25:07 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Wed, 21 Jan 2009 10:25:07 -0500
Subject: New Documentation for the Ensemble Dimension
In-Reply-To: <1232550898.17887.4.camel@eld188.desktop.frd.metoffice.com>
Message-ID: <WED.21.JAN.2009.102507.0500.>

Hi, Paul --
That is a very good suggestion, something I had overlooked. I will add
something about ensemble metadata to the sdfopen doc as well as the
ensemble doc. Thanks!
--Jennifer

On Jan 21, 2009, at 10:14 AM, Paul Earnshaw wrote:

> Jennifer,
>
> Thank you for the extra documentation, it has been a great help.
>
> There was one thing I think may need to be added, and I am sorry if
> this
> is in the documentation and I have missed it somewhere. Could you make
> it clear how GrADS identifies the ensemble dimension in SDF files?
>
> By this I mean having either the attribute "axis" or "grads_dim" set
> to
> "E" or "e" for the ensemble dimension.
>
> I did eventually work it out by digging in the source code, and I
> think
> it would be useful in the main documentation.
>
> Cheers,
> Paul
>
> On Thu, 2009-01-15 at 08:57 -0500, Jennifer Adams wrote:
>> Hi, Everyone --
>> I have added a new section to the User's Guide that explains how the
>> ensemble dimension has been implemented in GrADS 2.0. Additional
>> topics include how to create a descriptor file for an ensemble data
>> set, how to organize your data using templating, special
>> considerations for certain data formats, and some examples (with
>> illustrations!) using real data.
>>
>> http://iges.org/grads/gadoc/ensembles.html
>>
>> Real time NCEP GFS Ensemble forecasts are behind our GDS at http://monsoondata.org:9090/dods/gfsens/
>> Why not take a moment to fire up 'gradsdap', use 'sdfopen' to access
>> the latest forecast, and start playing with ensemble data. It's cool!
>>
>> --Jennifer
> --
> Dr. Paul Earnshaw   Global Model Development and Diagnostics
> Met Office  FitzRoy Road  Exeter  Devon  EX1 3PB  United Kingdom
> Tel: +44 (0)1392 884525  Fax: +44 (0)1392 885681
> email: paul.earnshaw at metoffice.gov.uk website: www.metoffice.gov.uk
>
> Met Office climate change predictions can now be viewed on Google
> Earth
> http://www.metoffice.gov.uk/research/hadleycentre/google/

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090121/c113991c/attachment.html 

From msponsler at COMCAST.NET  Wed Jan 21 18:54:18 2009
From: msponsler at COMCAST.NET (Mark Sponsler)
Date: Thu, 22 Jan 2009 00:54:18 +0100
Subject: Square Root Calculation
Message-ID: <THU.22.JAN.2009.005418.0100.>

Hello,
I'm trying to calculate the square root of a stack of variables in a Grads
script.  Is there any way to do this?

e.g.

answer  = sqrt(val1 + val2 + val3 + val4)

Any help would be greatly appreciated.


From romero619 at HOTMAIL.COM  Wed Jan 21 19:03:04 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Wed, 21 Jan 2009 16:03:04 -0800
Subject: Square Root Calculation
In-Reply-To: <20090121235515.D81A8217AD@mx2.cineca.it>
Message-ID: <WED.21.JAN.2009.160304.0800.>


hey mark,
I took a look at the source code
math_sqrt() is the function you're looking for.
 
P.Romero
> Date: Thu, 22 Jan 2009 00:54:18 +0100> From: msponsler at COMCAST.NET> Subject: Square Root Calculation> To: GRADSUSR at LIST.CINECA.IT> > Hello,> I'm trying to calculate the square root of a stack of variables in a Grads> script. Is there any way to do this?> > e.g.> > answer = sqrt(val1 + val2 + val3 + val4)> > Any help would be greatly appreciated.
_________________________________________________________________
Windows Live?: E-mail. Chat. Share. Get more ways to connect. 
http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t2_allup_howitworks_012009
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090121/d60856a6/attachment.html 

From ccjrl at MFCSV.CWB.GOV.TW  Wed Jan 21 19:39:04 2009
From: ccjrl at MFCSV.CWB.GOV.TW (Hsin-Hsing,Chia)
Date: Thu, 22 Jan 2009 08:39:04 +0800
Subject: the whiskers
In-Reply-To: <F01423DC-C25E-48AB-8216-8230B96FC5BF@cola.iges.org>
Message-ID: <THU.22.JAN.2009.083904.0800.>

Hi Jennifer:
 Thanks for your effort on the documentation about the ensemble dataset.
How to plot the figure 4, which has the blue whiskers in this documentation ?

 Thanks

Hsin-Hsing


From romero619 at HOTMAIL.COM  Wed Jan 21 19:24:11 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Wed, 21 Jan 2009 16:24:11 -0800
Subject: GRADS user list rejecting posts
Message-ID: <WED.21.JAN.2009.162411.0800.>

I keep receiving messages from the listserver saying that my posting or replies to postings have been rejected.
Is anyone else experiencing this problem?
P.Romero
_________________________________________________________________
Windows Live? Hotmail?:?more than just e-mail.
http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t2_hm_justgotbetter_explore_012009


From robock at ENVSCI.RUTGERS.EDU  Wed Jan 21 19:31:19 2009
From: robock at ENVSCI.RUTGERS.EDU (Alan Robock)
Date: Wed, 21 Jan 2009 19:31:19 -0500
Subject: GRADS user list rejecting posts
In-Reply-To: <BAY117-W2475B7524038D104E6659AF9CE0@phx.gbl>
Message-ID: <WED.21.JAN.2009.193119.0500.>

Dear Pablo,

I got the same message last week, but the message went through anyway.
Could whoever runs the listserver please fix this?

Alan

Alan Robock, Professor II
  Director, Meteorology Undergraduate Program
  Associate Director, Center for Environmental Prediction
Department of Environmental Sciences        Phone: +1-732-932-9800 x6222
Rutgers University                                  Fax: +1-732-932-8644
14 College Farm Road                   E-mail: robock at envsci.rutgers.edu
New Brunswick, NJ 08901-8551  USA      http://envsci.rutgers.edu/~robock



Pablo Romero wrote:

>I keep receiving messages from the listserver saying that my posting or replies to postings have been rejected.
>Is anyone else experiencing this problem?
>P.Romero
>_________________________________________________________________
>Windows Live? Hotmail?:?more than just e-mail.
>http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t2_hm_justgotbetter_explore_012009
>
>


From jma at COLA.IGES.ORG  Thu Jan 22 14:01:17 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Thu, 22 Jan 2009 14:01:17 -0500
Subject: the whiskers
In-Reply-To: <20090122003524.M15169@mfcsv.cwb.gov.tw>
Message-ID: <THU.22.JAN.2009.140117.0500.>

To draw the box and whiskers, use 'set gxout bar' and 'set gxout
errbar'. Then give the display command two variables, separated by a a
semi-colon. Here's the script that I used to draw the first four plots
in the ensemble documentation:

'open gfsens.2009010100.ctl'
'set lon -60'
'set lat 45'
'set lev 500'
'set parea 1 9.7 2.5 6'
'set t 0.5 65.5'
'set e 1'
'clear'
'set grads off'
'set tlsupp year'
'set cmark 0'
'set vrange 4950 5600'
'set ylint 100'
'd hgt'
'printim edemo1.png '

* Remaining members
e=2
cols='2 3 4 5 6 8 9 10 11 12 13 14 2 3 4 5 6  8 9 10 11 12 13 14 2 3 4
5 6 7 8 9 10 11 12 13 14 '
while (e<=21)
   'set e 'e
   c=subwrd(cols,e)
   'set ccolor 'c
   'set cmark 0'
   'd hgt'
   e=e+1
endwhile
'printim edemo2.png '

* all 21 members drawn as shaded boxes
'clear'
'set gxout grfill'
'set e 0.5 21.5'
'set clevs  5100 5150 5200 5250 5300 5350 5400 5450 '
'set ccols 14   4    11   5    3    7    12   8    2    '
'set grads off'
'set tlsupp year'
'd hgt'
'cbarn 0.8 1'
'printim edemo3.png'

* Calculate the ensemble mean
'set e 1'
'set t 1 last'
'define ensmean=ave(hgt,e=1,e=21)'

* Calculate the variance
diffsq = 'pow(hgt-ensmean,2)'
variance = 'ave('diffsq',e=1,e=21)'
'define stddev=sqrt('variance')'

* Calculate the min/max
'define ensmin=tloop(min(hgt,e=1,e=21))'
'define ensmax=tloop(max(hgt,e=1,e=21))'

* Plot the results
'clear'
'set t 0.5 65.5'
'set vrange 4950 5600'
'set ylint 100'
'set tlsupp year'
'set grads off'

* Draw error bars for min/max
'set gxout errbar'
'set ccolor 4'
'd ensmin;ensmax'

* Draw bars for +/- standard deviation
plus  = '(ensmean+stddev)'
minus = '(ensmean-stddev)'
'set gxout bar'
'set bargap 50'
'set baropts outline'
'set ccolor 3'
'd 'minus';'plus

* Draw line for Ensemble mean
'set gxout line'
'set cmark 0'
'set cthick 6'
'set digsiz 0.05'
'set ccolor 2'
'd ensmean'

'printim edemo4.png'



On Jan 21, 2009, at 7:39 PM, Hsin-Hsing,Chia wrote:

> Hi Jennifer:
> Thanks for your effort on the documentation about the ensemble
> dataset.
> How to plot the figure 4, which has the blue whiskers in this
> documentation ?
>
> Thanks
>
> Hsin-Hsing

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090122/daff9600/attachment.html 

From maldmm5 at GMAIL.COM  Thu Jan 22 17:28:09 2009
From: maldmm5 at GMAIL.COM (Mald MM5)
Date: Fri, 23 Jan 2009 03:28:09 +0500
Subject: GRADS user list rejecting posts
In-Reply-To: <4977BE57.2040001@envsci.rutgers.edu>
Message-ID: <FRI.23.JAN.2009.032809.0500.>

Me too had the same problem...But the messge seems to have gone through...

On Thu, Jan 22, 2009 at 5:31 AM, Alan Robock <robock at envsci.rutgers.edu>wrote:

> Dear Pablo,
>
> I got the same message last week, but the message went through anyway.
> Could whoever runs the listserver please fix this?
>
> Alan
>
> Alan Robock, Professor II
>  Director, Meteorology Undergraduate Program
>  Associate Director, Center for Environmental Prediction
> Department of Environmental Sciences        Phone: +1-732-932-9800 x6222
> Rutgers University                                  Fax: +1-732-932-8644
> 14 College Farm Road                   E-mail: robock at envsci.rutgers.edu
> New Brunswick, NJ 08901-8551  USA      http://envsci.rutgers.edu/~robock<http://envsci.rutgers.edu/%7Erobock>
>
>
>
>
> Pablo Romero wrote:
>
>  I keep receiving messages from the listserver saying that my posting or
>> replies to postings have been rejected.
>> Is anyone else experiencing this problem?
>> P.Romero
>> _________________________________________________________________
>> Windows Live? Hotmail(R):?more than just e-mail.
>>
>> http://windowslive.com/explore?ocid=TXT_TAGLM_WL_t2_hm_justgotbetter_explore_012009
>>
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090123/b696c0cc/attachment.html 

From fil at METEOPT.COM  Fri Jan 23 09:00:06 2009
From: fil at METEOPT.COM (Fil)
Date: Fri, 23 Jan 2009 15:00:06 +0100
Subject: Precipitation 3/6 h
Message-ID: <FRI.23.JAN.2009.150006.0100.>

Hi everyone,

I'm trying to plot a precipitation meteogram from GFS 3h data but, for
example, +33h timestep has precipitation from +30h to +33h but +36h timestep
also has precipitation from +30h to +33h and +33h to +36h. So what I wanted
to know is how can I subtract precipitation like 30h = 30h - 27h, 36h = 36h
- 33h, 42h = 42h - 39h, etc, so that +30h/+36h/+42h/etc timestep doesn't
include precipitation from +27h/+33h/+39h/etc timestep.

Also how can I sum precipitation from +3h to 180h? I'm trying to draw the
total precipitation above the bars graph and I tried something like:

totalprec=sum(APCPsfc,t=2,t=180)
'draw string x y 'totalprec

But it doesn't work at all.

Hope someone can help me.

Cheers.


From hmjbarbosa at GMAIL.COM  Fri Jan 23 09:43:06 2009
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Fri, 23 Jan 2009 12:43:06 -0200
Subject: grads script function return values and more
Message-ID: <FRI.23.JAN.2009.124306.0200.>

Dear all,

Some grads' intrinsic functions have return values with are longer
than a single line of text.
For example, when we do 'q dims' inside a gs script a variable called
result is defined
and we can read lines and words from it using sublin() and subwrd().

The question is: is it possible to do that with a grads script function?
I another words: is it possible to define a string variable which contains more
than one line of text? how? If so, can we give it as a function return value?

Second question:

When I call a script function grads will look for a file named like it in all
directories defined in the system variable $GASCRP and then executes
that script.... While the script function is being executed, is there anyway
to know the full file name of the script?
In KSH, for example, one can access the script name with $0

Thanks,
Henrique


From romero619 at HOTMAIL.COM  Fri Jan 23 10:39:52 2009
From: romero619 at HOTMAIL.COM (Pablo Romero)
Date: Fri, 23 Jan 2009 07:39:52 -0800
Subject: pygrads plot configuration
Message-ID: <FRI.23.JAN.2009.073952.0800.>

How can plot/map options be passed or adjusted in the galab object?
I cant seem to override the drawing of axes or the color bar, nor can I change the plot dimensions when using basemap through galab.

Does anyone have much experience using pygrads/galab, and customizing the projection plot area?

P.Romero
_________________________________________________________________
Windows Live?: E-mail. Chat. Share. Get more ways to connect.
http://windowslive.com/howitworks?ocid=TXT_TAGLM_WL_t2_allup_howitworks_012009


From doty at COLA.IGES.ORG  Fri Jan 23 11:40:22 2009
From: doty at COLA.IGES.ORG (Brian Doty)
Date: Fri, 23 Jan 2009 11:40:22 -0500
Subject: WeirdX
Message-ID: <FRI.23.JAN.2009.114022.0500.>

Hello everyone, I am curious if anyone has used WeirdX, which appears
to be an Xserver implemented in the Java language.  Particularly, if
anyone has used it with GrADS?  Thanks!  ....Brian


From arlindo.dasilva at GMAIL.COM  Fri Jan 23 12:36:16 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Fri, 23 Jan 2009 12:36:16 -0500
Subject: WeirdX
In-Reply-To: <EF91A9D6-B5BB-44BE-BB8B-3A16962AB3AC@cola.iges.org>
Message-ID: <FRI.23.JAN.2009.123616.0500.>

On Fri, Jan 23, 2009 at 11:40 AM, Brian Doty <doty at cola.iges.org> wrote:

> Hello everyone, I am curious if anyone has used WeirdX, which appears
> to be an Xserver implemented in the Java language.  Particularly, if
> anyone has used it with GrADS?  Thanks!  ....Brian
>

Yes, I have experimented with it and even considered using it for Win32
GrADS. I ended up going with Xming,

                 http://www.straightrunning.com/XmingNotes/

a native WIn32 build based on the x.org sources; currently Xming is bundled
with Win32 GrADS running off display :7. The main reason was that Xming was
relocatable (can be run from a USB mem stick), much faster, smaller and did
not require the overhead of the JVM. Otherwise, WeirdX seemed to work OK
with GrADS.

    Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090123/e1786aa2/attachment.html 

From arlindo.dasilva at GMAIL.COM  Fri Jan 23 13:04:46 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Fri, 23 Jan 2009 13:04:46 -0500
Subject: pygrads plot configuration
In-Reply-To: <BAY117-W55F71E151A6298E80B283EF9CF0@phx.gbl>
Message-ID: <FRI.23.JAN.2009.130446.0500.>

On Fri, Jan 23, 2009 at 10:39 AM, Pablo Romero <romero619 at hotmail.com>wrote:

> How can plot/map options be passed or adjusted in the galab object?
> I cant seem to override the drawing of axes or the color bar, nor can I
> change the plot dimensions when using basemap through galab.
>
> Does anyone have much experience using pygrads/galab, and customizing the
> projection plot area?
>


At this point galab provides syntatic sugar around Matplotlib and the
Basemap toolkit so that someone familiar with GrADS could get some basic
plots done with minimum effort. To effectively fine tune your plots,
though,  you need to understand how Matrplotlib/Basemap work; look at the
wiki for references on Matplotlib documentation.  Here are a few pointers:

   - The method ga.basemap() provides a front end to the Basemap() class,
   where the current GrADS dimension environment is used to specify the domain
   for the Map. The Basemap object is saved as an attribute of the grads
   object: *ga.map*
   - Likewise, when using ga.imshow(), ga.contour() the axis and maps for
   the plot are derived from the dimension enviroment. Indeed, these are
   wrappers around Basemap.imshow() and Basemap.contour().   Many of the
   options you give to ga.imshow() are passed down to Basemap.imshow(), so you
   need to look at the documentation for this; for the axis labeling you can
   specify dlon/dlat for the spacing in degrees between lat/lon labels.

If you tell me exactly what you are trying to accomplish I may be able to
give you more specific pointers or even write a recipe for it.

        Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090123/1af21121/attachment.html 

From dasilva at ALUM.MIT.EDU  Fri Jan 23 13:55:44 2009
From: dasilva at ALUM.MIT.EDU (Arlindo da Silva)
Date: Fri, 23 Jan 2009 13:55:44 -0500
Subject: pygrads plot configuration
In-Reply-To: <BAY117-W177901EB75711137B6D60BF9CF0@phx.gbl>
Message-ID: <FRI.23.JAN.2009.135544.0500.>

On Fri, Jan 23, 2009 at 1:08 PM, Pablo Romero <romero619 at hotmail.com> wrote:

>
> ok, thank you.
> Ive started digging through the galab.py source code, to see how its
> interfacing the matplotlib/basemap objects,and doing the plotting through
> matplotlib.


> question: if I want to avoid using galab, and simply instantiate my own
> matplotlib basemap object, do I simply need to export a variable from galab
> and pass it to matplotlib/basemap/contourf? or, do I need to do some sort of
> translation/interpolation/preparation???
>

No, simply export the array:

ts = ga.expr('ts')

The "ts" object will be a "GrADSField" object, which is a subclass of a
Masked Array class. A masked array has 3 main attributes

ts.data   - a numpy numeric array
ts.mask - a boolean array of same shape as ts.data which is used for
handling undefs.
ts.grid - information about the coordinates, dimension environment

So, ts can be used anywhere you could use a masked array.  Some functions
may not like a masked array, in which case you can use "ts.data" instead.

Notice that you can use either ga.exp() or ga.expr() with v1.9, but only
ga.expr() with v2.0. In principle ga.expr is more flexible, you can get back
x-y sections, for example. However, 'expr' is new, let me know if you run
into any snags.

   Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090123/d0609865/attachment.html 

From shimon at CYCLONE.TAU.AC.IL  Sun Jan 25 08:23:55 2009
From: shimon at CYCLONE.TAU.AC.IL (Simon Krichak)
Date: Sun, 25 Jan 2009 15:23:55 +0200
Subject: Calculation of south-Asian monsoon index
Message-ID: <SUN.25.JAN.2009.152355.0200.>

Dear Grads users,

I want to calculate index of south-Asian monsoon based on a global data
archive using GrADS.
Does anybody have (and hopefully ready to share) a ready script for
performing it ?
Will appreciate any help or recommendation.

Best,
 Simon Krichak


From dasilva at OPENGRADS.ORG  Mon Jan 26 11:33:39 2009
From: dasilva at OPENGRADS.ORG (Arlindo da Silva)
Date: Mon, 26 Jan 2009 11:33:39 -0500
Subject: duvida sobre o grads
In-Reply-To: <846824.14898.qm@web36205.mail.mud.yahoo.com>
Message-ID: <MON.26.JAN.2009.113339.0500.>

Dear Claudia,

  Please send questions like this to gradsusr, there are a lot more users
there that would be able to help you as well. In regards to your question,
you can put you data files in any directory of your choice. It is
recommended that you use Unix (posix) style syntax to refer to your data,
with forward slashes "/". For example, if your data file is in

   D:\PASTA\OCEANO\currents.nc

you can open this file in grads with the following command:

ga-> sdfopen d:/pasta/oceano/currents.nc

This should work as well,

ga-> sdfopen D:\PASTA\OCEANO\currents.nc

but the first syntax is preferred.

           Arlindo



2009/1/26 claudia klose parise <clau_oceano at yahoo.com.br>

> Ola Arlindo,
>
> Encontrei seu nome em muitos foruns sobre o GRADS e acho que vcoe pode me
> ajudar.
> Instalei o GRADS no meu computador que tem windows xp. So que os aqrquivos
> que o grads le precisam ser colocados no diretorio: C:\Arquivos de
> programas\GrADS20\data senao ele nao os reconhece. Porem gostaria colocar
> esse arquivo que serao lidos pelo GRADS dentro de uma pasta no D.
> Usei o GRADS (grads-2.0.a3.oga.1.win32_superpack)
> Como fa?o para mudar esse diretorio?
>
> Obrigada,
> Claudia Parise
>
> ------------------------------
> Veja quais s?o os assuntos do momento no Yahoo! + Buscados: Top 10<http://br.rd.yahoo.com/mail/taglines/mail/*http://br.maisbuscados.yahoo.com/>-
> Celebridades<http://br.rd.yahoo.com/mail/taglines/mail/*http://br.maisbuscados.yahoo.com/celebridades/>-
> M?sica<http://br.rd.yahoo.com/mail/taglines/mail/*http://br.maisbuscados.yahoo.com/m%C3%BAsica/>-
> Esportes<http://br.rd.yahoo.com/mail/taglines/mail/*http://br.maisbuscados.yahoo.com/esportes/>
>



-- 
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090126/36f1961a/attachment.html 

From hmjbarbosa at GMAIL.COM  Mon Jan 26 12:45:41 2009
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Mon, 26 Jan 2009 15:45:41 -0200
Subject: problem with 365_day calender under grads v2.0.a4.oga.1
Message-ID: <MON.26.JAN.2009.154541.0200.>

Dear All,

I have  a few netCDF files which use a 365day calendar. I used to open these
files using 'sdfopen' directly without any problem. Yesterday I
installed Grads2.0a4
(actually the open grads version) and now it complains about the calendar!

I wonder why does the new version reject the 365day calendar netCDF files,
while the previous versions accepted.... Is this a bug in the new grads or this
change was meant to be?

Of course I could write a ctl... but its so much easier to use the self
contained netCDF that I would rather go back to 1.9b4....

[]'s
Henrique

Scanning self-describing file:
original/cccma/pcmdi.ipcc4.cccma_cgcm3_1_t63.20c3m.run1.monthly.tos_o1_20c3m_1_cgcm3.1_t63_1850_2000.nc
SDF Error: 365 day calendars are not COARDS compliant.
  You can still open this file with GrADS if you write a descriptor file that
  contains a complete TDEF entry and OPTIONS 365_day_calendar. Please read the
  documentation at http://grads.iges.org/grads/gadoc/SDFdescriptorfile.html
ga->


From jma at COLA.IGES.ORG  Mon Jan 26 14:40:56 2009
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Mon, 26 Jan 2009 14:40:56 -0500
Subject: problem with 365_day calender under grads v2.0.a4.oga.1
In-Reply-To: <e71668fa0901260945l5c8d75b0vf09089ba6cdadd30@mail.gmail.com>
Message-ID: <MON.26.JAN.2009.144056.0500.>

The reason why you can't use sdfopen to open 365_day_calendar files is
because the udunits library (used by the sdfopen code) does not
support 365_day_calendars and so the initial time of your file may be
interpreted incorrectly. This error comes up often when using nc files
from the IPCC -- you must provide all the time axis information to
GrADS to make sure it is correct.
Jennifer


On Jan 26, 2009, at 12:45 PM, Henrique Barbosa wrote:

> Dear All,
>
> I have  a few netCDF files which use a 365day calendar. I used to
> open these
> files using 'sdfopen' directly without any problem. Yesterday I
> installed Grads2.0a4
> (actually the open grads version) and now it complains about the
> calendar!
>
> I wonder why does the new version reject the 365day calendar netCDF
> files,
> while the previous versions accepted.... Is this a bug in the new
> grads or this
> change was meant to be?
>
> Of course I could write a ctl... but its so much easier to use the
> self
> contained netCDF that I would rather go back to 1.9b4....
>
> []'s
> Henrique
>
> Scanning self-describing file:
> original/cccma/
> pcmdi
> .ipcc4
> .cccma_cgcm3_1_t63.20c3m
> .run1.monthly.tos_o1_20c3m_1_cgcm3.1_t63_1850_2000.nc
> SDF Error: 365 day calendars are not COARDS compliant.
>  You can still open this file with GrADS if you write a descriptor
> file that
>  contains a complete TDEF entry and OPTIONS 365_day_calendar. Please
> read the
>  documentation at http://grads.iges.org/grads/gadoc/SDFdescriptorfile.html
> ga->

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090126/afa60109/attachment.html 

From cbutenhoff at YAHOO.COM  Mon Jan 26 17:27:17 2009
From: cbutenhoff at YAHOO.COM (Chris Butenhoff)
Date: Mon, 26 Jan 2009 23:27:17 +0100
Subject: Memory Allocation Error: grid storage
Message-ID: <MON.26.JAN.2009.232717.0100.>

On my attempt to display a variable in Grads v1.9b4 I receive the following
message:

>>>Memory Allocation Error:  grid storage
>>>Data Request Error:  Error for variable 'lc'

The open file is about 600 Mb in size and has a lon/lat grid of 40320x16353.

My machine has over 3 Gb of RAM so I assume the memory limit must be set in
Grads.

Searching the archives I find the following info:

>>> "The limit of the cache size is (for now) hard-coded in gaio around
>>> line 76:
>>> #define MAXG2CACHE 500100100
>>> You can change this and recompile if you want to limit or expand your
>>> memory usage."

Before attempting, does anyone have experience with this problem, and is
this the right fix?

Thanks!

_Chris


From arlindo.dasilva at GMAIL.COM  Mon Jan 26 19:50:01 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Mon, 26 Jan 2009 19:50:01 -0500
Subject: Memory Allocation Error: grid storage
In-Reply-To: <20090126222815.7104221605@mx2.cineca.it>
Message-ID: <MON.26.JAN.2009.195001.0500.>

On Mon, Jan 26, 2009 at 5:27 PM, Chris Butenhoff <cbutenhoff at yahoo.com>wrote:

> On my attempt to display a variable in Grads v1.9b4 I receive the following
> message:
>
> >>>Memory Allocation Error:  grid storage
> >>>Data Request Error:  Error for variable 'lc'
>
> The open file is about 600 Mb in size and has a lon/lat grid of
> 40320x16353.
>

Let's do the math. The size of the file hardly matters as it may include
compression. Grads v1.9b4 holds gridded data in 4 byte floats.  So, a
singlevariable at your grid resolution would require:

40320*16353*4/(1024*1024) ~ 2.5 GB

So, you could be very close to your memory limit of 3GB. The particular
format also matters; there is a known memory leak with some form of
netcdf/hdf file.


>
> My machine has over 3 Gb of RAM so I assume the memory limit must be set in
> Grads.
>
> Searching the archives I find the following info:


> >>> "The limit of the cache size is (for now) hard-coded in gaio around
> >>> line 76:
> >>> #define MAXG2CACHE 500100100
> >>> You can change this and recompile if you want to limit or expand your
> >>> memory usage."
>
> Before attempting, does anyone have experience with this problem, and is
> this the right fix?
>

 This is a Grib-2 specific issue, certainly not an issue with v1.9b4 which
does not support Grib-2.

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090126/ac2e1183/attachment.html 

From aidab at PDX.EDU  Mon Jan 26 22:11:26 2009
From: aidab at PDX.EDU (Aida Biberic)
Date: Mon, 26 Jan 2009 19:11:26 -0800
Subject: land cover
In-Reply-To: <20090126222815.7104221605@mx2.cineca.it>
Message-ID: <MON.26.JAN.2009.191126.0800.>

Chris,

you may be familiar with this website already but just in case you are
not, check this out:

http://glcf.umiacs.umd.edu/data/landcover/

I haven't looked closely to be sure, but I think these are also
available as GIS files.


Aida


From furue at HAWAII.EDU  Tue Jan 27 02:38:19 2009
From: furue at HAWAII.EDU (Ryo Furue)
Date: Mon, 26 Jan 2009 21:38:19 -1000
Subject: problem with 365_day calender under grads v2.0.a4.oga.1
In-Reply-To: <F041FA5E-4CD0-4102-95A1-F24F9B703471@cola.iges.org>
Message-ID: <MON.26.JAN.2009.213819.1000.>

Hi

>     SDF Error: 365 day calendars are not COARDS compliant.

I think this error message is misleading.  The COARDS convention
doesn't talk about calendars, if I'm not mistaken.
For example,

  http://ferret.wrc.noaa.gov/noaa_coop/coop_cdf_profile.html

| The reason why you can't use sdfopen to open 365_day_calendar files
| is because the udunits library (used by the sdfopen code) does not
| support 365_day_calendars and so the initial time of your file may
| be interpreted incorrectly. This error comes up often when using nc
| files from the IPCC

For your information, the "365_day" or "noleap" calendar is defined
in the CF convention:

   http://cf-pcmdi.llnl.gov/

Both the COARDS and CF conventions are listed in

  http://www.unidata.ucar.edu/software/netcdf/conventions.html

among other conventions.

It seems that the udunits developers don't have any plans to include
these extra calendars:

  http://www.unidata.ucar.edu/support/help/MailArchives/netcdf/msg02262.html

Cheers,
Ryo


From Daniele.Gandini at ARPA.PIEMONTE.IT  Tue Jan 27 09:53:06 2009
From: Daniele.Gandini at ARPA.PIEMONTE.IT (Daniele Gandini)
Date: Tue, 27 Jan 2009 15:53:06 +0100
Subject: How to set xlabels in a date/time axis
Message-ID: <TUE.27.JAN.2009.155306.0100.>

Dear GrADS users,

I have a forecast file whose analysis starts at 00 UTC
of 27th January 2009 and consists of 41 time intervals 
ranged every 6 hours.

When I want to see the complete time series of a variable on a particular
grid point, i.e. setting

'set lat 45'
'set lon 8'
'set t 1 last'
'd var'

the x labels which GrADS shows are:
27Jan, 29Jan, 1Feb, 3Feb and 5Feb.

I want to see at least one label for each day, i.e. 27Jan, 28Jan, 29Jan, 30 Jan and so on.

I tried 'set xlint' but it doesn't apply on date/time axis.
'set xaxis 1 41 4' issues numbers not dates on x axis.
'set xlevs 27Jan 28Jan ...' doesn't work.

Is there any possibility?

Thanks to everyone in advance.

Daniele Gandini

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090127/5cf415d4/attachment.html 

From Daniele.Gandini at ARPA.PIEMONTE.IT  Tue Jan 27 10:25:29 2009
From: Daniele.Gandini at ARPA.PIEMONTE.IT (Daniele Gandini)
Date: Tue, 27 Jan 2009 16:25:29 +0100
Subject: date-time axis
Message-ID: <TUE.27.JAN.2009.162529.0100.>

Dear GrADS users,

I run into the same problem.
But when the time series crosses two months, i.e.
from 27 January to 6 February, how can I do?

It is a pity that 'set xlint' doesn't work on date/time axis (so it is in
GrADS 1.9,
I don't know on GrADS 2.0).

Daniele Gandini


----- Original Message -----
From: "Matt Reszka" <matt at ATMOSP.PHYSICS.UTORONTO.CA>
To: <GRADSUSR at LIST.CINECA.IT>
Sent: Monday, September 10, 2007 7:01 PM
Subject: Re: date-time axis


> You can try something like
>
>        'set xaxis 1 30 5'
>        'set xlab %2.0fJAN'
>
> Of course xlab will have to be redefined for other months.
>
> Cheers.
>
> Matt


From dasilva at ALUM.MIT.EDU  Tue Jan 27 12:17:07 2009
From: dasilva at ALUM.MIT.EDU (Arlindo da Silva)
Date: Tue, 27 Jan 2009 12:17:07 -0500
Subject: Lats4D new wiki page
Message-ID: <TUE.27.JAN.2009.121707.0500.>

All,

   Lats4D has now its own wiki page:

           http://opengrads.org/wiki/index.php?title=Lats4D

  where you will find download instructions and usage notes for both GrADS
v1.9 and v2.0.  A manual page can be found here:

          http://opengrads.org/doc/scripts/lats4d/

  You can also search the cookbooks for lats4d related recipes:

        http://cookbooks.opengrads.org/index.php?title=Table_of_Contents

   Enjoy,

        Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090127/2ac8a7b0/attachment.html 

From Ron.Goodson at EC.GC.CA  Tue Jan 27 18:12:27 2009
From: Ron.Goodson at EC.GC.CA (Goodson,Ron [Edm])
Date: Tue, 27 Jan 2009 16:12:27 -0700
Subject: suggestion for transparent colour in shaded-plots
In-Reply-To: A<45B51100-1885-4AA6-AEF7-E58CD32E75BE@cola.iges.org>
Message-ID: <TUE.27.JAN.2009.161227.0700.>

Just a suggestion -- any chance of adding the ability to designate a
colour number as transparent when doing "gxout shaded" plots?
 
I know you can set the lowest shade to the background colour so it
"looks" like nothing is being drawn (even though it is).  But that
doesn't help when you want to have two separately shaded parameters
(each between some limits so that neither cover the entire graphic)...
as drawing the 2nd parameter overwrite the first.
 
I currently do this by using maskout() to control what is valid data to
draw.  But that leaves big-jagged edges based upon the resolution of my
grid so is neither accurate, nor visually appealing (understandable as I
see maskout as being more for calculations/analysis than for display)
 
If a colour could be marked as transparent - then I could use that for
the colour of my lower/upper limits of my shaded area and be able to
smoothly overlay 2 different shaded areas.  Perhaps it would simply set
the operation to XOR or something like that?
 
ron
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090127/6c84e55c/attachment.html 

From matt.alonso at GMAIL.COM  Tue Jan 27 19:03:54 2009
From: matt.alonso at GMAIL.COM (Matt Alonso)
Date: Tue, 27 Jan 2009 19:03:54 -0500
Subject: suggestion for transparent colour in shaded-plots
In-Reply-To: <4A6AB38B55B49C44A22E021A83CBEDDB03412526@sr-pnr-exch3.prairie.int.ec.gc.ca>
Message-ID: <TUE.27.JAN.2009.190354.0500.>

Hi Ron,

I will second that suggestion.

In the meantime though what I typically end up doing is creating a rather
roundabout solution to the situation and making use of ImageMagick's (
http://www.imagemagick.org) transparency feature (
http://www.imagemagick.org/script/command-line-options.php#transparent).
You can draw the first variable and set the shading for the region you want
to be transparent to a specific RGB value and then print that image to a
file; clear the buffer and do the same for the other variable.  You can then
do something like the following:

mogrify -transparent rgb(<red value>,<green value>,<blue value>)
first_image.gif
mogrify -transparent rgb(<red value>,<green value>,<blue value>)
second_image.gif
composite atop first_image.gif second_image.gif combined_image.gif

I don't have any code samples in front of me but if anyone is in need of
them let me know and I can send them out.  Hope this helps.

Cheers,
Matt

On Tue, Jan 27, 2009 at 6:12 PM, Goodson,Ron [Edm] <Ron.Goodson at ec.gc.ca>wrote:

>  Just a suggestion -- any chance of adding the ability to designate a
> colour number as transparent when doing "gxout shaded" plots?
>
> I know you can set the lowest shade to the background colour so it "looks"
> like nothing is being drawn (even though it is).  But that doesn't help when
> you want to have two separately shaded parameters (each between some limits
> so that neither cover the entire graphic)... as drawing the
> 2nd parameter overwrite the first.
>
> I currently do this by using maskout() to control what is valid data to
> draw.  But that leaves big-jagged edges based upon the resolution of my grid
> so is neither accurate, nor visually appealing (understandable as I see
> maskout as being more for calculations/analysis than for display)
>
> If a colour could be marked as transparent - then I could use that for the
> colour of my lower/upper limits of my shaded area and be able to smoothly
> overlay 2 different shaded areas.  Perhaps it would simply set the operation
> to XOR or something like that?
>
> ron
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090127/c9502e5c/attachment.html 

From doty at COLA.IGES.ORG  Tue Jan 27 19:08:11 2009
From: doty at COLA.IGES.ORG (Brian Doty)
Date: Tue, 27 Jan 2009 19:08:11 -0500
Subject: suggestion for transparent colour in shaded-plots
In-Reply-To: <4A6AB38B55B49C44A22E021A83CBEDDB03412526@sr-pnr-exch3.prairie.int.ec.gc.ca>
Message-ID: <TUE.27.JAN.2009.190811.0500.>

I plan on having a solution for this in the next six months or so....
Brian

On Jan 27, 2009, at 6:12 PM, Goodson,Ron [Edm] wrote:

> Just a suggestion -- any chance of adding the ability to designate
> a colour number as transparent when doing "gxout shaded" plots?
>
> I know you can set the lowest shade to the background colour so it
> "looks" like nothing is being drawn (even though it is).  But that
> doesn't help when you want to have two separately shaded parameters
> (each between some limits so that neither cover the entire
> graphic)... as drawing the 2nd parameter overwrite the first.
>
> I currently do this by using maskout() to control what is valid
> data to draw.  But that leaves big-jagged edges based upon the
> resolution of my grid so is neither accurate, nor visually
> appealing (understandable as I see maskout as being more for
> calculations/analysis than for display)
>
> If a colour could be marked as transparent - then I could use that
> for the colour of my lower/upper limits of my shaded area and be
> able to smoothly overlay 2 different shaded areas.  Perhaps it
> would simply set the operation to XOR or something like that?
>
> ron


From cbutenhoff at YAHOO.COM  Wed Jan 28 21:33:18 2009
From: cbutenhoff at YAHOO.COM (Chris Butenhoff)
Date: Thu, 29 Jan 2009 03:33:18 +0100
Subject: Memory Allocation Error: grid storage
Message-ID: <THU.29.JAN.2009.033318.0100.>

Dear Arlindo:

Thanks for the response.

You might be right. My original thought was that since the file contains 1
byte char values, it would only use ~600 Mb of memory when read into Grads.
So Grads converts all input data into 4 byte floats? Which as you said, puts
me close to my memory limit. If that is correct, then there is no way around
this except to divide the input file into smaller chunks.

Thanks.


_Chris


My 600 Mb binary file is composed of 1 byte values. I use the tag "-1,40,1"
>> On my attempt to display a variable in Grads v1.9b4 I receive the following
>> message:
>>
>> >>>Memory Allocation Error:  grid storage
>> >>>Data Request Error:  Error for variable 'lc'
>>
>> The open file is about 600 Mb in size and has a lon/lat grid of
>> 40320x16353.
>>
>
>Let's do the math. The size of the file hardly matters as it may include
>compression. Grads v1.9b4 holds gridded data in 4 byte floats.  So, a
>singlevariable at your grid resolution would require:
>
>40320*16353*4/(1024*1024) ~ 2.5 GB
>
>So, you could be very close to your memory limit of 3GB. The particular
>format also matters; there is a known memory leak with some form of
>netcdf/hdf file.
>
>
>>
>> My machine has over 3 Gb of RAM so I assume the memory limit must be set in
>> Grads.
>>
>> Searching the archives I find the following info:
>
>
>> >>> "The limit of the cache size is (for now) hard-coded in gaio around
>> >>> line 76:
>> >>> #define MAXG2CACHE 500100100
>> >>> You can change this and recompile if you want to limit or expand your
>> >>> memory usage."
>>
>> Before attempting, does anyone have experience with this problem, and is
>> this the right fix?
>>
>
> This is a Grib-2 specific issue, certainly not an issue with v1.9b4 which
>does not support Grib-2.
>
>--
>Arlindo da Silva
>dasilva at alum.mit.edu
>


From Eric.Wolvovsky at NOAA.GOV  Thu Jan 29 09:25:44 2009
From: Eric.Wolvovsky at NOAA.GOV (Eric Wolvovsky)
Date: Thu, 29 Jan 2009 09:25:44 -0500
Subject: Windows GrADS
In-Reply-To: <20090129023416.F0D1920332@mx2.cineca.it>
Message-ID: <THU.29.JAN.2009.092544.0500.>

I have the attached GrADS files running daily in a batch job. For years
they were ok, this morning however they began failing, as do files from
earlier days that had once worked (like the sample file attached).  I
haven't changed anything within the files, nor the batch job. Any
thoughts or suggestions would be greatly appreciated as I cannot seem to
narrow down the problem.

The version of GrADS I am running is: GrADS1.8sl9.win32

Thank you for your help,

Eric
-------------- next part --------------
A non-text attachment was scrubbed...
Name: grads_problem_20090129.zip
Type: application/x-zip-compressed
Size: 1804248 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090129/ebb9c2a0/attachment.bin 

From arlindo.dasilva at GMAIL.COM  Thu Jan 29 12:14:24 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Thu, 29 Jan 2009 12:14:24 -0500
Subject: Windows GrADS
In-Reply-To: <4981BC68.4020407@noaa.gov>
Message-ID: <THU.29.JAN.2009.121424.0500.>

On Thu, Jan 29, 2009 at 9:25 AM, Eric Wolvovsky <Eric.Wolvovsky at noaa.gov>wrote:

> I have the attached GrADS files running daily in a batch job. For years
> they were ok, this morning however they began failing, as do files from
> earlier days that had once worked (like the sample file attached).  I
> haven't changed anything within the files, nor the batch job. Any
> thoughts or suggestions would be greatly appreciated as I cannot seem to
> narrow down the problem.
>
> The version of GrADS I am running is: GrADS1.8sl9.win32
>

Recent Windows updates have started to break these 10 year old builds of
Win32 GrADS. You best bet is to upgrade to either v2.0 (if it has all the
features you need) or v1.9.-rc1, see for example:

*v1.9*:
http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows
*v2.0*:
http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows

These are the unofficial OpenGrADS builds; I suggest you install the Win32
superpack. You can also obtain an official Win32 build for v2.0 from COLA:

*v2.0*: http://grads.iges.org/grads/downloads.html

    Good Luck,

           Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090129/40222337/attachment.html 

From Eric.Wolvovsky at NOAA.GOV  Thu Jan 29 13:37:51 2009
From: Eric.Wolvovsky at NOAA.GOV (Eric Wolvovsky)
Date: Thu, 29 Jan 2009 13:37:51 -0500
Subject: Windows GrADS
In-Reply-To: <77fcd6b20901290914m1918a244n42b26ded95581b64@mail.gmail.com>
Message-ID: <THU.29.JAN.2009.133751.0500.>

Arlindo,

I went ahead and uninstalled my 1.8 version, installed the 2.0 version
(the one from COLA) and ended up with the same results. A blank output.
Any other suggestions?

Eric


Arlindo da Silva wrote:
> On Thu, Jan 29, 2009 at 9:25 AM, Eric Wolvovsky
> <Eric.Wolvovsky at noaa.gov <mailto:Eric.Wolvovsky at noaa.gov>> wrote:
>
>     I have the attached GrADS files running daily in a batch job. For
>     years
>     they were ok, this morning however they began failing, as do files
>     from
>     earlier days that had once worked (like the sample file attached).  I
>     haven't changed anything within the files, nor the batch job. Any
>     thoughts or suggestions would be greatly appreciated as I cannot
>     seem to
>     narrow down the problem.
>
>     The version of GrADS I am running is: GrADS1.8sl9.win32
>
>
> Recent Windows updates have started to break these 10 year old builds
> of Win32 GrADS. You best bet is to upgrade to either v2.0 (if it has
> all the features you need) or v1.9.-rc1, see for example:
>
> *v1.9*:
> http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows
> *v2.0*:
> http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows
>
> These are the unofficial OpenGrADS builds; I suggest you install the
> Win32 superpack. You can also obtain an official Win32 build for v2.0
> from COLA:
>
> *v2.0*: http://grads.iges.org/grads/downloads.html
>
>     Good Luck,
>
>            Arlindo
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu <mailto:dasilva at alum.mit.edu>


--
=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=
Eric Wolvovsky
5200 Auth Road
Camp Springs, Maryland, USA 20746
PH: 301-763-8000 x7752
FX: 301-763-8125
eric.wolvovsky at noaa.gov
=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=


From mhess at OCEANI.COM  Thu Jan 29 17:52:47 2009
From: mhess at OCEANI.COM (Mark Hess)
Date: Thu, 29 Jan 2009 15:52:47 -0700
Subject: Grads with Windows Vista
Message-ID: <THU.29.JAN.2009.155247.0700.>

I installed GrADS 1.8 on my new Vista machine.  When I tried to run the
program it did not open the normal GrADS window, but just froze.  No error
message, nothing.  So I installed V2.0 from the web site
http://www.iges.org/grads/downloads.html (unzipped files from
grads-2.0.a3-cygwin-dap.zip into the C:\Program Files\PCGrADS\win32e folder
replacing the existing files).  Now I get two errors.

1) Windows brings up a User Account Control (UAC) window forcing me to click
on "Allow" in order to run the program.  On my old system GrADS was set up
to run automatically, so this will not work.
grads
2) Before GrADS has a chance to open I get an error window which reads:
"This application has failed to start because cygXext-6.dll was not found.
Re-installing the application may fis this probloem." The Windows "looks for
a solution" but fails to find one and closes.  The V2.0 .zip file does not
contain this file ( as far as I can tell).  Do I need to install via a
different method or can I find this file somewhere?

Thanks for your help.

Mark


From arlindo.dasilva at GMAIL.COM  Thu Jan 29 21:17:44 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Thu, 29 Jan 2009 21:17:44 -0500
Subject: Grads with Windows Vista
In-Reply-To: <0FD32F877EE24993A1714DA74298D36F@sailfish>
Message-ID: <THU.29.JAN.2009.211744.0500.>

On Thu, Jan 29, 2009 at 5:52 PM, Mark Hess <mhess at oceani.com> wrote:

> I installed GrADS 1.8 on my new Vista machine.


I as mentioned on an earlier posting, recent WIndows updates have been
reported to break these old v1.8 Win32 builds.


> When I tried to run the
> program it did not open the normal GrADS window, but just froze.  No error
> message, nothing.  So I installed V2.0 from the web site
> http://www.iges.org/grads/downloads.html (unzipped files from
> grads-2.0.a3-cygwin-dap.zip into the C:\Program Files\PCGrADS\win32e folder
> replacing the existing files).  Now I get two errors.
>
> 1) Windows brings up a User Account Control (UAC) window forcing me to
> click
> on "Allow" in order to run the program.  On my old system GrADS was set up
> to run automatically, so this will not work.
> grads
> 2) Before GrADS has a chance to open I get an error window which reads:
> "This application has failed to start because cygXext-6.dll was not found.
> Re-installing the application may fis this probloem." The Windows "looks
> for
> a solution" but fails to find one and closes.  The V2.0 .zip file does not
> contain this file ( as far as I can tell).  Do I need to install via a
> different method or can I find this file somewhere?
>

I'd suggest trying one of the OpenGrADS Win32 superpacks from either:

*v1.9*:
http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows
*
*
*v2.0*:
http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows

These are complete packages that should run out of the box; see links above
for more details. They are fairly easy to install and are relocatable (you
can copy it to a USB mem stick and run it from there, even on another
machine with no installation necessary.)  Although these are not supported
by COLA, v2.0.a4.oga.2 should have all the functionality that you find in
COLA's v2.0.a4. These are the continuation of the Win32 builds I have been
doing for the last 10 years.

    Arlindo

--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090129/b405f20b/attachment.html 

From arlindo.dasilva at GMAIL.COM  Thu Jan 29 21:57:00 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Thu, 29 Jan 2009 21:57:00 -0500
Subject: Windows GrADS
In-Reply-To: <4981F77F.9050503@noaa.gov>
Message-ID: <THU.29.JAN.2009.215700.0500.>

On Thu, Jan 29, 2009 at 1:37 PM, Eric Wolvovsky <Eric.Wolvovsky at noaa.gov>wrote:

> Arlindo,
>
> I went ahead and uninstalled my 1.8 version, installed the 2.0 version
> (the one from COLA) and ended up with the same results. A blank output.
> Any other suggestions?
>

I had no problem reading your sample data with the OpenGrADS build:

*v2.0*:
http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows

See attached plot. However, I had to alter the DSET record in the ctl you
supplied because it pointed to a file that did not exist. Could that be the
cause of your problem? I am attaching the revised ctl.

   Arlindo



>
> Eric
>
>
> Arlindo da Silva wrote:
>
>> On Thu, Jan 29, 2009 at 9:25 AM, Eric Wolvovsky
>> <Eric.Wolvovsky at noaa.gov <mailto:Eric.Wolvovsky at noaa.gov>> wrote:
>>
>>    I have the attached GrADS files running daily in a batch job. For
>>    years
>>    they were ok, this morning however they began failing, as do files
>>    from
>>    earlier days that had once worked (like the sample file attached).  I
>>    haven't changed anything within the files, nor the batch job. Any
>>    thoughts or suggestions would be greatly appreciated as I cannot
>>    seem to
>>    narrow down the problem.
>>
>>    The version of GrADS I am running is: GrADS1.8sl9.win32
>>
>>
>> Recent Windows updates have started to break these 10 year old builds
>> of Win32 GrADS. You best bet is to upgrade to either v2.0 (if it has
>> all the features you need) or v1.9.-rc1, see for example:
>>
>> *v1.9*:
>>
>> http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows
>> *v2.0*:
>>
>> http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows
>>
>> These are the unofficial OpenGrADS builds; I suggest you install the
>> Win32 superpack. You can also obtain an official Win32 build for v2.0
>> from COLA:
>>
>> *v2.0*: http://grads.iges.org/grads/downloads.html
>>
>>    Good Luck,
>>
>>           Arlindo
>>
>> --
>> Arlindo da Silva
>> dasilva at alum.mit.edu <mailto:dasilva at alum.mit.edu>
>>
>
>
> --
> =+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=
> Eric Wolvovsky
> 5200 Auth Road
> Camp Springs, Maryland, USA 20746
> PH: 301-763-8000 x7752
> FX: 301-763-8125
> eric.wolvovsky at noaa.gov
> =+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=
>



--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090129/ce0d3636/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cpc_rfe2_print_0511.ctl
Type: application/octet-stream
Size: 711 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090129/ce0d3636/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pmer3.png
Type: image/png
Size: 14591 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090129/ce0d3636/attachment.png 

From Eric.Wolvovsky at NOAA.GOV  Fri Jan 30 07:23:07 2009
From: Eric.Wolvovsky at NOAA.GOV (Eric Wolvovsky)
Date: Fri, 30 Jan 2009 07:23:07 -0500
Subject: Windows GrADS
In-Reply-To: <77fcd6b20901291856o43b5d436xb8e50397f18bf5b@mail.gmail.com>
Message-ID: <FRI.30.JAN.2009.072307.0500.>

No, that isn't the issue, but thank you. I have it now.

Eric

Arlindo da Silva wrote:
> On Thu, Jan 29, 2009 at 1:37 PM, Eric Wolvovsky
> <Eric.Wolvovsky at noaa.gov <mailto:Eric.Wolvovsky at noaa.gov>> wrote:
>
>     Arlindo,
>
>     I went ahead and uninstalled my 1.8 version, installed the 2.0 version
>     (the one from COLA) and ended up with the same results. A blank
>     output.
>     Any other suggestions?
>
>
> I had no problem reading your sample data with the OpenGrADS build:
>
> *v2.0*:
> http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows
>
> See attached plot. However, I had to alter the DSET record in the ctl
> you supplied because it pointed to a file that did not exist. Could
> that be the cause of your problem? I am attaching the revised ctl.
>
>    Arlindo
>
>
>
>
>     Eric
>
>
>     Arlindo da Silva wrote:
>
>         On Thu, Jan 29, 2009 at 9:25 AM, Eric Wolvovsky
>         <Eric.Wolvovsky at noaa.gov <mailto:Eric.Wolvovsky at noaa.gov>
>         <mailto:Eric.Wolvovsky at noaa.gov
>         <mailto:Eric.Wolvovsky at noaa.gov>>> wrote:
>
>            I have the attached GrADS files running daily in a batch
>         job. For
>            years
>            they were ok, this morning however they began failing, as
>         do files
>            from
>            earlier days that had once worked (like the sample file
>         attached).  I
>            haven't changed anything within the files, nor the batch
>         job. Any
>            thoughts or suggestions would be greatly appreciated as I
>         cannot
>            seem to
>            narrow down the problem.
>
>            The version of GrADS I am running is: GrADS1.8sl9.win32
>
>
>         Recent Windows updates have started to break these 10 year old
>         builds
>         of Win32 GrADS. You best bet is to upgrade to either v2.0 (if
>         it has
>         all the features you need) or v1.9.-rc1, see for example:
>
>         *v1.9*:
>         http://opengrads.org/wiki/index.php?title=Installing_GrADS_on_Microsoft_Windows
>         *v2.0*:
>         http://opengrads.org/wiki/index.php?title=Installing_GrADS_v2.0_on_Microsoft_Windows
>
>         These are the unofficial OpenGrADS builds; I suggest you
>         install the
>         Win32 superpack. You can also obtain an official Win32 build
>         for v2.0
>         from COLA:
>
>         *v2.0*: http://grads.iges.org/grads/downloads.html
>
>            Good Luck,
>
>                   Arlindo
>
>         --
>         Arlindo da Silva
>         dasilva at alum.mit.edu <mailto:dasilva at alum.mit.edu>
>         <mailto:dasilva at alum.mit.edu <mailto:dasilva at alum.mit.edu>>
>
>
>
>     --
>     =+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=
>     Eric Wolvovsky
>     5200 Auth Road
>     Camp Springs, Maryland, USA 20746
>     PH: 301-763-8000 x7752
>     FX: 301-763-8125
>     eric.wolvovsky at noaa.gov <mailto:eric.wolvovsky at noaa.gov>
>     =+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=
>
>
>
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu <mailto:dasilva at alum.mit.edu>
>
> ------------------------------------------------------------------------
>


--
=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=
Eric Wolvovsky
5200 Auth Road
Camp Springs, Maryland, USA 20746
PH: 301-763-8000 x7752
FX: 301-763-8125
eric.wolvovsky at noaa.gov
=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=+=


From nieves at LTP.UPR.CLU.EDU  Fri Jan 30 07:26:50 2009
From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves)
Date: Fri, 30 Jan 2009 08:26:50 -0400
Subject: build for freebsd-7/amd64
Message-ID: <FRI.30.JAN.2009.082650.0400.>

Hello

Just to inform that the

     GrADS v2.0.a4.oga.1 build for FreeBSD 7.0/AMD64

is available at

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032&release_id=651493

while supplibs-2.0.1 can be downloaded from:

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=241681

Jose


From nieves at LTP.UPR.CLU.EDU  Fri Jan 30 07:51:24 2009
From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves)
Date: Fri, 30 Jan 2009 08:51:24 -0400
Subject: build freebsd-7.0/amd64
Message-ID: <FRI.30.JAN.2009.085124.0400.>

Hello

Just to let you know that that the

     GrADS v2.0.a4.oga.1 build for FreeBSD 7.0/AMD64

is available at

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=305032&release_id=651493

while supplibs-2.0.1 can be downloaded from:

https://sourceforge.net/project/showfiles.php?group_id=161773&package_id=241681

Jose


From axl419 at YAHOO.COM  Fri Jan 30 09:28:10 2009
From: axl419 at YAHOO.COM (See Hai Ooi)
Date: Fri, 30 Jan 2009 06:28:10 -0800
Subject: Error in creating station data using FORTRAN
Message-ID: <FRI.30.JAN.2009.062810.0800.>

Hi,

     I used Digital Visual FORTRAN (window version) to write a station data set in GrADS format.  On using the stnmap utility from GrADS v1.9.-rc1 (window) to create the station map file, the following error appeared:-

-------------------------------------------------------------

D:\StnPlot>stnmap -i aug2003.ctl  [under DOS command]

  Name of binary data set: d:/stnplot/aug2003.dat
  Number of times in the data set: 1
  Number of surface variables: 4
  Number of level dependent variables: 0

Starting scan of station data binary file.
Binary data file open: d:/stnplot/aug2003.dat
Processing time = 1
  Invalid station hdr found in station binary file
  Possible causes:  Invalid level count in hdr
                    Descriptor file mismatch
                    File not station data
                    Invalid relative time
    levs = 0  flag = 1  time = 166.98

-----------------------------------------------------------------

     I can't detect my mistakes.  Please help.

     Thank you.


OOI See-hai


      Get your new Email address!
Grab the Email name you&#39;ve always wanted before someone else does!
http://mail.promotions.yahoo.com/newdomains/aa/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aug2003.ctl
Type: application/octet-stream
Size: 329 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090130/a3420c42/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aug2003.dat
Type: application/octet-stream
Size: 1236 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090130/a3420c42/attachment-0001.obj 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: aug2003.txt
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090130/a3420c42/attachment.txt 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: stndata.for
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090130/a3420c42/attachment.pl 

From axl419 at YAHOO.COM  Fri Jan 30 09:51:27 2009
From: axl419 at YAHOO.COM (See Hai Ooi)
Date: Fri, 30 Jan 2009 06:51:27 -0800
Subject: Programming/Structure Error in creating station data using FORTRAN
Message-ID: <FRI.30.JAN.2009.065127.0800.>

Hi,

    My message with the original title was rejected by the server.  I therefore altered my subject title with the hope of reaching someone for help in identifying the mistakes made in my work.

    Thank you.

Ooi


--- On Fri, 30/1/09, See Hai Ooi <axl419 at yahoo.com> wrote:

> From: See Hai Ooi <axl419 at yahoo.com>
> Subject: Error in creating station data using FORTRAN
> To: GRADSUSR at LIST.CINECA.IT
> Date: Friday, 30 January, 2009, 10:28 PM
> Hi,
> 
>      I used Digital Visual FORTRAN (window version) to
> write a station data set in GrADS format.  On using the
> stnmap utility from GrADS v1.9.-rc1 (window) to create the
> station map file, the following error appeared:-
> 
> -------------------------------------------------------------
> 
> D:\StnPlot>stnmap -i aug2003.ctl  [under DOS
> command]
> 
>   Name of binary data set: d:/stnplot/aug2003.dat
>   Number of times in the data set: 1
>   Number of surface variables: 4
>   Number of level dependent variables: 0
> 
> Starting scan of station data binary file.
> Binary data file open: d:/stnplot/aug2003.dat
> Processing time = 1
>   Invalid station hdr found in station binary file
>   Possible causes:  Invalid level count in hdr
>                     Descriptor file mismatch
>                     File not station data
>                     Invalid relative time
>     levs = 0  flag = 1  time = 166.98
> 
> -----------------------------------------------------------------
> 
>      I can't detect my mistakes.  Please help.
> 
>      Thank you.
> 
> 
> OOI See-hai


      Get your preferred Email name!
Now you can @ymail.com and @rocketmail.com. 
http://mail.promotions.yahoo.com/newdomains/aa/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aug2003.ctl
Type: application/octet-stream
Size: 329 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090130/11affb03/attachment.obj 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aug2003.dat
Type: application/octet-stream
Size: 1236 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20090130/11affb03/attachment-0001.obj 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: aug2003.txt
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090130/11affb03/attachment.txt 
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: stndata.for
Url: http://gradsusr.org/pipermail/gradsusr/attachments/20090130/11affb03/attachment.pl 

From meteo at SPORT123.IT  Fri Jan 30 11:04:37 2009
From: meteo at SPORT123.IT (Valerio Guerra)
Date: Fri, 30 Jan 2009 17:04:37 +0100
Subject: build for freebsd-7/amd64
Message-ID: <FRI.30.JAN.2009.170437.0100.>

Jose F. Nieves wrote:
> Just to inform that the
> GrADS v2.0.a4.oga.1 build for FreeBSD 7.0/AMD64
> is available at

Thank you very much for your build.
I tried it for a little while and it seems ok.

No more problems like:
> gradsdap -bl
> /libexec/ld-elf.so.1: /usr/X11R6/lib/libXext.so.6: unsupported file layout
or similar ones.

Anyway next week I'll try some more advanced tests.
In the meantime... thank you very much again!

Valerio Guerra


From wmjj47 at YAHOO.COM  Fri Jan 30 12:22:18 2009
From: wmjj47 at YAHOO.COM (wilson jones)
Date: Fri, 30 Jan 2009 09:22:18 -0800
Subject: Darkening boundary Lines
Message-ID: <FRI.30.JAN.2009.092218.0800.>

I would like to know if someone could help me in darkening/thickening my boundary lines on my map. The pictures I have are high resolution and concentrated in a small area so when I analyzed the pictures, the boundary lines are very hard to see. I tried thickening the contour lines but it did not work.


Thank You,
Wilson Jones
NOAA-ISET  Research Assistant

It's lack of faith that makes people afraid of meeting challenges, and I believed in myself.
- Muhammad Ali -


-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090130/1bcc2bb8/attachment.html 

From nieves at LTP.UPR.CLU.EDU  Fri Jan 30 13:55:58 2009
From: nieves at LTP.UPR.CLU.EDU (Jose F. Nieves)
Date: Fri, 30 Jan 2009 14:55:58 -0400
Subject: build for freebsd-7/amd64
In-Reply-To: <PC1870200901301704370066ef54f467@Poseidon>
Message-ID: <FRI.30.JAN.2009.145558.0400.>

> Thank you very much for your build.
> I tried it for a little while and it seems ok.
>
> No more problems like:
> > gradsdap -bl
> > /libexec/ld-elf.so.1: /usr/X11R6/lib/libXext.so.6: unsupported file lay=
> out
> or similar ones.
>
> Anyway next week I'll try some more advanced tests.
> In the meantime... thank you very much again!
>
> Valerio Guerra

You're welcome.

Jose


From mequitomz at YAHOO.COM.BR  Fri Jan 30 17:59:11 2009
From: mequitomz at YAHOO.COM.BR (Arlindo Meque)
Date: Fri, 30 Jan 2009 14:59:11 -0800
Subject: Res: Darkening boundary Lines
Message-ID: <FRI.30.JAN.2009.145911.0800.>


Hi wilson,

Try to use the "set map" command:

For example: set map 1 1 6


regards,

Arlindo


      Veja quais s?o os assuntos do momento no Yahoo! +Buscados
http://br.maisbuscados.yahoo.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090130/57a34782/attachment.html 

From scpiper at UCSD.EDU  Fri Jan 30 18:29:07 2009
From: scpiper at UCSD.EDU (Stephen Piper)
Date: Fri, 30 Jan 2009 15:29:07 -0800
Subject: Need help with grads ctl file for reading 3 dim instead of std 5
 dim array
Message-ID: <FRI.30.JAN.2009.152907.0800.>

I am trying to read a fortran binary file that is written with the
following Fortran code:

       integer, parameter                       :: statcnt=228,
       ^                                            tractot=26,tlen=2190
       real,dimension(statcnt,tractot,tlen)     :: statco2
     ie., 1 record with no. of stations, no. of tracers, and no. of days

Can I read this in grads by setting the dimensions of YDEF and ZDEF
to 1 with the
following ctl file?   At present, grads is coughing about the
variable definitions.

The values in the file are really just daily values for 26 different
tracers at 228 stations.
Current draft of my ctl file is below.

Thanks!
Steve


my ctl file - doesn't quite work:

DSET     /Users/stephenpiper/Transcom3/CSU.gurney.L1.statco2.le.dat
OPTIONS  sequential
OPTIONS  little_endian
OPTIONS   365_day_calendar

      * .le. in filename probably means little endian ...

TITLE     CSU gurney hi frequency station data
UNDEF    1.e36

   * fortran statements
   *    integer, parameter                       :: statcnt=228,
   *   ^                                            tractot=26,tlen=2190
   *    real,dimension(statcnt,tractot,tlen)     :: statco2

   *  fortran sequential file - 1 record of 228,26,2190
   *     ie, fortran binary with rec lengths in each record

XDEF     228 linear   1 1
    *  no. of items scaletype startcoord increment
YDEF     1 linear 1 1
ZDEF     1 levels 1000
TDEF     2190 linear 1jan91 1dy
VARS      26
           fos90  0 -99
      *  note that the format is standard x,y,z,var,t
           fos95  0 -99
           bioa   0 -99
           dum    0 -99
           bio1   0 -99
           bio2   0 -99
           bio3   0 -99
           bio4   0 -99
           bio5   0 -99
           bio6   0 -99
           bio7   0 -99
           bio8   0 -99
           bio9   0 -99
           bio10  0 -99
           bio11  0 -99
           oce1   0 -99
           oce2   0 -99
           oce3   0 -99
           oce4   0 -99
           oce5   0 -99
           oce6   0 -99
           oce7   0 -99
           oce8   0 -99
           oce9   0 -99
           oce10  0 -99
           oce11  0 -99
       ENDVARS


From Charles.Seman at NOAA.GOV  Fri Jan 30 20:24:00 2009
From: Charles.Seman at NOAA.GOV (Charles Seman)
Date: Fri, 30 Jan 2009 20:24:00 -0500
Subject: Need help with grads ctl file for reading 3 dim instead of
 std 5  dim array
In-Reply-To: <1091C0B5-9E2B-4371-81EC-878B7BC1A114@ucsd.edu>
Message-ID: <FRI.30.JAN.2009.202400.0500.>

Stephen,

Have you checked out:
http://grads.iges.org/grads/gadoc/aboutstationdata.html
contains a lot of info on station data, including a Fortran program to
write out data for GrADS, and a GrADS station data descriptor file...

Hope this helps,
Chuck

Stephen Piper wrote:
> I am trying to read a fortran binary file that is written with the
> following Fortran code:
>
>       integer, parameter                       :: statcnt=228,
>       ^                                            tractot=26,tlen=2190
>       real,dimension(statcnt,tractot,tlen)     :: statco2
>     ie., 1 record with no. of stations, no. of tracers, and no. of days
>
> Can I read this in grads by setting the dimensions of YDEF and ZDEF
> to 1 with the
> following ctl file?   At present, grads is coughing about the
> variable definitions.
>
> The values in the file are really just daily values for 26 different
> tracers at 228 stations.
> Current draft of my ctl file is below.
>
> Thanks!
> Steve
>
>
> my ctl file - doesn't quite work:
>
> DSET     /Users/stephenpiper/Transcom3/CSU.gurney.L1.statco2.le.dat
> OPTIONS  sequential
> OPTIONS  little_endian
> OPTIONS   365_day_calendar
>
>      * .le. in filename probably means little endian ...
>
> TITLE     CSU gurney hi frequency station data
> UNDEF    1.e36
>
>   * fortran statements
>   *    integer, parameter                       :: statcnt=228,
>   *   ^                                            tractot=26,tlen=2190
>   *    real,dimension(statcnt,tractot,tlen)     :: statco2
>
>   *  fortran sequential file - 1 record of 228,26,2190
>   *     ie, fortran binary with rec lengths in each record
>
> XDEF     228 linear   1 1
>    *  no. of items scaletype startcoord increment
> YDEF     1 linear 1 1
> ZDEF     1 levels 1000
> TDEF     2190 linear 1jan91 1dy
> VARS      26
>           fos90  0 -99
>      *  note that the format is standard x,y,z,var,t
>           fos95  0 -99
>           bioa   0 -99
>           dum    0 -99
>           bio1   0 -99
>           bio2   0 -99
>           bio3   0 -99
>           bio4   0 -99
>           bio5   0 -99
>           bio6   0 -99
>           bio7   0 -99
>           bio8   0 -99
>           bio9   0 -99
>           bio10  0 -99
>           bio11  0 -99
>           oce1   0 -99
>           oce2   0 -99
>           oce3   0 -99
>           oce4   0 -99
>           oce5   0 -99
>           oce6   0 -99
>           oce7   0 -99
>           oce8   0 -99
>           oce9   0 -99
>           oce10  0 -99
>           oce11  0 -99
>       ENDVARS

--

Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.

********************************************************************
 Charles Seman                                Charles.Seman at noaa.gov
 U.S. Department of Commerce / NOAA / OAR
 Geophysical Fluid Dynamics Laboratory         voice: (609) 452-6547
 201 Forrestal Road                              fax: (609) 987-5063
 Princeton, NJ  08540-6649            http://www.gfdl.noaa.gov/~cjs/
********************************************************************

"The contents of this message are mine personally and do not necessarily
reflect any position of the Government or NOAA."


From arlindo.dasilva at GMAIL.COM  Sat Jan 31 12:51:26 2009
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Sat, 31 Jan 2009 12:51:26 -0500
Subject: Need help with grads ctl file for reading 3 dim instead of
 std 5 dim array
In-Reply-To: <4983A830.5090207@noaa.gov>
Message-ID: <SAT.31.JAN.2009.125126.0500.>

On Fri, Jan 30, 2009 at 8:24 PM, Charles Seman <Charles.Seman at noaa.gov>wrote:

> Stephen,
>
> Have you checked out:
> http://grads.iges.org/grads/gadoc/aboutstationdata.html
> contains a lot of info on station data, including a Fortran program to
> write out data for GrADS, and a GrADS station data descriptor file...
>

There are also a couple of recipes about station data:

http://cookbooks.opengrads.org/index.php?title=Table_of_Contents#SEQUENTIAL_INDEX_OF_RECIPES

   Arlindo



--
Arlindo da Silva
dasilva at alum.mit.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20090131/118e737a/attachment.html