Thanks Arindam. It worked. Is it possible to do such time mean
for a 3D data. I wrote this script, but wrote only at one
time. Am I doing something wrong? Thanks in
advance. <br>
--------------------------------------------<br>
'reinit'<br>
'c'<br>
'set gxout fwrite'<br>
'set fwrite u_daily_3D.dat'<br>
'open u.ctl'<br>
'set lon 40 180'<br>
'set lat -10 10'<br>
t1=1<br>
z1=1<br>
while(z1<5)<br>
'set z 'z1<br>
while(t1<977)<br>
'set t 't1<br>
'd ave(u,t='t1',t='t1+7')'<br>
t1=t1+8<br>
endwhile<br>
z1=z1+1<br>
endwhile<br>
------------------------------------------------------------------------<br>
<br><br><div><span class="gmail_quote">On 2/22/06, <b class="gmail_sendername">Arindam Chakraborty</b> <<a href="mailto:arch@io.met.fsu.edu">arch@io.met.fsu.edu</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi,<br> Please change 'd ave(u,t=t1,t=t1+7)' to<br>'d ave(u,t='t1',t='t1+7')'<br><br>because t1 is a script variable and should be outside the single<br>quotes.<br><br>--<br><br>On Wed, 22 Feb 2006, John Guhin wrote:<br><br>
> Hi,<br>> I was wondering if someone can find the mistake. I am trying to<br>> convert a 3 hourly data to daily data.<br>> The original data is 3D. any help will be appreciated.<br>> ----------------------------------------------------------------------
<br>> * Program to write 3 hourly data to daily data<br>> 'reinit'<br>> 'c'<br>> 'set gxout fwrite'<br>> 'set fwrite u_daily.dat'<br>> 'open u.ctl'<br>> 'set lon 40 180'<br>> 'set lat -10 10'<br>> 'set lev 850'
<br>> t1=1<br>> while(t1<977)<br>> 'set t 't1<br>> 'd ave(u,t=t1,t=t1+7)'<br>> t1=t1+8<br>> endwhile<br>> --------------------------------------------------------------------------------------------------
<br>><br></blockquote></div><br>