[Thread Prev][Thread Next][Index]
Re: [ferret_users] how to compute the 5-day averaged climatology
Hi Peter,
It is a 142-day daily dataset.
set data imet_epic_pentad.sflx.mc
yes? show data/var
currently SET data sets:
1> ./imet_epic_pentad.sflx.mc (default)
Bouy Net Surface Heat Flux at the IMET Site
name title I J K L
SFLX Net Surface Heat Flux 1:1 1:1 1:1 1:142
W/m^2 on grid GMK1 with 1.000000E+35 for missing data
X=85.7W:84.7W Y=20.7S:19.7S Z=-0.7:0.3
time range: 03-JAN-2003 12:00 to 08-DEC-2004 12:00
I try other datasets which cover 5 years. it still does not work.
The major purpose is that I want to calculate the 5-day averaged anomalies.
I do not know if anyone else has a good idea to calculate the 5-day average climatology for several years 5-day datasets.
By the way, I try regrid 5-day data into daily dataset, like what the FAQ:
set data imet_epic_pentad.sflx.mc
define axis/t=03-JAN-2003:08-DEC-2004:1/units=days tsflxdaily
let sflx_daily = sflx[gt=tsflxdaily]
let sflx_clim = sflx_daily[gt=tdaily@mod]
but it does not work, either.
Anyone can help me out? Thanks!
Yangxing
--- On Mon, 7/7/08, Peter Szabo <szabpet83@xxxxxxxxx> wrote:
> From: Peter Szabo <szabpet83@xxxxxxxxx>
> Subject: Re: [ferret_users] how to compute the 5-day averaged climatology
> To: yangxingz@xxxxxxxxx
> Date: Monday, July 7, 2008, 8:58 PM
> hi Yangxingz,
> i see now. i have different data set (daily and i thought
> you wanted a 5-day
> averaged data set). anysleay wrote only two lines what i
> included in my last
> message.
>
> if you followed the daily climatology in the faq, then it
> should work.
> maybe you only have one year thats why it lists the first
> year as the only
> one to be taken into the "climatology". and it
> seems for me you only have
> one year.
> "SFLX[T=01-JAN-2003 00:00:01-JAN-2004 00:00]"
>
> another problem could be when you have your data on
> non-calendar axis.
>
> Peter
>
>
> On Fri, Jul 4, 2008 at 6:21 PM, yangxing zheng
> <yangxingz@xxxxxxxxx> wrote:
>
> > Hi Peter,
> >
> > Happy holiday!
> > Thank for your reply. It seems that your definition
> of sflx_clim is just
> > 5-day average from daily dataset, not really 5-day
> climatology. Also
> > actually I have 5-day average datasets, so I do not
> need to do conversion.
> >
> > It looks to me that regriding with time axis does
> not work properly. I do
> > not know why. For example, I try regrid 5-day average
> data to daily dataset,
> > as I define the daily axis "tdaily", then
> use
> > let sflx_daily = sflx[gt=tdaily]
> > as I list sflx_daily
> > the data is missing for all grid points.
> >
> > Could you send me the script that Ansley did it?
> >
> > Thanks.
> >
> > Yangxing
> >
> >
> > --- On Fri, 7/4/08, Peter Szabo
> <szabpet83@xxxxxxxxx> wrote:
> >
> > > From: Peter Szabo <szabpet83@xxxxxxxxx>
> > > Subject: Re: [ferret_users] how to compute the
> 5-day averaged climatology
> > > To: yangxingz@xxxxxxxxx
> > > Cc: ferret_users@xxxxxxxx
> > > Date: Friday, July 4, 2008, 4:12 PM
> > > Hi Yangxing,
> > > I also came across this problem, but someone
> helped me out
> > > (Ansley).
> > >
> > > Assume you have 1-day Net Surface Heat Flux and
> you want
> > > 5-day average of
> > > it, then you should use Boxcar-smoother.
> > > yourvar[L=@SBX:5] starts at L=3. and this is the
> average of
> > > the data at L=1
> > > through L=5.
> > > So you only need to do the following:
> > >
> > > let sflx_clim = sflx[gt=@sbx:5]
> > > !though the two of the start and end points will
> become
> > > missing data, this
> > > is not an axis-changing function
> > > list sflx_clim[l=3:365:5]
> > > !every 5 point will represent the 5-day average.
> and it
> > > will produce 73
> > > points.
> > >
> > > Hope this helps,
> > > Peter Szabo
> > >
> > >
> > >
> > > On Fri, Jul 4, 2008 at 12:13 AM, yangxing zheng
> > > <yangxingz@xxxxxxxxx> wrote:
> > >
> > > > Hi,
> > > >
> > > > This may not be a good and new question,
> but I am
> > > still wondering how to
> > > > do it.
> > > >
> > > > I read FAQ about how to create daily
> climatology, so
> > > I attempt to modify
> > > > it following the method. Below is what I
> change
> > > >
> > > > ---------------------------
> > > > !Get the length of a climatological
> year:365.2425
> > > > use climatological_axes
> > > > cancel data climatological_axes
> > > >
> > > > !use the length of the "climatological
> year"
> > > to define a pentad axis
> > > > define
> > >
> axis/t=0:365.2425:5/edges/npoints=73/T0=1-JAN-0001/Units=days/Modulo
> > > > t_pentad
> > > > !define
> > >
> axis/t=0:365.2425:5/edges/npoints=73/Units=days/Modulo
> > > t_pentad
> > > >
> > > > !for example:
> > > > set data imet_epic_pentad.sflx.mc
> > > >
> > > > !define the pentad climatology
> > > >
> > > > let sflx_clim = sflx[gt=t_pentad@mod]
> > > >
> > > > list sflx_clim,sflx[l=1:73]
> > > > ---------------------
> > > > The data imet_epic_pentad.sflx.mc is 5-day
> averaged
> > > surface heat flux that
> > > > I use as an example.
> > > >
> > > > I run the script, it does not any error
> message, but
> > > the result is not
> > > > correct.The first column is to list
> climatology, the
> > > second column is the
> > > > first year original datasets.
> > > > ---------
> > > > list sflx_clim,sflx[l=1:73]
> > > > WARNING: Listed variables have ambiguous
> coordinates
> > > on axes: T
> > > > DATA SET:
> ./imet_epic_pentad.sflx.mc
> > > > Bouy Net Surface Heat Flux at
> the IMET
> > > Site
> > > > LONGITUDE: 85.2W
> > > > LATITUDE: 20.2S
> > > > DEPTH (m): -0.2
> > > > Column 1: SFLX_CLIM[T=01-JAN
> 00:00:01-JAN-0002
> > > 05:49] is
> > > > SFLX[GT=T_PENTAD@MOD]
> > > > Column 2: SFLX[T=01-JAN-2003
> 00:00:01-JAN-2004
> > > 00:00] is Net Surface Heat
> > > > Flux (W/m^2)
> > > > SFLX_CLI SFLX
> > > > L / 1: .... 117.7
> > > > L / 2: .... 109.5
> > > > L / 3: .... 106.3
> > > > L / 4: .... 129.3
> > > > L / 5: .... 150.8
> > > > L / 6: .... 150.6
> > > > L / 7: .... 129.6
> > > > L / 8: .... 124.7
> > > > L / 9: .... 94.9
> > > > L / 10: .... 45.9
> > > > L / 11: .... 70.1
> > > > L / 12: .... 64.0
> > > > L / 13: .... 70.8
> > > > L / 14: .... 74.9
> > > > L / 15: .... 96.2
> > > > L / 16: .... 109.6
> > > > L / 17: .... 86.3
> > > > L / 18: .... 32.1
> > > > L / 19: .... 33.8
> > > > L / 20: .... 17.5
> > > > L / 21: .... 38.2
> > > > L / 22: .... -8.8
> > > > L / 23: .... 34.6
> > > > L / 24: .... -16.8
> > > > L / 25: .... -0.6
> > > > L / 26: .... 14.2
> > > > L / 27: .... -23.1
> > > > L / 28: .... -44.0
> > > > L / 29: .... 5.7
> > > > L / 30: .... 0.0
> > > > L / 31: .... -12.4
> > > > L / 32: .... -16.2
> > > > L / 33: .... 61.3
> > > > L / 34: .... -67.6
> > > > L / 35: .... -13.5
> > > > L / 36: .... -26.4
> > > > L / 37: .... 8.1
> > > > L / 38: .... -41.7
> > > > L / 39: .... -79.3
> > > > L / 40: .... -74.3
> > > > L / 41: .... -106.8
> > > > L / 42: .... -89.4
> > > > L / 43: .... -22.7
> > > > L / 44: .... -69.4
> > > > L / 45: .... -14.2
> > > > L / 46: .... 28.7
> > > > L / 47: .... -48.0
> > > > L / 48: .... -40.4
> > > > L / 49: .... 27.8
> > > > L / 50: .... 62.6
> > > > L / 51: .... -76.2
> > > > L / 52: .... 6.0
> > > > L / 53: .... 42.4
> > > > L / 54: .... -19.9
> > > > L / 55: .... 73.7
> > > > L / 56: .... 77.9
> > > > L / 57: .... 61.7
> > > > L / 58: .... 123.5
> > > > L / 59: .... 50.8
> > > > L / 60: .... 71.1
> > > > L / 61: .... 82.4
> > > > L / 62: .... 55.4
> > > > L / 63: .... 46.9
> > > > L / 64: .... 155.8
> > > > L / 65: .... 166.7
> > > > L / 66: .... 120.6
> > > > L / 67: .... 114.6
> > > > L / 68: .... 109.0
> > > > L / 69: .... 96.4
> > > > L / 70: .... 120.2
> > > > L / 71: .... 59.5
> > > > L / 72: .... 73.1
> > > > L / 73: 36.63 50.2
> > > > -------------
> > > > Can anyone help me out?
> > > >
> > > > Happy Independence Day!
> > > >
> > > > Yangxing
> > > >
> > > >
> > > >
> > > >
> > > >
> >
> >
> >
> >
[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Privacy Policy | Disclaimer | Accessibility Statement