On Jan 30, 2008 11:50 AM, Ranyére Nóbrega <<a href="mailto:ranyere.nobrega@sipam.gov.br">ranyere.nobrega@sipam.gov.br</a>> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Dear All,<br><br>How can calculate the precipitation accumulated diary in 3B42 TRMM<br>data, once do they belong to 3 and 3 hours?<br><br>Does anybody have a script for that? I need for the whole data period.<br></blockquote>
<div><br>Here is an example using command line lats4d for one month. The input ctl is attached below.<br><br>lats4d -i 3b42_TRMM200407.ctl -o test -format grads_grib -time 00Z01JUL2004 00Z31JUL2004 8 -func "sum(@,t+0,t+7)" -vars p<br>
<br>I would suggest wrapping it in a csh loop for each month of the period. However, it would also be possible to put your begin and end time into the -time option, and get all of the data in one file. That is not practical, if you plan on updating your new daily data in the future. <br>
<br>The result should be 31 times in the output file, a sum over each day.<br><br>Also, please check the 3b42 documentation and do a unit analysis, so you know what you get from this. This assume the input is mm/3hr, and the accumulation is mm in a day. you can add scale factors inside the -func option to suit your need.<br>
<br> Good Luck! (ctl follows)<br><br>Mike<br><br>dset 3B42/data/3B42.01%m2%d2.%h1.6.HDF<br>dtype hdfsds<br>options template<br>undef -9999.9<br>xdef 1440 linear 180 0.25<br>ydef 400 linear -50 0.25<br>zdef 1 levels 1<br>tdef 248 linear 00Z01jul2001 3hr<br>
vars 2<br>precipitation=>p 0 t,x,y precipitation<br>relativeError=>re 0 t,x,y relative error<br>endvars<br><br></div><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>Thanks.<br><font color="#888888"><br>--<br>__________________________________________________________<br>Ranyére Silva Nóbrega<br>Analista Intelectual<br>Divisão de Meteorologia e Climatologia - SIPAM<br>Centro Técnico e Operacional de Porto Velho<br>
</font></blockquote></div><br>