seasonal average
praise praise
praise_city1 at HOTMAIL.COM
Wed Nov 8 20:25:37 EST 2006
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
More information about the gradsusr
mailing list