For some reason, the metadata for your output has _FillValue and missing_value attributes for all of the grid definition variables (TLAT, TLONG, HTE, dz, etc.). I don't know why this happened yet, but it does cause problems in the ocean diagnostics. Here's a fix:
Inside gen_TS_meandiffrms_timeseries.csh, after the lines
cp ${TCLIMFILE} obs.nc
cp ${SCLIMFILE} SALT_obs.nc
ncks -A -v SALT SALT_obs.nc obs.nc
ncks -A -v TAREA,REGION_MASK ${CASE}.pop.h.${year}.nc obs.nc
add these lines
ncatted -a missing_value,TLAT,d,, obs.nc
ncatted -a missing_value,TLONG,d,, obs.nc
ncatted -a _FillValue,TLAT,d,, obs.nc
ncatted -a _FillValue,TLONG,d,, obs.nc
This will remove these two attributes from the obs.nc file and should allow the ncwa command to work.
Susan
