Hi George, Ok, it looks like your original grid is actually on a gregorian calendar (or not specified). To find out what type of calendar precip is on try list (`precip,ret=taxis`).calendar http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/variables-xpressions/VARIABLES#_VPID_81 (see the second example). You need to use the MONTH_IRREG time axis from climatological_axes.nc. You will also have to define your new monthly time axis as gregorian. Try go def_monthaxis_days.jnl gregorian 1979 2015 new_time ! Looks like you have a 365 day calendar let prec_on_monthly = precip[gt=new_time@nrst] ! get values closest let prec_clim_tmp = prec_on_monthly[gt=month_irreg@MOD] ! let prec_clim = prec_clim_tmp[gt=new_time@asn] ! Regrid to monthly axis let prec_anom_tmp = prec_on_monthly - prec_clim ! still got strange values on months not of interest let prec_anom = prec_anom_tmp[gt=precip@nrst] ! think this should get you back on the original grid. plot/X=80/Y=1.0 precip,prec_clim[gt=precip@nrst] ! plot original and climatology message plot/X=80/Y=1.0 prec_anom ! Anomaly For more info on time regridding see http://ferret.pmel.noaa.gov/Ferret/documentation/users-guide/Grids-Regions/GRIDS#_VPID_192 Cheers, Russ On 21/09/15 18:54, George Otieno wrote:
|