[Thread Prev][Thread Next][Index]
[ferret_users] how to compute the 5-day averaged climatology
- To: ferret_users@xxxxxxxx
- Subject: [ferret_users] how to compute the 5-day averaged climatology
- From: yangxing zheng <yangxingz@xxxxxxxxx>
- Date: Thu, 03 Jul 2008 15:13:08 -0700 (PDT)
- Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=s1024; d=yahoo.com; h=Received:X-Mailer:Date:From:Reply-To:Subject:To:MIME-Version:Content-Type:Message-ID; b=Cumr9W1UpRNnV2KqbVlzEfVWcBdIjpXJgw9jrP571goDUfY9IBVYD1fNQFbyRdkTwZueypq0KkyyVeB2mMISZWICu+coyRKXZRIUXZsM6kD5WQqL5lepfsLFbdYcgqojSd13M/NrM/qkc+RL1KiNNNyJ3Oh1eOOSeHG1Kx78vDE=;
- Reply-to: yangxingz@xxxxxxxxx
- Sender: owner-ferret_users@xxxxxxxx
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