;; Example HOPS script #8
;; interp: Linearly interpolating data from one grid to another

;; Open 1 year of prognostic ocean data (24 samples at half-month intervals)
;; for a 3-degree NCAR Ocean Model integration
;; (converted to CSM netCDF format using Frank Bryan's filters)

hopen, "/data/large1/CSM/3x3/G03-4.24/Y0989_0990_p.nc"

;; Get 3-D temperature field for first sample
T= hget("T",t=1,subscript=1)

;; Interpolate to 6-level uniformly spaced depth grid, with extrapolation
TINT= hinterp(T, 'z', grid=[0,1000,2000,3000,4000,5000], extrapolate=1)

;; Plot zonal-mean temperature on interpolated levels
hplot, hsub(TINT,x="avg"), fill=1

end