[gradsusr] How to convert 6hourly data to monthly data in ECMWF data file

[MY]

Dear Mubashar,

I don't have a script. May be "tmave()" can help you, or, you can just
use ave() in while loops such as:

mon.1='Jan'
mon.2='Feb'
.....
mon.12='Dec'

mdy.1=31
mdy.2=28
.....
mdy.12=31

'set gxout fwrite'
'set fwrite test.dat'

iy=1958
while( iy<=2001)
   im=1
   while(im<=12)

        if( math_mod(iy,4)=0)
            mdy.2=29
        else
            mdy.2=28
        endif

        id=1
        while(id<=mdy.im)
            'd ave(slp, time=00z1' % mon.im % iy % ',time=18z' %
mdy.im % mon.im % iy % ')'
            id=id+1
        endwhile

        im=im+1
   endwhile
   iy=iy+1
endwhile

'disable fwrite'

Please try. Good luck!

MY

2010/7/12 Mubashar Dogar <mubashardogar at gmail.com>:
> Dear GrADs users,
>
> I have to compute seasonal means from the ECMWF data set. The data set
> available is in 6 hourly format. I want to convert it into monthly format so
> that I can compute the seasonal and annual means and then anomalies from
> that data set. Could you please guide/provide me with a script to convert
> hourly data to monthly format. Your suggestion will be helpful for me.
>
>
> Regards,
>
> Mubashar
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>



-- 
++++++++++++++++++++++++++++++++++++++++++++++++++++
YING, Ming
Shanghai Typhoon Institute of China Meteorological Administration
166 Puxi Road, Xuhui District, Shanghai, 200030
PR China
++++++++++++++++++++++++++++++++++++++++++++++++++++