[gradsusr] lats4d for sum

[Lee Byerle]

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