[gradsusr] lats4d for sum

Mon Jan 21 10:33:00 EST 2013

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20130121/0620e3b2/attachment-0003.html 