[gradsusr] lats4d for sum
Estefanía Jiménez
estefaniajimenez.r at gmail.com
Mon Jan 21 15:25:10 EST 2013
Never mind! Made it work!! finally ;) Here is the final version:
*#!/usr/bin/csh*
*
*
*set mese = (jan feb mar apr may jun jul aug sep oct nov dec)*
*
*
*@ ano = 1961*
*while($ano <= 1990)*
*
*
* @ mes = 1*
* while($mes <= 12)*
* if($mes < 10) then*
* set mes = '0'$mes*
* endif*
*
*
* @ dia = 1*
* while($dia <= 30)*
* if($dia < 10) then*
* set dia = '0'$dia*
* endif*
*
*
*
*
* /opt/grads30d/bin/gradsc -bpc "run /opt/grads30d/scripts/lats4d.gs -i
/data/CR_presente_1960-1990/binctl/Eta_exp110km1960010100.ctl -o
../climat_6190sumaprec_diaria/sumaprec_$dia$mes$ano -ftype ctl -format
grads_grib -vars prec -time 00z$dia$mese[$mes]$ano 18z$dia$mese[$mes]$ano
-func sum(@*1000,time=00z$dia$mese[$mes]$ano,time=18z$dia$mese[$mes]$ano)
-levs 1000 925 850 500 250 -table cptec.table"*
*
*
*@ dia = $dia + 1*
*end*
*
*
*@ mes = $mes + 1*
*end*
*@ ano = $ano + 1*
*end*
*
*
If you still find an error please let me know.
Thank you!
2013/1/21 Estefanía Jiménez <estefaniajimenez.r at gmail.com>
> thank you Lee for all your help and patience!
>
> I got a new problem: I only get the script to do the first year... what am
> I missing??
>
> And I had to put
> *sum(@*1000,time=00z$dia$mese[$mes]$ano,time=18z$dia$mese[$mes]$ano) *
> instead of *'**
> sum(@*1000,time=00z$dia$mese[$mes]$ano,time=18z$dia$mese[$mes]$ano)'*
> 'cause if not I get another error
>
> I attach the newest version:
>
>
>
> *#!/usr/bin/csh*
> *
> *
> *set mese = (jan feb mar apr may jun jul aug sep oct nov dec)*
> *
> *
> *@ ano = 1961*
> *while($ano <= 1990)*
> *
> *
> * @ mes = 1*
> * while($mes <= 12)*
> * if($mes < 10) then*
> * set mes = '0'$mes*
> * endif*
> *
> *
> * @ dia = 1*
> * while($dia <= 30)*
> * if($dia < 10) then*
> * set dia = '0'$dia*
> * endif*
> *
> *
> * echo 00z$dia$mese[$mes]$ano 18z$dia$mese[$mes]$ano*
> *
> *
> * /opt/grads30d/bin/gradsc -bpc "run /opt/grads30d/scripts/lats4d.gs -i
> /data/CR_presente_1960-1990/binctl/Eta_exp110km1960010100.ctl -o
> ../climat_6190sumaprec_diaria/sumaprec_$dia$mes$ano -ftype ctl -format
> grads_grib -vars prec -time 00z$dia$mese[$mes]$ano 18z$dia$mese[$mes]$ano
> -func sum(@*1000,time=00z$dia$mese[$mes]$ano,time=18z$dia$mese[$mes]$ano)
> -levs 1000 925 850 500 250 -table cptec.table"*
> *
> *
> *@ dia = $dia + 1*
> *end*
> *
> *
> *@ mes = $mes + 1*
> *end*
>
> Regards and thanks again!
>
>
> 2013/1/18 Lee Byerle <lbyerle at yahoo.com>
>
>> Estefania,
>> Here is a sample lats4d command w/sum() that's working w/NOAA-CIRES 20th
>> Cent reanalysis. I added back in -time (selects single time with t1 and t2
>> the same). I moved -func {expression} to sum 00-18Z prate after -vars
>> expression, instead of before it.
>>
>> lats4d.sh -i
>> http://www.esrl.noaa.gov/psd/thredds/dodsC/Datasets/20thC_ReanV2/gaussian/monolevel/prate.1983.nc-vars prate -time 00z18jul1983 00z18jul1983 -func
>> 'sum(@,time=00z18jul1983,time=18z18jul1983)' -format grads_grib -center
>> ncep -o try.grb -v
>>
>> Lee
>>
>>
>> On Jan 18, 2013, at 11:56 AM, Estefanía Jiménez wrote:
>>
>> Tried all that... didn't work.
>>
>>
>>
>>
>> 2013/1/17 Lee Byerle <lbyerle at yahoo.com>
>>
>>> Hi Estefania,
>>> Instead of t=, try time=
>>>
>>> You may consider testing with a specific set of times without the
>>> strings in the sum expression if the error persists.
>>>
>>> I hope this works.
>>> Lee
>>>
>>>
>>> On Jan 17, 2013, at 1:30 PM, Estefanía Jiménez <
>>> estefaniajimenez.r at gmail.com> wrote:
>>>
>>> I got this error:
>>>
>>> GrADS-LATS Error: (lats_write) failed...
>>> Data Request Warning: Request beyond file limits
>>> LATS error (lats_write): LATS (lats.c) --> For GrADS convention, cannot
>>> skip timepoints; current time = 271716.000000 hours since 1960-1-1 0,
>>> previous time = 714.000000 hours since 1960-1-1 0, delta = 6, file =
>>> ../climat_6190sumaprec_diaria/sumaprec_12121961.grb, variable = prec
>>>
>>> any ideas? I attached the full script:
>>>
>>> #!/usr/bin/csh
>>>
>>> set mese = (jan feb mar apr may jun jul aug sep oct nov dec)
>>>
>>> @ ano = 1961
>>>
>>> while($ano <= 1990)
>>>
>>> @ mes = 1
>>> while($mes <= 12)
>>> if($mes < 10) then
>>> set mes = '0'$mes
>>> endif
>>>
>>> @ dia = 1
>>> while($dia <= 30)
>>> if($dia < 10) then
>>> set dia = '0'$dia
>>> endif
>>>
>>>
>>> /opt/grads30d/bin/gradsc -bpc "run /opt/grads30d/scripts/lats4d.gs-i /data/CR_presente_1960-1990/binctl/Eta_exp110km1960010100.ctl -o
>>> ../climat_6190sumaprec_diaria/sumaprec_$dia$mes$ano -ftype ctl -format
>>> grads_grib -func sum(@,t=00z$dia$mese[$mes]$ano,t=18z$dia$mese[$mes]$ano)
>>> -vars prec -levs 1000 925 850 500 250 -table cptec.table"
>>>
>>> @ dia = $dia + 1
>>> @ mes = $mes + 1
>>> end
>>> * *
>>> 25,5 Bot
>>>
>>>
>>>
>>> 2013/1/17 Estefanía Jiménez <estefaniajimenez.r at gmail.com>
>>>
>>>> yuju! thanks a lot Lee!
>>>>
>>>>
>>>> 2013/1/16 Lee Byerle <lbyerle at yahoo.com>
>>>>
>>>>> Estefania,
>>>>> One suggestion is to remove the options -time and -mean and add -func
>>>>> 'sum(@,time=00z01...,time=18z30...)' to the Lats command (with your month,
>>>>> year and day strings included as before).
>>>>>
>>>>> Good luck,
>>>>> Lee
>>>>>
>>>>> On Jan 16, 2013, at 12:08 PM, Estefanía Jiménez <
>>>>> estefaniajimenez.r at gmail.com> wrote:
>>>>>
>>>>> Hi all! Is there a way that I can transform this script to get the sum
>>>>> of the variables instead of the average??
>>>>>
>>>>>
>>>>> * /opt/grads30d/bin/gradsc -bpc "run /opt/grads30d/scripts/lats4d.gs-i /data/CR_presente_1960-1990/binctl/Eta_exp110km1960010100.ctl -o
>>>>> ../climat_6190/media_$mes$ano -ftype ctl -format grads_grib -mean -vars
>>>>> pslm pslc tp2m dp2m u10m v10m prec prcv neve clsf cssf ghfl tsfc qsfc tgsc
>>>>> tgrz ussl uzrs smav rnof evpp lwnv mdnv hinv cbnt ocis olis oces oles roce
>>>>> role albe agpl zgeo uvel vvel temp umrl omeg umes -time
>>>>> 00z01$mese[$mes]$ano 18z30$mese[$mes]$ano -levs 1000 925 850 500 250 -table
>>>>> cptec.table"*
>>>>> *
>>>>> *
>>>>>
>>>>> *
>>>>> *
>>>>>
>>>>> Thank you for your help!
>>>>>
>>>>> --
>>>>> *Estefanía Jiménez R.*
>>>>> Meteoróloga
>>>>> Departamento de Climatología e Investigaciones Aplicadas
>>>>> Instituto Meteorológico Nacional
>>>>> MINAET
>>>>> Tel: +506 2222 5616 ext 131
>>>>>
>>>>> _______________________________________________
>>>>> 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
>>>>>
>>>>>
>>>>
>>>>
>>>> --
>>>> *Estefanía Jiménez R.*
>>>> Meteoróloga
>>>> Departamento de Climatología e Investigaciones Aplicadas
>>>> Instituto Meteorológico Nacional
>>>> MINAET
>>>> Tel: +506 2222 5616 ext 131
>>>>
>>>
>>>
>>>
>>> --
>>> *Estefanía Jiménez R.*
>>> Meteoróloga
>>> Departamento de Climatología e Investigaciones Aplicadas
>>> Instituto Meteorológico Nacional
>>> MINAET
>>> Tel: +506 2222 5616 ext 131
>>>
>>> _______________________________________________
>>> 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
>>>
>>>
>>
>>
>> --
>> *Estefanía Jiménez R.*
>> Meteoróloga
>> Departamento de Climatología e Investigaciones Aplicadas
>> Instituto Meteorológico Nacional
>> MINAET
>> Tel: +506 2222 5616 ext 131
>> _______________________________________________
>> 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
>>
>>
>
>
> --
> *Estefanía Jiménez R.*
> Meteoróloga
> Departamento de Climatología e Investigaciones Aplicadas
> Instituto Meteorológico Nacional
> MINAET
> Tel: +506 2222 5616 ext 131
>
--
*Estefanía Jiménez R.*
Meteoróloga
Departamento de Climatología e Investigaciones Aplicadas
Instituto Meteorológico Nacional
MINAET
Tel: +506 2222 5616 ext 131
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130121/00d4087e/attachment-0003.html
More information about the gradsusr
mailing list