[Thread Prev][Thread Next][Index]

Re: [ferret_users] Monthly mean




Hi,

You're chopping off the end of your data with all the filtering. Is there any reason to do this? Have a look at uflsm and you'll see what I mean. Also defining your monthly axis isn't really optimal. I'd also use the @AVE transform to do your remapping to the monthly axis. This gets all the values within the relevant time window and weights them appropriately..

Try something like the following

! Nifty jnl script in contrib makes a monthly axis from 1993 to 2016

go def_monthaxis_days.jnl STANDARD 1993 2016 mytaxis

let ufmon=uf[gt=mytaxis@AVE]

Russ



On 28/02/17 16:12, Arijeet Dutta wrote:
Hi all
I have data in every 5 days interval, from 1993 to 2016, I want to
create monthly mean from it. I did the following

define axis/t=1:8766:1/t0="01-jan-1993:00:00"/unit=day tday ! from 5day to daily
let ufl = uf[gt=tday] ! interpolate (linearly by default)
let uflsm = uf[l=@spz:31] ! a filter with half-power point of about 37 days
define axis/t=15-jan-1993:15-dec-2016:1/t0="15-jan-1993"/ unit=month tmon
let ufmon=uflsm[gt=tmon] ! subsample to monthly

But when I am listing ufmon, I get values from l=4, no values from l=1 to l=3
Please let me know the mistake

Thanks & reg
Arijeet




[Thread Prev][Thread Next][Index]
Contact Us
Dept of Commerce / NOAA / OAR / PMEL / Ferret

Privacy Policy | Disclaimer | Accessibility Statement