seasonal average

praise praise praise_city1 at HOTMAIL.COM
Fri Nov 10 00:17:43 EST 2006 

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

More information about the gradsusr mailing list