I am sharing
what I did to save all the models into one file as suggested
by Russ, Satyesh and Ansley and its worked fine for me. Hope
this will help you as well.
!!!!!!!!!!!!!!
suppose you have 11 models !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!! all the models are having same variable
!!!!!!!!!!!!!!!!!
use
model1.nc
use
model2.nc
use "
use "
use "
use "
use "
use "
use "
use "
use model11
!!!!!!!!!!!!!!!! change the variable name for ensembling them
!!!!!!!!!!!!!!
set var/name=histbcc HISTZD[d=2]
set var/name=histccs HISTZD[d=3]
set var/name=histcmc HISTZD[d=4]
set var/name=histcnr HISTZD[d=5]
set var/name=histcsi HISTZD[d=6]
set var/name=histece HISTZD[d=7]
set var/name=histgfd HISTZD[d=8]
set var/name=histmpi HISTZD[d=9]
set var/name=histmri HISTZD[d=10]
set var/name=histnor HISTZD[d=11]
!!!!!!!!!!!!!!!!!!! regrid them for the same lat, lon, depth
and time !!!!!!!!!!!!!
!!!!!!! in my case only depth was different so i regrid it to
depth only !!!!!!!!!!
!!!!!!!!!!! I regrid with respect to my first data set you can
choose your own data set for regredding !!!!!!
let/d=2 histzd=histbcc[gz=HISTZD[d=1]@nrst]
let/d=3 histzd=histccs[gz=HISTZD[d=1]@nrst]
let/d=4 histzd=histcmc[gz=HISTZD[d=1]@nrst]
let/d=5 histzd=histcnr[gz=HISTZD[d=1]@nrst]
let/d=6 histzd=histcsi[gz=HISTZD[d=1]@nrst]
let/d=7 histzd=histece[gz=HISTZD[d=1]@nrst]
let/d=8 histzd=histgfd[gz=HISTZD[d=1]@nrst]
let/d=9 histzd=histmpi[gz=HISTZD[d=1]@nrst]
let/d=10 histzd=histmri[gz=HISTZD[d=1]@nrst]
let/d=11 histzd=histnor[gz=HISTZD[d=1]@nrst]
!!!!!!!!!!!!!!! make them into ensemble !!!!!!!!!!!!!!!!
!!!!!!!!!! numbers represents number of ensemble members i.e.
11 in my case !!!!!!!!!!!!!!
!!!!!!!!!! but the name of ensemble variable will we same as
original variable to save into netcdf file !!!!!!!!!!
ensemble histens=1,2,3,4,5,6,7,8,9,10,11
!!!!!!!!!!!!!!!
it will show you ensemble variable attribute
!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!! chunck them to improve efficieny as
discussed by russ in email !!!!!!!!!!!!!!!
!!!!!!!!!! I chuncked them into 4 depth levels and 1 year
!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!! 1 year because I was haivng monthly time series so
I chuncked into 12 months !!!!!!!!!
!!!!!!!!!! I was having 312 monthly time steps (26 years) and
36 vertical depth levels !!!!!!!!!!!!!!!!!!
set
list/xchunk=360/ychunk=126/zchunk=9/tchunk=12/echunk=1/deflate=1/shuffle/ncformat=4
!!!!!!!!!!!!!!!!!!!!!! define limits of each axes
!!!!!!!!!!!!!!!!!!!!!!!!!!
save/i=1/j=1/k=1/l=1/m=1/ilimits=1:360/jlimits=1:126/klimits=1:36/llimits=1:312/mlimits=1:11/file=
hist_russ.nc/clob
histzd
!!!!!!!!!!!!!!!!!!!!!!!!! save the ensemble data by a nested
loop !!!!!!!!!!!!!!!!!!!!!!!!!!
repeat/range=1:312:12/name=ll (repeat/range=1:11/name=mm
(repeat/range=1:36:9/name=kk save/append/file=
hist_russ.nc/l=`ll`:`ll+11`/m=`mm`/k=`kk`:`kk+8`
histzd;set mem/size=5000);set mem/size=9000)
!!!!!!!!!!!!! that is done !!!!!!!!!!!!!!!!!!!!!!!
exit
Hope this
will help you.