;; Example HOPS script #4
;; save_file: Save two hyperslabs to a netCDF file

;; Open winter (DJF) mean ECMWF analysis data (1986-1992)
;; interpolated to the T42, 18-level hybrid sigma-pressure CCM2 grid
;; (converted to CSM netCDF format using Brian Eaton's ccm2nc)

hopen,"/data/svn/ccm3/ET42DJF8692.nc"

;; Get temperature (T) and zonal wind (U) for first sample
T= hget("T",t=1,subscript=1)

;; Compute zonal-average temperature and rename variable
TBAR= hsub(T,x="avg")
TBAR.name= "TBAR"

;; Save the 3-d temperature field and the 2-d zonal-mean temperature field
;; together as a hyperslab netCDF file
hsave,"t_and_tbar.nc", T, TBAR

end