;; 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)
hplot,tmtser

end