seasonal average

praise praise praise_city1 at HOTMAIL.COM
Sun Nov 12 07:52:51 EST 2006 

Dear Chuck

Thanks again

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: Fri, 10 Nov 2006 17:59:12 -0500
>
>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."

_________________________________________________________________
MSN Tool Bar 幫你刪除惱人的廣告 ! http://toolbar.msn.com.hk

More information about the gradsusr mailing list