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