seasonal average
Charles Seman
Charles.Seman at NOAA.GOV
Fri Nov 10 17:59:12 EST 2006
Dear MC,
One idea would be to write out the results to a GrADS fwrite file
'set gxout fwrite' ; 'set fwrite tseries_tave.ieee'
'set x '...
'set y '...
'set z '...
ts=1
while ( ts <= 46 )
'set t 'ts*12
'define a = (air(t-1)+air(t+0)+air(t+1))/3' * unweighted time average
say 'ts 'ts
'd a'
'undefine a'
ts=ts+1
endwhile
'disable fwrite'
I hope this helps,
Chuck
praise praise wrote:
> Dear Chuck,
>
> Yes, I need a time series with 46 time levels, each time level
> corresponding to a DJF seasonal average for one year (46 years in the
> time
> series).
>
> Thank you very much.
> MC
>
>
>> From: Charles Seman <Charles.Seman at NOAA.GOV>
>> Reply-To: GRADSUSR at LIST.CINECA.IT
>> To: GRADSUSR at LIST.CINECA.IT
>> Subject: Re: seasonal average
>> Date: Thu, 9 Nov 2006 19:09:35 -0500
>>
>> MC,
>>
>> Do you need a time series with 46 time levels, each time level
>> corresponding to a DJF seasonal average for one year (46 years in the
>> time series)?
>>
>> Thanks,
>> Chuck
>>
>> praise praise wrote:
>> > Thank you Charles,
>> >
>> > But I need a variable having a time dimension of 46 years, not the
>> > seasonal
>> > climatology.
>> > Thanks in advance.
>> >
>> > mc
>> >
>> >
>> >> From: Charles Seman <Charles.Seman at NOAA.GOV>
>> >> Reply-To: GRADSUSR at LIST.CINECA.IT
>> >> To: GRADSUSR at LIST.CINECA.IT
>> >> Subject: Re: seasonal average
>> >> Date: Wed, 8 Nov 2006 10:16:57 -0500
>> >>
>> >> Dear praise praise,
>> >>
>> >> Please find attached files containing email posts on this subject to
> the
>> >> Listserv... please see also
>> >> http://grads.iges.org/grads/gadoc/gradcomdmodify.html
>> >>
>> >> I hope this helps,
>> >> Chuck
>> >>
>> >> praise praise wrote:
>> >> > Dear All,
>> >> >
>> >> > I have monthly data say starts from 1958 Jan till 2005 Dec (ie. t=1
> to
>> >> > 576).
>> >> >
>> >> > I would like to do an EOF analysis for the 46 winters only (say the
>> >> > average
>> >> > of DJF from 1958/59 to 2004/05).
>> >> >
>> >> > So, how can I define a new variable representing the 46 winters
>> (Not
>> > the
>> >> > climatology) before I can run the EOF.
>> >> >
>> >> > I have tried this:
>> >> > set t 1 46
>> >> > define a=ave(air,t+11,t+13,12)
>> >> >
>> >> > But clearly, its result was not I want.
>> >> >
>> >> > thanks
>> >> >
>> >> > mc
>> >> >
>> >> > _________________________________________________________________
>> >> > MSN Tool Bar 幫你刪除惱人的廣告 ! http://toolbar.msn.com.hk
>> >>
>> >> --
>> >>
>> >> Please note that Charles.Seman at noaa.gov should be considered my NOAA
>> >> email address, not cjs at gfdl.noaa.gov.
>> >>
>> >> ********************************************************************
>> >> Charles Seman Charles.Seman at noaa.gov
>> >> U.S. Department of Commerce / NOAA / OAR
>> >> Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
>> >> 201 Forrestal Road fax: (609) 987-5063
>> >> Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
>> >> ********************************************************************
>> >>
>> >> "The contents of this message are mine personally and do not reflect
>> >> any position of the Government or NOAA."
>> >>
>> >
>> >
>> >> From owner-gradsusr at LIST.CINECA.IT Tue Jan 16 13:23:32 2001
>> >> Received: from marge.cineca.it (marge1.cineca.it [130.186.1.9])
>> >> by GFDL.NOAA.GOV (8.9.3/8.9.3) with ESMTP id NAA04852;
>> >> Tue, 16 Jan 2001 13:23:28 -0500 (EST)
>> >> Received: from arianna (arianna.cineca.it [130.186.1.53])
>> >> by marge.cineca.it (8.11.1/8.11.1/CINECA 4.0) with ESMTP id
>> >> f0GI6Dg00074;
>> >> Tue, 16 Jan 2001 19:06:13 +0100 (MET)
>> >> Received: from LIST.CINECA.IT by LIST.CINECA.IT (LISTSERV-TCP/IP
> release
>> > 1.8d)
>> >> with spool id 1156776 for GRADSUSR at LIST.CINECA.IT; Tue, 16 Jan
>> > 2001
>> >> 19:06:12 +0100
>> >> Received: from cola.iges.org (cola.iges.org [192.239.84.10]) by
>> >> arianna.cineca.it (8.11.0.Beta1/8.11.0.Beta1/CINECA 4.0) with
>> > ESMTP
>> >> id f0GI6Bx21116 for <GRADSUSR at LIST.CINECA.IT>; Tue, 16 Jan 2001
>> >> 19:06:11 +0100 (MET)
>> >> Received: by cola.iges.org (8.8.8/1.1.22.3/16Jan00-0454PM) id
>> > NAA0000013641;
>> >> Tue, 16 Jan 2001 13:06:11 -0500 (EST)
>> >> Message-ID: <200101161806.NAA0000013641 at cola.iges.org>
>> >> Date: Tue, 16 Jan 2001 13:06:11 -0500
>> >> Reply-To: GRADSUSR at LIST.CINECA.IT
>> >> Sender: GRADSUSR at LIST.CINECA.IT
>> >> From: "Brian E. Doty" <doty at COLA.IGES.ORG>
>> >> Subject: Re: plotting question
>> >> To: GRADSUSR at LIST.CINECA.IT
>> >> Status: RO
>> >>
>> >> Seasonal averages can be done in two steps:
>> >>
>> >> First creat a 12 month seasonal average, and save it
>> >> as a defined variable:
>> >>
>> >> set lon 0 360
>> >> set lat -90 90
>> >> set z 1
>> >> set t 1 12
>> >> define ozave = ave(ozone,t+0,t=144,12)
>> >>
>> >> Now 'ozave' contains the seasonal mean for each month (the above
> assumes
>> >> the data starts with January).
>> >>
>> >> Then, to display a seasonal mean -- for example, DJF:
>> >>
>> >> display (ozave(t=12)+ozave(t=1)+ozave(t=2))/3
>> >>
>> >> or if you want to account for the lengths of the months (note that
>> ave
>> >> does not do this by default when averaging monthly data):
>> >>
>> >> display (ozave(t=12)*31+ozave(t=1)*31+ozave(t=2)*28)/90
>> >>
>> >> ...Brian
>> >
>> >
>> >> From arch at CAOS.IISC.ERNET.IN Fri Mar 19 13:12:29 2004
>> >> Date: Fri, 19 Mar 2004 15:18:22 +0530
>> >> From: Arindam Chakraborty <arch at CAOS.IISC.ERNET.IN>
>> >> Reply-To: GRADSUSR at LIST.CINECA.IT
>> >> To: GRADSUSR at LIST.CINECA.IT
>> >> Subject: Re: need confirmation on monthly average
>> >>
>> >> hi,
>> >> I think Sulchan requires seasonal average, that is average
>> >> over all jan's of the available years and so on for the other
>> >> months. For that do like this:
>> >>
>> >> 'sdfopen your file'
>> >> 'set the region in latlon'
>> >> 'set t 1 12'
>> >> 'define mnsn = ave(var,t+0,t=132,12)'
>> >> 'define mnsn seasonal'
>> >> 'set t 1'
>> >> * Now the variable mnsn is seasonal average of all the months over
>> the
>> >> * years available. It has no year associated with it. And it has 12
>> >> * values. t=1 corresponts to Jan and so on.
>> >> * Notice that in the period jan1991 - dec2001 there are 11 years,
>> >> * and so there are 132 months
>> >>
>> >> * Now to display, say feb do
>> >> 'd mnsn(t=2)'
>> >>
>> >>
>> >> Hope this helps.
>> >>
>> >> ~ Arindam Chakraborty
>> >>
>> >> On Fri, 19 Mar 2004, Don Hooper wrote:
>> >>
>> >> > > From owner-gradsusr at LIST.CINECA.IT Fri Mar 19 01:22:39 2004
>> >> > > Reply-To: M Sulchan Darmawan <msdarmawan at softhome.net>
>> >> > >
>> >> > > You mean I've done something wrong with the script ?
>> >> > Not as far as I know, but I'm no expert on GrADS scripting.
>> >> > > Please show me the right way ?
>> >> > See above. I'm good under the hood, but not so hot behind the
>> wheel.
>> >> > > I didn't see which discussion is not relevant, since I use grads
>> > to
>> >> > > open it.
>> >> > My point is, whether the data file is NetCDF, or HDF-SDS, or GRIB,
> or
>> >> > flat binary, the script wouldn't be any different after the first
>> >> line.
>> >> > The first line might change (to xdfopen or open maybe), but that's
> it.
>> >> > To be concerned about the fact that your data is in NetCDF is a red
>> > herring
>> >> > as far as this problem is concerned. It just doesn't matter.
>> >> > >
>> >> > > Best regards,
>> >> > > Sulchan
>> >> > >
>> >> > > Don Hooper, on Friday, March 19, 2004, 9:54:32 AM, you wrote:
>> >> > > DH> I suspect I'll get some criticism for this post. One thing I
>> > know for
>> >> > > DH> sure: the fact the file is NetCDF isn't relevant to this
>> > discussion.
>> >> > > DH> Once the file is open, the input format type is irrelevant.
>> >> > >
>> >> > > DH> -Hoop
>> >> > > >> From owner-gradsusr at LIST.CINECA.IT Thu Mar 18 19:40:57 2004
>> >> > > >> Reply-To: M Sulchan Darmawan <msdarmawan at softhome.net>
>> >> > > >>
>> >> > > >> Dear GRADSUSR...
>> >> > > >>
>> >> > > >> I have monthly data from NCEP realtime marine data. the
>> > format is
>> >> > > >> NetCDF. I'm trying to do a monthly average over 10 years
>> > data. I
>> >> > > >> need your confirmation to avoid mistaken in my analysis,
>> > does my
>> >> > > >> script is right or not ?, if its wrong, please show me the
>> > right
>> >> > > >> way...
>> >> > > >>
>> >> > > >> Data description :
>> >> > > >> monthly data : jan1991 - dec2001
>> >> > > >>
>> >> > > >> purpose :
>> >> > > >> Monthly average : Jan, Feb, ..., Dec = 12 picture
>> >> > > >>
>> >> > > >> Script :
>> >> > > >> sdfopen sst.mean.nc
>> >> > > >> set lon 80 110
>> >> > > >> set lat -10 10
>> >> > > >> set t 1
>> >> > > >> define zave=ave(sst, t+0, t=120, 1yr)
>> >> > > >> d zave
>> >> > > >>
>> >> > > >> does this script means January average ?
>> >> > > >> what is the different if I set t to 13 (set t 13) ?
>> >> > > >>
>> >> > > >> Just to make sure that I make the analysis right ;) Thank you
>> > very
>> >> > > >> much.
>> >> > > >>
>> >> > > >> Cheers,
>> >> > > >> M Sulchan Darmawan
>> >> > > >>
>> >> > > >> --
>> >> > > >> GEOSTECH -- Sub lab Teknologi Sistem Iklim
>> >> > > >> Laboratorium Teknologi Sistem Kebumian dan Mitigasi Bencana
> Alam
>> >> > > >> PUSPIPTEK -- http://www.clivarindo.org
>> >> >
>> >> > -Hoop
>> >> >
>> >>
>> >> From msdarmawan at softhome.net Thu Mar 25 13:10:57 2004
>> >> Date: Mon, 22 Mar 2004 20:11:25 +0700
>> >> From: M Sulchan Darmawan <msdarmawan at softhome.net>
>> >> To: GRADSUSR at LIST.CINECA.IT
>> >> Subject: Re: need confirmation on monthly average
>> >>
>> >> Many thanks to you. Your explanation is valuable for me as a beginner
>> >> in GrADS. Great explanation, it help very much.
>> >>
>> >> Warm regards from Indonesia,
>> >> Sulchan
>> >>
>> >> Arindam Chakraborty, on Friday, March 19, 2004, 4:48:22 PM, you
>> wrote:
>> >> AC> hi,
>> >> AC> I think Sulchan requires seasonal average, that is average
>> >> AC> over all jan's of the available years and so on for the other
>> >> AC> months. For that do like this:
>> >>
>> >> AC> 'sdfopen your file'
>> >> AC> 'set the region in latlon'
>> >> AC> 'set t 1 12'
>> >> AC> 'define mnsn = ave(var,t+0,t=132,12)'
>> >> AC> 'define mnsn seasonal'
>> >> AC> 'set t 1'
>> >> AC> * Now the variable mnsn is seasonal average of all the months
>> >> over the
>> >> AC> * years available. It has no year associated with it. And it has
> 12
>> >> AC> * values. t=1 corresponts to Jan and so on.
>> >> AC> * Notice that in the period jan1991 - dec2001 there are 11 years,
>> >> AC> * and so there are 132 months
>> >>
>> >> AC> * Now to display, say feb do
>> >> AC> 'd mnsn(t=2)'
>> >>
>> >>
>> >> AC> Hope this helps.
>> >>
>> >> AC> ~ Arindam Chakraborty
>> >>
>> >> <deleted>
>> >>
>> >> From arch at CAOS.IISC.ERNET.IN Thu Mar 25 13:11:22 2004
>> >> Date: Mon, 22 Mar 2004 19:59:04 +0530
>> >> From: Arindam Chakraborty <arch at CAOS.IISC.ERNET.IN>
>> >> Reply-To: GRADSUSR at LIST.CINECA.IT
>> >> To: GRADSUSR at LIST.CINECA.IT
>> >> Subject: Re: need confirmation on monthly average
>> >>
>> >> hi Sulchan,
>> >> I made a mistake perviously. Instead of
>> >> 'define mnsn seasonal'
>> >>
>> >> it has to be
>> >>
>> >> 'modify mnsn seasonal'
>> >>
>> >> Sorry for that. Everything else is ok.
>> >>
>> >> regards,
>> >>
>> >> ~ Arindam Chakraborty
>> >>
>> >> On Mon, 22 Mar 2004, M Sulchan Darmawan wrote:
>> >>
>> >> > Many thanks to you. Your explanation is valuable for me as a
> beginner
>> >> > in GrADS. Great explanation, it help very much.
>> >> >
>> >> > Warm regards from Indonesia,
>> >> > Sulchan
>> >> >
>> >> > Arindam Chakraborty, on Friday, March 19, 2004, 4:48:22 PM, you
> wrote:
>> >> > AC> hi,
>> >> > AC> I think Sulchan requires seasonal average, that is average
>> >> > AC> over all jan's of the available years and so on for the other
>> >> > AC> months. For that do like this:
>> >> >
>> >> > AC> 'sdfopen your file'
>> >> > AC> 'set the region in latlon'
>> >> > AC> 'set t 1 12'
>> >> > AC> 'define mnsn = ave(var,t+0,t=132,12)'
>> >> > AC> 'define mnsn seasonal'
>> >> > AC> 'set t 1'
>> >> > AC> * Now the variable mnsn is seasonal average of all the months
> over
>> > the
>> >> > AC> * years available. It has no year associated with it. And it
>> >> has 12
>> >> > AC> * values. t=1 corresponts to Jan and so on.
>> >> > AC> * Notice that in the period jan1991 - dec2001 there are 11
> years,
>> >> > AC> * and so there are 132 months
>> >> >
>> >> > AC> * Now to display, say feb do
>> >> > AC> 'd mnsn(t=2)'
>> >> >
>> >> >
>> >> > AC> Hope this helps.
>> >> >
>> >> > AC> ~ Arindam Chakraborty
>> >> >
>> >> > <deleted>
>> >> >
>> >
>> > _________________________________________________________________
>> > MSN Tool Bar 幫你刪除惱人的廣告 ! http://toolbar.msn.com.hk
>>
>> --
>>
>> Please note that Charles.Seman at noaa.gov should be considered my NOAA
>> email address, not cjs at gfdl.noaa.gov.
>>
>> ********************************************************************
>> Charles Seman Charles.Seman at noaa.gov
>> U.S. Department of Commerce / NOAA / OAR
>> Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
>> 201 Forrestal Road fax: (609) 987-5063
>> Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
>> ********************************************************************
>>
>> "The contents of this message are mine personally and do not reflect
>> any position of the Government or NOAA."
>
> _________________________________________________________________
> Hotmail Extra Storage讓你獲得10MB 額外儲存空間，請即申請！
> http://join.msn.com/?pgmarket=zh-hk
--
Please note that Charles.Seman at noaa.gov should be considered my NOAA
email address, not cjs at gfdl.noaa.gov.
********************************************************************
Charles Seman Charles.Seman at noaa.gov
U.S. Department of Commerce / NOAA / OAR
Geophysical Fluid Dynamics Laboratory voice: (609) 452-6547
201 Forrestal Road fax: (609) 987-5063
Princeton, NJ 08540-6649 http://www.gfdl.noaa.gov/~cjs/
********************************************************************
"The contents of this message are mine personally and do not reflect
any position of the Government or NOAA."
More information about the gradsusr
mailing list