[gradsusr] How to create a new Grads file?

Andrew Friedman andfried at berkeley.edu
Sat Nov 23 03:07:25 EST 2013


Hi Celso,

I’m not able to reproduce the problem using your script. What is pcp_2000M01_NEW.nc missing?
Also, I try to use a different name for the defined variable so as to avoid confusion.

Andrew

On Nov 22, 2013, at 7:32 AM, Celso Bandeira de Melo Ribeiro <celso.bandeira at ufjf.edu.br> wrote:

> Hello Sujata Mandke, good morning,
> Thank you for your attention.
> I put de x and y dimensions, how you suggest (see below), but de problem
> continues.
> 
> 'xdfopen amazon_2000N.ctl'
> 'set gxout shaded'
> 'open mask_new.ctl'
> 'set t 1'
> 'set lat -21.1 6.1'
> 'set lon -80.1 -48.9'
> 'rec1 = maskout(pcp,mask.2(t=1))'
> 'd rec1'
> 'define pcp = rec1'
> 'set sdfwrite pcp_2000M01_NEW.nc'
> 'sdfwrite pcp'
> 
> Can anybody help me?
> Thank you.
> Celso
> 
> 
> On Fri, 22 Nov 2013 12:40:06 +0530 (IST), Sujata Mandke
> <amin at tropmet.res.in> wrote:
>> Hello Celso,
>>   I am not sure if this is the solution of your problem,
>> but just try following - 
>> in your 'SAVE_pcp.gs' after opening
>> 'mask_new.ctl' you have 'set t 1', in addition
>> set x and y dimensions also, before defining variable 'rec1'.
>> regards
>> Sujata Mandke
>> 
>> ----- Original Message -----
>> From: Andrew Friedman <andfried at berkeley.edu>
>> To: GrADS Users Forum <gradsusr at gradsusr.org>
>> Sent: Fri, 22 Nov 2013 12:05:44 +0530 (IST)
>> Subject: Re: [gradsusr] How to create a new Grads file?
>> 
>> Hi Celso,
>> 
>> I’m not sure if this is all of it, but SAVE_pcp.gs is missing a single
>> quote at the end of the line 'define pcp = rec1’
>> Can you try it now?
>> 
>> Andrew
>> 
>> On Nov 21, 2013, at 3:55 PM, Celso Bandeira de Melo Ribeiro
>> <celso.bandeira at ufjf.edu.br> wrote:
>> 
>>> Hi Andrews and Grads users,
>>> 
>>> Can anybody help me?
>>> I am trying to apply masking in NetCdf file (ORIGINAL file =
>>> pcp_2000M01.nc) and trying to save de results in the same format
> "NetCdf"
>>> (result file = pcp_2000M01NEW.nc).
>>> 
>>> The script (SAVE_pcp.gs) is:
>>> 'xdfopen amazon_2000N.ctl'
>>> 'set gxout shaded'
>>> 'open mask_new.ctl'
>>> 'set t 1'
>>> 'rec1 = maskout(pcp,mask.2(t=1))'
>>> 'd rec1'
>>> 'define pcp = rec1
>>> 'set sdfwrite pcp_2000M01_NEW.nc'
>>> 'sdfwrite pcp'
>>> 
>>> The results appear that it has been saved file but when I try to open
> it
>>> (using amazon_2000NEW.ctl) only show the limits (boundary), but not the
>>> data inside limits.
>>> I am attaching the files that I am working.
>>> I appreciate some help.
>>> Celso
>>> 
>>> 
>>> On Fri, 15 Nov 2013 10:10:28 -0800, Andrew Friedman
>>> <andfried at berkeley.edu>
>>> wrote:
>>>> Hi Celso,
>>>> 
>>>> For save.gs, you also have to download the dependency scripts
>>>> parsestr.gsf, parseopt.gsf, qdims.gsf and put them in your GrADS
> script
>>>> folder.
>>>> The program should work for gridded data, and create a .ctl control
> file
>>>> and .dat data file. You can open them by using the command ‘open’ on
>>>> the
>>>> control file.
>>>> 
>>>> One issue I’ve had with save.gs recently is that the line OPTIONS _cal
>>>> causes problems in the control file, so try commenting it out.
>>>> 
>>>> Alternatively, you could use the GrADS command sdfwrite to create a
>>> netcdf
>>>> file.
>>>> 
>>>> If your data aren’t gridded, you may want to try writing as ASCII
>>>> following the instructions in the User’s Guide.
>>>> 
>>>> Andrew
>>>> 
>>>> 
>>>> On Nov 14, 2013, at 1:24 PM, Celso Bandeira de Melo Ribeiro
>>>> <celso.bandeira at ufjf.edu.br> wrote:
>>>> 
>>>>> Hello Andrew,
>>>>> Thank you for your help.
>>>>> I opened the script but I didn't understand how can I use it in may
>>> case.
>>>>> Could you help me more?
>>>>> Thank you again.
>>>>> Best wishes.
>>>>> Celso
>>>>> 
>>>>> 
>>>>> On Tue, 12 Nov 2013 17:01:07 -0800, Andrew Friedman
>>>>> <andfried at berkeley.edu>
>>>>> wrote:
>>>>>> Hi Celso,
>>>>>> 
>>>>>> Check out save.gs by Bin Guan.
>>>>>> 
>>>>>> http://www.atmos.umd.edu/~bguan/grads/GrADS_Scripts.html#save.gs
>>>>>> 
>>>>>> Andrew
>>>>>> 
>>>>>> On Nov 12, 2013, at 4:01 PM, Celso Bandeira de Melo Ribeiro
>>>>>> <celso.bandeira at ufjf.edu.br> wrote:
>>>>>> 
>>>>>>> Dear Grads user,
>>>>>>> I need a help.
>>>>>>> I'm trying to generate monthly rainfall files by accumulating daily
>>>>>>> rainfall from TRMM-MERGE.
>>>>>>> I created an scripts in this Grads and it is adding precipitation
>>> well
>>>>>>> for
>>>>>>> each month, but don't know how can I save this result, the monthly
>>> sum
>>>>>>> in a
>>>>>>> new Grads file.
>>>>>>> Could someone help me to save the monthly rainfall data in a new
>>> Grads
>>>>>>> file?
>>>>>>> Here the code I'm using in Grads.
>>>>>>> 
>>>>>>> 'open pcp.ctl'
>>>>>>> 'jan = sum(prec,t=1,t=31)'
>>>>>>> 'fev = sum(prec,t=32,t=59)'
>>>>>>> 'mar = sum(prec,t=60,t=90)'
>>>>>>> 'abr = sum(prec,t=91,t=120)'
>>>>>>> 'mai = sum(prec,t=121,t=151)'
>>>>>>> 'jun = sum(prec,t=152,t=181)'
>>>>>>> 'jul = sum(prec,t=182,t=212)'
>>>>>>> 'ago = sum(prec,t=213,t=243)'
>>>>>>> 'sep = sum(prec,t=244,t=273)'
>>>>>>> 'oct = sum(prec,t=274,t=304)'
>>>>>>> 'nov = sum(prec,t=305,t=334)'
>>>>>>> 'dec = sum(prec,t=335,t=365)'
>>>>>>> 'd jan'
>>>>>>> 
>>>>>>> I would like to save each monthly file: jan, fev, mar,...,dec.
>>>>>>> 
>>>>>>> Best,
>>>>>>> Celso
>>>>>>> 
>>>>>>> 
>>>>>>> -- 
>>>>>>> ____________________________________________
>>>>>>> Prof. Celso Bandeira de Melo Ribeiro
>>>>>>> Universidade Federal de Juiz de Fora - UFJF
>>>>>>> Faculdade de Engenharia, Plataforma 4
>>>>>>> Departamento de Engenharia Ambiental e Sanitaria
>>>>>>> Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900
>>>>>>> tel. (32) 2102-3400
>>>>>>> Visiting Scholar in Texas A&M University
>>>>>>> _______________________________________________
>>>>>>> gradsusr mailing list
>>>>>>> gradsusr at gradsusr.org
>>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>>> 
>>>>>> 
>>>>>> _______________________________________________
>>>>>> gradsusr mailing list
>>>>>> gradsusr at gradsusr.org
>>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>>> 
>>>>> -- 
>>>>> ____________________________________________
>>>>> Prof. Celso Bandeira de Melo Ribeiro
>>>>> Universidade Federal de Juiz de Fora - UFJF
>>>>> Faculdade de Engenharia, Plataforma 4
>>>>> Departamento de Engenharia Ambiental e Sanitaria
>>>>> Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900
>>>>> tel. (32) 2102-3400
>>>>> Visiting Scholar in Texas A&M University
>>>>> _______________________________________________
>>>>> gradsusr mailing list
>>>>> gradsusr at gradsusr.org
>>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>>> 
>>>> 
>>>> _______________________________________________
>>>> gradsusr mailing list
>>>> gradsusr at gradsusr.org
>>>> http://gradsusr.org/mailman/listinfo/gradsusr
>>> 
>>> -- 
>>> ____________________________________________
>>> Prof. Celso Bandeira de Melo Ribeiro
>>> Universidade Federal de Juiz de Fora - UFJF
>>> Faculdade de Engenharia, Plataforma 4
>>> Departamento de Engenharia Ambiental e Sanitaria
>>> Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900
>>> tel. (32) 2102-3400
>>> Visiting Scholar in Texas A&M
>>> 
> University<amazon_2000N.ctl><amazon_2000NEW.ctl><mask_new.ctl><mask_new.dat><pcp_2000M01.nc><pcp_2000M01_NEW.nc><SAVE_pcp.gs>_______________________________________________
>>> gradsusr mailing list
>>> gradsusr at gradsusr.org
>>> http://gradsusr.org/mailman/listinfo/gradsusr
>> 
>> 
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>> 
>> 
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
> 
> -- 
> ____________________________________________
> Prof. Celso Bandeira de Melo Ribeiro
> Universidade Federal de Juiz de Fora - UFJF
> Faculdade de Engenharia, Plataforma 4
> Departamento de Engenharia Ambiental e Sanitaria
> Bairro Martelos, Juiz de Fora (MG) CEP. 36036-900
> tel. (32) 2102-3400
> Visiting Scholar in Texas A&M University
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr




More information about the gradsusr mailing list