// 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;