seasonal average
Charles Seman
Charles.Seman at NOAA.GOV
Thu Nov 9 19:09:35 EST 2006
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."
More information about the gradsusr
mailing list