file='/glade/p/work/foster/timegcm_waccmx/timegcm_dres.s_waccmx_001.nc'
id=ncdf_open(file)
ncdf_varget,id,'lon',lon
ncdf_varget,id,'lat',lat
ncdf_varget,id,'lev',lev_tgcm
ncdf_varget,id,'TN',tn_tgcm
ncdf_varget,id,'TI',ti_tgcm
ncdf_varget,id,'TE',te_tgcm
ncdf_varget,id,'OP',ndop_tgcm

ncdf_close,id

nz=42
ndopr = dblarr(nz)
file1='Op_ambi.dat'
openr,1,file1
readf,1,ndopr
close,1

ndop = reverse(ndopr)

zplev = lev_tgcm(97-42:96)
tn = reform(tn_tgcm(72,36,97-42:96,0))
ti = reform(ti_tgcm(72,36,97-42:96,0))
te = reform(te_tgcm(72,36,97-42:96,0))
tp = .5*(ti+te)

file1='TN.dat'
openr,1,file1
readf,1,tn
close,1

file1='TP.dat'
openr,1,file1
readf,1,tp
close,1

rairv = dblarr(42)

file1='rairv.dat'
openr,1,file1
readf,1,rairv
close,1

opscaleheight,nz,ndop,tp,zplev,tn,rairv,sch,sch0

end