// Example HOPS script #3
// ocn_plot2: Plotting data from ocean model diagnostic history file

// Open 1 year of diagnostic 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_d.nc";

// Extract horizontal mean temperature for the Atlantic basin and plot it
T_horza= hget("T_horz",slice=["basins","atlantic"],t=1,subscript=1);
hplot, T_horza;

// Plot meridional overturning streamfunction
tmt= hget("tmt",t=1,subscript=1);
hplot, tmt, fill=1;

// Extract time-series of spatially averaged meridional overturning streamfn
tmtser= hget("tmt",limy=[50,60],y="avg",limz=[2000,3000],z="avg",noweight=1);
fma;
hplot,tmtser;