Re: [ferret_users] how to do climatology

[prescilla kurien <prescy@xxxxxxxxxxxxxxxxx>]

Hello Aditi:
I am attaching the programs which I used to make climatology from last 5
years output of model simulation
1) Append the all data to a single file (append_5_year_data.jnl)
2) Use ferret transformation to make monthly climatology
(CLIMATOLOGY.jnl)
Hope this helps
Prescilla 

 
On Wed, 2009-09-16 at 08:50 +0530, aditi srivast wrote:
> Dear ferret users,
> I am having 40 years model output (1 year) (1967-2007) with 12 months
> in each file. I would like to make climatology by using these 40 years
> data, all the these files are in Netcdf, I want the final output with
> only 12 months , how do I do climatological averages using ferret.
> 
> 
> Thank you,
> 
> S Aditi

SET MEMORY/SIZE=500
use 2010_hori_V_mom.nc
DEFINE AXIS/T=5-jan-2010:31-dec-2010:15/UNITS=DAYS t1nr
DEFINE GRID/T=t1nr my_t1grid
save/file=VVEL_appnd.nc VVEL[d=1,GT=my_t1grid]
use 2011_hori_V_mom.nc
DEFINE AXIS/T=15-jan-2011:26-dec-2011:15/UNITS=DAYS t2nr
DEFINE GRID/T=t2nr my_t2grid
save/append/file=VVEL_appnd.nc VVEL[d=2,GT=my_t2grid]
use 2012_hori_V_mom.nc
DEFINE AXIS/T=10-jan-2012:20-dec-2012:15/UNITS=DAYS t3nr
DEFINE GRID/T=t3nr my_t3grid
save/append/file=VVEL_appnd.nc VVEL[d=3,GT=my_t3grid]
use 2013_hori_V_mom.nc
DEFINE AXIS/T=5-jan-2013:31-dec-2013:15/UNITS=DAYS t4nr
DEFINE GRID/T=t4nr my_t4grid
save/append/file=VVEL_appnd.nc VVEL[d=4,GT=my_t4grid]
use 2014_hori_V_mom.nc
DEFINE AXIS/T=15-jan-2014:26-dec-2014:15/UNITS=DAYS t5nr
DEFINE GRID/T=t5nr my_t5grid
save/append/file=VVEL_appnd.nc VVEL[d=5,GT=my_t5grid]
SHOW DATA/VARIABLES 
QUIT

!this program is to create climatology
SET MEMORY/SIZE=500
use VVEL_appnd.nc
USE climatological_axes
CANCEL DATA climatological_axes
LET VVEL_CLIM=VVEL[GT=month_reg@MOD]
!SET VARIABLE/BAD=-999.9 SAL_CLIM 
save/file=HIGH_CLIM_VVEL.nc VVEL_CLIM 
SHADE/K=1/L=1 VVEL_CLIM
QUIT