load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin

a=addfile("test1.nc","r")
b=addfile("test2.nc","r")
psl=a->PSL(0,21:50,38:68)
cloud=b->CLOUD(0,9:25,21:50,38:68)
pcp=(b->PRECL(0,21:50,38:68))+(b->PRECC(0,21:50,38:68))
pcp=pcp*60.*60.*24.*1000. ; mm/d
lev=doubletofloat(a->lev(9:25))
lev=lev(::-1)
copy_VarCoords(a->lev(9:25),lev)
lat=doubletofloat(a->lat(21:50))
lon=doubletofloat(a->lon(38:68))

lat=lat(::-1)

copy_VarCoords(a->lat(21:50),lat)
copy_VarCoords(a->lon(38:68),lon)

psl=psl(::-1,:)
pcp=pcp(::-1,:)
cloud=cloud(:,::-1,:)

psl!0="lat"
psl!1="lon"
psl&lat=lat
psl&lon=lon
pcp!0="lat"
pcp!1="lon"
pcp&lat=lat
pcp&lon=lon
cloud!0="lev"
cloud!1="lat"
cloud!2="lon"
cloud&lev=lev
cloud&lat=lat
cloud&lon=lon

numtimes = 1      
numvars = 3  
nr = 30 
nc = 31 
nl = 17
nlev=17

varnames=(/"data_psl","data_pcp","data_cloud"/)

timestamp = (/1/)
datestamp = (/1/)
compress = 1 
northlat = 4.736842 
latinc = 1.89473
westlon=-95.
loninc = 2.5 
projection = 1    
proj_args = (/northlat, westlon, latinc, loninc/)
vert = 3 
vert_args=lev

data_psl=new((/nlev,nc,nr/),float)
data_pcp=new((/nlev,nc,nr/),float)
data_cloud=new((/nlev,nc,nr/),float)

do i=0,nlev-1
  data_psl(i,:,:)=psl(lon|:,lat|:)
  data_pcp(i,:,:)=pcp(lon|:,lat|:)
end do
data_cloud=cloud(lev|:,lon|:,lat|:)

data_psl!0="lev"
data_psl!1="lon"
data_psl!2="lat"
data_psl&lev=lev
data_psl&lon=lon
data_psl&lat=lat

data_pcp!0="lev"
data_pcp!1="lon"
data_pcp!2="lat"
data_pcp&lev=lev
data_pcp&lon=lon
data_pcp&lat=lat

data_cloud!0="lev"
data_cloud!1="lon"
data_cloud!2="lat"
data_cloud&lev=lev
data_cloud&lon=lon
data_cloud&lat=lat

system("rm thin.v5d")
v5d_create("thin.v5d", numtimes, numvars, nr, nc, (/nl,nl,nl/), varnames, timestamp, datestamp, compress, projection, proj_args, vert, vert_args)
v5d_write(numtimes, 1, data_psl)
v5d_write(numtimes, 2, data_pcp)
v5d_write(numtimes, 3, data_cloud)
v5d_close()

end