From zjuyanlibin at 163.COM  Mon Feb  1 02:14:49 2010
From: zjuyanlibin at 163.COM (=?utf-8?B?WWFuIExpYmlu?=)
Date: Mon, 1 Feb 2010 15:14:49 +0800
Subject: opengrads function requesting
Message-ID: <MON.1.FEB.2010.151449.0800.>

Hi, everyone,

Plz do me a favor.
Can anyone transfer the function file to me?
re.gex - A 2D regridding function for GrADS

Thank you in advance.

Best regards,
Libin Yan


2010-02-01 



Yan Libin 



???? M.Mahakur 
????? 2010-01-31  21:27:40 
???? Yan Libin 
??? 
??? Re: Different time indices 
 
Hi Libin,
Here it is..
Regards,
Mahakur
--
M. Mahakur
INDIAN INSTITUTE OF TROPICAL METEOROLOGY
PUNE - 411 008, INDIA 
---------- Original Message ---- Yan Libin <zjuyanlibin at 163.COM>-------
From:
To: GRADSUSR at list.cineca.it
Sent: Sun, 31 Jan 2010 09:07:53 +0800
Subject: Re: Different time indices
> http://ams.allenpress.com/perlserv/?request=get-
> pdf&doi=10.1175%2F2009JAMC2167.1
> 
> who can send me a copy?
> Tks a million.
> 
> Best regards,
> Libin Yan
> 
> 2010-01-31
> 
> Yan Libin
> 
> [UTF-8?]???????????? Piotr_Djak?w 
> [UTF-8?]??????????????? 2010-01-29  19:51:46 
> [UTF-8?]???????????? GRADSUSR 
> [UTF-8?]????????? 
> [UTF-8?]????????? Different time indices
> 
> Hello
> 
> I've try to open two GRIB files with different time indices. When i try to 
> make a plot :
> 
> d tmp2m.1 - tmp2m.2
> 
> I've got an error message:
> 
> "Data Request Error: Invalid grid coordinates
>   World coordinates convert to non-integer  grid coordinates"
> 
> Is there any possibility to display two files with different time indices?
> 
> Best Regards
> Piotr Djakow
------- End of Original Message -------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/3907180b/attachment.html 

From zjuyanlibin at 163.COM  Mon Feb  1 02:43:35 2010
From: zjuyanlibin at 163.COM (Yan Libin)
Date: Mon, 1 Feb 2010 15:43:35 +0800
Subject: How to start re function?
Message-ID: <MON.1.FEB.2010.154335.0800.>

Hi, every one,

I want to regrid files into one-standard-grid files.
and install "grads-2.0.a7.oga.3-win32_superpack.exe " under windows xp operation.

I  click the shortcut of opengrads and input "re" command,
but it pops up as ''no such command"
How can i resolve this problem?

Thanks in advance.

Regards,
Libin Yan


2010-02-01 



Yan Libin 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/9a27a90a/attachment.html 

From jithu9907 at GMAIL.COM  Mon Feb  1 06:32:33 2010
From: jithu9907 at GMAIL.COM (sreejith k.s)
Date: Mon, 1 Feb 2010 17:02:33 +0530
Subject: making discriptor file
Message-ID: <MON.1.FEB.2010.170233.0530.>

hai,
     can anyone tell me how discriptor file can be made of a data, which is
other than netcdf data.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/539bb0a8/attachment.html 

From pietrowicz at ZEDXINC.COM  Mon Feb  1 10:00:04 2010
From: pietrowicz at ZEDXINC.COM (Joseph Pietrowicz)
Date: Mon, 1 Feb 2010 16:00:04 +0100
Subject: using two files with different time ranges
Message-ID: <MON.1.FEB.2010.160004.0100.>

I'm using opengrads V 2.0.a7.oga.3 and opening two files, gfs and cmorph 3hr
precipitation estimates. I can open and use them one at a time just fine.
But when I open them in as file 1 and file 2, I can't display variables from
 both data sets. I think it has to do with the time coordinate.

The GFS ctl data set has time as 

tdef 9 linear 00Z27jan2010 3hr

While the cmorph time in the ctl file is given as

TDEF 99000 LINEAR  0z01dec2002 3hr

and the cmoph *ctl file uses the template option.

If I open both files, say the GFS first, I can display any of the GFS
variables for any of the GFS time steps just fine. However, I can't display
any of the cmorph data. I believe that it has something to do with setting time.

For example in the GFS
set t 1
d tmp2m

will give me the results properly.
In the cmorph file time step 20920 is the same time as time step 1 in the
GFS data set. Then I enter the command to display cmorph data.

d cpc.2

and I get a message that 'all undefined values'
And the time information for cmorph is correct, 'T= 20920 20920'

Also the GFS data set id grib2, while the cmorph is binary.

What do I need to do in order to work witht both data sets?

Thanks,

Joe Pietrowicz


From michael.l.sestak at NAVY.MIL  Mon Feb  1 10:47:43 2010
From: michael.l.sestak at NAVY.MIL (Sestak, Michael L CIV 63134)
Date: Mon, 1 Feb 2010 07:47:43 -0800
Subject: using two files with different time ranges
In-Reply-To: A<LISTSERV%201002011600049847.1C34@LIST.CINECA.IT>
Message-ID: <MON.1.FEB.2010.074743.0800.>

If only need to display variables from both data sets, rather than do calculations with both, use set dfile 1 and set dfile 2 before you plot each variable, as appropriate.  Grads uses the current default file to set the time range and increment.

Michael Sestak
Fleet Numerical Meteorology and Oceanography Center
Monterey, CA

-----Original Message-----
From: GRADSUSR at LIST.CINECA.IT [mailto:GRADSUSR at LIST.CINECA.IT] On Behalf Of Joseph Pietrowicz
Sent: Monday, February 01, 2010 7:00
To: GRADSUSR at LIST.CINECA.IT
Subject: using two files with different time ranges

I'm using opengrads V 2.0.a7.oga.3 and opening two files, gfs and cmorph 3hr precipitation estimates. I can open and use them one at a time just fine.
But when I open them in as file 1 and file 2, I can't display variables from  both data sets. I think it has to do with the time coordinate.

The GFS ctl data set has time as 

tdef 9 linear 00Z27jan2010 3hr

While the cmorph time in the ctl file is given as

TDEF 99000 LINEAR  0z01dec2002 3hr

and the cmoph *ctl file uses the template option.

If I open both files, say the GFS first, I can display any of the GFS variables for any of the GFS time steps just fine. However, I can't display any of the cmorph data. I believe that it has something to do with setting time.

For example in the GFS
set t 1
d tmp2m

will give me the results properly.
In the cmorph file time step 20920 is the same time as time step 1 in the GFS data set. Then I enter the command to display cmorph data.

d cpc.2

and I get a message that 'all undefined values'
And the time information for cmorph is correct, 'T= 20920 20920'

Also the GFS data set id grib2, while the cmorph is binary.

What do I need to do in order to work witht both data sets?

Thanks,

Joe Pietrowicz
-------------- next part --------------
A non-text attachment was scrubbed...
Name: smime.p7s
Type: application/x-pkcs7-signature
Size: 5212 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20100201/bf463033/attachment.bin 

From pietrowicz at ZEDXINC.COM  Mon Feb  1 12:27:27 2010
From: pietrowicz at ZEDXINC.COM (Joe Pietrowicz)
Date: Mon, 1 Feb 2010 12:27:27 -0500
Subject: using two files with different time ranges
In-Reply-To: <B500BFEE614AB548873A2C4727E6F20001124798@nawemuguez01v.nadsuswe.nads.navy.mil>
Message-ID: <MON.1.FEB.2010.122727.0500.>

Thanks, that is helpful.

Joe 
----- Original Message -----
From: "Michael L CIV 63134 Sestak" <michael.l.sestak at NAVY.MIL>
To: GRADSUSR at LIST.CINECA.IT
Sent: Monday, February 1, 2010 10:47:43 AM GMT -05:00 US/Canada Eastern
Subject: Re: using two files with different time ranges

If only need to display variables from both data sets, rather than do calculations with both, use set dfile 1 and set dfile 2 before you plot each variable, as appropriate.  Grads uses the current default file to set the time range and increment.

Michael Sestak
Fleet Numerical Meteorology and Oceanography Center
Monterey, CA

-----Original Message-----
From: GRADSUSR at LIST.CINECA.IT [mailto:GRADSUSR at LIST.CINECA.IT] On Behalf Of Joseph Pietrowicz
Sent: Monday, February 01, 2010 7:00
To: GRADSUSR at LIST.CINECA.IT
Subject: using two files with different time ranges

I'm using opengrads V 2.0.a7.oga.3 and opening two files, gfs and cmorph 3hr precipitation estimates. I can open and use them one at a time just fine.
But when I open them in as file 1 and file 2, I can't display variables from  both data sets. I think it has to do with the time coordinate.

The GFS ctl data set has time as 

tdef 9 linear 00Z27jan2010 3hr

While the cmorph time in the ctl file is given as

TDEF 99000 LINEAR  0z01dec2002 3hr

and the cmoph *ctl file uses the template option.

If I open both files, say the GFS first, I can display any of the GFS variables for any of the GFS time steps just fine. However, I can't display any of the cmorph data. I believe that it has something to do with setting time.

For example in the GFS
set t 1
d tmp2m

will give me the results properly.
In the cmorph file time step 20920 is the same time as time step 1 in the GFS data set. Then I enter the command to display cmorph data.

d cpc.2

and I get a message that 'all undefined values'
And the time information for cmorph is correct, 'T= 20920 20920'

Also the GFS data set id grib2, while the cmorph is binary.

What do I need to do in order to work witht both data sets?

Thanks,

Joe Pietrowicz


From frangix at YAHOO.IT  Mon Feb  1 13:35:43 2010
From: frangix at YAHOO.IT (mario frangipane)
Date: Mon, 1 Feb 2010 18:35:43 +0000
Subject: use of function cdiff
Message-ID: <MON.1.FEB.2010.183543.0000.>

Dear Grads Users,
I do not understand the example on using the cdiff function at page 76 of the Grads manual.
For example, if (u,v) is the wind, to calculate du/dy I believe that the calculation should be:
define du=cdiff(u,y)define dy=cdiff(lat,y)*(3.1416/180)*(6.37e6) ?(in meters)define dudy=du/dy
Instead, in the Manual, at pag. 76, I find (example 1):
define du=cdiff(u*cos(lat*3.1416/180),y).
Anyone can help me to understand?




      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/91003f44/attachment.html 

From jmssnyder at UCDAVIS.EDU  Mon Feb  1 15:38:52 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Mon, 1 Feb 2010 12:38:52 -0800
Subject: use of function cdiff
In-Reply-To: <936963.51217.qm@web23506.mail.ird.yahoo.com>
Message-ID: <MON.1.FEB.2010.123852.0800.>

the cos(lat*3.14/180),y) is used to convert degrees to meters.  I know that
for x which is longitude you have to account for the fact the length of a
degree longitude changes with latitude and that is accounted for when
calculating du.

-Jason

On Mon, Feb 1, 2010 at 10:35 AM, mario frangipane <frangix at yahoo.it> wrote:

> Dear Grads Users,
>
> I do not understand the example on using the cdiff function at page 76 of
> the Grads manual.
>
> For example, if (u,v) is the wind, to calculate du/dy I believe that the
> calculation should be:
>
> define du=cdiff(u,y)
> define dy=cdiff(lat,y)*(3.1416/180)*(6.37e6)  (in meters)
> define dudy=du/dy
>
> Instead, in the Manual, at pag. 76, I find (example 1):
>
> define du=cdiff(u*cos(lat*3.1416/180),y).
>
> Anyone can help me to understand?
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/f5bcf3f6/attachment.html 

From jmssnyder at UCDAVIS.EDU  Mon Feb  1 19:55:55 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Mon, 1 Feb 2010 16:55:55 -0800
Subject: Displaying Compound Variables
Message-ID: <MON.1.FEB.2010.165555.0800.>

I have a script to produce compound variables:

'sdfopen uwnd.2001.nc'
'set gxout shaded'
'set t 1'
'define U=uwnd'
i=1
count=1
While(i<5)
    count=count+1
    say count
   'set t 'i
   U.count=uwnd
    i=i+1
endwhile
the compound variable is U.count.  How do I go about producing and
displaying compound variables?
What modifications must I make to the script.

What I want to make a formula that would be like this
U(t)=U(t)+U(t-1)+U(t-2).  Where t-1 is a time step earlier than t and t-2 is
two timesteps earlier than t.  I want to have this formula work for
steptimes from 1 to say 5 (set t 1 5) How do I go about doing this?  Are
compound variables event the answer?

Thanks,

Jason
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/e92a8740/attachment.html 

From jdduda at IASTATE.EDU  Mon Feb  1 23:39:32 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Mon, 1 Feb 2010 22:39:32 -0600
Subject: Displaying Compound Variables
In-Reply-To: <fe8efa01002011655w7cfa839eq245f024acbd51815@mail.gmail.com>
Message-ID: <MON.1.FEB.2010.223932.0600.>

No.  Compounds only work for internal variables and you can't display them.
there's a difference between a displayable product/field and an internal
variable that is used in calculations.

If you want to make a field that is the sum of that field in previous times,
you'll have to code it like that, but you can do it with a while loop in
such a way:

'define product = product(first time)'
while (time step counter variable condition is true)
 product '=' product '-' product'(t-'time step counter variable')'
endwhile

'd 'product

Each time the code goes through the loop, the previous time step of the
field will be taken away from the current one, starting with the latest in
time.

Jeff Duda

On Mon, Feb 1, 2010 at 6:55 PM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> I have a script to produce compound variables:
>
> 'sdfopen uwnd.2001.nc'
> 'set gxout shaded'
> 'set t 1'
> 'define U=uwnd'
> i=1
> count=1
> While(i<5)
>     count=count+1
>     say count
>    'set t 'i
>    U.count=uwnd
>     i=i+1
> endwhile
> the compound variable is U.count.  How do I go about producing and
> displaying compound variables?
> What modifications must I make to the script.
>
> What I want to make a formula that would be like this
> U(t)=U(t)+U(t-1)+U(t-2).  Where t-1 is a time step earlier than t and t-2 is
> two timesteps earlier than t.  I want to have this formula work for
> steptimes from 1 to say 5 (set t 1 5) How do I go about doing this?  Are
> compound variables event the answer?
>
> Thanks,
>
> Jason
>
>
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/cc1fbef9/attachment.html 

From jmssnyder at UCDAVIS.EDU  Tue Feb  2 02:40:51 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Mon, 1 Feb 2010 23:40:51 -0800
Subject: Displaying Compound Variables
In-Reply-To: <c0f1d7891002012039t53b457b0h78ce991c64f30169@mail.gmail.com>
Message-ID: <MON.1.FEB.2010.234051.0800.>

I tried the following with no success:

'sdfopen uwnd.2001.nc'
'sdfopen vwnd.2001.nc'
'set gxout shaded'
'set t 1'
'define uwnd=uwnd(t=1)'
i=2
count=1
While(i<5)
    count=count+1
    say count
   'set t 'i
   uwnd = uwnd - uwnd(i-1)
*    V.count=vwnd.2
*   'd U.count'
    i=i+1
endwhile
what am I doing wrong here?

Thanks,

Jason
On Mon, Feb 1, 2010 at 8:39 PM, Jeffrey Duda <jdduda at iastate.edu> wrote:

> No.  Compounds only work for internal variables and you can't display
> them.  there's a difference between a displayable product/field and an
> internal variable that is used in calculations.
>
> If you want to make a field that is the sum of that field in previous
> times, you'll have to code it like that, but you can do it with a while loop
> in such a way:
>
> 'define product = product(first time)'
> while (time step counter variable condition is true)
>  product '=' product '-' product'(t-'time step counter variable')'
> endwhile
>
> 'd 'product
>
> Each time the code goes through the loop, the previous time step of the
> field will be taken away from the current one, starting with the latest in
> time.
>
> Jeff Duda
>
>
> On Mon, Feb 1, 2010 at 6:55 PM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:
>
>> I have a script to produce compound variables:
>>
>> 'sdfopen uwnd.2001.nc'
>> 'set gxout shaded'
>> 'set t 1'
>> 'define U=uwnd'
>> i=1
>> count=1
>> While(i<5)
>>     count=count+1
>>     say count
>>    'set t 'i
>>    U.count=uwnd
>>     i=i+1
>> endwhile
>> the compound variable is U.count.  How do I go about producing and
>> displaying compound variables?
>> What modifications must I make to the script.
>>
>> What I want to make a formula that would be like this
>> U(t)=U(t)+U(t-1)+U(t-2).  Where t-1 is a time step earlier than t and t-2 is
>> two timesteps earlier than t.  I want to have this formula work for
>> steptimes from 1 to say 5 (set t 1 5) How do I go about doing this?  Are
>> compound variables event the answer?
>>
>> Thanks,
>>
>> Jason
>>
>>
>>
>
>
>
> --
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100201/7bd8f08a/attachment.html 

From dsouza at CPTEC.INPE.BR  Tue Feb  2 07:02:24 2010
From: dsouza at CPTEC.INPE.BR (Diego Souza)
Date: Tue, 2 Feb 2010 10:02:24 -0200
Subject: Question
Message-ID: <TUE.2.FEB.2010.100224.0200.>

Hi,

I am using MERRA datasets for make a climatology (MERRA Hourly data)
When I use the command: define
prec=ave(prectot,time=00Z01JAN1999,time=23Z31JAN1999) , the command line
show this error:

gancgrid error: nc_get_vara_double failed; Unknown error 1001
Data Request Error:  Error for variable 'prectot'
Error from AVE:  Error getting grids
Operation Error:  Error from ave function
  Error ocurred at column 1
DEFINE error:  Invalid expression.

Someone can esplain this error?

Thanks,
_________________________________________

Diego Oliveira de Souza
Grupo de Estudos Ambientais Integrados (GEAI)
Instituto Nacional de Pesquisas Espaciais - INPE
Cachoeira Paulista - SP
Fone: (12) 31869512
Cel: (12) 91219668
_________________________________________
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100202/ae48b29a/attachment.html 

From gzepka at DGE.INPE.BR  Tue Feb  2 08:26:24 2010
From: gzepka at DGE.INPE.BR (Gisele Zepka)
Date: Tue, 2 Feb 2010 11:26:24 -0200
Subject: difference between two files
Message-ID: <TUE.2.FEB.2010.112624.0200.>

Hello users!!

I have two files with different resolutions and I want to compare the 
differences between them. I used the regrid function to take both in the 
same resolution. The regrided file is just one hour, while in the other 
I have a range of 22 times. I want to compare just the hour of the 
first, so I have two problems when I try to do that: an error of 
resolution, because the two files don't have the same quantity of points 
in x and y and I don't know how to set the time just for the second and 
then calculate the difference.

Please any help will be nice!

Thanks,
Gisele


-- 
Msc. Gisele dos Santos Zepka
Atmospheric Electricity Group (ELAT)
National Institute for Space Research (INPE)
Av. dos Astronautas, 1758
S?o Jos? dos Campos/ SP/ Brazil CEP 12227-010
Phone ++55(12)39456841


From paulopomatos at UOL.COM.BR  Tue Feb  2 08:36:20 2010
From: paulopomatos at UOL.COM.BR (Paulo)
Date: Tue, 2 Feb 2010 11:36:20 -0200
Subject: difference between two files
In-Reply-To: <4B682800.5090802@dge.inpe.br>
Message-ID: <TUE.2.FEB.2010.113620.0200.>

An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100202/ec5ef8a7/attachment.html 

From gzepka at DGE.INPE.BR  Tue Feb  2 08:47:52 2010
From: gzepka at DGE.INPE.BR (Gisele Zepka)
Date: Tue, 2 Feb 2010 11:47:52 -0200
Subject: difference between two files
In-Reply-To: <4b682a549459e_2147bbace70191@weasel18.tmail>
Message-ID: <TUE.2.FEB.2010.114752.0200.>

.grb

Gisele



Em 2/2/2010 11:36, Paulo escreveu:
>
> Hi Gisele,
>
> You may try the following:
>
> 1 - Cut both files in lat/lon (same lat/lon for both) and if they have 
> the same spatial resolution (if the regrid function worked ok) they 
> will have the same quantities of points in x and y.
>
> What kind of data are you using? Netcdf, binary or any other?
>
> Paulo.
>
> /?Se algu?m ama uma flor da qual s? existe um exemplar em milh?es de 
> estrelas, isso basta para que seja feliz quando a contempla." (Antoine 
> de Saint-Exup?ry)/
>
>
>
> Em 02/02/2010 11:26, *Gisele Zepka < http://biologias.com/ >* escreveu:
>
>
>     Hello users!!
>
>     I have two files with different resolutions and I want to compare the
>     differences between them. I used the regrid function to take both
>     in the
>     same resolution. The regrided file is just one hour, while in the
>     other
>     I have a range of 22 times. I want to compare just the hour of the
>     first, so I have two problems when I try to do that: an error of
>     resolution, because the two files don't have the same quantity of
>     points
>     in x and y and I don't know how to set the time just for the
>     second and
>     then calculate the difference.
>
>     Please any help will be nice!
>
>     Thanks,
>     Gisele
>
>
>     -- 
>     Msc. Gisele dos Santos Zepka
>     Atmospheric Electricity Group (ELAT)
>     National Institute for Space Research (INPE)
>     Av. dos Astronautas, 1758
>     S?o Jos? dos Campos/ SP/ Brazil CEP 122 27-010
>     Phone ++55(12)39456841
>


-- 
Msc. Gisele dos Santos Zepka
Atmospheric Electricity Group (ELAT)
National Institute for Space Research (INPE)
Av. dos Astronautas, 1758
S?o Jos? dos Campos/ SP/ Brazil CEP 12227-010
Phone ++55(12)39456841

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

From jdduda at IASTATE.EDU  Tue Feb  2 10:51:53 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Tue, 2 Feb 2010 09:51:53 -0600
Subject: Displaying Compound Variables
In-Reply-To: <fe8efa01002012340t58e4dc55n6cd3e947f6885674@mail.gmail.com>
Message-ID: <TUE.2.FEB.2010.095153.0600.>

On Tue, Feb 2, 2010 at 1:40 AM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> I tried the following with no success:
>
> 'sdfopen uwnd.2001.nc'
>
> 'sdfopen vwnd.2001.nc'
> 'set gxout shaded'
> 'set t 1'
> 'define uwnd=uwnd(t=1)'
> i=2
>
> count=1
> While(i<5)
>     count=count+1
>     say count
>    'set t 'i
>    uwnd = uwnd - uwnd(i-1)
> *    V.count=vwnd.2
> *   'd U.count'
>     i=i+1
> endwhile
> what am I doing wrong here?
>
> Thanks,
>
> Jason
>
>
You have to be very picky with your quotes.  If any part of any line in your
script would be something you would type in at the console (in Grads) to do
outside of a script, it needs to be in quotes.  That means any define,
display, or set plotting environment commands.  You need to change the line

uwnd = uwnd - uwnd(i-1)    to

'uwnd = uwnd - uwnd(t-'i')'     ("i" needs to be outside the quotes as you
see it.  That way the Grads script will sub in an integer value for "i" and
you'll get your t-1, t-2, etc)

-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100202/5b01c474/attachment.html 

From msponsler at COMCAST.NET  Tue Feb  2 16:22:37 2010
From: msponsler at COMCAST.NET (Mark Sponsler)
Date: Tue, 2 Feb 2010 21:22:37 +0000
Subject: WIN GRADS 2.0.a7.oga.3 - Possible Bug
In-Reply-To: <77fcd6b21001272059j539b278xee944ae919e50063@mail.gmail.com>
Message-ID: <TUE.2.FEB.2010.212237.0000.>



Hi Arlindo, 

So I edited the path and removed all instances of grads and wgrib.? I shut the machine down cold, then?booted it.??I then plugged in the flash drive with Grads 2 on it and tried '! ls <cr>'.?? It hung.? 



I'm thinking maybe I have another application in the path that might be causing a problem (GMT, maybe Perl).? Tonight I'll try pulling them out of the path too and see if it make a difference.? 


Thanks, 
Mark 

----- Original Message ----- 
From: "Arlindo da Silva" <dasilva at alum.mit.edu> 
To: "Mark Sponsler" <msponsler at comcast.net> 
Cc: "gradsusr" <GRADSUSR at list.cineca.it> 
Sent: Wednesday, January 27, 2010 8:59:01 PM GMT -08:00 US/Canada Pacific 
Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug 

Mark, 


??Because I am always trying to get a relocatable solution I never make use of the Windows registry. The only system change I make (with the user's permission) is to add something to the path. 


??Now, when you start a version of grads cygwin sets aside a chunk of shared memory for the POSIX emulation. Since v1.9 and v2.0 uses different versions o cygwin, running both versions without an intervening reboot could cause problems as these shared memory chunks could be inconsistent. So, in the machine where you are having problems, do a reboot and then try running v2.0 without having run v1.9. Any luck? 


??Arlindo ?? 





On Wed, Jan 27, 2010 at 10:45 PM, Mark Sponsler < msponsler at comcast.net > wrote: 






Hi Arlindo, 

Another finding: 



I have both the 1.9 and the 2.0 Grads installtions on 2 different flash drives. 



I went to a computer that has never had Grads installed, plugged in the 1.9 flash drive and fired up grads, then ran 

! ls? 

and sucessfully received a directory listing.? 



Did the same thing with the Grads 2.0 flash drive (on the same computer) and also received a directory liating. 



To me that suggests that indeed there is some sort of conflict occurring?on the machine that has previous grads installations. 



I went back to the machine that is having problems?and pulled wgrib out of the path too (along with all grads path variables). 



Plugged in the 1.9 flashdrive and got a Directory list, but still no luck with the 2.0 flash drive. 



I also have version 1.8 installed on that machine and am thinking that maybe?when I installed it I used the installer version instead of the zip installtion, and that maybe?some remnants are still lurking in the registry that might be causing a conflict.? That's purely a guess though.? Not even sure that would make a difference.? 



Any other suggestions????? 



Thanks, 
Mark 

----- Original Message ----- 
From: "Arlindo da Silva" < dasilva at alum.mit.edu > 
To: "Mark Sponsler" < msponsler at comcast.net > 
Cc: "gradsusr" < GRADSUSR at list.cineca.it > 
Sent: Wednesday, January 27, 2010 10:42:35 AM GMT -08:00 US/Canada Pacific 
Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug 







On Mon, Jan 25, 2010 at 3:37 PM, Mark Sponsler < msponsler at comcast.net > wrote: 






Hi Arlindo, 

Interesting..... 



If I start Grads then type in: 

! ls? <cr> 



Grads just hangs.?No answer set is returned.? I presume a directory list should appear. But instead?nothing.... 




Hmm... It works OK on my Windows platform. Anybody out there having similar problems on Windows? 


I wonder if your GrADS v1.9 installation is not conflicting with your GrADS v2. Do you explicitly set any GrADS environment (GADDIR, GASCRP, etc.)? Is the ?v1.9 bin/ directory on your PATH? Did you use the .exe or .zip version of the superpack? ? 


?? ? Arlindo 

-- 
Arlindo da Silva 
dasilva at alum.mit.edu 



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

From sushantpuranik at GMAIL.COM  Wed Feb  3 00:21:05 2010
From: sushantpuranik at GMAIL.COM (sushant puranik)
Date: Wed, 3 Feb 2010 10:51:05 +0530
Subject: How to reverse the plot
Message-ID: <WED.3.FEB.2010.105105.0530.>

Hello
I have ECMWF data at 1.5 deg resolution which is arranged in reverse order
showing 50mb at surface and 1000 mb at top.
I want to plot the vertical cross section.
Is there any command by which i can reverse the plot showing 1000 mb at
surface and 50 mb at the top in GrADS.


Any suggestion is appreciated



Thank you



Sushant
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100203/7d005a1d/attachment.html 

From cmunozc at UDEC.CL  Wed Feb  3 11:30:50 2010
From: cmunozc at UDEC.CL (Cristian Gonzal =?iso-8859-1?Q?Mu=F1oz_Castillo?=)
Date: Wed, 3 Feb 2010 13:30:50 -0300
Subject: Wind speed direction
In-Reply-To: <72e8bd3d1002022121v37263755h67c4bcca2b1fc6fe@mail.gmail.com>
Message-ID: <WED.3.FEB.2010.133050.0300.>

Hi Grads users,

I have wind speed data and i would like to display it on grads as vectors
on each grid point. I also have on a different variable their direction
measured as sexagesimal degrees. Is there any way to do what i want? That
is associate each wind speed with its corresponding direction and display
that association on grads as a vector.

Thank you all.

-- 
Cristian G. Mu?oz Castillo
Geof?sico
Licenciado en Ciencias F?sicas
Oficina 442
Depto. Geof?sica (DGEO)
Fono (41) 2203155
Facultad Ciencias F?sicas y Matem?ticas
Universidad de Concepci?n
Concepci?n
Chile


From spc_meteor at YAHOO.COM.BR  Wed Feb  3 11:59:05 2010
From: spc_meteor at YAHOO.COM.BR (Saulo Carvalho)
Date: Wed, 3 Feb 2010 08:59:05 -0800
Subject: Problems: OpenGrADS on OpenSUSE
Message-ID: <WED.3.FEB.2010.085905.0800.>

Dear all,

I want to run OpenGrADS on Open SUSE. However, it returns like that:

"Xlib: connection to ":0.0" refused by server
Xlib: No protocol specified

Error in GXSTRT: Unable to connect to X server"

What could it be? What can I correct that problem.

Thanks



 ________________________________________
Saulo Carvalho


SEMA
Diretoria de Recursos H?dricos



      ____________________________________________________________________________________
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/20100203/4a4330ab/attachment.html 

From sulcaf5 at GMAIL.COM  Wed Feb  3 13:23:09 2010
From: sulcaf5 at GMAIL.COM (juan carlos sulca jota)
Date: Wed, 3 Feb 2010 13:23:09 -0500
Subject: que comando se utilizan para que los vectores que grafique en
 grads tengan a su vez sus colores.
Message-ID: <WED.3.FEB.2010.132309.0500.>

-- 
Juan C. Sulca Jota
FCF - UNMSM
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100203/434ce682/attachment.html 

From msponsler at COMCAST.NET  Wed Feb  3 14:40:46 2010
From: msponsler at COMCAST.NET (Mark Sponsler)
Date: Wed, 3 Feb 2010 19:40:46 +0000
Subject: WIN GRADS 2.0.a7.oga.3 - Possible Bug
In-Reply-To: <2048727122.2665761265225966483.JavaMail.root@sz0052a.emeryville.ca.mail.comcast.net>
Message-ID: <WED.3.FEB.2010.194046.0000.>



Hi Arlindo, 

I tried again?to get the 'ls' command to work in Grads 2 without sucess.?? 



I've beleive?your statement earlier about there being some conflict beteween multiple versions of Cygwin to be correct.? I just can't seem to figure out where the conflict is occurring. 



I have 3 version of Grads installed (1.8, 1.9 and 2.0) on the workstation.? I deleted all reference to Grads from the Path environment variable (and most any other applic ation other than the basic windows directories), and also renamed Autoexec.bat (it had a reference to the old regrid function). Then did a cold reboot.?Plugged in the flash drive with grads 2 on it and started Grads from there, but?got not response from??! ls <cr>. 



I then renamed the actual top directories of those 3 versions of grads so even if?there was some buried reference, the code would not be found. Restarted .? Still no luck with the ! ls command. 



Any other suggestions?????? 

Thanks, 
Mark 

----- Original Message ----- 
From: "Arlindo da Silva" <dasilva at alum.mit.edu> 
To: "Mark Sponsler" <msponsler at comcast.net> 
Cc: "gradsusr" <GRADSUSR at list.cineca.it> 
Sent: Wednesday, January 27, 2010 8:59:01 PM GMT -08:00 US/Canada Pacific 
Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug 

Mark, 


??Because I am always trying to get a relocatable solution I never make use of the Windows registry. The only system change I make (with the user's permission) is to add something to the path. 


??Now, when you start a version of grads cygwin sets aside a chunk of shared memory for the POSIX emulation. Since v1.9 and v2.0 uses different versions o cygwin, running both versions without an intervening reboot could cause problems as these shared memory chunks could be inconsistent. So, in the machine where you are having problems, do a reboot and then try running v2.0 without having run v1.9. Any luck? 


??Arlindo ?? 





On Wed, Jan 27, 2010 at 10:45 PM, Mark Sponsler < msponsler at comcast.net > wrote: 






Hi Arlindo, 

Another finding: 



I have both the 1.9 and the 2.0 Grads installtions on 2 different flash drives. 



I went to a computer that has never had Grads installed, plugged in the 1.9 flash drive and fired up grads, then ran 

! ls? 

and sucessfully received a directory listing.? 



Did the same thing with the Grads 2.0 flash drive (on the same computer) and also received a directory liating. 



To me that suggests that indeed there is some sort of conflict occurring?on the machine that has previous grads installations. 



I went back to the machine that is having problems?and pulled wgrib out of the path too (along with all grads path variables). 



Plugged in the 1.9 flashdrive and got a Directory list, but still no luck with the 2.0 flash drive. 



I also have version 1.8 installed on that machine and am thinking that maybe?when I installed it I used the installer version instead of the zip installtion, and that maybe?some remnants are still lurking in the registry that might be causing a conflict.? That's purely a guess though.? Not even sure that would make a difference.? 



Any other suggestions????? 



Thanks, 
Mark 

----- Original Message ----- 
From: "Arlindo da Silva" < dasilva at alum.mit.edu > 
To: "Mark Sponsler" < msponsler at comcast.net > 
Cc: "gradsusr" < GRADSUSR at list.cineca.it > 
Sent: Wednesday, January 27, 2010 10:42:35 AM GMT -08:00 US/Canada Pacific 
Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug 







On Mon, Jan 25, 2010 at 3:37 PM, Mark Sponsler < msponsler at comcast.net > wrote: 






Hi Arlindo, 

Interesting..... 



If I start Grads then type in: 

! ls? <cr> 



Grads just hangs.?No answer set is returned.? I presume a directory list should appear. But instead?nothing.... 




Hmm... It works OK on my Windows platform. Anybody out there having similar problems on Windows? 


I wonder if your GrADS v1.9 installation is not conflicting with your GrADS v2. Do you explicitly set any GrADS environment (GADDIR, GASCRP, etc.)? Is the ?v1.9 bin/ directory on your PATH? Did you use the .exe or .zip version of the superpack? ? 


?? ? Arlindo 

-- 
Arlindo da Silva 
dasilva at alum.mit.edu 



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

From arlindo.dasilva at GMAIL.COM  Wed Feb  3 22:49:01 2010
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Wed, 3 Feb 2010 22:49:01 -0500
Subject: Problems: OpenGrADS on OpenSUSE
In-Reply-To: <746575.39080.qm@web110701.mail.gq1.yahoo.com>
Message-ID: <WED.3.FEB.2010.224901.0500.>

On Wed, Feb 3, 2010 at 11:59 AM, Saulo Carvalho <spc_meteor at yahoo.com.br>wrote:

> Dear all,
>
> I want to run OpenGrADS on Open SUSE. However, it returns like that:
>
> "Xlib: connection to ":0.0" refused by server
> Xlib: No protocol specified
>
> Error in GXSTRT: Unable to connect to X server"
>
> What could it be? What can I correct that problem.
>
>
Can you start any other X11 clients, e.g, xterm? Is your DISPLAY variable
set?

   Arlindo


> Thanks
>
>
>
> ________________________________________
> Saulo Carvalho
>
> SEMA
> Diretoria de Recursos H?dricos
>
>
> ------------------------------
> 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/20100203/324c3a69/attachment.html 

From arlindo.dasilva at GMAIL.COM  Wed Feb  3 22:58:19 2010
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Wed, 3 Feb 2010 22:58:19 -0500
Subject: difference between two files
In-Reply-To: <4B682800.5090802@dge.inpe.br>
Message-ID: <WED.3.FEB.2010.225819.0500.>

On Tue, Feb 2, 2010 at 8:26 AM, Gisele Zepka <gzepka at dge.inpe.br> wrote:

> Hello users!!
>
> I have two files with different resolutions and I want to compare the
> differences between them. I used the regrid function to take both in the
> same resolution. The regrided file is just one hour, while in the other I
> have a range of 22 times. I want to compare just the hour of the first, so I
> have two problems when I try to do that: an error of resolution, because the
> two files don't have the same quantity of points in x and y and I don't know
> how to set the time just for the second and then calculate the difference.
>
>
I am not sure I understand exactly what your problem is. First step is to
get the 2 files in the same horizontal grid. This means the same number of
points and same origin. Same grid. You can do that with the re() function
that comes with the opengrads bundle:

    http://opengrads.org/doc/udxt/re/

For the sake of argument let's say that you have done the regridding and
prepared 2 files on the same grid. Then open the 2 files

ga-> open file1
ga-> open file2

To compare variables valid at different times explicitly specify the
t-index, e.g.,

ga-> d var.2(t=3) - var.1(t=1)

  Good Luck,

     Arlindo







> Please any help will be nice!
>
> Thanks,
> Gisele
>
>
> --
> Msc. Gisele dos Santos Zepka
> Atmospheric Electricity Group (ELAT)
> National Institute for Space Research (INPE)
> Av. dos Astronautas, 1758
> S?o Jos? dos Campos/ SP/ Brazil CEP 12227-010
> Phone ++55(12)39456841
>



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

From arlindo.dasilva at GMAIL.COM  Wed Feb  3 23:03:54 2010
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Wed, 3 Feb 2010 23:03:54 -0500
Subject: How to start re function?
In-Reply-To: <201002011543345539423@163.com>
Message-ID: <WED.3.FEB.2010.230354.0500.>

On Mon, Feb 1, 2010 at 2:43 AM, Yan Libin <zjuyanlibin at 163.com> wrote:

>  Hi, every one,
>
> I want to regrid files into one-standard-grid files.
> and install "grads-2.0.a7.oga.3-win32_superpack.exe " under windows xp
> operation.
>
> I  click the shortcut of opengrads and input "re" command,
> but it pops up as ''no such command"
> How can i resolve this problem?
>
>

re is a function, not a command. Take a look at the re() man page for some
usage example:

          http://opengrads.org/doc/udxt/re/

  Arlindo

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

From mmahakur at TROPMET.RES.IN  Thu Feb  4 08:27:32 2010
From: mmahakur at TROPMET.RES.IN (M.Mahakur)
Date: Thu, 4 Feb 2010 18:57:32 +0530
Subject: stnmap error
In-Reply-To: <170200524.2666391265226046913.JavaMail.root@sz0052a.emeryville.ca.mail.comcast.net>
Message-ID: <THU.4.FEB.2010.185732.0530.>

Respected All,

While runnig stnmap [ $stnmap -i stn.ctl ] I am getting the following error. What does
this mean?

----------------------------------------------------------------------------------
$ stnmap -i stn.ctl
  Name of binary data set: K1V_dat.bin
  Number of times in the data set: 1
  Number of surface variables: 3
  Number of level dependent variables: 0

Starting scan of station data binary file.
Binary data file open: K1V_dat.bin

Processing time step 1
Sequential Read Error: Record size greater than one station report: 40 :: 32
----------------------------------------------------------------------------------

But, when I used:

$stnmap -0 -i stn.ctl it creates the map file (i.e. as if no error). But, later when I
tried to display, I am getting the following errors:

------------------------------------------------------------
Low Level I/O Error:  Sequential read error
  Record size exceeds report size
  Data file name = K1V_dat.bin 
Data Request Error:  Variable is 'gc'
  Error ocurred at column 1
DISPLAY error:  Invalid expression 
  Expression = gc
------------------------------------------------------------

Thanking you.


With Regards,

Mahakur



 


From azarzarrin at YAHOO.COM  Thu Feb  4 10:52:53 2010
From: azarzarrin at YAHOO.COM (Azar Zarrin)
Date: Thu, 4 Feb 2010 07:52:53 -0800
Subject: How to reverse the plot
In-Reply-To: <72e8bd3d1002022121v37263755h67c4bcca2b1fc6fe@mail.gmail.com>
Message-ID: <THU.4.FEB.2010.075253.0800.>

Hi Sushant,

You can use? set yflip on 
?to flip the order of the vertical axis.

Also try? set zlog on
to set log scaling of the Z dimension

Azar

--- On Wed, 2/3/10, sushant puranik <sushantpuranik at GMAIL.COM> wrote:

From: sushant puranik <sushantpuranik at GMAIL.COM>
Subject: How to reverse the plot
To: GRADSUSR at LIST.CINECA.IT
Date: Wednesday, February 3, 2010, 8:21 AM

Hello
I have ECMWF data at 1.5 deg resolution which is arranged in reverse order showing 50mb at surface and 1000 mb at top.
I want to plot the vertical cross section. 
Is there any command by which i can reverse the plot showing 1000 mb at surface and 50 mb at the top in GrADS.



Any suggestion is appreciated



Thank you



Sushant





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

From dasilva at ALUM.MIT.EDU  Thu Feb  4 12:04:48 2010
From: dasilva at ALUM.MIT.EDU (Arlindo da Silva)
Date: Thu, 4 Feb 2010 12:04:48 -0500
Subject: How to start re function?
In-Reply-To: <201002042035049840138@163.com>
Message-ID: <THU.4.FEB.2010.120448.0500.>

On Thu, Feb 4, 2010 at 7:35 AM, Yan Libin <zjuyanlibin at 163.com> wrote:

>  after installation, i can't initiate the function.
> normal installation, used it as 'ave' function
> but it poped up as "no command found"
>

Send me exactly the command you are issuing. Try this:

ga-> q udf
ga-> open $GADSET/model
ga-> d re(ps,1)

BTW, always send questions like this to the list -- please do not write
directly to me.

   Arlindo



> weird~
>
>
> 2010-02-04
> ------------------------------
>  Yan Libin
> ------------------------------
> *????* Arlindo da Silva
> *?????* 2010-02-04  12:04:00
> *????* GRADSUSR
> *???*
> *???* Re: How to start re function?
>  On Mon, Feb 1, 2010 at 2:43 AM, Yan Libin <zjuyanlibin at 163.com> wrote:
>
>>  Hi, every one,
>>
>> I want to regrid files into one-standard-grid files.
>> and install "grads-2.0.a7.oga.3-win32_superpack.exe " under windows xp
>> operation.
>>
>> I  click the shortcut of opengrads and input "re" command,
>> but it pops up as ''no such command"
>> How can i resolve this problem?
>>
>>
>
> re is a function, not a command. Take a look at the re() man page for some
> usage example:
>
>           http://opengrads.org/doc/udxt/re/
>
>   Arlindo
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu
>



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

From fallspinach at GMAIL.COM  Thu Feb  4 17:02:29 2010
From: fallspinach at GMAIL.COM (Ming Pan)
Date: Thu, 4 Feb 2010 17:02:29 -0500
Subject: loop over x,y,z in user defined function (opengrads)
In-Reply-To: <FC782C0C-A183-4FA6-A4C8-A0CAA9F37E1E@cola.iges.org>
Message-ID: <THU.4.FEB.2010.170229.0500.>

Hi Brian,

Thanks again for your kind response. My last email was intended to say
what kind of features of GrADS I'd hope for, from a user's point of
view, and I know essentially nothing about the coding philosophy of
GrADS.

Still, I hope the GrADS development team can consider good features
more favorably than coding style. If there is a conflict between the
features and issues like long-term maintainability, that suggests we
should probably improve our software engineering and migrate to a more
flexible/maintainable design. For example, c++ would let you
safely/conveniently reload functions and prevent users from accessing
what they shouldn't touch. It's worthy breaking some of the
compatibility (some is already broken in version 2) if we can make
GrADS grow more and more powerful. Otherwise we just let GrADS be
limited by some old designs.

By no mean to say any old design is not good -- just to say we should
let things grow.

ming



On Fri, Jan 22, 2010 at 2:46 PM, Brian Doty <doty at cola.iges.org> wrote:
> Ming, this question gets into the philosophy of how one writes code that is
> stable long term, code that is easy to change, how to avoid bloatware, etc.
> ?So it depends on what the goals are for the code. ?If it is for personal
> use, then one may not care much.
>
> By providing defined call-back entry points, and by isolating data
> structures, it is possible to have a base code that can be changed easily,
> and also allow external functions that don't have to change when the base
> code changes. ?Even with this sort of interface, there is nothing to stop an
> external function from ignoring this interface and calling back into grads
> anywhere, or use any internal data structure. ?In that case, a responsible
> code writer would be sure to flag this behavior very prominently, since such
> a code would be more effort to maintain over time. ?One can forsee that if
> there ends up being a lot of external functions like this, and they
> circulate through the community, then when certain changes are made to the
> grads base code there will be many headaches (and complaints). ?Not really a
> good place to be.
>
> So do we try to provide an interface that isolates dimension-looping
> external functions, such that they can be protected from changes to the base
> code? ?This is a tough question. ?It is a harder thing to do. ?I would
> personally hesitate to commit to such a path without really thinking it
> through.... Brian
>
>
> On Jan 22, 2010, at 1:59 PM, Ming Pan wrote:
>
>> Hi Brian,
>>
>> I really appreciate your explanations on the design of grads and its
>> user defined function interface.
>>
>> While I feel it absolutely necessary to limit memory consumption, I
>> wish the new interface will consider allowing user program to pull
>> more data (e.g. loop over an additional dimension) some way. For
>> example, we can allow passing only one grid at a time, but also allow
>> re-using that memory block for passing more in a loop. It's indeed a
>> trade-off between how much you allow user program to do and the risk
>> of breaking grads itself. (Re-design the software with C++ may solve
>> some fundamental issues on API?) Arlindo's UDF design gives me really
>> good flexibility and fast speed, though I have to read a lot of the
>> grads source and go through many crash/debug cycles with gdb. I don't
>> really mind about that since it's eventually the user's responsibility
>> not to blow up the memory or cause segmentation fault. To use the less
>> efficient external files for data transfer will completely isolate
>> grads so that it never crashes on user's fault, and users do not have
>> to learn too much about how grads works internally. On the other hand,
>> it doesn't really make a difference for me whether it crashes inside
>> grads or my functions -- both are my own fault.
>>
>> I completely agree that the returned grid from a user function should
>> not have fundamentally different dimensions. At the same time I wish
>> the new interface would allow returning more than one outputs. I once
>> needed this desperately and ended up writing a user defined command
>> (Arlindo's interface) to replace the default "define" command so that
>> I can receive/define multiple grids from my customized function in one
>> shot.
>>
>> ming
>>
>> On Tue, Jan 19, 2010 at 12:15 PM, Brian Doty <doty at cola.iges.org> wrote:
>>>
>>> The expression handling part of grads is fundamentally limited to grids
>>> with
>>> a maximum of 2D varying dimensions. ?This goes back to the early 90s when
>>> I
>>> designed and implemented the language known as the "grads expression".
>>> ?This
>>> design has stood the test of time -- if anything, experience since then
>>> has
>>> shown this to be extremely important, as it limits the amount of memory
>>> that
>>> grads consumes. ?You can, of course, override this at any time by staging
>>> data in memory via the "define" feature. ?Define can use great amounts of
>>> memory very rapidly, but it is under your control.
>>>
>>> There are functions, such as ave, that loop over some dimension and
>>> evaluate
>>> a target expression that results in a 2D or less grid, and then does some
>>> operation on those grids (such as sum them and divide by n). ?Note that
>>> at
>>> no time is there a grid present that is more than 2D. ?The looping
>>> results
>>> in potentially a lot of I/O. ?If the data is going to be looped through
>>> more
>>> than once (or sliced in some non-efficient way), then putting the data
>>> into
>>> memory in advance, via define, can yield large performance improvements.
>>> ?(It is important to understand the way grads evaluates expressions and
>>> what
>>> implications this has for I/O if you are working with a lot of data --
>>> the
>>> way you do certain things can have huge performance differences.)
>>>
>>> The functions that loop through a dimension (eg, ave) operate
>>> recursively.
>>> ?They are thus not suitable as UDF functions, since I don't see how to
>>> easily provide a standardized interface without very high long term
>>> maintenance costs. ?Thus I would personally recommend that you implement
>>> such a function as a source code mod, this being much safer in the long
>>> term, as it will make clear to future users the nature of the function
>>> and
>>> the level of burden for maintaining it.
>>>
>>> As a side note, it is also important that no function outputs a grid that
>>> has fundamentally different dimensions than the input dimension
>>> environment.
>>> ?This is a necessary restriction to preserve the internal consistency of
>>> the
>>> expression handling. ? What I mean by this can be shown by an example:
>>> ?if
>>> the input dimension environment is X and Y varying, it would be invalid
>>> to
>>> output a grid with some dimension other than X and Y varying (eg, T).
>>> It
>>> is ok to output a grid that is a subset, ie, with just X varying.
>>>
>>> Another approach is to do what I do in this situation: ?I write a script
>>> and
>>> a C program, and my script outputs the needed data into a file, and the C
>>> program does the calculation, and the script reads in the data output by
>>> the
>>> C program. ?Somewhat kludgy, but not all functions are worth the time and
>>> effort to implement within grads and maintain over the longer term, ?IMO.
>>>
>>> I am working through the design of some sort of interface to allow tasks
>>> to
>>> be done directly on "defined" objects. ?This is an outgrowth of our
>>> recent
>>> implementation of netCDF file output. ?Jennifer and I had extensive
>>> design
>>> discussions on this feature, and it became clear that we needed to
>>> collect
>>> large amounts of data in memory to optimize the I/O, and define was the
>>> obvious way to do this. ?Now I am working on a way to do a "sort"
>>> operation
>>> on a "defined" data object, where the sort would operate though one of
>>> the
>>> dimensions. ?I am also thinking that this type of thing could be
>>> generalized
>>> and a "UDF" style capability could be provided to perform some arbitrary
>>> task on a defined object. ?This would not be part of the expression
>>> capability, but would be for one-shot types of things. ?I've been trying
>>> to
>>> figure out how to provide a sort feature and this is the only way I can
>>> see
>>> to do it cleanly... Brian
>>>
>>>
>>> On Jan 13, 2010, at 3:37 PM, Ming Pan wrote:
>>>
>>>> Dear Users/Developers,
>>>>
>>>> I'm using grads-2.0.a7.oga.3 and trying to write a customized function
>>>> to compute low-level cape. In this case, I need to loop over 3
>>>> dimension x, y, and z, but it seems the current udf interface only
>>>> allows passing 2-D grids (isiz, jsiz). Anyone ran into a similar
>>>> situation and came up with a solution/work-around?
>>>>
>>>> Thanks!
>>>>
>>>> ming
>>>
>


From javier.p at EREDA.COM  Fri Feb  5 05:16:02 2010
From: javier.p at EREDA.COM (=?ISO-8859-1?Q?Javier_Pe=F1a?=)
Date: Fri, 5 Feb 2010 11:16:02 +0100
Subject: Interpolate 3D variables on pressure levels or model levels to
 height levels
Message-ID: <FRI.5.FEB.2010.111602.0100.>

Dear users:

 I have a doubt which someone could solve it:

I have two data type from:
For one hand I have 16 levels of model levels (45,46,...60) wiht variables
u-wind, v-wind, temperature, humedity relative and  the surface level with
variables: geopotential and log of surface pressure.
For other hand I have 10 levels of pressure levels ( 1000, 975, 950,
....775) with variables u-wind, v-wind and geopotential.

They are grib files and, with gribctl2.pl, gribmap and grads, I can extract
all data in ASCII format but, the question is the following:

Is it possible to interpolate the data in height and obtain those variables
in the heigth that I want. I mean, Is it possible to obtain, from the model
levels or pressure levels, heigth levels with GRADS or other program?


Thank you very much at all

Kind regards:

Javier Pe?a


-- 
Javier Pe?a ?lvarez
Energ?as Renovables
javier.p at ereda.com
EREDA

Tel:  +34 91 5014755
Fax: +34 91 5014756
P?. Marqu?s de Monistrol 7
28011 MADRID
www.ereda.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100205/d9cfd29a/attachment.html 

From mequitomz at YAHOO.COM.BR  Fri Feb  5 06:46:22 2010
From: mequitomz at YAHOO.COM.BR (Arlindo Meque)
Date: Fri, 5 Feb 2010 03:46:22 -0800
Subject: Problem Using sublin and subwrd
Message-ID: <FRI.5.FEB.2010.034622.0800.>

Dear grads Users,
?
?
1.I am trying to use the sublin and subwrd commands/functions. Unfortunately, I am getting error messages as follow:
?
?
?
?a)??ga-> q dims
Default file number is: 1 
X is varying?? Lon = -100 to -20?? X = -39 to -7
Y is varying?? Lat = -60 to 20?? Y = 13 to 45
Z is fixed???? Lev = 1000? Z = 1
T is fixed???? Time = 00Z01JAN1979? T = 1
E is fixed???? Ens = 1? E = 1
ga-> ll=sublin(result,2)
Syntax Error:? Invalid Operand 
? 'sublin' not a variable or function name
? Error ocurred at column 1
DEFINE error:? Invalid expression. 
?
?
b)??ga-> q time
Time = 00Z01JAN1979 to 00Z01JAN1979? Mon to Mon
ga-> res=subwrd(result,3)
Syntax Error:? Invalid Operand 
? 'subwrd' not a variable or function name
? Error ocurred at column 1
DEFINE error:? Invalid expression. 
ga-> 
?
?
?
2. My version of GrADS is:
?
?
ga-> q config
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds hdf5 opendap-grids,stn athena geotiff
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
Copyright (c) 1988-2009 by Brian Doty and the
Institute for Global Environment and Society (IGES) 
This program is distributed WITHOUT ANY WARRANTY 
See file COPYRIGHT for more information. 
?
?
?
?
?
Help is needed
?
Arlindo 


      ____________________________________________________________________________________
Veja quais s?o os assuntos do momento no Yahoo! +Buscados
http://br.maisbuscados.yahoo.com


From madan at EES.HOKUDAI.AC.JP  Fri Feb  5 07:03:21 2010
From: madan at EES.HOKUDAI.AC.JP (madan sigdel)
Date: Fri, 5 Feb 2010 21:03:21 +0900
Subject: Problem Using sublin and subwrd
In-Reply-To: <316270.7968.qm@web52305.mail.re2.yahoo.com>
Message-ID: <FRI.5.FEB.2010.210321.0900.>

These functions only work with in GRADS scripts, not with command line.

-----Original Message-----
From: GRADSUSR at LIST.CINECA.IT [mailto:GRADSUSR at LIST.CINECA.IT] On Behalf Of
Arlindo Meque
Sent: Friday, February 05, 2010 8:46 PM
To: GRADSUSR at LIST.CINECA.IT
Subject: Problem Using sublin and subwrd

Dear grads Users,
?
?
1.I am trying to use the sublin and subwrd commands/functions.
Unfortunately, I am getting error messages as follow:
?
?
?
?a)??ga-> q dims
Default file number is: 1 
X is varying?? Lon = -100 to -20?? X = -39 to -7
Y is varying?? Lat = -60 to 20?? Y = 13 to 45
Z is fixed???? Lev = 1000? Z = 1
T is fixed???? Time = 00Z01JAN1979? T = 1
E is fixed???? Ens = 1? E = 1
ga-> ll=sublin(result,2)
Syntax Error:? Invalid Operand 
? 'sublin' not a variable or function name
? Error ocurred at column 1
DEFINE error:? Invalid expression. 
?
?
b)??ga-> q time
Time = 00Z01JAN1979 to 00Z01JAN1979? Mon to Mon
ga-> res=subwrd(result,3)
Syntax Error:? Invalid Operand 
? 'subwrd' not a variable or function name
? Error ocurred at column 1
DEFINE error:? Invalid expression. 
ga-> 
?
?
?
2. My version of GrADS is:
?
?
ga-> q config
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf hdf4-sds
hdf5 opendap-grids,stn athena geotiff
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
Copyright (c) 1988-2009 by Brian Doty and the
Institute for Global Environment and Society (IGES) 
This program is distributed WITHOUT ANY WARRANTY 
See file COPYRIGHT for more information. 
?
?
?
?
?
Help is needed
?
Arlindo 


 
____________________________________________________________________________
________
Veja quais s?o os assuntos do momento no Yahoo! +Buscados
http://br.maisbuscados.yahoo.com


From spc_meteor at YAHOO.COM.BR  Fri Feb  5 07:08:19 2010
From: spc_meteor at YAHOO.COM.BR (Saulo Carvalho)
Date: Fri, 5 Feb 2010 04:08:19 -0800
Subject: Res: Problems: OpenGrADS on OpenSUSE
In-Reply-To: <77fcd6b21002031949m72cffebbm967a5e958685b5a5@mail.gmail.com>
Message-ID: <FRI.5.FEB.2010.040819.0800.>


Arlindo,

I started the xterm, but it's not  working yet.

 ________________________________________
Saulo Carvalho


SEMA
Diretoria de Recursos H?dricos




________________________________
De: Arlindo da Silva <arlindo.dasilva at GMAIL.COM>
Para: GRADSUSR at LIST.CINECA.IT
Enviadas: Quinta-feira, 4 de Fevereiro de 2010 0:49:01
Assunto: Re: Problems: OpenGrADS on OpenSUSE

On Wed, Feb 3, 2010 at 11:59 AM, Saulo Carvalho <spc_meteor at yahoo.com.br> wrote:

Dear all,
>
>I want to run OpenGrADS on Open SUSE. However, it returns like that:
>
>"Xlib: connection to ":0.0" refused by server
>>Xlib: No protocol specified
>
>Error in GXSTRT: Unable to connect to X server"
>
>What could it be? What can I correct that problem.
>
>

Can you start any other X11 clients, e.g, xterm? Is your DISPLAY variable set?

   Arlindo
 
>Thanks
>
>
>
> ________________________________________
>Saulo Carvalho
>
>
>SEMA
>Diretoria de Recursos H?dricos
>
>
>>
>
>________________________________
Veja quais s?o os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - M?sica - Esportes


-- 
Arlindo da Silva
dasilva at alum.mit.edu



      ____________________________________________________________________________________
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/20100205/2fed8968/attachment.html 

From jdduda at IASTATE.EDU  Fri Feb  5 11:50:19 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Fri, 5 Feb 2010 10:50:19 -0600
Subject: Interpolate 3D variables on pressure levels or model levels
 to height levels
In-Reply-To: <607bb2141002050216h7c6fb0dbt2e10c1289beaa2ce@mail.gmail.com>
Message-ID: <FRI.5.FEB.2010.105019.0600.>

If you're running WRF, there's either the WRF Post Processor (
http://www.dtcenter.org/wrf-nmm/users/downloads/index.php) or ARWPost (
http://www.mmm.ucar.edu/wrf/users/download/get_sources.html).  Each one can
take model level data and interpret it to height levels, pressure levels, or
just keep it on model levels but add diagnostic variables.  You do need to
register an email address to download them.

Jeff Duda

On Fri, Feb 5, 2010 at 4:16 AM, Javier Pe?a <javier.p at ereda.com> wrote:

> Dear users:
>
>  I have a doubt which someone could solve it:
>
> I have two data type from:
> For one hand I have 16 levels of model levels (45,46,...60) wiht variables
> u-wind, v-wind, temperature, humedity relative and  the surface level with
> variables: geopotential and log of surface pressure.
> For other hand I have 10 levels of pressure levels ( 1000, 975, 950,
> ....775) with variables u-wind, v-wind and geopotential.
>
> They are grib files and, with gribctl2.pl, gribmap and grads, I can
> extract all data in ASCII format but, the question is the following:
>
> Is it possible to interpolate the data in height and obtain those variables
> in the heigth that I want. I mean, Is it possible to obtain, from the model
> levels or pressure levels, heigth levels with GRADS or other program?
>
>
> Thank you very much at all
>
> Kind regards:
>
> Javier Pe?a
>
>
> --
> Javier Pe?a ?lvarez
> Energ?as Renovables
> javier.p at ereda.com
> EREDA
>
> Tel:  +34 91 5014755
> Fax: +34 91 5014756
> P?. Marqu?s de Monistrol 7
> 28011 MADRID
> www.ereda.com
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100205/ad5f4557/attachment.html 

From arlindo.dasilva at GMAIL.COM  Fri Feb  5 13:14:46 2010
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Fri, 5 Feb 2010 13:14:46 -0500
Subject: Problems: OpenGrADS on OpenSUSE
In-Reply-To: <850326.25592.qm@web110716.mail.gq1.yahoo.com>
Message-ID: <FRI.5.FEB.2010.131446.0500.>

On Fri, Feb 5, 2010 at 7:08 AM, Saulo Carvalho <spc_meteor at yahoo.com.br>wrote:

>
> Arlindo,
>
> I started the xterm, but it's not  working yet.
>
>

I believe this is the same problem: your DISPLAY variable is not properly
set. If you get xterm to start, very likely grads would start as well.

   Arlindo.



> ________________________________________
> Saulo Carvalho
>
> SEMA
> Diretoria de Recursos H?dricos
>
>
> ------------------------------
> *De:* Arlindo da Silva <arlindo.dasilva at GMAIL.COM>
> *Para:* GRADSUSR at LIST.CINECA.IT
> *Enviadas:* Quinta-feira, 4 de Fevereiro de 2010 0:49:01
> *Assunto:* Re: Problems: OpenGrADS on OpenSUSE
>
> On Wed, Feb 3, 2010 at 11:59 AM, Saulo Carvalho <spc_meteor at yahoo.com.br>wrote:
>
>> Dear all,
>>
>> I want to run OpenGrADS on Open SUSE. However, it returns like that:
>>
>> "Xlib: connection to ":0.0" refused by server
>> Xlib: No protocol specified
>>
>> Error in GXSTRT: Unable to connect to X server"
>>
>> What could it be? What can I correct that problem.
>>
>>
> Can you start any other X11 clients, e.g, xterm? Is your DISPLAY variable
> set?
>
>    Arlindo
>
>
>> Thanks
>>
>>
>>
>> ________________________________________
>> Saulo Carvalho
>>
>> SEMA
>> Diretoria de Recursos H?dricos
>>
>>
>> ------------------------------
>> 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
>
> ------------------------------
> 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/20100205/03be3225/attachment.html 

From wharms at BFS.DE  Fri Feb  5 13:46:42 2010
From: wharms at BFS.DE (walter harms)
Date: Fri, 5 Feb 2010 19:46:42 +0100
Subject: Res: Problems: OpenGrADS on OpenSUSE
In-Reply-To: <850326.25592.qm@web110716.mail.gq1.yahoo.com>
Message-ID: <FRI.5.FEB.2010.194642.0100.>

Saulo Carvalho schrieb:
> Arlindo,
> 
> I started the xterm, but it's not  working yet.
> 



I would be important to know: how did you login ?

when using SSH you may need to enable X11 forwarding with the -X option.
(add -C for compression also).

Then try 'xterm' and it should open.

re,
 wh


From rsolanob at GMAIL.COM  Fri Feb  5 14:57:18 2010
From: rsolanob at GMAIL.COM (Ramon Solano)
Date: Fri, 5 Feb 2010 12:57:18 -0700
Subject: Res: Problems: OpenGrADS on OpenSUSE
In-Reply-To: <850326.25592.qm@web110716.mail.gq1.yahoo.com>
Message-ID: <FRI.5.FEB.2010.125718.0700.>

Hello Saulo,

Your problem is not related to GrADS, but with your computer's graphic system (i.e . your X-Window system). The best would be to ask to your local computer guru, but in the mean time, the problem is that you are not allowed to display in your screen.

This can be due to several reasons, depending on if you are working locally in your computer ore remotely on a server. If you are remote user, be sure to enable X-Window in your local computer, configure it to accept incoming X-11 connections and to forward the graphics  from the server to your local computer.  If you are working locally, be sure you have a X-Window system installed in your computer. If you have done all this, then your X-Window system may have some misconfiguration.

Once you are able to run xterm or xclock or some other graphic program, then GrADS should open with no problems as well.

Ramon

On Feb 5, 2010, at 5:08 AM, Saulo Carvalho wrote:

> 
> Arlindo,
> 
> I started the xterm, but it's not  working yet.
>  
> ________________________________________
> Saulo Carvalho
> 
> SEMA
> Diretoria de Recursos H?dricos
> 
> 
> De: Arlindo da Silva <arlindo.dasilva at GMAIL.COM>
> Para: GRADSUSR at LIST.CINECA.IT
> Enviadas: Quinta-feira, 4 de Fevereiro de 2010 0:49:01
> Assunto: Re: Problems: OpenGrADS on OpenSUSE
> 
> On Wed, Feb 3, 2010 at 11:59 AM, Saulo Carvalho <spc_meteor at yahoo.com.br> wrote:
> Dear all,
> 
> I want to run OpenGrADS on Open SUSE. However, it returns like that:
> 
> "Xlib: connection to ":0.0" refused by server
> Xlib: No protocol specified
> 
> Error in GXSTRT: Unable to connect to X server"
> 
> What could it be? What can I correct that problem.
> 
> 
> Can you start any other X11 clients, e.g, xterm? Is your DISPLAY variable set?
> 
>    Arlindo
>  
> Thanks
> 
> 
>  
> ________________________________________
> Saulo Carvalho
> 
> SEMA
> Diretoria de Recursos H?dricos
> 
> 
> Veja quais s?o os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - M?sica - Esportes
> 
> 
> 
> -- 
> Arlindo da Silva
> dasilva at alum.mit.edu
> 
> Veja quais s?o os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - M?sica - Esportes

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100205/2b167dfe/attachment.html 

From spc_meteor at YAHOO.COM.BR  Fri Feb  5 19:37:30 2010
From: spc_meteor at YAHOO.COM.BR (Saulo Carvalho)
Date: Fri, 5 Feb 2010 16:37:30 -0800
Subject: Res: Res: Problems: OpenGrADS on OpenSUSE
In-Reply-To: <7F3B3F6A-4E1E-43E2-9A1B-C54421F793A7@gmail.com>
Message-ID: <FRI.5.FEB.2010.163730.0800.>

Ok Ramon. I'll check this out.

 Thanks for the advice.________________________________________
Saulo Carvalho


SEMA
Diretoria de Recursos H?dricos




________________________________
De: Ramon Solano <rsolanob at GMAIL.COM>
Para: GRADSUSR at LIST.CINECA.IT
Enviadas: Sexta-feira, 5 de Fevereiro de 2010 16:57:18
Assunto: Re: Res: Problems: OpenGrADS on OpenSUSE

Hello Saulo,

Your problem is not related to GrADS, but with your computer's graphic system (i.e . your X-Window system). The best would be to ask to your local computer guru, but in the mean time, the problem is that you are not allowed to display in your screen.

This can be due to several reasons, depending on if you are working locally in your computer ore remotely on a server. If you are remote user, be sure to enable X-Window in your local computer, configure it to accept incoming X-11 connections and to forward the graphics  from the server to your local computer.  If you are working locally, be sure you have a X-Window system installed in your computer. If you have done all this, then your X-Window system may have some misconfiguration.

Once you are able to run xterm or xclock or some other graphic program, then GrADS should open with no problems as well.

Ramon


On Feb 5, 2010, at 5:08 AM, Saulo Carvalho wrote:


>Arlindo,
>
>I started the xterm, but it's not  working yet.
>
> ________________________________________
>Saulo Carvalho
>
>
>SEMA
>Diretoria de Recursos H?dricos
>
>
>
>
>
________________________________
De: Arlindo da Silva <arlindo.dasilva at GMAIL.COM>
>Para: GRADSUSR at LIST.CINECA.IT
>Enviadas: Quinta-feira, 4 de Fevereiro de 2010 0:49:01
>Assunto: Re: Problems: OpenGrADS on OpenSUSE
>
>On Wed, Feb 3, 2010 at 11:59 AM, Saulo Carvalho <spc_meteor at yahoo.com.br> wrote:
>
>Dear all,
>>
>>I want to run OpenGrADS on Open SUSE. However, it returns like that:
>>
>>"Xlib: connection to ":0.0" refused by server
>>Xlib: No protocol specified
>>
>>Error in GXSTRT: Unable to connect to X server"
>>
>>What could it be? What can I correct that problem.
>>
>>
>
>
>Can you start any other X11 clients, e.g, xterm? Is your DISPLAY variable set?
>
>
>   Arlindo
> 
>Thanks
>>
>>
>>
>> ________________________________________
>>Saulo Carvalho
>>
>>
>>SEMA
>>Diretoria de Recursos H?dricos
>>
>>
>>
>>________________________________
Veja quais s?o os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - M?sica - Esportes
>
>
>-- 
>Arlindo da Silva
>dasilva at alum.mit.edu
>
>________________________________
Veja quais s?o os assuntos do momento no Yahoo! + Buscados: Top 10 - Celebridades - M?sica - Esportes



      ____________________________________________________________________________________
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/20100205/d3e726e9/attachment.html 

From max.rythmos at YAHOO.COM  Sat Feb  6 04:52:20 2010
From: max.rythmos at YAHOO.COM (MARZIO VENEMAN)
Date: Sat, 6 Feb 2010 01:52:20 -0800
Subject: Open Earth
Message-ID: <SAT.6.FEB.2010.015220.0800.>

Dear subscribers Grads Users,,

With reference to the above, I take pleasure informing you as follows:


http://public.deltares.nl/display/OET/OpenEarth
?
OpenEarth is the open source initiative to archive, host and disseminate Data, Models and Tools for marine & coastal scientist and engineers. It aims to remedy the above-described inefficiencies by providing a project-superseding approach.
?OpenEarth is, amongst others, supported by the concerted effort of professionals from Deltares (formerly Delft Hydraulics), Delft University of Technology's Hydraulic Engineering and Environmental Fluid Mechanics sections and UNESCO-IHE. It is currently the central platform for data and knowledge management in the research programmes Delft Cluster - Northsea and Coast, Building with Nature and MICORE.

I trust this information is sufficient for your purposes, in case you require any additional details, please do not hesitate to contact the undersigned.


Yours sincerely,
Cordiali Saluti


Marzio Veneman
The Netherlands

?

Click here to visit my international professional profile and connect!
http://www.linkedin-ech3.com/in/Rythmomachy
http://en.wikipedia.org/wiki/User:Rythmomachy
http://www.thefacebook.com/Rythmomachy
Free?Mobile telephone worldwide http://www.jaxtr.com/max.rythmos
?

?
________________________________

This email message may contain privilliged information and is solely intended for the recipient(s) mentioned above.? To ensure that you continue receiving our emails, please add max.rythmos [at] yahoo_com to your address book or safe list.


________________________________

________________________________


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

From owandaotieno at GMAIL.COM  Sun Feb  7 16:53:05 2010
From: owandaotieno at GMAIL.COM (Owanda Otieno)
Date: Sun, 7 Feb 2010 16:53:05 -0500
Subject: ctl file
Message-ID: <SUN.7.FEB.2010.165305.0500.>

Hi all,
Can anyone help with a descriptor file that can open the data with ncdump
result as shown.

[otieno at cesehpcc NARCCAP]$ ncdump -h pr_CRCM_cgcm3_2038010103.nc
netcdf pr_CRCM_cgcm3_2038010103 {
dimensions:
        xc = 140 ;
        yc = 115 ;
        time = UNLIMITED ; // (8760 currently)
        bnds = 2 ;
variables:
        double yc(yc) ;
                yc:units = "m" ;
                yc:long_name = "y-coordinate of polar_stereographic
projection" ;
                yc:standard_name = "projection_y_coordinate" ;
                yc:axis = "Y" ;
                yc:coordinate_defines = "point" ;
                yc:actual_range = 0.f, 5700000.f ;
        double xc(xc) ;
                xc:units = "m" ;
                xc:long_name = "x-coordinate of polar_stereographic
projection" ;
                xc:standard_name = "projection_x_coordinate" ;
                xc:axis = "X" ;
                xc:coordinate_defines = "point" ;
                xc:actual_range = 0.f, 6950000.f ;
        double time(time) ;
                time:long_name = "time" ;
                time:calendar = "365_day" ;
                time:standard_name = "time" ;
                time:axis = "T" ;
                time:units = "days since 2038-01-01 00:00:0.0" ;
                time:delta_t = "0000-00-00 03:00:00" ;
                time:coordinate_defines = "point" ;
                time:actual_range = 0.125f, 1095.f ;
                time:bounds = "time_bnds" ;
        double lon(yc, xc) ;
                lon:units = "degrees_east" ;
                lon:long_name = "longitude" ;
                lon:standard_name = "longitude" ;
                lon:axis = "X" ;
                lon:actual_range = 199.8962f, 326.4349f ;
        double lat(yc, xc) ;
                lat:units = "degrees_north" ;
                lat:long_name = "latitude" ;
                lat:standard_name = "latitude" ;
                lat:axis = "Y" ;
                lat:actual_range = 20.60496f, 73.2522f ;
        char polar_stereographic ;
                polar_stereographic:straight_vertical_longitude_from_pole =
263. ;
                polar_stereographic:standard_parallel = 60. ;
                polar_stereographic:false_easting = 3450000. ;
                polar_stereographic:false_northing = 7450000. ;
                polar_stereographic:latitude_of_projection_origin = 90. ;
                polar_stereographic:grid_mapping_name =
"polar_stereographic" ;
                polar_stereographic:resolution_at_standard_parallel = 50000.
;
                polar_stereographic:earth_radius = 6371000. ;
        float pr(time, yc, xc) ;
                pr:units = "kg m-2 s-1" ;
                pr:type = "average" ;
                pr:long_name = "Precipitation" ;
                pr:standard_name = "precipitation_flux" ;
                pr:cell_methods = "time:mean(interval: 3 hours)" ;
                pr:missing_value = 1.e+20f ;
                pr:_FillValue = 1.e+20f ;
                pr:actual_range = 0.f, 0.003255121f ;
                pr:add_offset = 0.f ;
                pr:scale_factor = 1.f ;
                pr:coordinates = "lon lat" ;
                pr:grid_mapping = "polar_stereographic" ;
                pr:level_desc = "Surface" ;
                pr:grid_desc = "polar_stereographic" ;
        double time_bnds(time, bnds) ;
// global attributes:
                :Conventions = "CF-1.0" ;
                :institution = "Ouranos (Ouranos, Montreal, PQ, Canada)" ;
                :source =
"CRCM(2006):atmosphere:CRCMv4.2.0(50km,29levels):sea ice:AMASSv2.7" ;
                :project_id = "NARCCAP" ;
                :realization = "1" ;
                :experiment_id = "2038-2070 climate simulation using CGCM3
member #4 riment)" ;
                :experiment_id2 = "Ouranos operational simulation afs" ;
                :contact1 = "Sebastien Biner" ;
                :contact2 = "biner.sebastien at ouranos.ca" ;
                :contact3 = "550 Sherbrooke Ouest, 19e etage, Montreal, PQ,
H3A 1B9,
                :contact4 = "tel : (514)282.6464 poste 263" ;
                :history = "Tue Feb  3 18:17:22 2009: ncrename -a
conventions,Conventpr_CRCM_2038010103.nc\n",
                        "Mon Dec 15 16:31:26 2008: ncatted -a
experiment_id,global,m,climate simulation using CGCM3 member #4 (mc_abv
experiment) 4_NETCDF_FINAUX/pr_CRCM_c\n",
                        "created 24/04/2004 by Ton Nom (netCDF3.6.1
udunits1.12.4)" ;
                :table_id = "Table 2" ;
}
Thank you

-- 
Owanda Otieno
Department of Natural Resources and the Environment
University of Connecticut
1376 Storrs Road, Unit 4087
Storrs, CT 06269-4087
Phone: (860) 486-1876
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100207/b8bd68ca/attachment.html 

From claudiomet at GMAIL.COM  Mon Feb  8 11:11:28 2010
From: claudiomet at GMAIL.COM (claudio cortes)
Date: Mon, 8 Feb 2010 13:11:28 -0300
Subject: Grads custom maps
Message-ID: <MON.8.FEB.2010.131128.0300.>

Greetings

The maps included in GrADS only have international political boundaries for
my region.
Now I need to add regional political boundaries for my country.
How can I create custom maps for GrADS ?

Cheers !



-- 
Claudio Cortes
+56 (2) 2994121

Meteorologo
Laboratorio de Meteorologia y Pronostico de la Calidad del Aire (LMPCA)
Unidad de Modelacion y Gestion de la Calidad del Aire (UMGCA)
Centro Nacional del Medio Ambiente (CENMA)

--
Claudio Cortes
+56 (2) 2994121

Meteorologist
Meteorology and Air Quality Forecast Labs
Modeling and Air Quality Management Unit
National Enviroment Center, Chile (CENMA)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100208/bf3b259f/attachment.html 

From dasilva at ALUM.MIT.EDU  Tue Feb  9 11:11:44 2010
From: dasilva at ALUM.MIT.EDU (Arlindo da Silva)
Date: Tue, 9 Feb 2010 11:11:44 -0500
Subject: Question
In-Reply-To: <c70cd7eb1002040125n13d4b21boa848c4f4017ece36@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.111144.0500.>

On Thu, Feb 4, 2010 at 4:25 AM, Diego Souza <dsouza at cptec.inpe.br> wrote:

> Hi Arlindo,
>
> I am using grads script avaliable in opengrads home page to open the MERRA
> datasets. I am using the merra_hourly.gs with the option lnd and using the
> variable prectot.
> In the case of MERRA monthly mean files, this isn't my objective, but I
> will consider.
>
>
It could be a problem with the OpenDAP server. In this case, contact the
Goddard M-DISC.

                      http://goldsmr1.sci.gsfc.nasa.gov/

   Arlindo

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

From arlindo.dasilva at GMAIL.COM  Tue Feb  9 11:13:41 2010
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Tue, 9 Feb 2010 11:13:41 -0500
Subject: Grads custom maps
In-Reply-To: <5f3ba4e01002080811n69687d55rb1bdd87124815f5c@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.111341.0500.>

On Mon, Feb 8, 2010 at 11:11 AM, claudio cortes <claudiomet at gmail.com>wrote:

> Greetings
>
> The maps included in GrADS only have international political boundaries for
> my region.
> Now I need to add regional political boundaries for my country.
> How can I create custom maps for GrADS ?
>
>

One option is to use shapefiles. See this

    http://opengrads.org/doc/udxt/shape/

Also, search this list for "shapefiles".

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

From hmjbarbosa at GMAIL.COM  Tue Feb  9 12:12:41 2010
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Tue, 9 Feb 2010 15:12:41 -0200
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
Message-ID: <TUE.9.FEB.2010.151241.0200.>

Dear All,

I need help to access CFSR (NCEP reanalysis 3) over opendap with opengrads.
I get an horrendous error message (see below) when I try to:

ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2

I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
have opendap enabled (see output of 'q config' below).

The strange thing is that I can open other sources of data, for example

ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z

works perfectly and I can plot, subset, etc...

Any ideas?

[]s
Henrique


--------------------------------------------------------------------------------------------------------------------------------------------------
ga-> q config
Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
hdf4-sds hdf5 opendap-grids,stn athena geotiff
Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
Copyright (c) 1988-2009 by Brian Doty and the
Institute for Global Environment and Society (IGES)
This program is distributed WITHOUT ANY WARRANTY
See file COPYRIGHT for more information.

Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu

This version of GrADS has been configured with the following options:
  o Built on a LITTLE ENDIAN machine
  o Command line editing ENABLED
      http://tiswww.case.edu/php/chet/readline/rltop.html
  o printim command for image output ENABLED
      http://www.zlib.net
      http://www.libpng.org/pub/png/libpng.html
      http://www.libgd.org/Main_Page
  o GRIB2 interface ENABLED
      http://www.ijg.org
      http://www.ece.uvic.ca/~mdadams/jasper
      http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
      g2clib-1.0.5
  o NetCDF interface ENABLED
      http://www.opendap.org
      libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar  3 2009 14:13:33 $
  o HDF interface ENABLED
      http://hdfgroup.org
      HDF 4.2r3
      HDF5 1.8.2
  o Athena Widget GUI ENABLED
  o OPeNDAP gridded data interface ENABLED
      http://www.opendap.org
      libdap 3.7.10
  o OPeNDAP station data interface ENABLED
      http://iges.org/grads/gadoc/supplibs.html
      libgadap 2.0.oga.1
  o GeoTIFF and KML output ENABLED
      http://www.libtiff.org
      http://geotiff.osgeo.org

For additional information please consult http://iges.org/grads

--------------------------------------------------------------------------------------------------------------------------------------------------
ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2
Scanning self-describing file:
http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
*** glibc detected *** grads: free(): invalid pointer: 0x00000000018ff354 ***
======= Backtrace: =========
/lib/libc.so.6[0x2b6884e9fcb8]
/lib/libc.so.6(cfree+0x76)[0x2b6884ea2276]
grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea]
grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9)[0x833a05]
grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea]
grads(nc3_get_vara_double+0x9e)[0x80e3b8]
grads(nc_get_vara_double+0x1a)[0x60bcca]
grads(read_one_dimension+0x89)[0x541ac1]
grads(gadsdf+0x17d6)[0x540080]
grads(gasdfopen+0x1d4)[0x53e700]
grads(gacmd+0xfb2)[0x524504]
grads(Main+0x4e1)[0x496f01]
grads(main+0x2e)[0x496a16]
/lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6]
grads(__gxx_personality_v0+0x1a9)[0x496959]
======= Memory map: ========
00400000-00cca000 r-xp 00000000 08:01 604869
  /usr/local/grads-2.0.a7.oga.3/bin/grads
00dc9000-00dea000 rwxp 008c9000 08:01 604869
  /usr/local/grads-2.0.a7.oga.3/bin/grads
00dea000-010aa000 rwxp 00dea000 00:00 0
0180f000-01a04000 rwxp 0180f000 00:00 0                                  [heap]
2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056
  /lib/ld-2.9.so
2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0


From jimp at HAWAII.EDU  Tue Feb  9 13:29:30 2010
From: jimp at HAWAII.EDU (James T. Potemra)
Date: Tue, 9 Feb 2010 08:29:30 -1000
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <e71668fa1002090912g37c33699gfb17084afd734754@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.082930.1000.>

Henrique (and Arlindo):

This problem is related to my post back in May 2009.  I have similar 
issues with the
data set below.  If I try the "sdfopen" exactly as below, using both 
grads-2.0.a7.oga.2
and grads-2.0.a5.oga.5 I get the following (followed by grads exiting):

*** glibc detected *** 
/usr/local/grads-2.0.a7.oga.2/Contents/Linux/Versions/2.0.a7.oga.2/x86_64/grads: 
free(): invalid pointer: 0x0000000004dfdff4 ***
======= Backtrace: =========
.... <cut>

If instead I "sdfopen" with grads-1.9.0-rc1 (gradsdods), grads-2.0.a5 
(gradsdap), or
grads-2.0.a6 (gradsdap), all works well. 

Back in May when I raised this issue Arlindo suspected it had something 
to do with
integer size (in my case I was reading netCDF files, not grib).

Jim

PS I'm running on a 64-bit linux, redhat EL-5 (2.6.18-128).

Henrique Barbosa wrote:
> Dear All,
>
> I need help to access CFSR (NCEP reanalysis 3) over opendap with opengrads.
> I get an horrendous error message (see below) when I try to:
>
> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>
> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
> have opendap enabled (see output of 'q config' below).
>
> The strange thing is that I can open other sources of data, for example
>
> ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z
>
> works perfectly and I can plot, subset, etc...
>
> Any ideas?
>
> []s
> Henrique
>
>
> --------------------------------------------------------------------------------------------------------------------------------------------------
> ga-> q config
> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
> hdf4-sds hdf5 opendap-grids,stn athena geotiff
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2009 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> This program is distributed WITHOUT ANY WARRANTY
> See file COPYRIGHT for more information.
>
> Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu
>
> This version of GrADS has been configured with the following options:
>   o Built on a LITTLE ENDIAN machine
>   o Command line editing ENABLED
>       http://tiswww.case.edu/php/chet/readline/rltop.html
>   o printim command for image output ENABLED
>       http://www.zlib.net
>       http://www.libpng.org/pub/png/libpng.html
>       http://www.libgd.org/Main_Page
>   o GRIB2 interface ENABLED
>       http://www.ijg.org
>       http://www.ece.uvic.ca/~mdadams/jasper
>       http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
>       g2clib-1.0.5
>   o NetCDF interface ENABLED
>       http://www.opendap.org
>       libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar  3 2009 14:13:33 $
>   o HDF interface ENABLED
>       http://hdfgroup.org
>       HDF 4.2r3
>       HDF5 1.8.2
>   o Athena Widget GUI ENABLED
>   o OPeNDAP gridded data interface ENABLED
>       http://www.opendap.org
>       libdap 3.7.10
>   o OPeNDAP station data interface ENABLED
>       http://iges.org/grads/gadoc/supplibs.html
>       libgadap 2.0.oga.1
>   o GeoTIFF and KML output ENABLED
>       http://www.libtiff.org
>       http://geotiff.osgeo.org
>
> For additional information please consult http://iges.org/grads
>
> --------------------------------------------------------------------------------------------------------------------------------------------------
> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2
> Scanning self-describing file:
> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
> *** glibc detected *** grads: free(): invalid pointer: 0x00000000018ff354 ***
> ======= Backtrace: =========
> /lib/libc.so.6[0x2b6884e9fcb8]
> /lib/libc.so.6(cfree+0x76)[0x2b6884ea2276]
> grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea]
> grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9)[0x833a05]
> grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea]
> grads(nc3_get_vara_double+0x9e)[0x80e3b8]
> grads(nc_get_vara_double+0x1a)[0x60bcca]
> grads(read_one_dimension+0x89)[0x541ac1]
> grads(gadsdf+0x17d6)[0x540080]
> grads(gasdfopen+0x1d4)[0x53e700]
> grads(gacmd+0xfb2)[0x524504]
> grads(Main+0x4e1)[0x496f01]
> grads(main+0x2e)[0x496a16]
> /lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6]
> grads(__gxx_personality_v0+0x1a9)[0x496959]
> ======= Memory map: ========
> 00400000-00cca000 r-xp 00000000 08:01 604869
>   /usr/local/grads-2.0.a7.oga.3/bin/grads
> 00dc9000-00dea000 rwxp 008c9000 08:01 604869
>   /usr/local/grads-2.0.a7.oga.3/bin/grads
> 00dea000-010aa000 rwxp 00dea000 00:00 0
> 0180f000-01a04000 rwxp 0180f000 00:00 0                                  [heap]
> 2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056
>   /lib/ld-2.9.so
> 2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0
>   


From Dan.Swank at NOAA.GOV  Tue Feb  9 15:34:54 2010
From: Dan.Swank at NOAA.GOV (Dan Swank)
Date: Tue, 9 Feb 2010 15:34:54 -0500
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <e71668fa1002090912g37c33699gfb17084afd734754@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.153454.0500.>

Hello,

GrADS in general is not ready for the CFSR dataset.
I was able to use gradsdap alpha5 to sdfopen the source
  you listed ~ but this is a THREDDS server, which is
  presenting variable names which are too long for
  grads to handle.  [(q file) snipit]:

     relative_humidi  0  t,y,x  Relative humidity @ sigma
     relative_humidi  0  t,y,x  Relative humidity @ sigma_layer
     relative_humidi  0  t,y,x  Relative humidity @ zeroDegC_isotherm
     specific_humidi  0  ty,x  Specific humidity @ pressure
     specific_humidi  0  ty,x  Specific humidity @ pressure_difference_layer
     storm_relative_  0  ty,x  Storm relative helicity @
height_above_ground_layer
     stream_function  0  ty,x  Stream function @ pressure
     surface_lifted_  0  t,y,x  Surface Lifted Index @ surface
     temperature  0  ty,x  Temperature @ pressure
     temperature_alt  4  t,z,y,x  Temperature @ altitude_above_msl
     temperature_pot  0  ty,x  Temperature @ potential_vorticity_surface
     temperature_max  0  t,y,x  Temperature @ maximum_wind

So even if you managed to open it, you wouldn't get much further.

Connecting a GrADS client to a THREDDS server has many known issues.
It will flat out refuse to identify the x/y coordinates when opening
 non-lat/lon grids served from THREDDS, for instance.  Although the
 source you specified IS a lat/lon, some CFSR sets are on
 gaussian coordinates.  The problem there is that THREDDS presents
 the native grid x/y's without projecting it to a lat/lon
 like a GDS does.

Attempting to download a GRIB2 file and creating a control
  file will result in a gribmap2 not being able to MATCH the
  time coordinates because of the new local GRIB table the headers
  use.  This is especially evident in the monthly means.  This
  pretty much stops us from putting the CFSR behind a GDS2.

I hope the GrADS dev-team will take some time to get the
CFSR working in GrADS.

-Dan Swank
NCDC-NOMADS Support




Henrique Barbosa wrote the following on 2/9/2010 12:12 PM:
> Dear All,
> 
> I need help to access CFSR (NCEP reanalysis 3) over opendap with opengrads.
> I get an horrendous error message (see below) when I try to:
> 
> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
> 
> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
> have opendap enabled (see output of 'q config' below).
> 
> The strange thing is that I can open other sources of data, for example
> 
> ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z
> 
> works perfectly and I can plot, subset, etc...
> 
> Any ideas?
> 
> []s
> Henrique
> 
> 
> --------------------------------------------------------------------------------------------------------------------------------------------------
> ga-> q config
> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
> hdf4-sds hdf5 opendap-grids,stn athena geotiff
> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
> Copyright (c) 1988-2009 by Brian Doty and the
> Institute for Global Environment and Society (IGES)
> This program is distributed WITHOUT ANY WARRANTY
> See file COPYRIGHT for more information.
> 
> Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu
> 
> This version of GrADS has been configured with the following options:
>   o Built on a LITTLE ENDIAN machine
>   o Command line editing ENABLED
>       http://tiswww.case.edu/php/chet/readline/rltop.html
>   o printim command for image output ENABLED
>       http://www.zlib.net
>       http://www.libpng.org/pub/png/libpng.html
>       http://www.libgd.org/Main_Page
>   o GRIB2 interface ENABLED
>       http://www.ijg.org
>       http://www.ece.uvic.ca/~mdadams/jasper
>       http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
>       g2clib-1.0.5
>   o NetCDF interface ENABLED
>       http://www.opendap.org
>       libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar  3 2009 14:13:33 $
>   o HDF interface ENABLED
>       http://hdfgroup.org
>       HDF 4.2r3
>       HDF5 1.8.2
>   o Athena Widget GUI ENABLED
>   o OPeNDAP gridded data interface ENABLED
>       http://www.opendap.org
>       libdap 3.7.10
>   o OPeNDAP station data interface ENABLED
>       http://iges.org/grads/gadoc/supplibs.html
>       libgadap 2.0.oga.1
>   o GeoTIFF and KML output ENABLED
>       http://www.libtiff.org
>       http://geotiff.osgeo.org
> 
> For additional information please consult http://iges.org/grads
> 
> --------------------------------------------------------------------------------------------------------------------------------------------------
> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2
> Scanning self-describing file:
> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
> *** glibc detected *** grads: free(): invalid pointer: 0x00000000018ff354 ***
> ======= Backtrace: =========
> /lib/libc.so.6[0x2b6884e9fcb8]
> /lib/libc.so.6(cfree+0x76)[0x2b6884ea2276]
> grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea]
> grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9)[0x833a05]
> grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea]
> grads(nc3_get_vara_double+0x9e)[0x80e3b8]
> grads(nc_get_vara_double+0x1a)[0x60bcca]
> grads(read_one_dimension+0x89)[0x541ac1]
> grads(gadsdf+0x17d6)[0x540080]
> grads(gasdfopen+0x1d4)[0x53e700]
> grads(gacmd+0xfb2)[0x524504]
> grads(Main+0x4e1)[0x496f01]
> grads(main+0x2e)[0x496a16]
> /lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6]
> grads(__gxx_personality_v0+0x1a9)[0x496959]
> ======= Memory map: ========
> 00400000-00cca000 r-xp 00000000 08:01 604869
>   /usr/local/grads-2.0.a7.oga.3/bin/grads
> 00dc9000-00dea000 rwxp 008c9000 08:01 604869
>   /usr/local/grads-2.0.a7.oga.3/bin/grads
> 00dea000-010aa000 rwxp 00dea000 00:00 0
> 0180f000-01a04000 rwxp 0180f000 00:00 0                                  [heap]
> 2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056
>   /lib/ld-2.9.so
> 2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0

-- 
Dan Swank
STG, Incorporated - Government Contractor
NCDC-NOMADS Project:  Software & Data Management
Data Access Branch
National Climatic Data Center
Veach-Baley Federal Building
151 Patton Avenue
Asheville, NC 28801-5001
Email: Dan.Swank at noaa.gov
Phone: 828-271-4007


From claudiomet at GMAIL.COM  Tue Feb  9 16:37:53 2010
From: claudiomet at GMAIL.COM (claudio cortes)
Date: Tue, 9 Feb 2010 18:37:53 -0300
Subject: Grads custom maps
In-Reply-To: <77fcd6b21002090813q787d0271xb8e09a10dcdbb6ef@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.183753.0300.>

Thanks Arlindo.

But, this feature is only implemented in OpenGrADS?
Because I have GrADS v2.0.a7 (on Ubuntu 9.10) and I have "command not found"
when I apply the "shp_polyf" command.

2010/2/9 Arlindo da Silva <arlindo.dasilva at gmail.com>

>
>
> On Mon, Feb 8, 2010 at 11:11 AM, claudio cortes <claudiomet at gmail.com>wrote:
>
>> Greetings
>>
>> The maps included in GrADS only have international political boundaries
>> for my region.
>> Now I need to add regional political boundaries for my country.
>> How can I create custom maps for GrADS ?
>>
>>
>
> One option is to use shapefiles. See this
>
>     http://opengrads.org/doc/udxt/shape/
>
> Also, search this list for "shapefiles".
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu
>



-- 
Claudio Cortes
+56 (2) 2994121

Meteorologo
Laboratorio de Meteorologia y Pronostico de la Calidad del Aire (LMPCA)
Unidad de Modelacion y Gestion de la Calidad del Aire (UMGCA)
Centro Nacional del Medio Ambiente (CENMA)

--
Claudio Cortes
+56 (2) 2994121

Meteorologist
Meteorology and Air Quality Forecast Labs
Modeling and Air Quality Management Unit
National Enviroment Center, Chile (CENMA)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100209/e1571d5c/attachment.html 

From jdduda at IASTATE.EDU  Tue Feb  9 16:55:44 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Tue, 9 Feb 2010 15:55:44 -0600
Subject: Writing out to a file and cdiff
Message-ID: <TUE.9.FEB.2010.155544.0600.>

I have two questions.

1) How do you write data out to a file in text format?  I know there is an
fwrite command, but I don't really know how it works.

2) Can the cdiff command work when three environment dimensions are
varying?  I'm trying to compute EPV, but to do so, I need to do a cdiff
command in the x, y, and z dimensions, but I always get an error and usually
Grads crashes when I do it.

Jeff Duda

-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100209/35e1139e/attachment.html 

From jmssnyder at UCDAVIS.EDU  Tue Feb  9 17:43:45 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Tue, 9 Feb 2010 14:43:45 -0800
Subject: Writing out to a file and cdiff
In-Reply-To: <c0f1d7891002091355p6520c404n989f64981029e3ba@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.144345.0800.>

In regards to your second question would this work?

'sdfopen uwnd.2001.nc'
'define dx=cdiff(uwnd,x)'
'define dxdy=cdiff(dx,y)'
'define dxdydz=cdiff(dxdy,z)'


or something like this.

-Jason

On Tue, Feb 9, 2010 at 1:55 PM, Jeffrey Duda <jdduda at iastate.edu> wrote:

> I have two questions.
>
> 1) How do you write data out to a file in text format?  I know there is an
> fwrite command, but I don't really know how it works.
>
> 2) Can the cdiff command work when three environment dimensions are
> varying?  I'm trying to compute EPV, but to do so, I need to do a cdiff
> command in the x, y, and z dimensions, but I always get an error and usually
> Grads crashes when I do it.
>
> Jeff Duda
>
> --
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100209/650723c6/attachment.html 

From jdduda at IASTATE.EDU  Tue Feb  9 17:50:37 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Tue, 9 Feb 2010 16:50:37 -0600
Subject: Writing out to a file and cdiff
In-Reply-To: <fe8efa01002091443i3c3084do7db1b66d60c22fe@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.165037.0600.>

Jason,
That's not exactly what I'm asking about.  The exact formulation and
computations I'm trying to make are:

'OMEGA = 7.2722e-5'
'pi = 3.14159265358988'
'f = 2*sin(OMEGA*(pi/180))'
'dx = 6.371e6*cos(lat*(pi/180))*cdiff(lon,x)*(pi/180)'
'dy = 6.371e6*cdiff(lat,y)*(pi/180)'
'dz = cdiff((ph+phb)/9.81,z)'

*'icomp = (cdiff(w,y)/dy - cdiff(v,z)/dz + f)*(cdiff(thetae,x)/dx)'
'jcomp = (cdiff(u,z)/dz - cdiff(w,x)/dx + f)*(cdiff(thetae,y)/dy)'
'kcomp = (cdiff(v,x)/dx - cdiff(u,y)/dy + f)*(cdiff(thetae,z)/dz)'*
'EPV = icomp + jcomp + kcomp'

Jeff

On Tue, Feb 9, 2010 at 4:43 PM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> In regards to your second question would this work?
>
> 'sdfopen uwnd.2001.nc'
> 'define dx=cdiff(uwnd,x)'
> 'define dxdy=cdiff(dx,y)'
> 'define dxdydz=cdiff(dxdy,z)'
>
>
> or something like this.
>
> -Jason
>
> On Tue, Feb 9, 2010 at 1:55 PM, Jeffrey Duda <jdduda at iastate.edu> wrote:
>
>> I have two questions.
>>
>> 1) How do you write data out to a file in text format?  I know there is an
>> fwrite command, but I don't really know how it works.
>>
>> 2) Can the cdiff command work when three environment dimensions are
>> varying?  I'm trying to compute EPV, but to do so, I need to do a cdiff
>> command in the x, y, and z dimensions, but I always get an error and usually
>> Grads crashes when I do it.
>>
>> Jeff Duda
>>
>> --
>> Jeff Duda
>> Iowa State University
>> Meteorology Graduate Student
>> 3134 Agronomy Hall
>> www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>
>>
>
>


-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100209/44aec2ca/attachment.html 

From thi_veloso at YAHOO.COM.BR  Tue Feb  9 19:16:16 2010
From: thi_veloso at YAHOO.COM.BR (Thiago Veloso)
Date: Tue, 9 Feb 2010 16:16:16 -0800
Subject: Writing out to a file and cdiff
In-Reply-To: <c0f1d7891002091355p6520c404n989f64981029e3ba@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.161616.0800.>

??Hello, Jeffrey.
??Regardind doubt #2, I have a script which saves the "display" command into a file, instead of?exhibiting the contents in the monitor. Here you are:
'reinit''open /media/Dados/Academico/cmorph/data/2009/CMORPH_8km_hourly.ctl'
itime=1ftime=24while (itime <= 744)'set t 'itime'q time'data=subwrd(result,3)diaf=substr(data,4,2)mesf=substr(data,6,3)anof=substr(data,9,4)
lx=-796while (lx <= -673)ly=357while (ly <= 455)'set x 'lxlon=subwrd(result,4)'set y 'lylat=subwrd(result,4)'d sum(cmorph*0.2,t='itime',t='ftime')'rec=sublin(result,2)prec=subwrd(rec,4)prec1='%3.2f'precf=math_format(prec1,prec)
a=write('prec_'anof''mesf''diaf'.txt', lon';'lat';'precf';'diaf''mesf''anof)ly=ly+1endwhilelx=lx+1endwhileitime=itime+24ftime=ftime+24endwhilesay 'End of the process!'
??Basically there are three loops: longitude, latitude and time. And then just wrote the variables (including coordinates and time) into a ascii file. You can modify it according to your requirements.
??Best regards,
??Thiago.
--- On Tue, 9/2/10, Jeffrey Duda <jdduda at IASTATE.EDU> wrote:

From: Jeffrey Duda <jdduda at IASTATE.EDU>
Subject: Writing out to a file and cdiff
To: GRADSUSR at LIST.CINECA.IT
Date: Tuesday, 9 February, 2010, 19:55

I have two questions.

1) How do you write data out to a file in text format?? I know there is an fwrite command, but I don't really know how it works.

2) Can the cdiff command work when three environment dimensions are varying?? I'm trying to compute EPV, but to do so, I need to do a cdiff command in the x, y, and z dimensions, but I always get an error and usually Grads crashes when I do it.


Jeff Duda

-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda




      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100209/1139c487/attachment.html 

From hersala at GMAIL.COM  Tue Feb  9 20:32:06 2010
From: hersala at GMAIL.COM (hersala)
Date: Tue, 9 Feb 2010 22:32:06 -0300
Subject: Writing out to a file and cdiff
In-Reply-To: <c0f1d7891002091355p6520c404n989f64981029e3ba@mail.gmail.com>
Message-ID: <TUE.9.FEB.2010.223206.0300.>

Hi, Jeffrey
With regard your first question, please give a glance to the following
cookbook:

http://cookbooks.opengrads.org/index.php?title=Recipe-002:_Saving_GrADS_variable_data_to_a_text_file

Hern?n

2010/2/9 Jeffrey Duda <jdduda at iastate.edu>

> I have two questions.
>
> 1) How do you write data out to a file in text format?  I know there is an
> fwrite command, but I don't really know how it works.
>
> 2) Can the cdiff command work when three environment dimensions are
> varying?  I'm trying to compute EPV, but to do so, I need to do a cdiff
> command in the x, y, and z dimensions, but I always get an error and usually
> Grads crashes when I do it.
>
> Jeff Duda
>
> --
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100209/3f5855cc/attachment.html 

From zmumba at YAHOO.COM  Wed Feb 10 11:59:01 2010
From: zmumba at YAHOO.COM (zilore mumba)
Date: Wed, 10 Feb 2010 08:59:01 -0800
Subject: Using CYGWIN GrADS on vista and Windows 7
Message-ID: <WED.10.FEB.2010.085901.0800.>

Hello gentle Grads people,
is there anyone who is able to download cygwin correctly on Vista or Windows 7? I have tried several times, each time startxwin .bat is not picked and the software does not seem to work correctly.
i have tried to run grads with native server, it also does not work on Windows 7. Any information or solution to this (other than going back to XP) would be appreciated>
Zilore




      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100210/2326bb4f/attachment.html 

From fmashingia at YAHOO.COM  Wed Feb 10 12:18:57 2010
From: fmashingia at YAHOO.COM (Fredrick Mashingia)
Date: Wed, 10 Feb 2010 09:18:57 -0800
Subject: Using CYGWIN GrADS on vista and Windows 7
In-Reply-To: <687574.57581.qm@web57008.mail.re3.yahoo.com>
Message-ID: <WED.10.FEB.2010.091857.0800.>

Dear Zilore

You need to?change its setting to window xp compactibility mode.?You can do that while installing or even after installation. Both option works fine with windows vista.

Fredrick


Kupata si werevu, na kukosa si ujinga
(Swahili proverb: Getting something is not necessarily because of cleverness, and missing something is not necessarily because of ignorance)
------------------------------------------------------
Fredrick A.J. Mashingia
PhD?Scholar (catchment hydrology and cover/climate change)
G66,Centre for Water Resources Research
Departiment of Civil Engineering
University Colllege Dublin - UCD
Belfield, Dublin 4, IRELAND.
Mobile: ?+ ...??
E-mail: fmashingia at yahoo.com 
------------------------------------------------------




________________________________
From: zilore mumba <zmumba at YAHOO.COM>
To: GRADSUSR at LIST.CINECA.IT
Sent: Wed, February 10, 2010 4:59:01 PM
Subject: Using CYGWIN GrADS on vista and Windows 7


Hello gentle Grads people,
is there anyone who is able to download cygwin correctly on Vista or Windows 7? I have tried several times, each time startxwin .bat is not picked and the software does not seem to work correctly.
i have tried to run grads with native server, it also does not work on Windows 7. Any information or solution to this (other than going back to XP) would be appreciated>
Zilore


      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100210/794915d1/attachment.html 

From hmjbarbosa at GMAIL.COM  Wed Feb 10 13:59:02 2010
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Wed, 10 Feb 2010 16:59:02 -0200
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <4B71C6EE.9060700@noaa.gov>
Message-ID: <WED.10.FEB.2010.165902.0200.>

Dear Dan Swank,

Thank you for clarifying the problem. The difference between GDS and TDS
servers was not clear to me. I have not mentioned in my previous email,
but I already tried to download the data and write a ctl...
... and I got the time coordinate error you mentioned!

> I hope the GrADS dev-team will take some time to get the
> CFSR working in GrADS.

I can only enforce your request. +1 here.

[]'s
Henrique

On Tue, Feb 9, 2010 at 6:34 PM, Dan Swank <Dan.Swank at noaa.gov> wrote:
> Hello,
>
> GrADS in general is not ready for the CFSR dataset.
> I was able to use gradsdap alpha5 to sdfopen the source
> ?you listed ~ but this is a THREDDS server, which is
> ?presenting variable names which are too long for
> ?grads to handle. ?[(q file) snipit]:
>
> ? ? relative_humidi ?0 ?t,y,x ?Relative humidity @ sigma
> ? ? relative_humidi ?0 ?t,y,x ?Relative humidity @ sigma_layer
> ? ? relative_humidi ?0 ?t,y,x ?Relative humidity @ zeroDegC_isotherm
> ? ? specific_humidi ?0 ?ty,x ?Specific humidity @ pressure
> ? ? specific_humidi ?0 ?ty,x ?Specific humidity @ pressure_difference_layer
> ? ? storm_relative_ ?0 ?ty,x ?Storm relative helicity @
> height_above_ground_layer
> ? ? stream_function ?0 ?ty,x ?Stream function @ pressure
> ? ? surface_lifted_ ?0 ?t,y,x ?Surface Lifted Index @ surface
> ? ? temperature ?0 ?ty,x ?Temperature @ pressure
> ? ? temperature_alt ?4 ?t,z,y,x ?Temperature @ altitude_above_msl
> ? ? temperature_pot ?0 ?ty,x ?Temperature @ potential_vorticity_surface
> ? ? temperature_max ?0 ?t,y,x ?Temperature @ maximum_wind
>
> So even if you managed to open it, you wouldn't get much further.
>
> Connecting a GrADS client to a THREDDS server has many known issues.
> It will flat out refuse to identify the x/y coordinates when opening
> ?non-lat/lon grids served from THREDDS, for instance. ?Although the
> ?source you specified IS a lat/lon, some CFSR sets are on
> ?gaussian coordinates. ?The problem there is that THREDDS presents
> ?the native grid x/y's without projecting it to a lat/lon
> ?like a GDS does.
>
> Attempting to download a GRIB2 file and creating a control
> ?file will result in a gribmap2 not being able to MATCH the
> ?time coordinates because of the new local GRIB table the headers
> ?use. ?This is especially evident in the monthly means. ?This
> ?pretty much stops us from putting the CFSR behind a GDS2.
>
> I hope the GrADS dev-team will take some time to get the
> CFSR working in GrADS.
>
> -Dan Swank
> NCDC-NOMADS Support
>
>
>
>
> Henrique Barbosa wrote the following on 2/9/2010 12:12 PM:
>> Dear All,
>>
>> I need help to access CFSR (NCEP reanalysis 3) over opendap with opengrads.
>> I get an horrendous error message (see below) when I try to:
>>
>> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>>
>> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
>> have opendap enabled (see output of 'q config' below).
>>
>> The strange thing is that I can open other sources of data, for example
>>
>> ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z
>>
>> works perfectly and I can plot, subset, etc...
>>
>> Any ideas?
>>
>> []s
>> Henrique
>>
>>
>> --------------------------------------------------------------------------------------------------------------------------------------------------
>> ga-> q config
>> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
>> hdf4-sds hdf5 opendap-grids,stn athena geotiff
>> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
>> Copyright (c) 1988-2009 by Brian Doty and the
>> Institute for Global Environment and Society (IGES)
>> This program is distributed WITHOUT ANY WARRANTY
>> See file COPYRIGHT for more information.
>>
>> Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu
>>
>> This version of GrADS has been configured with the following options:
>> ? o Built on a LITTLE ENDIAN machine
>> ? o Command line editing ENABLED
>> ? ? ? http://tiswww.case.edu/php/chet/readline/rltop.html
>> ? o printim command for image output ENABLED
>> ? ? ? http://www.zlib.net
>> ? ? ? http://www.libpng.org/pub/png/libpng.html
>> ? ? ? http://www.libgd.org/Main_Page
>> ? o GRIB2 interface ENABLED
>> ? ? ? http://www.ijg.org
>> ? ? ? http://www.ece.uvic.ca/~mdadams/jasper
>> ? ? ? http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
>> ? ? ? g2clib-1.0.5
>> ? o NetCDF interface ENABLED
>> ? ? ? http://www.opendap.org
>> ? ? ? libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar ?3 2009 14:13:33 $
>> ? o HDF interface ENABLED
>> ? ? ? http://hdfgroup.org
>> ? ? ? HDF 4.2r3
>> ? ? ? HDF5 1.8.2
>> ? o Athena Widget GUI ENABLED
>> ? o OPeNDAP gridded data interface ENABLED
>> ? ? ? http://www.opendap.org
>> ? ? ? libdap 3.7.10
>> ? o OPeNDAP station data interface ENABLED
>> ? ? ? http://iges.org/grads/gadoc/supplibs.html
>> ? ? ? libgadap 2.0.oga.1
>> ? o GeoTIFF and KML output ENABLED
>> ? ? ? http://www.libtiff.org
>> ? ? ? http://geotiff.osgeo.org
>>
>> For additional information please consult http://iges.org/grads
>>
>> --------------------------------------------------------------------------------------------------------------------------------------------------
>> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2
>> Scanning self-describing file:
>> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>> *** glibc detected *** grads: free(): invalid pointer: 0x00000000018ff354 ***
>> ======= Backtrace: =========
>> /lib/libc.so.6[0x2b6884e9fcb8]
>> /lib/libc.so.6(cfree+0x76)[0x2b6884ea2276]
>> grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea]
>> grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9)[0x833a05]
>> grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea]
>> grads(nc3_get_vara_double+0x9e)[0x80e3b8]
>> grads(nc_get_vara_double+0x1a)[0x60bcca]
>> grads(read_one_dimension+0x89)[0x541ac1]
>> grads(gadsdf+0x17d6)[0x540080]
>> grads(gasdfopen+0x1d4)[0x53e700]
>> grads(gacmd+0xfb2)[0x524504]
>> grads(Main+0x4e1)[0x496f01]
>> grads(main+0x2e)[0x496a16]
>> /lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6]
>> grads(__gxx_personality_v0+0x1a9)[0x496959]
>> ======= Memory map: ========
>> 00400000-00cca000 r-xp 00000000 08:01 604869
>> ? /usr/local/grads-2.0.a7.oga.3/bin/grads
>> 00dc9000-00dea000 rwxp 008c9000 08:01 604869
>> ? /usr/local/grads-2.0.a7.oga.3/bin/grads
>> 00dea000-010aa000 rwxp 00dea000 00:00 0
>> 0180f000-01a04000 rwxp 0180f000 00:00 0 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?[heap]
>> 2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056
>> ? /lib/ld-2.9.so
>> 2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0
>
> --
> Dan Swank
> STG, Incorporated - Government Contractor
> NCDC-NOMADS Project: ?Software & Data Management
> Data Access Branch
> National Climatic Data Center
> Veach-Baley Federal Building
> 151 Patton Avenue
> Asheville, NC 28801-5001
> Email: Dan.Swank at noaa.gov
> Phone: 828-271-4007
>


From luis.blacutt at GMAIL.COM  Wed Feb 10 15:08:46 2010
From: luis.blacutt at GMAIL.COM (luis blacutt)
Date: Wed, 10 Feb 2010 17:08:46 -0300
Subject: Converting outputs to netcdf
Message-ID: <WED.10.FEB.2010.170846.0300.>

Dear grads users,

I would like to save the display directly to netcdf. I'm already aware that
there is a way to convert to a txt file, through the fprintf function, or
the set gxout print command.

Regards
Luis
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100210/f5a7f50f/attachment.html 

From krishmet at GMAIL.COM  Wed Feb 10 23:20:07 2010
From: krishmet at GMAIL.COM (krishnamohan)
Date: Thu, 11 Feb 2010 09:50:07 +0530
Subject: Converting outputs to netcdf
In-Reply-To: <5825ce3f1002101208l32509bccsa782b145db7fd2f@mail.gmail.com>
Message-ID: <THU.11.FEB.2010.095007.0530.>

Hai Luis,

       You can use grads sdfwrite function  for writing netcdf output. The
sdfwrite command works with GrADS version 2.0.a3 (or higher). The example as
given in GrADS documentation index is

open global_forecast.ctl
set lon -111.3 -103.8
set lat 40.9 45.0
set lev 1000 10
set t 1 last
set e 1 last
define precip = ptot
set sdfwrite wyoming_precip.nc
sdfwrite precip

this will write output to the netcdf file wyoming_precip.nc


Krishnamohan.K.S
Junior Research Fellow
Department Of Atmospheric Sciences
Cochin University of Science and Technology
Cochin,India


On Thu, Feb 11, 2010 at 1:38 AM, luis blacutt <luis.blacutt at gmail.com>wrote:

> Dear grads users,
>
> I would like to save the display directly to netcdf. I'm already aware that
> there is a way to convert to a txt file, through the fprintf function, or
> the set gxout print command.
>
> Regards
> Luis
>



--
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100211/2d70fb6b/attachment.html 

From zmumba at YAHOO.COM  Thu Feb 11 11:06:37 2010
From: zmumba at YAHOO.COM (zilore mumba)
Date: Thu, 11 Feb 2010 08:06:37 -0800
Subject: Using CYGWIN GrADS on vista and Windows 7
In-Reply-To: <389005.16675.qm@web51503.mail.re2.yahoo.com>
Message-ID: <THU.11.FEB.2010.080637.0800.>

Fredrick,
thanks very much for your assistance. But where exactly do I change the setting?
Zilore



--- On Wed, 2/10/10, Fredrick Mashingia <fmashingia at YAHOO.COM> wrote:

From: Fredrick Mashingia <fmashingia at YAHOO.COM>
Subject: Re: Using CYGWIN GrADS on vista and Windows 7
To: GRADSUSR at LIST.CINECA.IT
Date: Wednesday, February 10, 2010, 7:18 PM


Dear Zilore

You need to?change its setting to window xp compactibility mode.?You can do that while installing or even after installation. Both option works fine with windows vista.
?
Fredrick
?
?

Kupata si werevu, na kukosa si ujinga
(Swahili proverb: Getting something is not necessarily because of cleverness, and missing something is not necessarily because of ignorance)

------------------------------------------------------
Fredrick A.J. Mashingia
PhD?Scholar (catchment hydrology and cover/climate change)
G66,Centre for Water Resources Research
Departiment of Civil Engineering
University Colllege Dublin - UCD
Belfield, Dublin 4, IRELAND.
Mobile: ?+ ...??
E-mail: fmashingia at yahoo.com 

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






From: zilore mumba <zmumba at YAHOO.COM>
To: GRADSUSR at LIST.CINECA.IT
Sent: Wed, February 10, 2010 4:59:01 PM
Subject: Using CYGWIN GrADS on vista and Windows 7





Hello gentle Grads people,
is there anyone who is able to download cygwin correctly on Vista or Windows 7? I have tried several times, each time startxwin .bat is not picked and the software does not seem to work correctly.
i have tried to run grads with native server, it also does not work on Windows 7. Any information or solution to this (other than going back to XP) would be appreciated>
Zilore





      


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

From c.hayes at WEERONLINE.NL  Thu Feb 11 11:43:52 2010
From: c.hayes at WEERONLINE.NL (Chris Hayes)
Date: Thu, 11 Feb 2010 17:43:52 +0100
Subject: grib to grib
Message-ID: <THU.11.FEB.2010.174352.0100.>

Dear group

I would like to use grib files (grib 1) to create new grib (grib 1) files,
with hourly change values instead of accumulated values.

 

E.g.  grib_03:00 contains accumulative_precipitation value 12

And grib_04:00 contains accumulative_precipitation value 15

Then I would like to create a new grib file with hourly_precipitation value
3 (15-12).

 

The reason I do not just open them both  and do the math to create a map (I
asked that question earlier) is that we would like to use the new grib file
for other purposes as well.

 

An example of this would be very helpful.

 

Sincerely,

Chris

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100211/3bd27de0/attachment.html 

From jdduda at IASTATE.EDU  Thu Feb 11 13:43:38 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Thu, 11 Feb 2010 12:43:38 -0600
Subject: grib to grib
In-Reply-To: <237769565863322513@unknownmsgid>
Message-ID: <THU.11.FEB.2010.124338.0600.>

Look into using pcp_combine, a tool that comes with the MET package that
includes MODE. (
http://www.dtcenter.org/met/users/support/online_tutorial/METv2.0/pcp_combine/index.php)
I use it to change the run-total accumulated precip in my WRF runs into 1-hr
accumulated precip using the argument pcp_combine -subtract (WRF file 1)
(hour of the run) (WRF file 2) (hour of the run) (name of output file)

Jeff Duda

On Thu, Feb 11, 2010 at 10:43 AM, Chris Hayes <c.hayes at weeronline.nl> wrote:

>  Dear group
>
> I would like to use grib files (grib 1) to create new grib (grib 1) files,
> with hourly change values instead of accumulated values.
>
>
>
> E.g.  grib_03:00 contains accumulative_precipitation value 12
>
> And grib_04:00 contains accumulative_precipitation value 15
>
> Then I would like to create a new grib file with hourly_precipitation value
> 3 (15-12).
>
>
>
> The reason I do not just open them both  and do the math to create a map (I
> asked that question earlier) is that we would like to use the new grib file
> for other purposes as well.
>
>
>
> An example of this would be very helpful.
>
>
>
> Sincerely,
>
> Chris
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100211/add19477/attachment.html 

From remi.meynadier at AERO.JUSSIEU.FR  Thu Feb 11 14:27:40 2010
From: remi.meynadier at AERO.JUSSIEU.FR (RM)
Date: Thu, 11 Feb 2010 20:27:40 +0100
Subject: hdvig script
Message-ID: <THU.11.FEB.2010.202740.0100.>

Hi grads's users,
I have a question concerning the way hdivg function is done.
I looked at the script code and it seems that the finite difference method
does not exactly correspond to the horizontal divergence expression in
spherical coordinates. More specifically the term in tan(latitude) is ignored.
Is this term negligible ?

Thanks in advance
RM


From jm08j at FSU.EDU  Thu Feb 11 14:44:45 2010
From: jm08j at FSU.EDU (J-P Michael)
Date: Thu, 11 Feb 2010 14:44:45 -0500
Subject: Writing out to a file and cdiff
In-Reply-To: <c0f1d7891002091355p6520c404n989f64981029e3ba@mail.gmail.com>
Message-ID: <THU.11.FEB.2010.144445.0500.>

Jeff, 
This may be a foregone conclusion for most but there is some issue with cdiff in z, depending on the version of GrADS you are using. Are you able to cdiff in x and y and failing in the z? 

J-P Michael
Florida State University

----- Original Message -----
From: Jeffrey Duda <jdduda at IASTATE.EDU>
Date: Tuesday, February 9, 2010 4:55 pm
Subject: Writing out to a file and cdiff
To: GRADSUSR at LIST.CINECA.IT

> I have two questions.
> 
> 1) How do you write data out to a file in text format?  I know 
> there is an
> fwrite command, but I don't really know how it works.
> 
> 2) Can the cdiff command work when three environment dimensions are
> varying?  I'm trying to compute EPV, but to do so, I need to do a 
> cdiffcommand in the x, y, and z dimensions, but I always get an 
> error and usually
> Grads crashes when I do it.
> 
> Jeff Duda
> 
> -- 
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda
> 


From jdduda at IASTATE.EDU  Thu Feb 11 15:05:11 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Thu, 11 Feb 2010 14:05:11 -0600
Subject: Writing out to a file and cdiff
In-Reply-To: <f6f3c0124668.4b7417dd@fsu.edu>
Message-ID: <THU.11.FEB.2010.140511.0600.>

Yeah.  I can't define a product that involves the cdiff command in the
z-dimension.  However, for some odd reason, I can display a product using
the cdiff command in the z-dimension.  What I mean is, I can't have the
following INSIDE a grads script:

'define dz = cdiff(hght,z)'

However, If I wanted to display something like vertical wind shear, I could
do so by typing in the command window:

'display cdiff(u,z)/cdiff(hght,z)'

and it would display just fine.  I'm not sure why this is.  I'm running
version 2.0.a5 though, so maybe it's a problem with the version.

Jeff

On Thu, Feb 11, 2010 at 1:44 PM, J-P Michael <jm08j at fsu.edu> wrote:

> Jeff,
> This may be a foregone conclusion for most but there is some issue with
> cdiff in z, depending on the version of GrADS you are using. Are you able to
> cdiff in x and y and failing in the z?
>
> J-P Michael
> Florida State University
>
> ----- Original Message -----
> From: Jeffrey Duda <jdduda at IASTATE.EDU>
> Date: Tuesday, February 9, 2010 4:55 pm
> Subject: Writing out to a file and cdiff
> To: GRADSUSR at LIST.CINECA.IT
>
> > I have two questions.
> >
> > 1) How do you write data out to a file in text format?  I know
> > there is an
> > fwrite command, but I don't really know how it works.
> >
> > 2) Can the cdiff command work when three environment dimensions are
> > varying?  I'm trying to compute EPV, but to do so, I need to do a
> > cdiffcommand in the x, y, and z dimensions, but I always get an
> > error and usually
> > Grads crashes when I do it.
> >
> > Jeff Duda
> >
> > --
> > Jeff Duda
> > Iowa State University
> > Meteorology Graduate Student
> > 3134 Agronomy Hall
> > www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>
> >
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100211/043d7bf6/attachment.html 

From fmashingia at YAHOO.COM  Thu Feb 11 16:31:34 2010
From: fmashingia at YAHOO.COM (Fredrick Mashingia)
Date: Thu, 11 Feb 2010 13:31:34 -0800
Subject: Using CYGWIN GrADS on vista and Windows 7
In-Reply-To: <614919.79976.qm@web57005.mail.re3.yahoo.com>
Message-ID: <THU.11.FEB.2010.133134.0800.>

Zilore

If its installed,right-click the cygwin shortcut from the the desktop,go to properties,compactibility and select windows xp. It should work. Otherwise go to control panel..programs and features..run program made for previous versions..

Happy?working

Fredrick
?
Kupata si werevu, na kukosa si ujinga
(Swahili proverb: Getting something is not necessarily because of cleverness, and missing something is not necessarily because of ignorance)
------------------------------------------------------
Fredrick A.J. Mashingia
PhD?Scholar (catchment hydrology and cover/climate change)
G66,Centre for Water Resources Research
Departiment of Civil Engineering
University Colllege Dublin - UCD
Belfield, Dublin 4, IRELAND.
Mobile: ?+ ...??
E-mail: fmashingia at yahoo.com 
------------------------------------------------------




________________________________
From: zilore mumba <zmumba at YAHOO.COM>
To: GRADSUSR at LIST.CINECA.IT
Sent: Thu, February 11, 2010 4:06:37 PM
Subject: Re: Using CYGWIN GrADS on vista and Windows 7


Fredrick,
thanks very much for your assistance. But where exactly do I change the setting?
Zilore



--- On Wed, 2/10/10, Fredrick Mashingia <fmashingia at YAHOO.COM> wrote:


>From: Fredrick Mashingia <fmashingia at YAHOO.COM>
>Subject: Re: Using CYGWIN GrADS on vista and Windows 7
>To: GRADSUSR at LIST.CINECA.IT
>Date: Wednesday, February 10, 2010, 7:18 PM
>
>
>Dear Zilore
>
>You need to?change its setting to window xp compactibility mode.?You can do that while installing or even after installation. Both option works fine with windows vista.
>
>Fredrick
>
>
>Kupata si werevu, na kukosa si ujinga
>(Swahili proverb: Getting something is not necessarily because of cleverness, and missing something is not necessarily because of ignorance)
>------------------------------------------------------
>Fredrick A.J. Mashingia
>PhD?Scholar (catchment hydrology and cover/climate change)
>G66,Centre for Water Resources Research
>Departiment of Civil Engineering
>University Colllege Dublin - UCD
>Belfield, Dublin 4, IRELAND.
>Mobile: ?+ ...??
>E-mail: fmashingia at yahoo.com 
>------------------------------------------------------
>
>
>
>
>
________________________________
From: zilore mumba <zmumba at YAHOO.COM>
>To: GRADSUSR at LIST.CINECA.IT
>Sent: Wed, February 10, 2010 4:59:01 PM
>Subject: Using CYGWIN GrADS on vista and Windows 7
>
>
>Hello gentle Grads people,
>is there anyone who is able to download cygwin correctly on Vista or Windows 7? I have tried several times, each time startxwin .bat is not picked and the software does not seem to work correctly.
>i have tried to run grads with native server, it also does not work on Windows 7. Any information or solution to this (other than going back to XP) would be appreciated>
>Zilore
>
> 
>
> 



      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100211/50a62ef0/attachment.html 

From Charles.Seman at NOAA.GOV  Thu Feb 11 17:34:55 2010
From: Charles.Seman at NOAA.GOV (Charles Seman)
Date: Thu, 11 Feb 2010 17:34:55 -0500
Subject: Writing out to a file and cdiff
In-Reply-To: <c0f1d7891002111205x7f4a22dm9efad6cf7d11c374@mail.gmail.com>
Message-ID: <THU.11.FEB.2010.173455.0500.>

Following what Davide Sacchetti sent to the Listserv a number of years 
ago, try something like:

set z 2 nz-1
define dvardz = (var(z+1)-var(z-1))/(lev(z+1)-lev(z-1))

if your "zdef" vertical coordinate is "z"; or else if say pressure is 
your "zdef" vertical coordinate and you have zht(x,y,z) as a GrADS 
variable, then substitute (zht(z+1)-zht(z-1)) for centered vertical 
differences...

Hope this helps,
Chuck

Jeffrey Duda wrote:
> Yeah.  I can't define a product that involves the cdiff command in the 
> z-dimension.  However, for some odd reason, I can display a product 
> using the cdiff command in the z-dimension.  What I mean is, I can't 
> have the following INSIDE a grads script:
>
> 'define dz = cdiff(hght,z)'
>
> However, If I wanted to display something like vertical wind shear, I 
> could do so by typing in the command window:
>
> 'display cdiff(u,z)/cdiff(hght,z)'
>
> and it would display just fine.  I'm not sure why this is.  I'm 
> running version 2.0.a5 though, so maybe it's a problem with the version.
>
> Jeff
>
> On Thu, Feb 11, 2010 at 1:44 PM, J-P Michael <jm08j at fsu.edu 
> <mailto:jm08j at fsu.edu>> wrote:
>
>     Jeff,
>     This may be a foregone conclusion for most but there is some issue
>     with cdiff in z, depending on the version of GrADS you are using.
>     Are you able to cdiff in x and y and failing in the z?
>
>     J-P Michael
>     Florida State University
>
>     ----- Original Message -----
>     From: Jeffrey Duda <jdduda at IASTATE.EDU <mailto:jdduda at IASTATE.EDU>>
>     Date: Tuesday, February 9, 2010 4:55 pm
>     Subject: Writing out to a file and cdiff
>     To: GRADSUSR at LIST.CINECA.IT <mailto:GRADSUSR at LIST.CINECA.IT>
>
>     > I have two questions.
>     >
>     > 1) How do you write data out to a file in text format?  I know
>     > there is an
>     > fwrite command, but I don't really know how it works.
>     >
>     > 2) Can the cdiff command work when three environment dimensions are
>     > varying?  I'm trying to compute EPV, but to do so, I need to do a
>     > cdiffcommand in the x, y, and z dimensions, but I always get an
>     > error and usually
>     > Grads crashes when I do it.
>     >
>     > Jeff Duda
>     >
>     > --
>     > Jeff Duda
>     > Iowa State University
>     > Meteorology Graduate Student
>     > 3134 Agronomy Hall
>     > www.meteor.iastate.edu/~jdduda
>     <http://www.meteor.iastate.edu/%7Ejdduda>
>     >
>
>
>
>
> -- 
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>

-- 

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 Dan.Swank at NOAA.GOV  Fri Feb 12 11:52:30 2010
From: Dan.Swank at NOAA.GOV (Dan Swank)
Date: Fri, 12 Feb 2010 11:52:30 -0500
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <e71668fa1002101059o440a6b4fx57df2ec423f9ca38@mail.gmail.com>
Message-ID: <FRI.12.FEB.2010.115230.0500.>

Hello,

I need to take back a few statements I made in
my previous post.

I have gotten several CFSR/CFSRMON files to plot successfully.
Here are some samples:

http://nomads.ncdc.noaa.gov/data/cfsr/200001/wnd10m.gdas.200001.grb2
http://nomads.ncdc.noaa.gov/data/cfsr/200001/wnd10m.gdas.200001.ctl
http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2.idx

http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2
http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.ctl
http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2.idx

http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.ctl
http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.grb2.idx
http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.grb2

http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.grb2
http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.ctl
http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.grb2.idx

The -b should be used when invoking GrADS2 gribmap.
The -0 opt should be set when working with monthly means
  (cfsrmon) and diurnal fraction monthlies.

Here are the known gotchas:
wgrib2 v.1.8.2 or earlier may produce control files
with incorrect TDEF's when used with g2ctl.pl .

The original CFSR was produced in GRIB1 at NCEP.
NCEP converted this data to GRIB2 for transmission
and storage purposes.  There are possible issues
in the translation between GRIB1 and GRIB2 [indicator
of statistical process] header.  This affects monthly
means (which are all statistically processed in some way).
gribmap's -0 option nulls this *potential* issue by matching
using the reference time.  Again, this issue is unproven.

The GRIB2 headers contain local extensions
used by Table 4.10 (-fix_ncep wgrib2 option).
I had thought this might cause issues with gribmap
  but am now not certain.


-Dan




Henrique Barbosa wrote the following on 2/10/2010 1:59 PM:
> Dear Dan Swank,
> 
> Thank you for clarifying the problem. The difference between GDS and TDS
> servers was not clear to me. I have not mentioned in my previous email,
> but I already tried to download the data and write a ctl...
> ... and I got the time coordinate error you mentioned!
> 
>> I hope the GrADS dev-team will take some time to get the
>> CFSR working in GrADS.
> 
> I can only enforce your request. +1 here.
> 
> []'s
> Henrique
> 
> On Tue, Feb 9, 2010 at 6:34 PM, Dan Swank <Dan.Swank at noaa.gov> wrote:
>> Hello,
>>
>> GrADS in general is not ready for the CFSR dataset.
>> I was able to use gradsdap alpha5 to sdfopen the source
>>  you listed ~ but this is a THREDDS server, which is
>>  presenting variable names which are too long for
>>  grads to handle.  [(q file) snipit]:
>>
>>     relative_humidi  0  t,y,x  Relative humidity @ sigma
>>     relative_humidi  0  t,y,x  Relative humidity @ sigma_layer
>>     relative_humidi  0  t,y,x  Relative humidity @ zeroDegC_isotherm
>>     specific_humidi  0  ty,x  Specific humidity @ pressure
>>     specific_humidi  0  ty,x  Specific humidity @ pressure_difference_layer
>>     storm_relative_  0  ty,x  Storm relative helicity @
>> height_above_ground_layer
>>     stream_function  0  ty,x  Stream function @ pressure
>>     surface_lifted_  0  t,y,x  Surface Lifted Index @ surface
>>     temperature  0  ty,x  Temperature @ pressure
>>     temperature_alt  4  t,z,y,x  Temperature @ altitude_above_msl
>>     temperature_pot  0  ty,x  Temperature @ potential_vorticity_surface
>>     temperature_max  0  t,y,x  Temperature @ maximum_wind
>>
>> So even if you managed to open it, you wouldn't get much further.
>>
>> Connecting a GrADS client to a THREDDS server has many known issues.
>> It will flat out refuse to identify the x/y coordinates when opening
>>  non-lat/lon grids served from THREDDS, for instance.  Although the
>>  source you specified IS a lat/lon, some CFSR sets are on
>>  gaussian coordinates.  The problem there is that THREDDS presents
>>  the native grid x/y's without projecting it to a lat/lon
>>  like a GDS does.
>>
>> Attempting to download a GRIB2 file and creating a control
>>  file will result in a gribmap2 not being able to MATCH the
>>  time coordinates because of the new local GRIB table the headers
>>  use.  This is especially evident in the monthly means.  This
>>  pretty much stops us from putting the CFSR behind a GDS2.
>>
>> I hope the GrADS dev-team will take some time to get the
>> CFSR working in GrADS.
>>
>> -Dan Swank
>> NCDC-NOMADS Support
>>
>>
>>
>>
>> Henrique Barbosa wrote the following on 2/9/2010 12:12 PM:
>>> Dear All,
>>>
>>> I need help to access CFSR (NCEP reanalysis 3) over opendap with opengrads.
>>> I get an horrendous error message (see below) when I try to:
>>>
>>> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>>>
>>> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
>>> have opendap enabled (see output of 'q config' below).
>>>
>>> The strange thing is that I can open other sources of data, for example
>>>
>>> ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z
>>>
>>> works perfectly and I can plot, subset, etc...
>>>
>>> Any ideas?
>>>
>>> []s
>>> Henrique
>>>
>>>
>>> --------------------------------------------------------------------------------------------------------------------------------------------------
>>> ga-> q config
>>> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
>>> hdf4-sds hdf5 opendap-grids,stn athena geotiff
>>> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
>>> Copyright (c) 1988-2009 by Brian Doty and the
>>> Institute for Global Environment and Society (IGES)
>>> This program is distributed WITHOUT ANY WARRANTY
>>> See file COPYRIGHT for more information.
>>>
>>> Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu
>>>
>>> This version of GrADS has been configured with the following options:
>>>   o Built on a LITTLE ENDIAN machine
>>>   o Command line editing ENABLED
>>>       http://tiswww.case.edu/php/chet/readline/rltop.html
>>>   o printim command for image output ENABLED
>>>       http://www.zlib.net
>>>       http://www.libpng.org/pub/png/libpng.html
>>>       http://www.libgd.org/Main_Page
>>>   o GRIB2 interface ENABLED
>>>       http://www.ijg.org
>>>       http://www.ece.uvic.ca/~mdadams/jasper
>>>       http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
>>>       g2clib-1.0.5
>>>   o NetCDF interface ENABLED
>>>       http://www.opendap.org
>>>       libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar  3 2009 14:13:33 $
>>>   o HDF interface ENABLED
>>>       http://hdfgroup.org
>>>       HDF 4.2r3
>>>       HDF5 1.8.2
>>>   o Athena Widget GUI ENABLED
>>>   o OPeNDAP gridded data interface ENABLED
>>>       http://www.opendap.org
>>>       libdap 3.7.10
>>>   o OPeNDAP station data interface ENABLED
>>>       http://iges.org/grads/gadoc/supplibs.html
>>>       libgadap 2.0.oga.1
>>>   o GeoTIFF and KML output ENABLED
>>>       http://www.libtiff.org
>>>       http://geotiff.osgeo.org
>>>
>>> For additional information please consult http://iges.org/grads
>>>
>>> --------------------------------------------------------------------------------------------------------------------------------------------------
>>> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2
>>> Scanning self-describing file:
>>> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>>> *** glibc detected *** grads: free(): invalid pointer: 0x00000000018ff354 ***
>>> ======= Backtrace: =========
>>> /lib/libc.so.6[0x2b6884e9fcb8]
>>> /lib/libc.so.6(cfree+0x76)[0x2b6884ea2276]
>>> grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea]
>>> grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9)[0x833a05]
>>> grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea]
>>> grads(nc3_get_vara_double+0x9e)[0x80e3b8]
>>> grads(nc_get_vara_double+0x1a)[0x60bcca]
>>> grads(read_one_dimension+0x89)[0x541ac1]
>>> grads(gadsdf+0x17d6)[0x540080]
>>> grads(gasdfopen+0x1d4)[0x53e700]
>>> grads(gacmd+0xfb2)[0x524504]
>>> grads(Main+0x4e1)[0x496f01]
>>> grads(main+0x2e)[0x496a16]
>>> /lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6]
>>> grads(__gxx_personality_v0+0x1a9)[0x496959]
>>> ======= Memory map: ========
>>> 00400000-00cca000 r-xp 00000000 08:01 604869
>>>   /usr/local/grads-2.0.a7.oga.3/bin/grads
>>> 00dc9000-00dea000 rwxp 008c9000 08:01 604869
>>>   /usr/local/grads-2.0.a7.oga.3/bin/grads
>>> 00dea000-010aa000 rwxp 00dea000 00:00 0
>>> 0180f000-01a04000 rwxp 0180f000 00:00 0                                  [heap]
>>> 2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056
>>>   /lib/ld-2.9.so
>>> 2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0
>> --
>> Dan Swank
>> STG, Incorporated - Government Contractor
>> NCDC-NOMADS Project:  Software & Data Management
>> Data Access Branch
>> National Climatic Data Center
>> Veach-Baley Federal Building
>> 151 Patton Avenue
>> Asheville, NC 28801-5001
>> Email: Dan.Swank at noaa.gov
>> Phone: 828-271-4007
>>

-- 
Dan Swank
STG, Incorporated - Government Contractor
NCDC-NOMADS Project:  Software & Data Management
Data Access Branch
National Climatic Data Center
Veach-Baley Federal Building
151 Patton Avenue
Asheville, NC 28801-5001
Email: Dan.Swank at noaa.gov
Phone: 828-271-4007


From dasilva at ALUM.MIT.EDU  Fri Feb 12 12:31:09 2010
From: dasilva at ALUM.MIT.EDU (Arlindo da Silva)
Date: Fri, 12 Feb 2010 12:31:09 -0500
Subject: Fwd: WIN GRADS 2.0.a7.oga.3 - Possible Bug
In-Reply-To: <1347555851.2122391265820197981.JavaMail.root@sz0052a.emeryville.ca.mail.comcast.net>
Message-ID: <FRI.12.FEB.2010.123109.0500.>

Dear Win32 Users,

   If you find yourself in a (rare) situation when

ga-> ! ls

does not work, read the thread below for a way to repair your installation.
This problem happens when mixing both v1.9 and v2.0, and either version is
exited uncleanly.  I'll add this information to the Troubleshooting sections
of the wiki.

   Arlindo



---------- Forwarded message ----------
From: Mark Sponsler <msponsler at comcast.net>
Date: Wed, Feb 10, 2010 at 11:43 AM
Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug
To: Arlindo da Silva <dasilva at alum.mit.edu>


Arlindo,

Think I got it.  Took a bit of work but figured out I needed to stay in bash
AND use short file names and use the"./" before many commands:



sh-3.2$ cd C:/PROGRA~1/opengr~1/Contents/Cygwin/Versions/2.0.a7.oga.2/i686

*Once I was able to CD to the Grads2 directory I could not unmount the /bin
directory, But removing all mounts seemed to work  *



sh-3.2$ ./umount /bin
umount: /bin: No such file or directory

sh-3.2$ ./umount -U
sh-3.2$ ./umount -A



*Then when I did an 'ls' I was still in the directory where grads lived.
Making progress
*

sh-3.2$ ./ls
ADVAPI32.DLL                                   SHELL32.dll
cygpopt-0.dll             libmf.pod
API-MS-WIN-Service-Core-L1-1-0.dll             SHLWAPI.dll
cygreadline6.dll          ln.exe
API-MS-WIN-Service-Management-L1-1-0.dll       USER32.dll
cygssl-0.9.8.dll          ls.exe
API-MS-WIN-Service-Management-L2-1-0.dll       USP10.dll
cygwin1.dll               merra.exe
API-MS-WIN-Service-winsvc-L1-1-0.dll           VERSION
cygxcb-1.dll              mkdir.exe
API-MS-Win-Core-Console-L1-1-0.dll             WINSPOOL.DRV
cygxcb-render-0.dll       mount.exe
API-MS-Win-Core-DateTime-L1-1-0.dll            awk.exe
cygxcb-render-util-0.dll  msvcrt.dll
API-MS-Win-Core-Debug-L1-1-0.dll               bufrscan.exe
cygxcb-xlib-0.dll         mv.exe
API-MS-Win-Core-DelayLoad-L1-1-0.dll           cat.exe
cygz.dll                  nomads.exe
API-MS-Win-Core-ErrorHandling-L1-1-0.dll       comdlg32.dll
date.exe                  ntdll.dll
API-MS-Win-Core-Fibers-L1-1-0.dll              cp.exe
diff.exe                  opengrads.exe
API-MS-Win-Core-File-L1-1-0.dll                cut.exe
dir.exe                   pngquant.exe
API-MS-Win-Core-Handle-L1-1-0.dll              cygICE-6.dll
dos2unix.exe              pwd.exe
API-MS-Win-Core-Heap-L1-1-0.dll                cygSM-6.dll
env.pod                   re.pod
API-MS-Win-Core-IO-L1-1-0.dll                  cygX11-6.dll
fish.pod                  rm.exe
API-MS-Win-Core-Interlocked-L1-1-0.dll         cygXau-6.dll
gex                       rxvt.exe
API-MS-Win-Core-LibraryLoader-L1-1-0.dll       cygXaw-7.dll
grads.exe                 saakeskus.pod_
API-MS-Win-Core-LocalRegistry-L1-1-0.dll       cygXdmcp-6.dll
gradsdap.exe              sed.exe
API-MS-Win-Core-Localization-L1-1-0.dll        cygXext-6.dll
gradsgui.exe              sh.exe
API-MS-Win-Core-Memory-L1-1-0.dll              cygXmu-6.dll
grep.exe                  shape.pod
API-MS-Win-Core-Misc-L1-1-0.dll                cygXpm-4.dll
grib2scan.exe             shfilt.pod
API-MS-Win-Core-NamedPipe-L1-1-0.dll           cygXrender-1.dll
gribmap.exe               sleep.exe
API-MS-Win-Core-ProcessEnvironment-L1-1-0.dll  cygXt-6.dll
gribscan.exe              sort.exe
API-MS-Win-Core-ProcessThreads-L1-1-0.dll      cygcairo-2.dll
gsudf.pod                 stnmap.exe
API-MS-Win-Core-Profile-L1-1-0.dll             cygcrypto-0.9.8.dll
gv32.exe                  tail.exe
API-MS-Win-Core-RtlSupport-L1-1-0.dll          cygexpat-1.dll
gv32.hlp                  tar.exe
API-MS-Win-Core-String-L1-1-0.dll              cygfontconfig-1.dll
gxeps.exe                 tcsh.exe
API-MS-Win-Core-Synch-L1-1-0.dll               cygfreetype-6.dll
gxps.exe                  tee.exe
API-MS-Win-Core-SysInfo-L1-1-0.dll             cyggcc_s-1.dll
gxtran.exe                test.exe
API-MS-Win-Core-ThreadPool-L1-1-0.dll          cyggfortran-3.dll
gxyat.exe                 touch.exe
API-MS-Win-Core-Util-L1-1-0.dll                cygglitz-1.dll
gzip.exe                  umount.exe
API-MS-Win-Security-Base-L1-1-0.dll            cygiconv-2.dll
head.exe                  uname.exe
COMCTL32.dll                                   cygintl-2.dll
lats4d.sh                 uniq.exe
GDI32.dll                                      cygintl-3.dll
less.exe                  unix2dos.exe
KERNEL32.dll                                   cygintl-8.dll
libW11.dll                wc.exe
KERNELBASE.dll                                 cygncurses-9.dll
libbjt.pod                wget.exe
LPK.dll                                        cygpixman-1-0.dll
libgrads.dll              wgrib.exe
RPCRT4.dll                                     cygpng12.dll
libipc.pod                wgrib2.exe

*Next I started up Grads*



sh-3.2$ ./grads
Starting X server under c:\PROGRA~1\opengr~1\Contents\Resources\Xming
User-mode mounting c:\PROGRA~1\opengr~1\Contents as / and
c:\PROGRA~1\opengr~1\Contents\Cygwin\Versions\20A7OG~1.2\i686
as /bin ... done!
Starting grads under
c:\PROGRA~1\opengr~1\Contents\Cygwin\Versions\20A7OG~1.2\i686 ...

Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.2
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.a7.oga.2 little-endian readline printim grib2 netcdf hdf4-sds
hdf5 opendap-grids,stn athena geotiff
Issue 'q config' command for more information.
Loading User Defined Extensions table
</cygdrive/c/PROGRA~1/opengr~1/Contents/Cygwin/Versions/20A7OG~1.2/i686/gex/udxt>

... ok.
Landscape mode? ('n' for portrait):
GX Package Initialization: Size = 11 8.5


*Then I tried ls and it worked!*

ga-> !ls


ADVAPI32.DLL                                   SHELL32.dll
cygpopt-0.dll             libmf.pod
API-MS-WIN-Service-Core-L1-1-0.dll             SHLWAPI.dll
cygreadline6.dll          ln.exe
API-MS-WIN-Service-Management-L1-1-0.dll       USER32.dll
cygssl-0.9.8.dll          ls.exe
API-MS-WIN-Service-Management-L2-1-0.dll       USP10.dll
cygwin1.dll               merra.exe
API-MS-WIN-Service-winsvc-L1-1-0.dll           VERSION
cygxcb-1.dll              mkdir.exe
API-MS-Win-Core-Console-L1-1-0.dll             WINSPOOL.DRV
cygxcb-render-0.dll       mount.exe
API-MS-Win-Core-DateTime-L1-1-0.dll            awk.exe
cygxcb-render-util-0.dll  msvcrt.dll
API-MS-Win-Core-Debug-L1-1-0.dll               bufrscan.exe
cygxcb-xlib-0.dll         mv.exe
API-MS-Win-Core-DelayLoad-L1-1-0.dll           cat.exe
cygz.dll                  nomads.exe
API-MS-Win-Core-ErrorHandling-L1-1-0.dll       comdlg32.dll
date.exe                  ntdll.dll
API-MS-Win-Core-Fibers-L1-1-0.dll              cp.exe
diff.exe                  opengrads.exe
API-MS-Win-Core-File-L1-1-0.dll                cut.exe
dir.exe                   pngquant.exe
API-MS-Win-Core-Handle-L1-1-0.dll              cygICE-6.dll
dos2unix.exe              pwd.exe
API-MS-Win-Core-Heap-L1-1-0.dll                cygSM-6.dll
env.pod                   re.pod
API-MS-Win-Core-IO-L1-1-0.dll                  cygX11-6.dll
fish.pod                  rm.exe
API-MS-Win-Core-Interlocked-L1-1-0.dll         cygXau-6.dll
gex                       rxvt.exe
API-MS-Win-Core-LibraryLoader-L1-1-0.dll       cygXaw-7.dll
grads.exe                 saakeskus.pod_
API-MS-Win-Core-LocalRegistry-L1-1-0.dll       cygXdmcp-6.dll
gradsdap.exe              sed.exe
API-MS-Win-Core-Localization-L1-1-0.dll        cygXext-6.dll
gradsgui.exe              sh.exe
API-MS-Win-Core-Memory-L1-1-0.dll              cygXmu-6.dll
grep.exe                  shape.pod
API-MS-Win-Core-Misc-L1-1-0.dll                cygXpm-4.dll
grib2scan.exe             shfilt.pod
API-MS-Win-Core-NamedPipe-L1-1-0.dll           cygXrender-1.dll
gribmap.exe               sleep.exe
API-MS-Win-Core-ProcessEnvironment-L1-1-0.dll  cygXt-6.dll
gribscan.exe              sort.exe
API-MS-Win-Core-ProcessThreads-L1-1-0.dll      cygcairo-2.dll
gsudf.pod                 stnmap.exe
API-MS-Win-Core-Profile-L1-1-0.dll             cygcrypto-0.9.8.dll
gv32.exe                  tail.exe
API-MS-Win-Core-RtlSupport-L1-1-0.dll          cygexpat-1.dll
gv32.hlp                  tar.exe
API-MS-Win-Core-String-L1-1-0.dll              cygfontconfig-1.dll
gxeps.exe                 tcsh.exe
API-MS-Win-Core-Synch-L1-1-0.dll               cygfreetype-6.dll
gxps.exe                  tee.exe
API-MS-Win-Core-SysInfo-L1-1-0.dll             cyggcc_s-1.dll
gxtran.exe                test.exe
API-MS-Win-Core-ThreadPool-L1-1-0.dll          cyggfortran-3.dll
gxyat.exe                 touch.exe
API-MS-Win-Core-Util-L1-1-0.dll                cygglitz-1.dll
gzip.exe                  umount.exe
API-MS-Win-Security-Base-L1-1-0.dll            cygiconv-2.dll
head.exe                  uname.exe
COMCTL32.dll                                   cygintl-2.dll
lats4d.sh                 uniq.exe
GDI32.dll                                      cygintl-3.dll
less.exe                  unix2dos.exe
KERNEL32.dll                                   cygintl-8.dll
libW11.dll                wc.exe
KERNELBASE.dll                                 cygncurses-9.dll
libbjt.pod                wget.exe
LPK.dll                                        cygpixman-1-0.dll
libgrads.dll              wgrib.exe
RPCRT4.dll                                     cygpng12.dll
libipc.pod                wgrib2.exe
ga->




*Finally I shut down bash and tried running my script that has the calls to
imagemagik - and it worked!*

**

*Many thanks for all your patience and help.  I learned a bit along the way
too about bash!  Really appreciate your help. *


Thanks,
Mark


----- Original Message -----
From: "Arlindo da Silva" <dasilva at alum.mit.edu>
To: "Mark Sponsler" <msponsler at comcast.net>
Sent: Wednesday, February 10, 2010 7:10:39 AM GMT -08:00 US/Canada Pacific
Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug

On Tue, Feb 9, 2010 at 10:58 PM, Mark Sponsler <msponsler at comcast.net>wrote:

>  Hi Arlindo,
>
> I'm stumped.  First I tried it in the command prompt as you suggested:
>
>
>
> Microsoft Windows XP [Version 5.1.2600]
> (C) Copyright 1985-2001 Microsoft Corp.
>
> D:\Documents and Settings\Sponsler> ./umount /bin
>
>
> '.' is not recognized as an internal or external command,
>
> operable program or batch file.
>
>
>
> D:\Documents and Settings\Sponsler> dir
>
>  Volume in drive D has no label.
>  Volume Serial Number is B466-61C3
>
>   Directory of D:\Documents and Settings\Sponsler
>
> 02/02/2010  08:06 AM    <DIR>          .
> 02/02/2010  08:06 AM    <DIR>          ..
> 04/18/2008  07:12 AM                 0 !
> 02/21/2009  08:06 AM               666 .dodsrc
> 02/07/2010  07:27 PM    <DIR>          Desktop
> 02/07/2010  07:22 PM    <DIR>          Favorites
>
> 07/28/2009  07:18 PM             8,765 gsview32.ini
> 03/25/2006  02:04 PM    <DIR>          My Documents
> 01/09/2010  09:19 AM               600 PUTTY.RND
> 03/24/2006  02:33 PM    <DIR>          Start Menu
> 04/17/2008  07:59 PM            42,348 test.bin
> 03/25/2006  12:26 PM    <DIR>          WINDOWS
>                5 File(s)         52,379 bytes
>                7 Dir(s)   8,130,625,536 bytes free
>
>
>
> D:\Documents and Settings\Sponsler>
>
>
>
> Clearly "." is a directory.
>
>
>
> No luck there. So I tried it in sh
>
>
>
> sh-3.2$ ./mount
>
> c:\PROGRA~1\PCGRAD~2\win32 on /bin type user (binmode)
> c:\PROGRA~1\PCGRAD~2 on / type user (binmode)
> c: on /cygdrive/c type system (binmode,noumount)
> d: on /cygdrive/d type system (binmode,noumount)
>
>
> sh-3.2$ ./umount /bin
> umount: /bin: No such file or directory
>
>
>
> So I dug around and found there is no /bin directory in either Grads 1.8 or
> 1.9, but there is a /bin directory in 2.0.  So I tried to do that.
>

Make sure to use mount.exe shipped with v2.

>   sh-3.2$ ./umount c:/Program Files/opengrads/Classic/bin
>

This would not work...


>   Usage: umount [OPTION] [<posixpath>]
> Unmount filesystems
>
>   -A, --remove-all-mounts       remove all mounts
>   -c, --remove-cygdrive-prefix  remove cygdrive prefix
>   -h, --help                    output usage information and exit
>   -s, --system                  remove system mount (default)
>   -S, --remove-system-mounts    remove all system mounts
>   -u, --user                    remove user mount
>   -U, --remove-user-mounts      remove all user mounts
>   -v, --version                 output version information and exit
>
> sh-3.2$ ./umount -s c:/Program Files/opengrads/Classic/bin
> Usage: umount [OPTION] [<posixpath>]
> Unmount filesystems
>
>   -A, --remove-all-mounts       remove all mounts
>   -c, --remove-cygdrive-prefix  remove cygdrive prefix
>   -h, --help                    output usage information and exit
>   -s, --system                  remove system mount (default)
>   -S, --remove-system-mounts    remove all system mounts
>   -u, --user                    remove user mount
>   -U, --remove-user-mounts      remove all user mounts
>   -v, --version                 output version information and exit
>
> sh-3.2$
>
> But I couldn't figure out the command syntax and figured I was going down
> the wrong path anyway. Any other suggestions?
>

Try

./umount -U

or even

./umount -A

The point is, we need to have these mount points cleared.

   Arlindo


>   Thanks,
> Mark
>
>
> ----- Original Message -----
> From: "Arlindo da Silva" <dasilva at alum.mit.edu>
> To: "Mark Sponsler" <msponsler at comcast.net>
>  Sent: Tuesday, February 9, 2010 8:04:57 AM GMT -08:00 US/Canada Pacific
> Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug
>
> On Fri, Feb 5, 2010 at 10:27 PM, Mark Sponsler <msponsler at comcast.net>wrote:
>
>>  Hi Arlindo,
>>
>> I pulled everyting out of the path but the same windows directories as
>> before, did a cold reboot, then put the flash drive in and it became an F
>> drive.
>>
>>
>>
>> I think this is what you wanted:
>>
>>  Then I tried something else:
>>
>>
>>
>> I navigated to Grads 2.0 on the Flash drive and started sh. I issued the
>> ./mount command and then the ./ls command and things started looking
>> good.Could even get grads to start.
>>
>> sh-3.2$ ./mount
>> c:\PROGRA~1\PCGRAD~2\win32 on /bin type user (binmode)
>> c:\PROGRA~1\PCGRAD~2 on / type user (binmode)
>>
>
> Here is the problem! You most likely have these mounted from v1.9. It may
> have been left behind somehow... Try from the DOS cmd prompt:
>
> ./umount /bin
> ./umount /
>
> then
>
> ./mount
>
> The 2 first lines should disappear.  If so, start grads and try the "!ls".
>
>     Arlindo
>
>
>
>
>>  c: on /cygdrive/c type system (binmode,noumount)
>> d: on /cygdrive/d type system (binmode,noumount)
>> f: on /cygdrive/f type system (binmode,noumount)
>> sh-3.2$ ./ls
>> ADVAPI32.DLL                                   cygpopt-0.dll
>> API-MS-WIN-Service-Core-L1-1-0.dll             cygreadline6.dll
>> API-MS-WIN-Service-Management-L1-1-0.dll       cygssl-0.9.8.dll
>> API-MS-WIN-Service-Management-L2-1-0.dll       cygwin1.dll
>> API-MS-WIN-Service-winsvc-L1-1-0.dll           cygxcb-1.dll
>> API-MS-Win-Core-Console-L1-1-0.dll             cygxcb-render-0.dll
>> API-MS-Win-Core-DateTime-L1-1-0.dll            cygxcb-render-util-0.dll
>> API-MS-Win-Core-Debug-L1-1-0.dll               cygxcb-xlib-0.dll
>> API-MS-Win-Core-DelayLoad-L1-1-0.dll           cygz.dll
>> API-MS-Win-Core-ErrorHandling-L1-1-0.dll       date.exe
>> API-MS-Win-Core-Fibers-L1-1-0.dll              diff.exe
>> API-MS-Win-Core-File-L1-1-0.dll                dir.exe
>> API-MS-Win-Core-Handle-L1-1-0.dll              dos2unix.exe
>> API-MS-Win-Core-Heap-L1-1-0.dll                env.pod
>> API-MS-Win-Core-IO-L1-1-0.dll                  fish.pod
>> API-MS-Win-Core-Interlocked-L1-1-0.dll         gex
>> API-MS-Win-Core-LibraryLoader-L1-1-0.dll       grads.exe
>> API-MS-Win-Core-LocalRegistry-L1-1-0.dll       gradsdap.exe
>> API-MS-Win-Core-Localization-L1-1-0.dll        gradsgui.exe
>> API-MS-Win-Core-Memory-L1-1-0.dll              grep.exe
>> API-MS-Win-Core-Misc-L1-1-0.dll                grib2scan.exe
>> API-MS-Win-Core-NamedPipe-L1-1-0.dll           gribmap.exe
>> API-MS-Win-Core-ProcessEnvironment-L1-1-0.dll  gribscan.exe
>> API-MS-Win-Core-ProcessThreads-L1-1-0.dll      gsudf.pod
>> API-MS-Win-Core-Profile-L1-1-0.dll             gv32.exe
>> API-MS-Win-Core-RtlSupport-L1-1-0.dll          gv32.hlp
>> API-MS-Win-Core-String-L1-1-0.dll              gxeps.exe
>> API-MS-Win-Core-Synch-L1-1-0.dll               gxps.exe
>> API-MS-Win-Core-SysInfo-L1-1-0.dll             gxtran.exe
>> API-MS-Win-Core-ThreadPool-L1-1-0.dll          gxyat.exe
>> API-MS-Win-Core-Util-L1-1-0.dll                gzip.exe
>> API-MS-Win-Security-Base-L1-1-0.dll            head.exe
>> COMCTL32.dll                                   lats4d.sh
>> GDI32.dll                                      less.exe
>> KERNEL32.dll                                   libW11.dll
>> KERNELBASE.dll                                 libbjt.pod
>> LPK.dll                                        libgrads.dll
>> RPCRT4.dll                                     libipc.pod
>> SHELL32.dll                                    libmf.pod
>> SHLWAPI.dll                                    ln.exe
>> USER32.dll                                     ls.exe
>> USP10.dll                                      merra.exe
>> VERSION                                        mkdir.exe
>> WINSPOOL.DRV                                   mount.exe
>> awk.exe                                        msvcrt.dll
>> bufrscan.exe                                   mv.exe
>> cat.exe                                        nomads.exe
>> comdlg32.dll                                   ntdll.dll
>> cp.exe                                         opengrads.exe
>> cut.exe                                        pngquant.exe
>> cygICE-6.dll                                   pwd.exe
>> cygSM-6.dll                                    re.pod
>> cygX11-6.dll                                   rm.exe
>> cygXau-6.dll                                   rxvt.exe
>> cygXaw-7.dll                                   saakeskus.pod_
>> cygXdmcp-6.dll                                 sed.exe
>> cygXext-6.dll                                  sh.exe
>> cygXmu-6.dll                                   shape.pod
>> cygXpm-4.dll                                   shfilt.pod
>> cygXrender-1.dll                               sleep.exe
>> cygXt-6.dll                                    sort.exe
>> cygcairo-2.dll                                 stnmap.exe
>> cygcrypto-0.9.8.dll                            tail.exe
>> cygexpat-1.dll                                 tar.exe
>> cygfontconfig-1.dll                            tcsh.exe
>> cygfreetype-6.dll                              tee.exe
>> cyggcc_s-1.dll                                 test.exe
>> cyggfortran-3.dll                              touch.exe
>> cygglitz-1.dll                                 umount.exe
>> cygiconv-2.dll                                 uname.exe
>> cygintl-2.dll                                  uniq.exe
>> cygintl-3.dll                                  unix2dos.exe
>> cygintl-8.dll                                  wc.exe
>> cygncurses-9.dll                               wget.exe
>> cygpixman-1-0.dll                              wgrib.exe
>> cygpng12.dll                                   wgrib2.exe
>>
>>
>>
>> You can't see it here but the veriosn of grads in the list above is
>> version 2.0a07
>>
>>
>> sh-3.2$ ./grads
>>
>> Starting X server under f:\OPENGR~1\Contents\Resources\Xming
>> Starting grads under f:\OPENGR~1\Contents\Cygwin\Versions\20A7OG~1.2\i686
>> ...
>>
>>  Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.2
>> 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.a7.oga.2 little-endian readline printim grib2 netcdf hdf4-sds
>> hdf5
>> opendap-grids,stn athena geotiff
>> Issue 'q config' command for more information.
>> Loading User Defined Extensions table
>> </cygdrive/f/OPENGR~1/Contents/Cygwin/Vers
>> ions/20A7OG~1.2/i686/gex/udxt> ... ok.
>> Landscape mode? ('n' for portrait):
>> GX Package Initialization: Size = 11 8.5
>> ga-> ./!ls
>> Unknown command: ./!ls
>> ga->
>>
>>
>>
>> Couldn't get beyond that though....
>> Thanks,
>> Mark
>>
>>
>> ----- Original Message -----
>> From: "Arlindo da Silva" <dasilva at alum.mit.edu>
>> To: "Mark Sponsler" <msponsler at comcast.net>
>> Sent: Friday, February 5, 2010 10:22:36 AM GMT -08:00 US/Canada Pacific
>> Subject: Re: WIN GRADS 2.0.a7.oga.3 - Possible Bug
>>
>>
>>
>> On Fri, Feb 5, 2010 at 10:38 AM, Mark Sponsler <msponsler at comcast.net>wrote:
>>
>>>  Hi Arlindo,
>>>
>>> So I went to  Control Panel>System>Advanced>Env Variables>Path and pulled
>>> everything out of the path but the essential windows stuff. Also renamed the
>>> AUTOEXEC.bat file to something so it could not be found. Then did a cold
>>> reboot.  Then issuesd the path command after it restarted just to verify my
>>> path was clean:
>>>
>>>
>>>
>>> Microsoft Windows XP [Version 5.1.2600]
>>> (C) Copyright 1985-2001 Microsoft Corp.
>>>
>>> D:\Documents and Settings\Sponsler>path
>>> PATH=D:\Perl\bin\;D:\WINDOWS\system32;D:\WINDOWS;D:\WINDOWS\System32\Wbem
>>>
>>> D:\Documents and Settings\Sponsler>
>>>
>>>
>>>
>>> Then I plugged in the Flash Drive with Grads on it and navigated to the
>>> directory with the grads binaries. They were here:
>>>
>>>
>>>
>>> F:\opengrads\Contents\Cygwin\Versions\2.0.a7.oga.2\i686
>>>
>>> Next I did as you requested, fired up bash and tried to mount and ls a
>>> tried to start grads:
>>>
>>>
>>>
>>   sh-3.2$ mount
>>>
>>
>> It is likely that "." may not be in your path. Try this instead:
>>
>> C:\> cd F:\Opengrads\Contents\Cygwin\Versions\2.0.a7.oga.3\i686
>> C:\> dir
>> C:\> sh
>> sh-3.2$ ./mount
>> sh-3.2$ ./ls
>>
>> Make sure mount.exe and ls.exe are in your current directory.
>>
>>
>>     Arlindo
>>
>>
>> --
>> Arlindo da Silva
>> dasilva at alum.mit.edu
>>
>
>
>
> --
> Arlindo da Silva
> dasilva at alum.mit.edu
>



-- 
Arlindo da Silva
dasilva at alum.mit.edu



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

From jma at COLA.IGES.ORG  Fri Feb 12 13:48:13 2010
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Fri, 12 Feb 2010 13:48:13 -0500
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <4B75874E.9000403@noaa.gov>
Message-ID: <FRI.12.FEB.2010.134813.0500.>

Dear All,
Everyone at COLA is interested in looking at the new CFS Reanalysis,  
so you can be sure that GrADS will be able handle it. The -b option  
for gribmap has not been implemented for grib2 yet, so that will have  
to be patched quickly. We're working on some other changes to gribmap  
to handle large (>2GB) grib files, so the timing will be just about  
right for a new release with these features. Please don't ask when, I  
don't know.

As for the CFSR data behind TDS at NCDC, GrADS works fine as a client,  
but you will likely have to create a descriptor file, especially if  
the data are not on a regular lat/lon grid and need a PDEF statement.  
The variable names are very long and need to be aliased. The time axis  
metadata from TDS may also be incompatible with 'sdfopen'. I was able  
to open the URL posted earlier in this thread and display several  
variables.

I have not had a chance to grab CFSR samples and see where the  
problems are (as far as GrADS is concerned). It sounds like it will  
depend a lot on which file/product you are looking at. I know the data  
wranglers at NCEP and NCDC are working hard to make this data  
available and usable and it's not an easy job. If you are eager to  
look at this data, make sure you are using the most recent versions of  
all the available tools. And if you post questions to the list, be  
specific about the version you're using and the operating system  
you're running on as well as the data file (or URL) you're trying to  
read.

--Jennifer

On Feb 12, 2010, at 11:52 AM, Dan Swank wrote:

> Hello,
>
> I need to take back a few statements I made in
> my previous post.
>
> I have gotten several CFSR/CFSRMON files to plot successfully.
> Here are some samples:
>
> http://nomads.ncdc.noaa.gov/data/cfsr/200001/wnd10m.gdas.200001.grb2
> http://nomads.ncdc.noaa.gov/data/cfsr/200001/wnd10m.gdas.200001.ctl
> http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2.idx
>
> http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2
> http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.ctl
> http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2.idx
>
> http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.ctl
> http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.grb2.idx
> http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas. 
> 200001.grb2
>
> http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.grb2
> http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.ctl
> http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.grb2.idx
>
> The -b should be used when invoking GrADS2 gribmap.
> The -0 opt should be set when working with monthly means
>  (cfsrmon) and diurnal fraction monthlies.
>
> Here are the known gotchas:
> wgrib2 v.1.8.2 or earlier may produce control files
> with incorrect TDEF's when used with g2ctl.pl .
>
> The original CFSR was produced in GRIB1 at NCEP.
> NCEP converted this data to GRIB2 for transmission
> and storage purposes.  There are possible issues
> in the translation between GRIB1 and GRIB2 [indicator
> of statistical process] header.  This affects monthly
> means (which are all statistically processed in some way).
> gribmap's -0 option nulls this *potential* issue by matching
> using the reference time.  Again, this issue is unproven.
>
> The GRIB2 headers contain local extensions
> used by Table 4.10 (-fix_ncep wgrib2 option).
> I had thought this might cause issues with gribmap
>  but am now not certain.
>
>
> -Dan
>
>
>
>
> Henrique Barbosa wrote the following on 2/10/2010 1:59 PM:
>> Dear Dan Swank,
>>
>> Thank you for clarifying the problem. The difference between GDS  
>> and TDS
>> servers was not clear to me. I have not mentioned in my previous  
>> email,
>> but I already tried to download the data and write a ctl...
>> ... and I got the time coordinate error you mentioned!
>>
>>> I hope the GrADS dev-team will take some time to get the
>>> CFSR working in GrADS.
>>
>> I can only enforce your request. +1 here.
>>
>> []'s
>> Henrique
>>
>> On Tue, Feb 9, 2010 at 6:34 PM, Dan Swank <Dan.Swank at noaa.gov> wrote:
>>> Hello,
>>>
>>> GrADS in general is not ready for the CFSR dataset.
>>> I was able to use gradsdap alpha5 to sdfopen the source
>>> you listed ~ but this is a THREDDS server, which is
>>> presenting variable names which are too long for
>>> grads to handle.  [(q file) snipit]:
>>>
>>>    relative_humidi  0  t,y,x  Relative humidity @ sigma
>>>    relative_humidi  0  t,y,x  Relative humidity @ sigma_layer
>>>    relative_humidi  0  t,y,x  Relative humidity @ zeroDegC_isotherm
>>>    specific_humidi  0  ty,x  Specific humidity @ pressure
>>>    specific_humidi  0  ty,x  Specific humidity @  
>>> pressure_difference_layer
>>>    storm_relative_  0  ty,x  Storm relative helicity @
>>> height_above_ground_layer
>>>    stream_function  0  ty,x  Stream function @ pressure
>>>    surface_lifted_  0  t,y,x  Surface Lifted Index @ surface
>>>    temperature  0  ty,x  Temperature @ pressure
>>>    temperature_alt  4  t,z,y,x  Temperature @ altitude_above_msl
>>>    temperature_pot  0  ty,x  Temperature @  
>>> potential_vorticity_surface
>>>    temperature_max  0  t,y,x  Temperature @ maximum_wind
>>>
>>> So even if you managed to open it, you wouldn't get much further.
>>>
>>> Connecting a GrADS client to a THREDDS server has many known issues.
>>> It will flat out refuse to identify the x/y coordinates when opening
>>> non-lat/lon grids served from THREDDS, for instance.  Although the
>>> source you specified IS a lat/lon, some CFSR sets are on
>>> gaussian coordinates.  The problem there is that THREDDS presents
>>> the native grid x/y's without projecting it to a lat/lon
>>> like a GDS does.
>>>
>>> Attempting to download a GRIB2 file and creating a control
>>> file will result in a gribmap2 not being able to MATCH the
>>> time coordinates because of the new local GRIB table the headers
>>> use.  This is especially evident in the monthly means.  This
>>> pretty much stops us from putting the CFSR behind a GDS2.
>>>
>>> I hope the GrADS dev-team will take some time to get the
>>> CFSR working in GrADS.
>>>
>>> -Dan Swank
>>> NCDC-NOMADS Support
>>>
>>>
>>>
>>>
>>> Henrique Barbosa wrote the following on 2/9/2010 12:12 PM:
>>>> Dear All,
>>>>
>>>> I need help to access CFSR (NCEP reanalysis 3) over opendap with  
>>>> opengrads.
>>>> I get an horrendous error message (see below) when I try to:
>>>>
>>>> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>>>>
>>>> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
>>>> have opendap enabled (see output of 'q config' below).
>>>>
>>>> The strange thing is that I can open other sources of data, for  
>>>> example
>>>>
>>>> ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z
>>>>
>>>> works perfectly and I can plot, subset, etc...
>>>>
>>>> Any ideas?
>>>>
>>>> []s
>>>> Henrique
>>>>
>>>>
>>>> --------------------------------------------------------------------------------------------------------------------------------------------------
>>>> ga-> q config
>>>> Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf
>>>> hdf4-sds hdf5 opendap-grids,stn athena geotiff
>>>> Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3
>>>> Copyright (c) 1988-2009 by Brian Doty and the
>>>> Institute for Global Environment and Society (IGES)
>>>> This program is distributed WITHOUT ANY WARRANTY
>>>> See file COPYRIGHT for more information.
>>>>
>>>> Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu
>>>>
>>>> This version of GrADS has been configured with the following  
>>>> options:
>>>>  o Built on a LITTLE ENDIAN machine
>>>>  o Command line editing ENABLED
>>>>      http://tiswww.case.edu/php/chet/readline/rltop.html
>>>>  o printim command for image output ENABLED
>>>>      http://www.zlib.net
>>>>      http://www.libpng.org/pub/png/libpng.html
>>>>      http://www.libgd.org/Main_Page
>>>>  o GRIB2 interface ENABLED
>>>>      http://www.ijg.org
>>>>      http://www.ece.uvic.ca/~mdadams/jasper
>>>>      http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2
>>>>      g2clib-1.0.5
>>>>  o NetCDF interface ENABLED
>>>>      http://www.opendap.org
>>>>      libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar  3 2009  
>>>> 14:13:33 $
>>>>  o HDF interface ENABLED
>>>>      http://hdfgroup.org
>>>>      HDF 4.2r3
>>>>      HDF5 1.8.2
>>>>  o Athena Widget GUI ENABLED
>>>>  o OPeNDAP gridded data interface ENABLED
>>>>      http://www.opendap.org
>>>>      libdap 3.7.10
>>>>  o OPeNDAP station data interface ENABLED
>>>>      http://iges.org/grads/gadoc/supplibs.html
>>>>      libgadap 2.0.oga.1
>>>>  o GeoTIFF and KML output ENABLED
>>>>      http://www.libtiff.org
>>>>      http://geotiff.osgeo.org
>>>>
>>>> For additional information please consult http://iges.org/grads
>>>>
>>>> --------------------------------------------------------------------------------------------------------------------------------------------------
>>>> ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2
>>>> Scanning self-describing file:
>>>> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>>>> *** glibc detected *** grads: free(): invalid pointer:  
>>>> 0x00000000018ff354 ***
>>>> ======= Backtrace: =========
>>>> /lib/libc.so.6[0x2b6884e9fcb8]
>>>> /lib/libc.so.6(cfree+0x76)[0x2b6884ea2276]
>>>> grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea]
>>>> grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9) 
>>>> [0x833a05]
>>>> grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea]
>>>> grads(nc3_get_vara_double+0x9e)[0x80e3b8]
>>>> grads(nc_get_vara_double+0x1a)[0x60bcca]
>>>> grads(read_one_dimension+0x89)[0x541ac1]
>>>> grads(gadsdf+0x17d6)[0x540080]
>>>> grads(gasdfopen+0x1d4)[0x53e700]
>>>> grads(gacmd+0xfb2)[0x524504]
>>>> grads(Main+0x4e1)[0x496f01]
>>>> grads(main+0x2e)[0x496a16]
>>>> /lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6]
>>>> grads(__gxx_personality_v0+0x1a9)[0x496959]
>>>> ======= Memory map: ========
>>>> 00400000-00cca000 r-xp 00000000 08:01 604869
>>>>  /usr/local/grads-2.0.a7.oga.3/bin/grads
>>>> 00dc9000-00dea000 rwxp 008c9000 08:01 604869
>>>>  /usr/local/grads-2.0.a7.oga.3/bin/grads
>>>> 00dea000-010aa000 rwxp 00dea000 00:00 0
>>>> 0180f000-01a04000 rwxp 0180f000 00:00  
>>>> 0                                  [heap]
>>>> 2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056
>>>>  /lib/ld-2.9.so
>>>> 2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0
>>> --
>>> Dan Swank
>>> STG, Incorporated - Government Contractor
>>> NCDC-NOMADS Project:  Software & Data Management
>>> Data Access Branch
>>> National Climatic Data Center
>>> Veach-Baley Federal Building
>>> 151 Patton Avenue
>>> Asheville, NC 28801-5001
>>> Email: Dan.Swank at noaa.gov
>>> Phone: 828-271-4007
>>>
>
> -- 
> Dan Swank
> STG, Incorporated - Government Contractor
> NCDC-NOMADS Project:  Software & Data Management
> Data Access Branch
> National Climatic Data Center
> Veach-Baley Federal Building
> 151 Patton Avenue
> Asheville, NC 28801-5001
> Email: Dan.Swank at noaa.gov
> Phone: 828-271-4007

--
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/20100212/e79f1334/attachment.html 

From pietrowicz at ZEDXINC.COM  Fri Feb 12 14:11:42 2010
From: pietrowicz at ZEDXINC.COM (Joe Pietrowicz)
Date: Fri, 12 Feb 2010 14:11:42 -0500
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <11461271.213781266001679038.JavaMail.root@shark.zedxinc.com>
Message-ID: <FRI.12.FEB.2010.141142.0500.>

All,

I don't know if this is a related problem. I downloaded model and observation climate files to use with the CFS climate model to adjust the model for bias. Most of the time the monthly climate files work fine. However, it looks like the date in the grib1 files are wrong at times. This example I have below is taken from file available on the 30th of January. The February data, in grib1 clearly has the date of 19810230, an impossible date. If the code used to generate the new CFSR monthly data sets is unchanged then the same error will continue. My guess is that the code increments the month, not looking at the day. 

wgrib -verf -4yr dd.xxxx03_pa.modclim_tl.sflux | more

1:0:d=1981023000:PRATE:kpds5=59:kpds6=1:kpds7=0:TR=133:P1=1:P2=2:TimeU=3:sfc:time?:NAve=22
2:33924:d=1981023000:TMP:kpds5=11:kpds6=1:kpds7=0:TR=133:P1=1:P2=2:TimeU=3:sfc:time?:NAve=22
3:70104:d=1981023000:TMP:kpds5=11:kpds6=105:kpds7=2:TR=133:P1=1:P2=2:TimeU=3:2 m above gnd:time?:NAve=22
4:106284:d=1981023000:DLWRF:kpds5=205:kpds6=1:kpds7=0:TR=133:P1=1:P2=2:TimeU=3:sfc:time?:NAve=22
5:142464:d=1981023000:DSWRF:kpds5=204:kpds6=1:kpds7=0:TR=133:P1=1:P2=2:TimeU=3:sfc:time?:NAve=22
6:176388:d=1981023000:GFLUX:kpds5=155:kpds6=1:kpds7=0:TR=133:P1=1:P2=2:TimeU=3:sfc:time?:NAve=22
7:192936:d=1981023000:HPBL:kpds5=221:kpds6=1:kpds7=0:TR=133:P1=1:P2=2:TimeU=3:sfc:time?:NAve=22
...


As you can see the grib file contains an invalid date of 1981023000. It's this date that causes the problem. When I change the date in the control *ctl file it will not map properly because I get 'Entire Grid Undefined' because gribmap reads the grib file and generates an error that time is out of bounds, no matter what I put into the *ctl file. The problem is the invalid time of 30Feb in the grib file.

I hope that I'm on the right track and not adding to confusion.

Joe Pietrowicz, Ph.D.
Senior Research Scientist
ZedX, Inc.


----- Original Message -----
From: "Jennifer Adams" <jma at COLA.IGES.ORG>
To: GRADSUSR at LIST.CINECA.IT
Sent: Friday, February 12, 2010 1:48:13 PM GMT -05:00 US/Canada Eastern
Subject: Re: Problem to open grib2 Nomad/CFSR data with opengrads

Dear All, 
Everyone at COLA is interested in looking at the new CFS Reanalysis, so you can be sure that GrADS will be able handle it. The -b option for gribmap has not been implemented for grib2 yet, so that will have to be patched quickly. We're working on some other changes to gribmap to handle large (>2GB) grib files, so the timing will be just about right for a new release with these features. Please don't ask when, I don't know. 



As for the CFSR data behind TDS at NCDC, GrADS works fine as a client, but you will likely have to create a descriptor file, especially if the data are not on a regular lat/lon grid and need a PDEF statement. The variable names are very long and need to be aliased. The time axis metadata from TDS may also be incompatible with 'sdfopen'. I was able to open the URL posted earlier in this thread and display several variables. 


I have not had a chance to grab CFSR samples and see where the problems are (as far as GrADS is concerned). It sounds like it will depend a lot on which file/product you are looking at. I know the data wranglers at NCEP and NCDC are working hard to make this data available and usable and it's not an easy job. If you are eager to look at this data, make sure you are using the most recent versions of all the available tools. And if you post questions to the list, be specific about the version you're using and the operating system you're running on as well as the data file (or URL) you're trying to read. 


--Jennifer 


On Feb 12, 2010, at 11:52 AM, Dan Swank wrote: 



Hello, 

I need to take back a few statements I made in 
my previous post. 

I have gotten several CFSR/CFSRMON files to plot successfully. 
Here are some samples: 

http://nomads.ncdc.noaa.gov/data/cfsr/200001/wnd10m.gdas.200001.grb2 
http://nomads.ncdc.noaa.gov/data/cfsr/200001/wnd10m.gdas.200001.ctl 
http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2.idx 

http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2 
http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.ctl 
http://nomads.ncdc.noaa.gov/data/cfsr/200001/pwat.gdas.200001.grb2.idx 

http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.ctl 
http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.grb2.idx 
http://nomads.ncdc.noaa.gov/data/cfsrmon/200001/ipvh02.gdas.200001.grb2 

http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.grb2 
http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.ctl 
http://nomads.ncdc.noaa.gov/data/cfsrmon/200909/ocnf04.gdas.200909.06Z.grb2.idx 

The -b should be used when invoking GrADS2 gribmap. 
The -0 opt should be set when working with monthly means 
(cfsrmon) and diurnal fraction monthlies. 

Here are the known gotchas: 
wgrib2 v.1.8.2 or earlier may produce control files 
with incorrect TDEF's when used with g2ctl.pl . 

The original CFSR was produced in GRIB1 at NCEP. 
NCEP converted this data to GRIB2 for transmission 
and storage purposes. There are possible issues 
in the translation between GRIB1 and GRIB2 [indicator 
of statistical process] header. This affects monthly 
means (which are all statistically processed in some way). 
gribmap's -0 option nulls this *potential* issue by matching 
using the reference time. Again, this issue is unproven. 

The GRIB2 headers contain local extensions 
used by Table 4.10 (-fix_ncep wgrib2 option). 
I had thought this might cause issues with gribmap 
but am now not certain. 


-Dan 




Henrique Barbosa wrote the following on 2/10/2010 1:59 PM: 


Dear Dan Swank, 





Thank you for clarifying the problem. The difference between GDS and TDS 


servers was not clear to me. I have not mentioned in my previous email, 


but I already tried to download the data and write a ctl... 


... and I got the time coordinate error you mentioned! 







I hope the GrADS dev-team will take some time to get the 




CFSR working in GrADS. 





I can only enforce your request. +1 here. 





[]'s 


Henrique 





On Tue, Feb 9, 2010 at 6:34 PM, Dan Swank <Dan.Swank at noaa.gov> wrote: 




Hello, 









GrADS in general is not ready for the CFSR dataset. 




I was able to use gradsdap alpha5 to sdfopen the source 




you listed ~ but this is a THREDDS server, which is 




presenting variable names which are too long for 




grads to handle. [(q file) snipit]: 









relative_humidi 0 t,y,x Relative humidity @ sigma 




relative_humidi 0 t,y,x Relative humidity @ sigma_layer 




relative_humidi 0 t,y,x Relative humidity @ zeroDegC_isotherm 




specific_humidi 0 ty,x Specific humidity @ pressure 




specific_humidi 0 ty,x Specific humidity @ pressure_difference_layer 




storm_relative_ 0 ty,x Storm relative helicity @ 




height_above_ground_layer 




stream_function 0 ty,x Stream function @ pressure 




surface_lifted_ 0 t,y,x Surface Lifted Index @ surface 




temperature 0 ty,x Temperature @ pressure 




temperature_alt 4 t,z,y,x Temperature @ altitude_above_msl 




temperature_pot 0 ty,x Temperature @ potential_vorticity_surface 




temperature_max 0 t,y,x Temperature @ maximum_wind 









So even if you managed to open it, you wouldn't get much further. 









Connecting a GrADS client to a THREDDS server has many known issues. 




It will flat out refuse to identify the x/y coordinates when opening 




non-lat/lon grids served from THREDDS, for instance. Although the 




source you specified IS a lat/lon, some CFSR sets are on 




gaussian coordinates. The problem there is that THREDDS presents 




the native grid x/y's without projecting it to a lat/lon 




like a GDS does. 









Attempting to download a GRIB2 file and creating a control 




file will result in a gribmap2 not being able to MATCH the 




time coordinates because of the new local GRIB table the headers 




use. This is especially evident in the monthly means. This 




pretty much stops us from putting the CFSR behind a GDS2. 









I hope the GrADS dev-team will take some time to get the 




CFSR working in GrADS. 









-Dan Swank 




NCDC-NOMADS Support 
























Henrique Barbosa wrote the following on 2/9/2010 12:12 PM: 






Dear All, 













I need help to access CFSR (NCEP reanalysis 3) over opendap with opengrads. 






I get an horrendous error message (see below) when I try to: 













ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2 













I am using v2.0.a7.oga.3 downloaded from opengrads website. I do 






have opendap enabled (see output of 'q config' below). 













The strange thing is that I can open other sources of data, for example 













ga-> sdfopen http://nomads.ncep.noaa.gov:9090/dods/gfs/gfs20100111/gfs_00z 













works perfectly and I can plot, subset, etc... 













Any ideas? 













[]s 






Henrique 




















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






ga-> q config 






Config: v2.0.a7.oga.3 little-endian readline printim grib2 netcdf 






hdf4-sds hdf5 opendap-grids,stn athena geotiff 






Grid Analysis and Display System (GrADS) Version 2.0.a7.oga.3 






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






Institute for Global Environment and Society (IGES) 






This program is distributed WITHOUT ANY WARRANTY 






See file COPYRIGHT for more information. 













Built Thu Oct 29 17:57:08 EDT 2009 for x86_64-unknown-linux-gnu 













This version of GrADS has been configured with the following options: 






o Built on a LITTLE ENDIAN machine 






o Command line editing ENABLED 






http://tiswww.case.edu/php/chet/readline/rltop.html 






o printim command for image output ENABLED 






http://www.zlib.net 






http://www.libpng.org/pub/png/libpng.html 






http://www.libgd.org/Main_Page 






o GRIB2 interface ENABLED 






http://www.ijg.org 






http://www.ece.uvic.ca/~mdadams/jasper 






http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2 






g2clib-1.0.5 






o NetCDF interface ENABLED 






http://www.opendap.org 






libnc-dap 4.0.1-beta3-snapshot2009021712 of Mar 3 2009 14:13:33 $ 






o HDF interface ENABLED 






http://hdfgroup.org 






HDF 4.2r3 






HDF5 1.8.2 






o Athena Widget GUI ENABLED 






o OPeNDAP gridded data interface ENABLED 






http://www.opendap.org 






libdap 3.7.10 






o OPeNDAP station data interface ENABLED 






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






libgadap 2.0.oga.1 






o GeoTIFF and KML output ENABLED 






http://www.libtiff.org 






http://geotiff.osgeo.org 













For additional information please consult http://iges.org/grads 













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






ga-> sdfopen http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/2001912/pgblnl.gdas.200912.grb2 






Scanning self-describing file: 






http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2 






*** glibc detected *** grads: free(): invalid pointer: 0x00000000018ff354 *** 






======= Backtrace: ========= 






/lib/libc.so.6[0x2b6884e9fcb8] 






/lib/libc.so.6(cfree+0x76)[0x2b6884ea2276] 






grads(_Z15convert_nc_type7nc_typeiiPvS0_+0x2ca)[0x8011ea] 






grads(_ZN7NCArray14extract_valuesEPviiRK12ClientParams+0x2e9)[0x833a05] 






grads(_Z8DODvarioiiPKmS0_PKlPvi+0xc6e)[0x803aea] 






grads(nc3_get_vara_double+0x9e)[0x80e3b8] 






grads(nc_get_vara_double+0x1a)[0x60bcca] 






grads(read_one_dimension+0x89)[0x541ac1] 






grads(gadsdf+0x17d6)[0x540080] 






grads(gasdfopen+0x1d4)[0x53e700] 






grads(gacmd+0xfb2)[0x524504] 






grads(Main+0x4e1)[0x496f01] 






grads(main+0x2e)[0x496a16] 






/lib/libc.so.6(__libc_start_main+0xe6)[0x2b6884e465a6] 






grads(__gxx_personality_v0+0x1a9)[0x496959] 






======= Memory map: ======== 






00400000-00cca000 r-xp 00000000 08:01 604869 






/usr/local/grads-2.0.a7.oga.3/bin/grads 






00dc9000-00dea000 rwxp 008c9000 08:01 604869 






/usr/local/grads-2.0.a7.oga.3/bin/grads 






00dea000-010aa000 rwxp 00dea000 00:00 0 






0180f000-01a04000 rwxp 0180f000 00:00 0 [heap] 






2b6882c7a000-2b6882c9a000 r-xp 00000000 08:01 83056 






/lib/ld-2.9.so 






2b6882c9a000-2b6882c9d000 rwxp 2b6882c9a000 00:00 0 




-- 




Dan Swank 




STG, Incorporated - Government Contractor 




NCDC-NOMADS Project: Software & Data Management 




Data Access Branch 




National Climatic Data Center 




Veach-Baley Federal Building 




151 Patton Avenue 




Asheville, NC 28801-5001 




Email: Dan.Swank at noaa.gov 




Phone: 828-271-4007 






-- 
Dan Swank 
STG, Incorporated - Government Contractor 
NCDC-NOMADS Project: Software & Data Management 
Data Access Branch 
National Climatic Data Center 
Veach-Baley Federal Building 
151 Patton Avenue 
Asheville, NC 28801-5001 
Email: Dan.Swank at noaa.gov 
Phone: 828-271-4007 




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


From jmssnyder at UCDAVIS.EDU  Fri Feb 12 16:20:50 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Fri, 12 Feb 2010 13:20:50 -0800
Subject: Two Questions about NETCDF and GRIB
Message-ID: <FRI.12.FEB.2010.132050.0800.>

I am running an aerosol puff model that can only use GRIB data sets and also
needs the data to be in hourly time steps.  The data I currently have is in
daily time steps and I was wondering if I could interpolate this data to
have to display in hourly time steps.

Also how do I convert netcdf to grib?

Thanks,

Jason
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100212/ff04b50b/attachment.html 

From jdduda at IASTATE.EDU  Fri Feb 12 16:52:17 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Fri, 12 Feb 2010 15:52:17 -0600
Subject: Two Questions about NETCDF and GRIB
In-Reply-To: <fe8efa01002121320o4f2817f2v56465fe1c6003e4a@mail.gmail.com>
Message-ID: <FRI.12.FEB.2010.155217.0600.>

One way I know how to convert netcdf to grib is by using the WRF Post
Processor (generally only works if the model you're running is the WRF
though).

Regarding the linear interpolation of time between sets, I don't think you
can create new files that have interpolated data in them, but you can always
compute the interpolated values yourself in Grads and just display them.
For example, if you have wind speed 24 hours apart, you could compute a
linearly interpolated value of wind speed at a given hour in between by
defining a variable this way:

define dwinddt = windspeed(t=t1) - windspeed(t=t0)/(24 hours, 1440 minutes,
86400 s)
display windspeed(t=t0) + dwinddt*(1 hr, 60 min, or 3600 s depending on what
unit you want to use)

Hope this helps.

Jeff Duda

On Fri, Feb 12, 2010 at 3:20 PM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> I am running an aerosol puff model that can only use GRIB data sets and
> also needs the data to be in hourly time steps.  The data I currently have
> is in daily time steps and I was wondering if I could interpolate this data
> to have to display in hourly time steps.
>
> Also how do I convert netcdf to grib?
>
> Thanks,
>
> Jason
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100212/79692e68/attachment.html 

From sigmond at ATMOSP.PHYSICS.UTORONTO.CA  Fri Feb 12 17:15:19 2010
From: sigmond at ATMOSP.PHYSICS.UTORONTO.CA (Michael Sigmond)
Date: Fri, 12 Feb 2010 23:15:19 +0100
Subject: masked continents on polar stereographic map with circular frame?
Message-ID: <FRI.12.FEB.2010.231519.0100.>

Hello,
I'm trying to plot a variable on a polar stereographic map, filling the the
continents with a gray color, and with a circular frame (I know, that's a
lot to ask..)

I'm using:

set mproj sps
set lat -90 -50
set mpvals -180 180 -90 -50
set frame circle
d variable

I tried using basemap, but then the frame turns square, which I don't really
like. Any suggestions are much appreciated!

Thank you,

Michael 


From schan at FSU.EDU  Sat Feb 13 17:01:03 2010
From: schan at FSU.EDU (Steven Chan)
Date: Sat, 13 Feb 2010 17:01:03 -0500
Subject: Two Questions about NETCDF and GRIB
In-Reply-To: <c0f1d7891002121352if43048ej73b6818506e15e5d@mail.gmail.com>
Message-ID: <SAT.13.FEB.2010.170103.0500.>

http://www.nco.ncep.noaa.gov/pmb/codes/GRIB2/

If the grib format you want is grib2 then you can consider the NCEP 
wgrib2 tool. If it is grib1 than it is quite a bit more tricky. NCEP has 
a Fortran library called w3lib that you can try for doing I/O with 
grib1; it may take some time and research to understand how it works -- 
there is a code example inside the Readme in w3lib tar file, but dealing 
with the Grib header is tricky. A trick is to read from a similar grib1 
file header, make necessary changes (like the date), and just change the 
output data.

Steven Chan


On 2/12/2010 4:52 PM, Jeffrey Duda wrote:
> One way I know how to convert netcdf to grib is by using the WRF Post 
> Processor (generally only works if the model you're running is the WRF 
> though).
>
> Regarding the linear interpolation of time between sets, I don't think 
> you can create new files that have interpolated data in them, but you 
> can always compute the interpolated values yourself in Grads and just 
> display them.  For example, if you have wind speed 24 hours apart, you 
> could compute a linearly interpolated value of wind speed at a given 
> hour in between by defining a variable this way:
>
> define dwinddt = windspeed(t=t1) - windspeed(t=t0)/(24 hours, 1440 
> minutes, 86400 s)
> display windspeed(t=t0) + dwinddt*(1 hr, 60 min, or 3600 s depending 
> on what unit you want to use)
>
> Hope this helps.
>
> Jeff Duda
>
> On Fri, Feb 12, 2010 at 3:20 PM, Jason Snyder <jmssnyder at ucdavis.edu 
> <mailto:jmssnyder at ucdavis.edu>> wrote:
>
>     I am running an aerosol puff model that can only use GRIB data
>     sets and also needs the data to be in hourly time steps.  The data
>     I currently have is in daily time steps and I was wondering if I
>     could interpolate this data to have to display in hourly time steps.
>     Also how do I convert netcdf to grib?
>     Thanks,
>     Jason
>
>
>
>
> -- 
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda <http://www.meteor.iastate.edu/%7Ejdduda>

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

From arlindo.dasilva at GMAIL.COM  Sat Feb 13 18:31:06 2010
From: arlindo.dasilva at GMAIL.COM (Arlindo da Silva)
Date: Sat, 13 Feb 2010 18:31:06 -0500
Subject: Two Questions about NETCDF and GRIB
In-Reply-To: <fe8efa01002121320o4f2817f2v56465fe1c6003e4a@mail.gmail.com>
Message-ID: <SAT.13.FEB.2010.183106.0500.>

On Fri, Feb 12, 2010 at 4:20 PM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> I am running an aerosol puff model that can only use GRIB data sets and
> also needs the data to be in hourly time steps.  The data I currently have
> is in daily time steps and I was wondering if I could interpolate this data
> to have to display in hourly time steps.
>
> Also how do I convert netcdf to grib?
>
>

Try Lats4d:

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

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

From seemanth.m at GMAIL.COM  Mon Feb 15 07:40:58 2010
From: seemanth.m at GMAIL.COM (seemanth m)
Date: Mon, 15 Feb 2010 18:10:58 +0530
Subject: How to plot cyclone track ?
Message-ID: <MON.15.FEB.2010.181058.0530.>

Dear grads users,
I am new to grads software.
Can anybody please point out a solution for the following problem?

I have an atmospheric data set for an area (u wind, v wind, atmospheric
pressure in different levels) during a cyclone period, in grads usable
format (*.ctl and *.dat files) as well as in netcdf (*.nc file) format. .How
can I plot the cyclone track (with a line) accurately, by analysing this
data set?

Thanking you in advance
seemanth m
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100215/3c8e3d3b/attachment.html 

From taranaphy at YAHOO.COM  Mon Feb 15 14:20:07 2010
From: taranaphy at YAHOO.COM (tarana mahzabin)
Date: Mon, 15 Feb 2010 11:20:07 -0800
Subject: How to plot cyclone track ?
In-Reply-To: <ea96d58b1002150440i546bb6eds6fb746206e56415a@mail.gmail.com>
Message-ID: <MON.15.FEB.2010.112007.0800.>

Dear Seemanth,

Here is a gs file for cyclone track and you have  to change in some places. 
Hope it will help you.

 Tarana Mahzabin
Graduate Student,
Dept. of Earth and Atmospheric Science,
U of A, Edmonton, Canada.




________________________________
From: seemanth m <seemanth.m at GMAIL.COM>
To: GRADSUSR at LIST.CINECA.IT
Sent: Monday, February 15, 2010 5:40:58
Subject: How to plot cyclone track ?

Dear grads users,
I am new to grads software.
Can anybody please point out a solution for the following problem?

I have an atmospheric data set for
an area (u wind, v wind, atmospheric pressure in different levels)
during a cyclone period, in grads usable format (*.ctl and *.dat files)
as well as in netcdf (*.nc file) format. .How can I plot the cyclone track (with a line) accurately, by analysing this data set? 

Thanking you in advance
seemanth m



      Get your preferred Email name!
Now you can @ymail.com and @rocketmail.com. 
http://mail.promotions.yahoo.com/newdomains/aa/
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100215/43c3ed3f/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot_hurricane_center3.gs
Type: application/octet-stream
Size: 9941 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20100215/43c3ed3f/attachment.obj 

From jithu9907 at GMAIL.COM  Tue Feb 16 09:23:46 2010
From: jithu9907 at GMAIL.COM (sreejith k.s)
Date: Tue, 16 Feb 2010 19:53:46 +0530
Subject: aggregation of discriptor files
Message-ID: <TUE.16.FEB.2010.195346.0530.>

hai,
     i have 12 data sets and 12 discriptor files to read this 12 data files.
i want to make one discriptor file that can read all these 12 data sets. i
have combined 12 data set, now i want to combine these 12 discriptor files
into one which can read all the 12 data files .The sequence of data file is
like t001_mon.unf,t002_mon.unf ........... t012_mon.unf


with hope,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100216/b71a48cd/attachment.html 

From hersala at GMAIL.COM  Tue Feb 16 09:43:12 2010
From: hersala at GMAIL.COM (hersala)
Date: Tue, 16 Feb 2010 11:43:12 -0300
Subject: aggregation of discriptor files
In-Reply-To: <971af1ba1002160623s2b5da8c8q269ba6e59af1ecac@mail.gmail.com>
Message-ID: <TUE.16.FEB.2010.114312.0300.>

Hello Sreejith
You should use *options template* in your ctl file to read all the files as
if they were only one.
Please read this reference:

http://www.iges.org/grads/gadoc/gadocindex.html

I also send you an example.
Cheers

Hern?n


2010/2/16 sreejith k.s <jithu9907 at gmail.com>

> hai,
>      i have 12 data sets and 12 discriptor files to read this 12 data
> files. i want to make one discriptor file that can read all these 12 data
> sets. i have combined 12 data set, now i want to combine these 12 discriptor
> files into one which can read all the 12 data files .The sequence of data
> file is like t001_mon.unf,t002_mon.unf ........... t012_mon.unf
>
>
> with hope,
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100216/024d786a/attachment.html 
-------------- next part --------------
dset air.sig995.%y4.nc 
options template
tdef time 11323 linear 1jan1979 1dy 

*The number 11534 indicates the number of days from 01jan1979. You should adapt this number at your convenience.



From jithu9907 at GMAIL.COM  Tue Feb 16 09:57:35 2010
From: jithu9907 at GMAIL.COM (sreejith k.s)
Date: Tue, 16 Feb 2010 20:27:35 +0530
Subject: aggregation of discriptor files
In-Reply-To: <5b40c91002160643t63e087d0we53b71e8aea625b0@mail.gmail.com>
Message-ID: <TUE.16.FEB.2010.202735.0530.>

thanks,
           i would like to know if you have done this. because iam not able
ti display the variable. it  is giving a wrong display


On Tue, Feb 16, 2010 at 8:13 PM, hersala <hersala at gmail.com> wrote:

> Hello Sreejith
> You should use *options template* in your ctl file to read all the files
> as if they were only one.
> Please read this reference:
>
> http://www.iges.org/grads/gadoc/gadocindex.html
>
> I also send you an example.
> Cheers
>
> Hern?n
>
>
> 2010/2/16 sreejith k.s <jithu9907 at gmail.com>
>
> hai,
>>      i have 12 data sets and 12 discriptor files to read this 12 data
>> files. i want to make one discriptor file that can read all these 12 data
>> sets. i have combined 12 data set, now i want to combine these 12 discriptor
>> files into one which can read all the 12 data files .The sequence of data
>> file is like t001_mon.unf,t002_mon.unf ........... t012_mon.unf
>>
>>
>> with hope,
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100216/4e23abfd/attachment.html 

From hersala at GMAIL.COM  Tue Feb 16 10:53:57 2010
From: hersala at GMAIL.COM (hersala)
Date: Tue, 16 Feb 2010 12:53:57 -0300
Subject: aggregation of discriptor files
In-Reply-To: <971af1ba1002160657o5fe9a057sc4000229614dbace@mail.gmail.com>
Message-ID: <TUE.16.FEB.2010.125357.0300.>

Yes, I used with no problem, of course after some errors and some help.
Hernan



2010/2/16 sreejith k.s <jithu9907 at gmail.com>

> thanks,
>            i would like to know if you have done this. because iam not able
> ti display the variable. it  is giving a wrong display
>
>
>
> On Tue, Feb 16, 2010 at 8:13 PM, hersala <hersala at gmail.com> wrote:
>
>> Hello Sreejith
>> You should use *options template* in your ctl file to read all the files
>> as if they were only one.
>> Please read this reference:
>>
>> http://www.iges.org/grads/gadoc/gadocindex.html
>>
>> I also send you an example.
>> Cheers
>>
>> Hern?n
>>
>>
>> 2010/2/16 sreejith k.s <jithu9907 at gmail.com>
>>
>> hai,
>>>      i have 12 data sets and 12 discriptor files to read this 12 data
>>> files. i want to make one discriptor file that can read all these 12 data
>>> sets. i have combined 12 data set, now i want to combine these 12 discriptor
>>> files into one which can read all the 12 data files .The sequence of data
>>> file is like t001_mon.unf,t002_mon.unf ........... t012_mon.unf
>>>
>>>
>>> with hope,
>>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100216/9480f2fa/attachment.html 

From jithu9907 at GMAIL.COM  Tue Feb 16 11:05:16 2010
From: jithu9907 at GMAIL.COM (sreejith k.s)
Date: Tue, 16 Feb 2010 21:35:16 +0530
Subject: aggregation of discriptor files
In-Reply-To: <5b40c91002160753h18d0dbccj8fd91c980edc414b@mail.gmail.com>
Message-ID: <TUE.16.FEB.2010.213516.0530.>

i will attach the discriptor file could you please tell me if it is right or
wrong. The data file is t001_mat_33.unf , t002_mat_33.unf and till
t012_mat_33.unf. 1 month interval.


On Tue, Feb 16, 2010 at 9:23 PM, hersala <hersala at gmail.com> wrote:

> Yes, I used with no problem, of course after some errors and some help.
> Hernan
>
>
>
>
> 2010/2/16 sreejith k.s <jithu9907 at gmail.com>
>
>> thanks,
>>            i would like to know if you have done this. because iam not
>> able ti display the variable. it  is giving a wrong display
>>
>>
>>
>> On Tue, Feb 16, 2010 at 8:13 PM, hersala <hersala at gmail.com> wrote:
>>
>>> Hello Sreejith
>>> You should use *options template* in your ctl file to read all the files
>>> as if they were only one.
>>> Please read this reference:
>>>
>>> http://www.iges.org/grads/gadoc/gadocindex.html
>>>
>>> I also send you an example.
>>> Cheers
>>>
>>> Hern?n
>>>
>>>
>>> 2010/2/16 sreejith k.s <jithu9907 at gmail.com>
>>>
>>> hai,
>>>>      i have 12 data sets and 12 discriptor files to read this 12 data
>>>> files. i want to make one discriptor file that can read all these 12 data
>>>> sets. i have combined 12 data set, now i want to combine these 12 discriptor
>>>> files into one which can read all the 12 data files .The sequence of data
>>>> file is like t001_mon.unf,t002_mon.unf ........... t012_mon.unf
>>>>
>>>>
>>>> with hope,
>>>>
>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100216/7633385a/attachment.html 

From jithu9907 at GMAIL.COM  Tue Feb 16 11:06:52 2010
From: jithu9907 at GMAIL.COM (sreejith k.s)
Date: Tue, 16 Feb 2010 21:36:52 +0530
Subject: aggregation of discriptor files
In-Reply-To: <971af1ba1002160805l231d7427pbccded27fbefb9bb@mail.gmail.com>
Message-ID: <TUE.16.FEB.2010.213652.0530.>

On Tue, Feb 16, 2010 at 9:35 PM, sreejith k.s <jithu9907 at gmail.com> wrote:

> i will attach the discriptor file could you please tell me if it is right
> or wrong. The data file is t001_mat_33.unf , t002_mat_33.unf and till
> t012_mat_33.unf. 1 month interval.
>
>
> On Tue, Feb 16, 2010 at 9:23 PM, hersala <hersala at gmail.com> wrote:
>
>> Yes, I used with no problem, of course after some errors and some help.
>> Hernan
>>
>>
>>
>>
>> 2010/2/16 sreejith k.s <jithu9907 at gmail.com>
>>
>>> thanks,
>>>            i would like to know if you have done this. because iam not
>>> able ti display the variable. it  is giving a wrong display
>>>
>>>
>>>
>>> On Tue, Feb 16, 2010 at 8:13 PM, hersala <hersala at gmail.com> wrote:
>>>
>>>> Hello Sreejith
>>>> You should use *options template* in your ctl file to read all the
>>>> files as if they were only one.
>>>> Please read this reference:
>>>>
>>>> http://www.iges.org/grads/gadoc/gadocindex.html
>>>>
>>>> I also send you an example.
>>>> Cheers
>>>>
>>>> Hern?n
>>>>
>>>>
>>>> 2010/2/16 sreejith k.s <jithu9907 at gmail.com>
>>>>
>>>> hai,
>>>>>      i have 12 data sets and 12 discriptor files to read this 12 data
>>>>> files. i want to make one discriptor file that can read all these 12 data
>>>>> sets. i have combined 12 data set, now i want to combine these 12 discriptor
>>>>> files into one which can read all the 12 data files .The sequence of data
>>>>> file is like t001_mon.unf,t002_mon.unf ........... t012_mon.unf
>>>>>
>>>>>
>>>>> with hope,
>>>>>
>>>>
>>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100216/bcc38d47/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: DDF.LEVMON_0.25_T01_IND
Type: application/octet-stream
Size: 881 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20100216/bcc38d47/attachment.obj 

From hmjbarbosa at GMAIL.COM  Tue Feb 16 15:03:20 2010
From: hmjbarbosa at GMAIL.COM (Henrique Barbosa)
Date: Tue, 16 Feb 2010 18:03:20 -0200
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <7D3751A6-23B1-4213-A6DB-92E459AF9F56@cola.iges.org>
Message-ID: <TUE.16.FEB.2010.180320.0200.>

Dear Jennifer,

> And if you post questions to the
> list, be specific about the version you're using and the operating system
> you're running on as well as the data file (or URL) you're trying to read.
> --Jennifer

In my original message I have said I am using OpenGrads v2.0.a7.oga.3,
but ok, I forgot to say its on linux 64bits, ubuntu 8.10 server :-)

>> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do

I also posted the link to the data file:

>> ga-> sdfopen
>> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2

Could you try to open it?

The output of 'q config' and the crash message I get are attached to
my original message.

Thanks in advance,
Henrique


From jma at COLA.IGES.ORG  Tue Feb 16 16:32:38 2010
From: jma at COLA.IGES.ORG (Jennifer Adams)
Date: Tue, 16 Feb 2010 16:32:38 -0500
Subject: Problem to open grib2 Nomad/CFSR data with opengrads
In-Reply-To: <e71668fa1002161203h4bdfbefayb865834bb3238fa1@mail.gmail.com>
Message-ID: <TUE.16.FEB.2010.163238.0500.>


On Feb 16, 2010, at 3:03 PM, Henrique Barbosa wrote:

> Dear Jennifer,
>
>> And if you post questions to the
>> list, be specific about the version you're using and the operating  
>> system
>> you're running on as well as the data file (or URL) you're trying  
>> to read.
>> --Jennifer
>
> In my original message I have said I am using OpenGrads v2.0.a7.oga.3,
> but ok, I forgot to say its on linux 64bits, ubuntu 8.10 server :-)

Very good. It never hurts to remind all users about this -- this  
information really matters for debugging problems.

>
>>> I am using v2.0.a7.oga.3 downloaded from opengrads website. I do
>
> I also posted the link to the data file:
>
>>> ga-> sdfopen
>>> http://nomads.ncdc.noaa.gov/thredds/dodsC/cfsrmon/200912/pgblnl.gdas.200912.grb2
>
> Could you try to open it?

I can open this URL and display variables (even though the names are  
shortened awkwardly) with the latest COLA version of GrADS.
--Jennifer

>
> The output of 'q config' and the crash message I get are attached to
> my original message.
>
> Thanks in advance,
> Henrique

--
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/20100216/21db2c57/attachment.html 

From faisal.saeed at GCISC.ORG.PK  Tue Feb 16 09:25:06 2010
From: faisal.saeed at GCISC.ORG.PK (faisal.saeed)
Date: Tue, 16 Feb 2010 19:25:06 +0500
Subject: automated response
Message-ID: <TUE.16.FEB.2010.192506.0500.>

Please contact me at faisal at misu.su.se ,or  faisal.met at gmail.com


From antoyc at YAHOO.COM  Wed Feb 17 21:04:28 2010
From: antoyc at YAHOO.COM (Antoy Chang)
Date: Wed, 17 Feb 2010 18:04:28 -0800
Subject: Convert binary to Grib 1 FM-92 format
Message-ID: <WED.17.FEB.2010.180428.0800.>

Hi there,
Could someone guide on how this can be done...
Your assistance is highly appreciated
Antoy


      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100217/18c653aa/attachment.html 

From luis.blacutt at GMAIL.COM  Thu Feb 18 11:27:17 2010
From: luis.blacutt at GMAIL.COM (luis blacutt)
Date: Thu, 18 Feb 2010 13:27:17 -0300
Subject: problem opening netcdf file
Message-ID: <THU.18.FEB.2010.132717.0300.>

Dear Grads users

When I try to open a netcdf file I have the following error,
udunits(3): Couldn't open units database "/usr/local/lib/grads/udunits.dat":
No such file or directory
gadsdf: UDUNITS package initialization failure.

I'm using grads  Version 2.0.a7.1

Can you help me please?

Luis
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100218/f4622403/attachment.html 

From jithu9907 at GMAIL.COM  Fri Feb 19 04:49:09 2010
From: jithu9907 at GMAIL.COM (sreejith k.s)
Date: Fri, 19 Feb 2010 15:19:09 +0530
Subject: Fwd: aggregation of discriptor files
In-Reply-To: <971af1ba1002160623s2b5da8c8q269ba6e59af1ecac@mail.gmail.com>
Message-ID: <FRI.19.FEB.2010.151909.0530.>

---------- Forwarded message ----------
From: sreejith k.s <jithu9907 at gmail.com>
Date: Tue, Feb 16, 2010 at 7:53 PM
Subject: aggregation of discriptor files
To: GRADSUSR at list.cineca.it


hai,
     I have 12 data sets and 12 discriptor files to read this 12 data files.
I want to make one discriptor file that can read all these 12 data sets. I
have combined 12 data set, now i want to combine these 12 discriptor files
into one which can read all the 12 data files .The sequence of data file is
like t001_mon.unf,t002_mon.unf ........... t012_mon.unf . If anyone could
help me.


with hope,
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100219/e28e6d04/attachment.html 

From gzepka at DGE.INPE.BR  Sat Feb 20 12:38:34 2010
From: gzepka at DGE.INPE.BR (Gisele Zepka)
Date: Sat, 20 Feb 2010 15:38:34 -0200
Subject: theta-e and cape
Message-ID: <SAT.20.FEB.2010.153834.0200.>

Hello users,

I've simulated surface theta-e and cape to thunderstorm events in 
Brazil. The values of theta-e vary from 300 to 312 K and of cape, from 
500 to 2500 J/kg. What is the variation range to these variables 
according thunderstorm and convection occurrence? Can someone indicate 
references in literature about this?

Thanks in advance.

Gisele

-- 
Msc. Gisele dos Santos Zepka
Atmospheric Electricity Group (ELAT)
National Institute for Space Research (INPE)
Av. dos Astronautas, 1758
S?o Jos? dos Campos/ SP/ Brazil CEP 12227-010
Phone ++55(12)39456841


From barbara.turato at ARPAL.ORG  Sat Feb 20 12:48:31 2010
From: barbara.turato at ARPAL.ORG (barbara.turato)
Date: Sat, 20 Feb 2010 18:48:31 +0100
Subject: automated response
Message-ID: <SAT.20.FEB.2010.184831.0100.>

sar? assente dal 20 al 28 febbraio 2010. Riprender? la lettura della mail il 1 marzo.
grazie


From ivtoman at INET.HR  Sun Feb 21 05:54:33 2010
From: ivtoman at INET.HR (Ivan Toman)
Date: Sun, 21 Feb 2010 11:54:33 +0100
Subject: Wind Direction Change
Message-ID: <SUN.21.FEB.2010.115433.0100.>

Dear GrADS users,

I need some help with scripting. I have to calculate change in wind
direction in degrees between two different times. I do this this way:

'define hwd0 = 57.3*atan2(ugrd10m(t='i'),vgrd10m(t='j'))'
'define hwd1 = 57.3*atan2(ugrd10m(t='j'),vgrd10m(t='j'))'
'define hwdc = hwd0-hwd1'

where hwd0 is present wind direction, and hwd1 is wind direction in past time.

The problem arises when there is, for example, hwd0=10deg, hwd1=350deg, the
result of this equation is of course -340 degrees. This is not good, because
I have to get something like +20 degrees as result (by words: wind changed
direction for 20 degrees to the right).

In fact, when wind is backing result should be negative value (from 0 to
-180), and when wind is veering the result must be positive value (from 0 to
180), and never < -180 or > 180 degrees.

I hope you got the point? Can somebody help me please here? I feel this is a
trivial problem but I'm stuck unfortunately. Thanks a lot.


From jdduda at IASTATE.EDU  Sun Feb 21 09:18:11 2010
From: jdduda at IASTATE.EDU (Jeffrey Duda)
Date: Sun, 21 Feb 2010 08:18:11 -0600
Subject: Wind Direction Change
In-Reply-To: <LISTSERV%201002211154332284.0C09@LIST.CINECA.IT>
Message-ID: <SUN.21.FEB.2010.081811.0600.>

Ivan,
I've tried to do this in the past and have encountered the same problem.  I
have not yet found a solution to it.  You'd really need to do those
calculations in a few DO loops (in FORTRAN or something like it) using an IF
statement in the calculation.  Grads does not do that kind of thing
unfortunately.  You might be able to use the maskout function, however, but
that is just speculation.  Good luck.

Jeff Duda

On Sun, Feb 21, 2010 at 4:54 AM, Ivan Toman <ivtoman at inet.hr> wrote:

> Dear GrADS users,
>
> I need some help with scripting. I have to calculate change in wind
> direction in degrees between two different times. I do this this way:
>
> 'define hwd0 = 57.3*atan2(ugrd10m(t='i'),vgrd10m(t='j'))'
> 'define hwd1 = 57.3*atan2(ugrd10m(t='j'),vgrd10m(t='j'))'
> 'define hwdc = hwd0-hwd1'
>
> where hwd0 is present wind direction, and hwd1 is wind direction in past
> time.
>
> The problem arises when there is, for example, hwd0=10deg, hwd1=350deg, the
> result of this equation is of course -340 degrees. This is not good,
> because
> I have to get something like +20 degrees as result (by words: wind changed
> direction for 20 degrees to the right).
>
> In fact, when wind is backing result should be negative value (from 0 to
> -180), and when wind is veering the result must be positive value (from 0
> to
> 180), and never < -180 or > 180 degrees.
>
> I hope you got the point? Can somebody help me please here? I feel this is
> a
> trivial problem but I'm stuck unfortunately. Thanks a lot.
>



-- 
Jeff Duda
Iowa State University
Meteorology Graduate Student
3134 Agronomy Hall
www.meteor.iastate.edu/~jdduda
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100221/4f127290/attachment.html 

From dipakmath at GMAIL.COM  Mon Feb 22 10:41:18 2010
From: dipakmath at GMAIL.COM (dipak sahu)
Date: Mon, 22 Feb 2010 21:11:18 +0530
Subject: error with while loop
Message-ID: <MON.22.FEB.2010.211118.0530.>

Dear frnds
     here i have some difficulty with the following script, i want to write
the maximum value of daily temperature of 4-day forecast, but it is only
writing a single value (i.e the maximum value of last day only). I didin't
found where is the mistake. Thanks in advance if anybody has any solution to
share on it.The  script is given below.

'reinit'
'open singldom-Nov.ctl'
'set lat 20.25'
'set lon 85.83'
tt = 1
while (tt<97)
'set t 'tt
   'define maxt= max((t2-273.15),t+0,t+23)'
    'set gxout print'
   'fprintf maxt maxt2-bbsr-DA.txt';
    tt=tt+24
  endwhile
~
with regards
Dipak
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100222/9ea7a71f/attachment.html 

From Charles.Seman at NOAA.GOV  Mon Feb 22 16:28:36 2010
From: Charles.Seman at NOAA.GOV (Charles Seman)
Date: Mon, 22 Feb 2010 16:28:36 -0500
Subject: Wind Direction Change
In-Reply-To: <LISTSERV%201002211154332284.0C09@LIST.CINECA.IT>
Message-ID: <MON.22.FEB.2010.162836.0500.>

Ivan,

http://en.wikipedia.org/wiki/Law_of_cosines#Vector_formulation
...try the vector formulation for two wind vectors

Hope this helps,
Chuck

Ivan Toman wrote:
> Dear GrADS users,
>
> I need some help with scripting. I have to calculate change in wind
> direction in degrees between two different times. I do this this way:
>
> 'define hwd0 = 57.3*atan2(ugrd10m(t='i'),vgrd10m(t='j'))'
> 'define hwd1 = 57.3*atan2(ugrd10m(t='j'),vgrd10m(t='j'))'
> 'define hwdc = hwd0-hwd1'
>
> where hwd0 is present wind direction, and hwd1 is wind direction in past time.
>
> The problem arises when there is, for example, hwd0=10deg, hwd1=350deg, the
> result of this equation is of course -340 degrees. This is not good, because
> I have to get something like +20 degrees as result (by words: wind changed
> direction for 20 degrees to the right).
>
> In fact, when wind is backing result should be negative value (from 0 to
> -180), and when wind is veering the result must be positive value (from 0 to
> 180), and never < -180 or > 180 degrees.
>
> I hope you got the point? Can somebody help me please here? I feel this is a
> trivial problem but I'm stuck unfortunately. Thanks a lot.
>   

-- 

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 Charles.Seman at NOAA.GOV  Mon Feb 22 17:02:10 2010
From: Charles.Seman at NOAA.GOV (Charles Seman)
Date: Mon, 22 Feb 2010 17:02:10 -0500
Subject: error with while loop
In-Reply-To: <8562142e1002220741s4bb02697j13e6bcf462540b2d@mail.gmail.com>
Message-ID: <MON.22.FEB.2010.170210.0500.>

Dipak,

Try using the "write()" function: 
http://grads.iges.org/grads/gadoc/script.html#intrinsic and no "set 
gxout print" or if you do, then the "subwrd(...)" statement would have 
to be changed...

'define maxt = ...'
'd maxt' ; maxt = subwrd(result,4)
rc = write('maxt2-bbsr-DA.txt','t = 'tt' maxt = 'maxt)
* or
rc =  write('maxt2-bbsr-DA.txt','t = 'tt' maxt = 'maxt,append)

* ",append" : if you want to append to any existing file on the first 
write; see documentation; try with/without ",append" to experiment...

Hope this helps,
Chuck

dipak sahu wrote:
> Dear frnds
>      here i have some difficulty with the following script, i want to 
> write the maximum value of daily temperature of 4-day forecast, but it 
> is only writing a single value (i.e the maximum value of last day 
> only). I didin't found where is the mistake. Thanks in advance if 
> anybody has any solution to share on it.The  script is given below.
>  
> 'reinit'
> 'open singldom-Nov.ctl'
> 'set lat 20.25'
> 'set lon 85.83'
> tt = 1
> while (tt<97)
> 'set t 'tt
>    'define maxt= max((t2-273.15),t+0,t+23)'
>     'set gxout print'
>    'fprintf maxt maxt2-bbsr-DA.txt';
>     tt=tt+24
>   endwhile
> ~
> with regards
> Dipak
>  

-- 

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 thi_veloso at YAHOO.COM.BR  Tue Feb 23 08:01:44 2010
From: thi_veloso at YAHOO.COM.BR (Thiago Veloso)
Date: Tue, 23 Feb 2010 05:01:44 -0800
Subject: How to modify netCDF header?
Message-ID: <TUE.23.FEB.2010.050144.0800.>

??Dear all,
??I have a precipitation dataset which was distributed in binary files which I converted to netCDF by employing opengrads' lats4d.
??However, some information of the data is missing/incorrect. Please see the attached picture. I'd like to edit the "comments" and "units" fields. Is it possible? Using what software?
??Thanks in advance for any direction,
??Thiago.


      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100223/be59956f/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: netCDF header.png
Type: image/png
Size: 40861 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20100223/be59956f/attachment.png 

From rodrigo at DGF.UCHILE.CL  Tue Feb 23 08:25:40 2010
From: rodrigo at DGF.UCHILE.CL (Rodrigo Sanchez)
Date: Tue, 23 Feb 2010 10:25:40 -0300
Subject: How to modify netCDF header?
In-Reply-To: <503610.67434.qm@web54301.mail.re2.yahoo.com>
Message-ID: <TUE.23.FEB.2010.102540.0300.>

Try,
 http://nco.sourceforge.net/

 look for..
 ncatted, ncrename ...

 Rodrigo Sanchez
 DGF, U. de Chile



> ??Dear all,
> ??I have a precipitation dataset which was distributed in binary files
> which I converted to netCDF by employing opengrads' lats4d.
> ??However, some information of the data is missing/incorrect. Please see
> the attached picture. I'd like to edit the "comments" and "units" fields.
> Is it possible? Using what software?
> ??Thanks in advance for any direction,
> ??Thiago.
>
>
>


From appopson at YAHOO.FR  Tue Feb 23 10:07:28 2010
From: appopson at YAHOO.FR (louis Vonder)
Date: Tue, 23 Feb 2010 15:07:28 +0000
Subject: How to modify netCDF header?
In-Reply-To: <503610.67434.qm@web54301.mail.re2.yahoo.com>
Message-ID: <TUE.23.FEB.2010.150728.0000.>

look to 

http://www.mpimet.mpg.de/~cdo



--- En date de?: Mar 23.2.10, Thiago Veloso <thi_veloso at YAHOO.COM.BR> a ?crit?:

De: Thiago Veloso <thi_veloso at YAHOO.COM.BR>
Objet: How to modify netCDF header?
?: GRADSUSR at LIST.CINECA.IT
Date: Mardi 23 f?vrier 2010, 14h01

??Dear all,
??I have a precipitation dataset which was distributed in binary files which I converted to netCDF by employing opengrads' lats4d.
??However, some information of the data is missing/incorrect. Please see the attached picture. I'd like to edit the "comments" and "units" fields. Is it possible? Using what software?
??Thanks in advance for any direction,
??Thiago.




      


      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100223/883ac0b2/attachment.html 

From mhsanabary9 at YAHOO.COM  Tue Feb 23 16:45:01 2010
From: mhsanabary9 at YAHOO.COM (MMME)
Date: Tue, 23 Feb 2010 13:45:01 -0800
Subject: Reading Group of files
Message-ID: <TUE.23.FEB.2010.134501.0800.>

Hi
?? I have group of files and I want to read all of them by one script. Is it possible?

I have:

i=1
while (i < 9) 'sdfopen BlueNILEB-6i-ECMWF.nc'
'set prnopts %10.8
'set lat 
'set lon 
'set t 1
?
?
'd 
ascii = result
asciifile = 'precipitation.asc'
numbers = sublin(result,2)
if (numbers != '')
rc = write (asciifile,numbers)
else
say 'no ascii data to write'
return
endif
* write out remaining ascii data
i=3
while(1)
numbers = sublin(result,i)
if (numbers != '')
rc = write (asciifile,numbers,append)
else
break
endif
i = i+1
endwhile
rc = close(asciifile)
say 'return code from close is:'rc'set gxout print'e1 1'7.5'35'1461'(cp+lsp+sf)*1000/60/60/6''set gxout contour'
'd cp'i = i + 1 
endwhile 


      


From ivtoman at INET.HR  Wed Feb 24 04:52:14 2010
From: ivtoman at INET.HR (Ivan Toman)
Date: Wed, 24 Feb 2010 10:52:14 +0100
Subject: Reading Group of files
In-Reply-To: <369938.93199.qm@web112501.mail.gq1.yahoo.com>
Message-ID: <WED.24.FEB.2010.105214.0100.>

MMME wrote:
> Hi
>    I have group of files and I want to read all of them by one script. Is it possible?
>
> I have:
>
> i=1
> while (i < 9) 'sdfopen BlueNILEB-6i-ECMWF.nc'
>   
Hello,

this couldn't work, but you may try something like:

i=1
while (i < 9) 'sdfopen BlueNILEB-6'i'-ECMWF.nc'


Not sure though...
Ivan


From hersala at GMAIL.COM  Wed Feb 24 09:23:45 2010
From: hersala at GMAIL.COM (hersala)
Date: Wed, 24 Feb 2010 11:23:45 -0300
Subject: Reading Group of files
In-Reply-To: <369938.93199.qm@web112501.mail.gq1.yahoo.com>
Message-ID: <WED.24.FEB.2010.112345.0300.>

Hi MMME:

As far as I know, it is possible to open several files (at least .nc files)
using Templates. Please see
http://www.iges.org/grads/gadoc/gadocindex.htmlfor details.
In my case I used two files to do that, one control file and a txt file with
some command lines. I have attached both files, perhaps it serves you as
example.
Cheers,

Hern?n Sala
DCT-IAA


2010/2/23 MMME <mhsanabary9 at yahoo.com>

> Hi
>    I have group of files and I want to read all of them by one script. Is
> it possible?
>
> I have:
>
> i=1
> while (i < 9) 'sdfopen BlueNILEB-6i-ECMWF.nc'
> 'set prnopts %10.8
> 'set lat
> 'set lon
> 'set t 1
>
>
> 'd
> ascii = result
> asciifile = 'precipitation.asc'
> numbers = sublin(result,2)
> if (numbers != '')
> rc = write (asciifile,numbers)
> else
> say 'no ascii data to write'
> return
> endif
> * write out remaining ascii data
> i=3
> while(1)
> numbers = sublin(result,i)
> if (numbers != '')
> rc = write (asciifile,numbers,append)
> else
> break
> endif
> i = i+1
> endwhile
> rc = close(asciifile)
> say 'return code from close is:'rc'set gxout print'e1
> 1'7.5'35'1461'(cp+lsp+sf)*1000/60/60/6''set gxout contour'
> 'd cp'i = i + 1
> endwhile
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100224/3cd03aa9/attachment.html 
-------------- next part --------------
reinit

*This file allow to open simultaneously all the following files: from air.sig995.1975.nc to air.sig995.2005.nc
*It works associated with the control file daily.ctl
*In the way they have been written, both must be placed in the same directory with the data files.
*Then it creates a temporal series using aave function.
*It works properly in gradsNC-v1.8 and newer versions

xdfopen daily.ctl

*The content of daily.ctl file is:
*dset air.sig995.%y4.nc 
*options template
*tdef time 11323 linear 1jan1975 1dy

set lat 0
set lon 0

*Establish marker size
*set digsiz 0.02

*Sets color and line thickness
*set annot 1 1

*set the style of the marker for line plots
*set cmark 0


*%y4 replace the part of the name file which has 4 digit year


*This number should be set in function of the t units which are desired to analyse.
*In this case t unit is day. So, 11323 days beginning in 01ene1975 will be shown.
*this means 31dec2005

set t 1 11323

tt= aave(air, lon=275, lon=295, lat=-73, lat=-68) - 273.15

d tt
-------------- next part --------------
A non-text attachment was scrubbed...
Name: daily.ctl
Type: application/octet-stream
Size: 145 bytes
Desc: not available
Url : http://gradsusr.org/pipermail/gradsusr/attachments/20100224/3cd03aa9/attachment.obj 

From jgmsantos at GMAIL.COM  Wed Feb 24 10:06:01 2010
From: jgmsantos at GMAIL.COM (Guilherme Martins)
Date: Wed, 24 Feb 2010 13:06:01 -0200
Subject: Reading Group of files
In-Reply-To: <5b40c91002240623l6611c042u81946a45fa66d7e0@mail.gmail.com>
Message-ID: <WED.24.FEB.2010.130601.0200.>

Try something like

For example, you have file this way:

u.1990.nc
u.1991.nc
u.1992.nc

Just do it:

'sdfopen u.1990.nc u.%y4.nc 100'
'd variable'

u.1990.nc = your fist nc file
u.%y4.nc = your template
100          = numbers times you want display

Guilherme.



2010/2/24 hersala <hersala at gmail.com>

> Hi MMME:
>
> As far as I know, it is possible to open several files (at least .nc files)
> using Templates. Please see
> http://www.iges.org/grads/gadoc/gadocindex.html for details.
> In my case I used two files to do that, one control file and a txt file
> with some command lines. I have attached both files, perhaps it serves you
> as example.
> Cheers,
>
> Hern?n Sala
> DCT-IAA
>
>
> 2010/2/23 MMME <mhsanabary9 at yahoo.com>
>
> Hi
>>    I have group of files and I want to read all of them by one script. Is
>> it possible?
>>
>> I have:
>>
>> i=1
>> while (i < 9) 'sdfopen BlueNILEB-6i-ECMWF.nc'
>> 'set prnopts %10.8
>> 'set lat
>> 'set lon
>> 'set t 1
>>
>>
>> 'd
>> ascii = result
>> asciifile = 'precipitation.asc'
>> numbers = sublin(result,2)
>> if (numbers != '')
>> rc = write (asciifile,numbers)
>> else
>> say 'no ascii data to write'
>> return
>> endif
>> * write out remaining ascii data
>> i=3
>> while(1)
>> numbers = sublin(result,i)
>> if (numbers != '')
>> rc = write (asciifile,numbers,append)
>> else
>> break
>> endif
>> i = i+1
>> endwhile
>> rc = close(asciifile)
>> say 'return code from close is:'rc'set gxout print'e1
>> 1'7.5'35'1461'(cp+lsp+sf)*1000/60/60/6''set gxout contour'
>> 'd cp'i = i + 1
>> endwhile
>>
>>
>>
>>
>


-- 
-----------------------------------------------------------
Jos? Guilherme Martins dos Santos
Grupo de Modelagem da Atmosfera e Interfaces - GMAI
Divis?o de Modelagem e Desenvolvimento - DMD
Centro de Previs?o de Tempo e Estudos Clim?ticos - INPE
Rodovia Presidente Dutra, km 39 - CEP 12630-000
Cachoeira Paulista, SP, Brasil
Telefone(CP): +55 12 3186 8502
Celular(SJC): +55 12 91678770 / +55 12 81528580
E-mail: guilherme.martins at cptec.inpe.br
E-mail alternativo: jgmsantos at gmail.com
Homepage: http://guilhermemartins.hd1.com.br/
-----------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100224/1b010c4c/attachment.html 

From roger.torres at CPTEC.INPE.BR  Wed Feb 24 13:59:44 2010
From: roger.torres at CPTEC.INPE.BR (Roger Rodrigues Torres)
Date: Wed, 24 Feb 2010 15:59:44 -0300
Subject: Netcdf problem
Message-ID: <WED.24.FEB.2010.155944.0300.>

Dear Grads users,

I have some problems in opening some CMIP3 netcdf data. The problem is
attached below.
I tried to find some CDO command to convert the data to grib or binary
format but I wasn't succeeded.
Someone could help me please?

Cheers,

Roger Torres

ga-> sdfopen pcmdi.ipcc4.csiro_mk3_0.20c3m.run1.monthly.tas_A1.nc
Scanning self-describing file:
pcmdi.ipcc4.csiro_mk3_0.20c3m.run1.monthly.tas_A1.nc
SDF Error: 365 day calendars are no longer supported by sdfopen.
  To open this file with GrADS, use a descriptor file with
  a complete TDEF entry and OPTIONS 365_day_calendar.
  Documentation is at http://iges.org/grads/gadoc/SDFdescriptorfile.html

-- 
Roger Rodrigues Torres
Ph.D. Student
National Institute for Space Research (INPE)
Center for Weather Forecast and Climate Studies (CPTEC)
Phone: 55 12 3186-9541
Personal Webpage: http://rtorres.webnode.com.br
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100224/a879721e/attachment.html 

From jgmsantos at GMAIL.COM  Wed Feb 24 14:25:52 2010
From: jgmsantos at GMAIL.COM (Guilherme Martins)
Date: Wed, 24 Feb 2010 17:25:52 -0200
Subject: Netcdf problem
In-Reply-To: <ec7a8b9b1002241059g7e878b72n17a0f6fee864e698@mail.gmail.com>
Message-ID: <WED.24.FEB.2010.172552.0200.>

J? usaste o lats4d.gs para gerar em outro formato???

GUilherme

2010/2/24 Roger Rodrigues Torres <roger.torres at cptec.inpe.br>

> Dear Grads users,
>
> I have some problems in opening some CMIP3 netcdf data. The problem is
> attached below.
> I tried to find some CDO command to convert the data to grib or binary
> format but I wasn't succeeded.
> Someone could help me please?
>
> Cheers,
>
> Roger Torres
>
> ga-> sdfopen pcmdi.ipcc4.csiro_mk3_0.20c3m.run1.monthly.tas_A1.nc
> Scanning self-describing file:
> pcmdi.ipcc4.csiro_mk3_0.20c3m.run1.monthly.tas_A1.nc
> SDF Error: 365 day calendars are no longer supported by sdfopen.
>   To open this file with GrADS, use a descriptor file with
>   a complete TDEF entry and OPTIONS 365_day_calendar.
>   Documentation is at http://iges.org/grads/gadoc/SDFdescriptorfile.html
>
> --
> Roger Rodrigues Torres
> Ph.D. Student
> National Institute for Space Research (INPE)
> Center for Weather Forecast and Climate Studies (CPTEC)
> Phone: 55 12 3186-9541
> Personal Webpage: http://rtorres.webnode.com.br
>



-- 
-----------------------------------------------------------
Jos? Guilherme Martins dos Santos
Grupo de Modelagem da Atmosfera e Interfaces - GMAI
Divis?o de Modelagem e Desenvolvimento - DMD
Centro de Previs?o de Tempo e Estudos Clim?ticos - INPE
Rodovia Presidente Dutra, km 39 - CEP 12630-000
Cachoeira Paulista, SP, Brasil
Telefone(CP): +55 12 3186 8502
Celular(SJC): +55 12 91678770 / +55 12 81528580
E-mail: guilherme.martins at cptec.inpe.br
E-mail alternativo: jgmsantos at gmail.com
Homepage: http://guilhermemartins.hd1.com.br/
-----------------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100224/aced9059/attachment.html 

From roger.torres at CPTEC.INPE.BR  Wed Feb 24 14:36:10 2010
From: roger.torres at CPTEC.INPE.BR (Roger Rodrigues Torres)
Date: Wed, 24 Feb 2010 16:36:10 -0300
Subject: Netcdf problem
In-Reply-To: <ae12c391002241125v6ef4dd20m330cae3b4cf5a753@mail.gmail.com>
Message-ID: <WED.24.FEB.2010.163610.0300.>

E ai Guilherme, blz?

Ainda n?o usei. Na verdade eu nunca usei o lats4d para trocar o formato, Eu
acabo usando aquele software do Max Plank, o CDO. Vc sabe onde encontro os
comandos do lats4d?

Abra??o,

Roger

2010/2/24 Guilherme Martins <jgmsantos at gmail.com>

> J? usaste o lats4d.gs para gerar em outro formato???
>
> GUilherme
>
> 2010/2/24 Roger Rodrigues Torres <roger.torres at cptec.inpe.br>
>
> Dear Grads users,
>>
>> I have some problems in opening some CMIP3 netcdf data. The problem is
>> attached below.
>> I tried to find some CDO command to convert the data to grib or binary
>> format but I wasn't succeeded.
>> Someone could help me please?
>>
>> Cheers,
>>
>> Roger Torres
>>
>> ga-> sdfopen pcmdi.ipcc4.csiro_mk3_0.20c3m.run1.monthly.tas_A1.nc
>> Scanning self-describing file:
>> pcmdi.ipcc4.csiro_mk3_0.20c3m.run1.monthly.tas_A1.nc
>> SDF Error: 365 day calendars are no longer supported by sdfopen.
>>   To open this file with GrADS, use a descriptor file with
>>   a complete TDEF entry and OPTIONS 365_day_calendar.
>>   Documentation is at http://iges.org/grads/gadoc/SDFdescriptorfile.html
>>
>> --
>> Roger Rodrigues Torres
>> Ph.D. Student
>> National Institute for Space Research (INPE)
>> Center for Weather Forecast and Climate Studies (CPTEC)
>> Phone: 55 12 3186-9541
>> Personal Webpage: http://rtorres.webnode.com.br
>>
>
>
>
> --
> -----------------------------------------------------------
> Jos? Guilherme Martins dos Santos
> Grupo de Modelagem da Atmosfera e Interfaces - GMAI
> Divis?o de Modelagem e Desenvolvimento - DMD
> Centro de Previs?o de Tempo e Estudos Clim?ticos - INPE
> Rodovia Presidente Dutra, km 39 - CEP 12630-000
> Cachoeira Paulista, SP, Brasil
> Telefone(CP): +55 12 3186 8502
> Celular(SJC): +55 12 91678770 / +55 12 81528580
> E-mail: guilherme.martins at cptec.inpe.br
> E-mail alternativo: jgmsantos at gmail.com
> Homepage: http://guilhermemartins.hd1.com.br/
> -----------------------------------------------------------
>



-- 
Roger Rodrigues Torres
Ph.D. Student
National Institute for Space Research (INPE)
Center for Weather Forecast and Climate Studies (CPTEC)
Phone: 55 12 3186-9541
Personal Webpage: http://rtorres.webnode.com.br
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100224/9f0e1dd7/attachment.html 

From lalvarez at IDEMOC.COM  Thu Feb 25 14:39:21 2010
From: lalvarez at IDEMOC.COM (Leopoldo Alvarez)
Date: Thu, 25 Feb 2010 20:39:21 +0100
Subject: Air Density script -
Message-ID: <THU.25.FEB.2010.203921.0100.>

Dear GrADS-users,

Anyone have any script that calculates only Air Density (kg/m3) ?? for a point
given (i,j) and for any length of time series...

cheers, 

Leopoldo ??lvarez

--
Leopoldo ??lvarez
Director de Explotaci??n T??cnica.
Instituto de Estudios Meteorol??gicos y Oceanogr??ficos S.L. (IDEMOC).

www.idemoc.com


From jmssnyder at UCDAVIS.EDU  Fri Feb 26 13:27:09 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Fri, 26 Feb 2010 10:27:09 -0800
Subject: Grads for Ubuntu Problem
Message-ID: <FRI.26.FEB.2010.102709.0800.>

I have grads on Ubuntu 9.10 but have been having problems with opening up
netcdf files on it.  The program opens all right but I get the following
message when I open a netcdf file:

Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds hdf5
geotiff
Issue 'q config' command for more information.
Landscape mode? ('n' for portrait):  no
GX Package Initialization: Size = 8.5 11
ga-> sdfopen uwnd.1998.nc
Scanning self-describing file:  uwnd.1998.nc
udunits(3): Couldn't open units database "/usr/local/lib/grads/udunits.dat":
No such file or directory
gadsdf: UDUNITS package initialization failure.
ga->

How do I resolve this problem?

Thanks,

Jason
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100226/23ec58cf/attachment.html 

From luis.blacutt at GMAIL.COM  Sat Feb 27 08:49:52 2010
From: luis.blacutt at GMAIL.COM (luis blacutt)
Date: Sat, 27 Feb 2010 09:49:52 -0400
Subject: Grads for Ubuntu Problem
In-Reply-To: <fe8efa01002261027w4794de89w30e070b8b0c7ba82@mail.gmail.com>
Message-ID: <SAT.27.FEB.2010.094952.0400.>

Hi Jason,

I think you are running on the same problem I run. You need to download the
data file

check this from http://www.iges.org/grads/downloads.html

Supplementary Files
The font and map files are supplementary data sets that are required in
order to run GrADS. Their default location is: /usr/local/lib/grads/. If you
do not have write permission for this directory you can place the files
elsewhere, but you must also change the environment variable GADDIR so the
GrADS executables will know where to find these files. You can download the
data files separately by clicking here:
*data.tar.Z*<ftp://grads.iges.org/grads/data.tar.Z>
 ** <ftp://grads.iges.org/grads/data.tar.Z>.

cd *<dirname>*
uncompress data.tar.Z
tar xvf data.tar

Regards
Luis

On Fri, Feb 26, 2010 at 2:27 PM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> I have grads on Ubuntu 9.10 but have been having problems with opening up
> netcdf files on it.  The program opens all right but I get the following
> message when I open a netcdf file:
>
> Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds hdf5
> geotiff
> Issue 'q config' command for more information.
> Landscape mode? ('n' for portrait):  no
> GX Package Initialization: Size = 8.5 11
> ga-> sdfopen uwnd.1998.nc
> Scanning self-describing file:  uwnd.1998.nc
> udunits(3): Couldn't open units database
> "/usr/local/lib/grads/udunits.dat": No such file or directory
> gadsdf: UDUNITS package initialization failure.
> ga->
>
> How do I resolve this problem?
>
> Thanks,
>
> Jason
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100227/20359b47/attachment.html 

From jmssnyder at UCDAVIS.EDU  Sun Feb 28 14:37:11 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Sun, 28 Feb 2010 11:37:11 -0800
Subject: Grads for Ubuntu Problem
In-Reply-To: <5825ce3f1002270549n6ec86ab3p331e6ecb0a92e74b@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.113711.0800.>

I am still ahving the same problems even after downloading data.tar.  Where
do I find the variable GADDIR?

Thanks,

jason

On Sat, Feb 27, 2010 at 5:49 AM, luis blacutt <luis.blacutt at gmail.com>wrote:

> Hi Jason,
>
> I think you are running on the same problem I run. You need to download the
> data file
>
> check this from http://www.iges.org/grads/downloads.html
>
> Supplementary Files
> The font and map files are supplementary data sets that are required in
> order to run GrADS. Their default location is: /usr/local/lib/grads/. If
> you do not have write permission for this directory you can place the files
> elsewhere, but you must also change the environment variable GADDIR so the
> GrADS executables will know where to find these files. You can download the
> data files separately by clicking here: *data.tar.Z*<ftp://grads.iges.org/grads/data.tar.Z>
>  ** <ftp://grads.iges.org/grads/data.tar.Z>.
>
> cd *<dirname>*
> uncompress data.tar.Z
> tar xvf data.tar
>
> Regards
> Luis
>
> On Fri, Feb 26, 2010 at 2:27 PM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:
>
>> I have grads on Ubuntu 9.10 but have been having problems with opening up
>> netcdf files on it.  The program opens all right but I get the following
>> message when I open a netcdf file:
>>
>> Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds hdf5
>> geotiff
>> Issue 'q config' command for more information.
>> Landscape mode? ('n' for portrait):  no
>> GX Package Initialization: Size = 8.5 11
>> ga-> sdfopen uwnd.1998.nc
>> Scanning self-describing file:  uwnd.1998.nc
>> udunits(3): Couldn't open units database
>> "/usr/local/lib/grads/udunits.dat": No such file or directory
>> gadsdf: UDUNITS package initialization failure.
>> ga->
>>
>> How do I resolve this problem?
>>
>> Thanks,
>>
>> Jason
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/3de8b271/attachment.html 

From luis.blacutt at GMAIL.COM  Sun Feb 28 15:21:56 2010
From: luis.blacutt at GMAIL.COM (luis blacutt)
Date: Sun, 28 Feb 2010 16:21:56 -0400
Subject: Grads for Ubuntu Problem
In-Reply-To: <fe8efa01002281137g497b7f54k1773e8d2295a50b0@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.162156.0400.>

Have you unziped the file and installed following the instructions?
About the GADDIR file don't worry, I haven't made anything with it and still
works fine

Luis

On Sun, Feb 28, 2010 at 3:37 PM, Jason Snyder <jmssnyder at ucdavis.edu> wrote:

> I am still ahving the same problems even after downloading data.tar.  Where
> do I find the variable GADDIR?
>
> Thanks,
>
> jason
>
>
> On Sat, Feb 27, 2010 at 5:49 AM, luis blacutt <luis.blacutt at gmail.com>wrote:
>
>> Hi Jason,
>>
>> I think you are running on the same problem I run. You need to download
>> the data file
>>
>> check this from http://www.iges.org/grads/downloads.html
>>
>> Supplementary Files
>> The font and map files are supplementary data sets that are required in
>> order to run GrADS. Their default location is: /usr/local/lib/grads/. If
>> you do not have write permission for this directory you can place the files
>> elsewhere, but you must also change the environment variable GADDIR so
>> the GrADS executables will know where to find these files. You can download
>> the data files separately by clicking here: *data.tar.Z*<ftp://grads.iges.org/grads/data.tar.Z>
>>  ** <ftp://grads.iges.org/grads/data.tar.Z>.
>>
>> cd *<dirname>*
>> uncompress data.tar.Z
>> tar xvf data.tar
>>
>> Regards
>> Luis
>>
>> On Fri, Feb 26, 2010 at 2:27 PM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:
>>
>>> I have grads on Ubuntu 9.10 but have been having problems with opening up
>>> netcdf files on it.  The program opens all right but I get the following
>>> message when I open a netcdf file:
>>>
>>> Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds hdf5
>>> geotiff
>>> Issue 'q config' command for more information.
>>> Landscape mode? ('n' for portrait):  no
>>> GX Package Initialization: Size = 8.5 11
>>> ga-> sdfopen uwnd.1998.nc
>>> Scanning self-describing file:  uwnd.1998.nc
>>> udunits(3): Couldn't open units database
>>> "/usr/local/lib/grads/udunits.dat": No such file or directory
>>> gadsdf: UDUNITS package initialization failure.
>>> ga->
>>>
>>> How do I resolve this problem?
>>>
>>> Thanks,
>>>
>>> Jason
>>>
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/f7add27e/attachment.html 

From jmssnyder at UCDAVIS.EDU  Sun Feb 28 15:24:07 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Sun, 28 Feb 2010 12:24:07 -0800
Subject: Grads for Ubuntu Problem
In-Reply-To: <fe8efa01002261027w4794de89w30e070b8b0c7ba82@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.122407.0800.>

still having the same problem and the same messages as below.  How do I
modify the paths of GADDIR ,etc.:

on my computer GrADS in the following directory:

/home/snyderjm/Desktop/Grads-2.0.a7.1/bin

I have the following files in my Grads folder after installing data.tar

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ ls
bufrscan   font1.dat  font5.dat  gribmap   gxtran   mres
COPYRIGHT  font2.dat  grads      gribscan  hires    stnmap
data.tar   font3.dat  gradsdap   gxeps     INSTALL  udunits.dat
font0.dat  font4.dat  grib2scan  gxps      lowres   wgrib

Based on this can you help me out with setting up GrADS to run properly.

I tried doing setenv in the csh.  Am I doing this correctly?

-Jason



On Fri, Feb 26, 2010 at 10:27 AM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:

> I have grads on Ubuntu 9.10 but have been having problems with opening up
> netcdf files on it.  The program opens all right but I get the following
> message when I open a netcdf file:
>
> Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds hdf5
> geotiff
> Issue 'q config' command for more information.
> Landscape mode? ('n' for portrait):  no
> GX Package Initialization: Size = 8.5 11
> ga-> sdfopen uwnd.1998.nc
> Scanning self-describing file:  uwnd.1998.nc
> udunits(3): Couldn't open units database
> "/usr/local/lib/grads/udunits.dat": No such file or directory
> gadsdf: UDUNITS package initialization failure.
> ga->
>
> How do I resolve this problem?
>
> Thanks,
>
> Jason
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/8409fba6/attachment.html 

From jmssnyder at UCDAVIS.EDU  Sun Feb 28 15:24:31 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Sun, 28 Feb 2010 12:24:31 -0800
Subject: Grads for Ubuntu Problem
In-Reply-To: <5825ce3f1002270549n6ec86ab3p331e6ecb0a92e74b@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.122431.0800.>

still having the same problem and the same messages as below.  How do I
modify the paths of GADDIR ,etc.:

on my computer GrADS in the following directory:

/home/snyderjm/Desktop/Grads-2.0.a7.1/bin

I have the following files in my Grads folder after installing data.tar

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ ls
bufrscan   font1.dat  font5.dat  gribmap   gxtran   mres
COPYRIGHT  font2.dat  grads      gribscan  hires    stnmap
data.tar   font3.dat  gradsdap   gxeps     INSTALL  udunits.dat
font0.dat  font4.dat  grib2scan  gxps      lowres   wgrib

Based on this can you help me out with setting up GrADS to run properly.

I tried doing setenv in the csh.  Am I doing this correctly?

-Jason

On Sat, Feb 27, 2010 at 5:49 AM, luis blacutt <luis.blacutt at gmail.com>wrote:

> Hi Jason,
>
> I think you are running on the same problem I run. You need to download the
> data file
>
> check this from http://www.iges.org/grads/downloads.html
>
> Supplementary Files
> The font and map files are supplementary data sets that are required in
> order to run GrADS. Their default location is: /usr/local/lib/grads/. If
> you do not have write permission for this directory you can place the files
> elsewhere, but you must also change the environment variable GADDIR so the
> GrADS executables will know where to find these files. You can download the
> data files separately by clicking here: *data.tar.Z*<ftp://grads.iges.org/grads/data.tar.Z>
>  ** <ftp://grads.iges.org/grads/data.tar.Z>.
>
> cd *<dirname>*
> uncompress data.tar.Z
> tar xvf data.tar
>
> Regards
> Luis
>
> On Fri, Feb 26, 2010 at 2:27 PM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:
>
>> I have grads on Ubuntu 9.10 but have been having problems with opening up
>> netcdf files on it.  The program opens all right but I get the following
>> message when I open a netcdf file:
>>
>> Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds hdf5
>> geotiff
>> Issue 'q config' command for more information.
>> Landscape mode? ('n' for portrait):  no
>> GX Package Initialization: Size = 8.5 11
>> ga-> sdfopen uwnd.1998.nc
>> Scanning self-describing file:  uwnd.1998.nc
>> udunits(3): Couldn't open units database
>> "/usr/local/lib/grads/udunits.dat": No such file or directory
>> gadsdf: UDUNITS package initialization failure.
>> ga->
>>
>> How do I resolve this problem?
>>
>> Thanks,
>>
>> Jason
>>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/ad61d53e/attachment.html 

From jmssnyder at UCDAVIS.EDU  Sun Feb 28 15:35:18 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Sun, 28 Feb 2010 12:35:18 -0800
Subject: Grads for Ubuntu Problem
In-Reply-To: <5825ce3f1002281221g54343e10q838b4365c9736c6@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.123518.0800.>

I tried doing this according to the instructions but to no avail.
The problem is is when I do setenv in the unbuntu terminal I get the
following message:

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ setenv GADDIR
/home/snyderjm/Desktop/grads-2.0.a7.1/bin
No command 'setenv' found, did you mean:
 Command 'netenv' from package 'netenv' (universe)
setenv: command not found


Therefore, I tried the following:

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ csh
% setenv GADDIR /home/snyderjm/Desktop/grads-2.0/a7.1/bin
% exit

This is my Grads directory now as it stands with data.tar.gz installed,
unzipped and tarred:

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ ls
bufrscan   font1.dat  font5.dat  gribmap   gxtran   mres          wgrib
COPYRIGHT  font2.dat  grads      gribscan  hires    stnmap
data.tar   font3.dat  gradsdap   gxeps     INSTALL  udunits.dat
font0.dat  font4.dat  grib2scan  gxps      lowres   uwnd.1998.nc

I am continuing to have the same problems with Grads.  How do I resolve
this  and my setenv problems?

Thanks,

Jason



On Sun, Feb 28, 2010 at 12:21 PM, luis blacutt <luis.blacutt at gmail.com>wrote:

> Have you unziped the file and installed following the instructions?
> About the GADDIR file don't worry, I haven't made anything with it and
> still works fine
>
> Luis
>
>
> On Sun, Feb 28, 2010 at 3:37 PM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:
>
>> I am still ahving the same problems even after downloading data.tar.
>> Where do I find the variable GADDIR?
>>
>> Thanks,
>>
>> jason
>>
>>
>> On Sat, Feb 27, 2010 at 5:49 AM, luis blacutt <luis.blacutt at gmail.com>wrote:
>>
>>> Hi Jason,
>>>
>>> I think you are running on the same problem I run. You need to download
>>> the data file
>>>
>>> check this from http://www.iges.org/grads/downloads.html
>>>
>>> Supplementary Files
>>> The font and map files are supplementary data sets that are required in
>>> order to run GrADS. Their default location is: /usr/local/lib/grads/. If
>>> you do not have write permission for this directory you can place the files
>>> elsewhere, but you must also change the environment variable GADDIR so
>>> the GrADS executables will know where to find these files. You can download
>>> the data files separately by clicking here: *data.tar.Z*<ftp://grads.iges.org/grads/data.tar.Z>
>>>  ** <ftp://grads.iges.org/grads/data.tar.Z>.
>>>
>>> cd *<dirname>*
>>> uncompress data.tar.Z
>>> tar xvf data.tar
>>>
>>> Regards
>>> Luis
>>>
>>> On Fri, Feb 26, 2010 at 2:27 PM, Jason Snyder <jmssnyder at ucdavis.edu>wrote:
>>>
>>>> I have grads on Ubuntu 9.10 but have been having problems with opening
>>>> up netcdf files on it.  The program opens all right but I get the following
>>>> message when I open a netcdf file:
>>>>
>>>> Config: v2.0.a7 little-endian readline printim grib2 netcdf hdf4-sds
>>>> hdf5 geotiff
>>>> Issue 'q config' command for more information.
>>>> Landscape mode? ('n' for portrait):  no
>>>> GX Package Initialization: Size = 8.5 11
>>>> ga-> sdfopen uwnd.1998.nc
>>>> Scanning self-describing file:  uwnd.1998.nc
>>>> udunits(3): Couldn't open units database
>>>> "/usr/local/lib/grads/udunits.dat": No such file or directory
>>>> gadsdf: UDUNITS package initialization failure.
>>>> ga->
>>>>
>>>> How do I resolve this problem?
>>>>
>>>> Thanks,
>>>>
>>>> Jason
>>>>
>>>
>>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/636ea9a4/attachment.html 

From Ron.Goodson at EC.GC.CA  Sun Feb 28 16:01:18 2010
From: Ron.Goodson at EC.GC.CA (Goodson,Ron [Edm])
Date: Sun, 28 Feb 2010 14:01:18 -0700
Subject: Grads for Ubuntu Problem
In-Reply-To: A<fe8efa01002281235y78443de9yef5a48dba7a2e81c@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.140118.0700.>

I am reading only this latest message .. so may have missed something
 
setenv is to set the environment variable for C-shell or similar.
 
If you are running Ubuntu - you are probably using bash (look in your
/etc/passwd to see what shell you are using).
 
If so .. instead of "setenv" try
 
export GADDIR=/home/snyderjm/whatever.
 
ron
 


 

 

________________________________

From: GRADSUSR at LIST.CINECA.IT [mailto:GRADSUSR at LIST.CINECA.IT] On Behalf
Of Jason Snyder
Sent: February 28, 2010 1:35 PM
To: GRADSUSR at LIST.CINECA.IT
Subject: Re: Grads for Ubuntu Problem


I tried doing this according to the instructions but to no avail.  
The problem is is when I do setenv in the unbuntu terminal I get the
following message:

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ setenv GADDIR
/home/snyderjm/Desktop/grads-2.0.a7.1/bin
No command 'setenv' found, did you mean:
 Command 'netenv' from package 'netenv' (universe)
setenv: command not found


Therefore, I tried the following:

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ csh
% setenv GADDIR /home/snyderjm/Desktop/grads-2.0/a7.1/bin 
% exit

This is my Grads directory now as it stands with data.tar.gz installed,
unzipped and tarred:

snyderjm at snyderjm-laptop:~/Desktop/grads-2.0.a7.1/bin$ ls
bufrscan   font1.dat  font5.dat  gribmap   gxtran   mres          wgrib
COPYRIGHT  font2.dat  grads      gribscan  hires    stnmap
data.tar   font3.dat  gradsdap   gxeps     INSTALL  udunits.dat
font0.dat  font4.dat  grib2scan  gxps      lowres   uwnd.1998.nc

I am continuing to have the same problems with Grads.  How do I resolve
this  and my setenv problems?

Thanks,

Jason




On Sun, Feb 28, 2010 at 12:21 PM, luis blacutt <luis.blacutt at gmail.com>
wrote:


	Have you unziped the file and installed following the
instructions? 
	About the GADDIR file don't worry, I haven't made anything with
it and still works fine

	Luis 


	On Sun, Feb 28, 2010 at 3:37 PM, Jason Snyder
<jmssnyder at ucdavis.edu> wrote:
	

		I am still ahving the same problems even after
downloading data.tar.  Where do I find the variable GADDIR?
		
		Thanks,
		
		jason 


		On Sat, Feb 27, 2010 at 5:49 AM, luis blacutt
<luis.blacutt at gmail.com> wrote:
		

			Hi Jason, 

			I think you are running on the same problem I
run. You need to download the data file

			check this from
http://www.iges.org/grads/downloads.html

			Supplementary Files
			The font and map files are supplementary data
sets that are required in order to run GrADS. Their default location is:
/usr/local/lib/grads/. If you do not have write permission for this
directory you can place the files elsewhere, but you must also change
the environment variable GADDIR so the GrADS executables will know where
to find these files. You can download the data files separately by
clicking here: data.tar.Z <ftp://grads.iges.org/grads/data.tar.Z>
<ftp://grads.iges.org/grads/data.tar.Z> .

				cd <dirname>
				uncompress data.tar.Z
				tar xvf data.tar
				

			Regards
			Luis

			On Fri, Feb 26, 2010 at 2:27 PM, Jason Snyder
<jmssnyder at ucdavis.edu> wrote:
			

				I have grads on Ubuntu 9.10 but have
been having problems with opening up netcdf files on it.  The program
opens all right but I get the following message when I open a netcdf
file:
				
				Config: v2.0.a7 little-endian readline
printim grib2 netcdf hdf4-sds hdf5 geotiff
				Issue 'q config' command for more
information.
				Landscape mode? ('n' for portrait):  no
				GX Package Initialization: Size = 8.5 11

				ga-> sdfopen uwnd.1998.nc
				Scanning self-describing file:
uwnd.1998.nc
				udunits(3): Couldn't open units database
"/usr/local/lib/grads/udunits.dat": No such file or directory
				gadsdf: UDUNITS package initialization
failure.
				ga-> 
				
				How do I resolve this problem?
				
				Thanks,
				
				Jason
				





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/382bc83a/attachment.html 

From jmssnyder at UCDAVIS.EDU  Sun Feb 28 16:16:09 2010
From: jmssnyder at UCDAVIS.EDU (Jason Snyder)
Date: Sun, 28 Feb 2010 13:16:09 -0800
Subject: Installing GrADS on Ubuntu 9.10
Message-ID: <SUN.28.FEB.2010.131609.0800.>

Can you help give me information on how to install GrADS on Ubuntu?  I tried
doing it on my own but am encountering numerous problems.  Can anyone who
had installed GrADS on Ubuntu 9.10 give me detailed instructions on how to
go about doing this including setting environments, etc.

Thanks,

Jason
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/3b2f5d28/attachment.html 

From thi_veloso at YAHOO.COM.BR  Sun Feb 28 17:42:31 2010
From: thi_veloso at YAHOO.COM.BR (Thiago Veloso)
Date: Sun, 28 Feb 2010 14:42:31 -0800
Subject: Installing GrADS on Ubuntu 9.10
In-Reply-To: <fe8efa01002281316i4df0cc2cu73bde37ee5ea4603@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.144231.0800.>

??Jason, I'd suggest you to follow these directions:
??1) Download and untar grads to a directory of your preference, for example /usr/local:> cd /usr/local> sudo tar -zxf grads-bin-x.xxxxxx.tar.gz (x "varies" according to GrADS version)
  2) Create a directory to hold "data.tar.Z":> cd /usr/local/grads-x.xxx> sudo mkdir data> cd data> sudo tar -zxf data.tar.Z
 3) Add the following lines to your "/etc/.bashrc" or "/etc/profile" if you want all users to access GrADS by typing "grads" on terminal:#--- GRADS ----export GADDIR=/usr/local/grads-x.xxx/dataexport GASCRP=/usr/local/grads-x.xxx/lib (this is optional, if you prefer to have a directory to run scripts from - create it just as described on Step 2)export GRADS=/usr/local/grads-x.xxx/binexport PATH=$PATH:$GRADS
  Adapt the commands to your system when necessary and it should work.
  Best regards,
  Thiago.
--- On Sun, 28/2/10, Jason Snyder <jmssnyder at UCDAVIS.EDU> wrote:

From: Jason Snyder <jmssnyder at UCDAVIS.EDU>
Subject: Installing GrADS on Ubuntu 9.10
To: GRADSUSR at LIST.CINECA.IT
Date: Sunday, 28 February, 2010, 18:16

Can you help give me information on how to install GrADS on Ubuntu?? I tried doing it on my own but am encountering numerous problems.? Can anyone who had installed GrADS on Ubuntu 9.10 give me detailed instructions on how to go about doing this including setting environments, etc.


Thanks,

Jason




      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/35b95edd/attachment.html 

From thi_veloso at YAHOO.COM.BR  Sun Feb 28 22:09:33 2010
From: thi_veloso at YAHOO.COM.BR (Thiago Veloso)
Date: Sun, 28 Feb 2010 19:09:33 -0800
Subject: Installing GrADS on Ubuntu 9.10
In-Reply-To: <fe8efa01002281316i4df0cc2cu73bde37ee5ea4603@mail.gmail.com>
Message-ID: <SUN.28.FEB.2010.190933.0800.>

??Jason, that's something I forgot to tell on the previous e-mail. Currently there's a more practical way to install grads on Ubuntu. The package is now available on the official channels (or will be, when Ubuntu 10.04 - Lucid- is released). To download the latest version visit the folowing website:?http://packages.ubuntu.com/lucid/grads.
??You will be asked to install one or two dependencies, which can be found at the same link...
??Best regards,
??Thiago.

--- On Sun, 28/2/10, Jason Snyder <jmssnyder at UCDAVIS.EDU> wrote:

From: Jason Snyder <jmssnyder at UCDAVIS.EDU>
Subject: Installing GrADS on Ubuntu 9.10
To: GRADSUSR at LIST.CINECA.IT
Date: Sunday, 28 February, 2010, 18:16

Can you help give me information on how to install GrADS on Ubuntu?? I tried doing it on my own but am encountering numerous problems.? Can anyone who had installed GrADS on Ubuntu 9.10 give me detailed instructions on how to go about doing this including setting environments, etc.


Thanks,

Jason




      
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100228/52e9ebcf/attachment.html