# Diagnostic Package Choices # # 6-hourly 2-D fields: # ==================== # # ETAN Perturbation of Surface (pressure, height) (Pa,m) # PHIBOT ocean bottom pressure / top. atmos geo-Potential # oceTAUX zonal surface wind stress, >0 increases uVel (N/m^2) # oceTAUY meridional surf. wind stress, >0 increases vVel (N/m^2) # surForcT model surface forcing for Temperature, >0 increases theta (W/m^2) # oceQsw net Short-Wave radiation (+=down), >0 increases theta (W/m^2) # surForcS model surface forcing for Salinity, >0 increases salinity (g/m^2/s) # KPPhbl KPP boundary layer depth, bulk Ri criterion # KPPmld Mixed layer depth, dT=.8degC density criterion # SSS Sea Surface Salinity (g/kg) # SST Sea Surface Temperature (degC,K) # UVEL_k2 Zonal Component of Velocity at level 2 (m/s) # VVEL_k2 Meridional Component of Velocity at level 2 (m/s) # # daily 2-D fields: # ================= # # SIarea SEAICE fractional ice-covered area [0 to 1] # SIheff SEAICE effective ice thickness (m) # SIuice SEAICE zonal ice velocity, >0 from West to East (m/s) # SIvice SEAICE merid. ice velocity, >0 from South to North (m/s) # SIhsnow SEAICE snow thickness (m) # # monthly 2-D fields: # =================== # # ETANSQ Square of Perturbation of Sfc (Pa2,m2) # EXFhs Sensible heat flux into ocean, >0 increases theta (W/m^2) # EXFhl Latent heat flux into ocean, >0 increases theta (W/m^2) # EXFlwnet Net upward longwave radiation, >0 decreases theta (W/m^2) # oceFWflx net surface Fresh-Water flux into ocean, >0 decreases salinity (kg/m^2/s) # oceSflux net surface Salt flux into the ocean, >0 increases salinity (g/m^2/s) # oceQnet net surface heat flux into the ocean, >0 increases theta (W/m^2) # SRELAX surface salinity relaxation, >0 increases salt (g/m^2/s) # TFLUX total heat flux (match heat-content variations), >0 increases theta (W/m^2) # SFLUX total salt flux (match salt-content variations), >0 increases salt (g/m^2/s) # # monthly 3-D fields: # =================== # # SALTanom Salt anomaly (=SALT-35; g/kg) # THETA Potential Temperature (degC,K) # UVEL Zonal Component of Velocity (m/s) # VVEL Meridional Component of Velocity (m/s) # UVELMASS Zonal Mass-Weighted Comp of Velocity (m/s) # VVELMASS Meridional Mass-Weighted Comp of Velocity (m/s) # WVELMASS Vertical Mass-Weighted Comp of Velocity (m/s) # # SALTSQan Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2) # THETASQ Square of Potential Temperature (K2) # UVELSQ Square of Zonal Comp of Velocity (m2/s2) # VVELSQ Square of Meridional Comp of Velocity (m2/s2) # WVELSQ Square of Vertical Comp of Velocity (m2/s2) # UV_VEL_Z Meridional Transport of Zonal Momentum (m2/s2) # WU_VEL Vertical Transport of Zonal Momentum (m^2/s^2) # WV_VEL Vertical Transport of Meridional Momentum (m^2/s^2) # # UTHMASS Zonal Mass-Weight Transp of Pot Temp (K.m/s) # VTHMASS Meridional Mass-Weight Transp of Pot Temp (K.m/s) # WTHMASS Vertical Mass-Weight Transp of Pot Temp (K.m/s) # USLTMASS Zonal Mass-Weight Transp of Salt (g/kg.m/s) # VSLTMASS Meridional Mass-Weight Transp of Salt (g/kg.m/s) # WSLTMASS Vertical Mass-Weight Transp of Salt (g/kg.m/s) # # RHOAnoma Density Anomaly (=Rho-rhoConst; kg/m^3) # DRHODR Stratification: d.Sigma/dr (kg/m3/r_unit; kg/m^4) # RHOANOSQ Square of Density Anomaly (=(Rho-rhoConst)^2; kg^2/m^6) # URHOMASS Zonal Transport of Density (kg/m^2/s) # VRHOMASS Meridional Transport of Density (kg/m^2/s) # WRHOMASS Vertical Transport of Potential Density (kg/m^2/s) ===================================================================== Nov 19, 2006 (after tag checkpoint58r_post) I) Some diagnostics have been renamed (essentially, to better match the content): 1) name: description (oceanic set-up): PRESSURE Cell-Center Height replaced by: RCENTER Cell-Center Height 2) name: description (oceanic set-up): TICE heat from melt/freeze of sea-ice, >0 increases theta replaced by: oceFreez heating from freezing of sea-water (allowFreezing=T) 3) name: description (oceanic set-up): TAUX zonal surface wind stress, >0 increases uVel TAUY meridional surf. wind stress, >0 increases vVel replaced by: oceTAUX zonal surface wind stress, >0 increases uVel oceTAUY meridional surf. wind stress, >0 increases vVel 4) name: description (oceanic set-up): SWFLUX net upward SW radiation, >0 increases theta replaced by: oceQsw net Short-Wave radiation (+=down), >0 increases theta 5) name: description (oceanic set-up): DIFx_TH Zonal Diffusive Flux of Pot.Temperature DIFy_TH Meridional Diffusive Flux of Pot.Temperature replaced by: DFxE_TH Zonal Diffusive Flux of Pot.Temperature DFyE_TH Meridional Diffusive Flux of Pot.Temperature 6) name: description (oceanic set-up): DIFx_SLT Zonal Diffusive Flux of Salinity DIFy_SLT Meridional Diffusive Flux of Salinity replaced by: DFxE_SLT Zonal Diffusive Flux of Salinity DFyE_SLT Meridional Diffusive Flux of Salinity ----------------------------------------------- II) Change description of existing diagnostics: name: old description (oceanic set-up): TFLUX net surface heat flux, >0 increases theta SFLUX net surface salt flux, >0 increases salt name: new description (oceanic set-up): TFLUX total heat flux (match heat-content variations), >0 increases theta SFLUX total salt flux (match salt-content variations), >0 increases salt ----------------------------------------------- III) New diagnostics have been added: name: description (oceanic set-up): atmPload Atmospheric pressure loading sIceLoad sea-ice loading (in Mass of ice+snow / area unit) oceQnet net surface heat flux into the ocean (+=down), >0 increases theta oceFWflx net surface Fresh-Water flux into the ocean (+=down), >0 decreases salinity oceSflux net surface Salt flux into the ocean (+=down), >0 increases salinity surForcT model surface forcing for Temperature, >0 increases theta surForcS model surface forcing for Salinity, >0 increases salinity -------------------------------------------------------------------------- Relation between surForcT,surForcS and others surface forcing diagnostics: [x] = average of model variable "x" over the diagnostic time interval a) if useRealFreshWaterFlux=F or (nonlinFreeSurf=0 & usingZCoords=T) surForcT = oceQnet + TRELAX - oceQsw surForcS = oceSflux + SRELAX - [PmEpR*So] (with So = local Sea-Surface Salinity (SSS) if convertFW2Salt=-1 and So = convertFW2Salt otherwise) oceFWflx = [PmEpR] TFLUX = surForcT + oceQsw + oceFreez SFLUX = surForcS b) if useRealFreshWaterFlux=T & (nonlinFreeSurf>0 or usingPCoords=T), In general: surForcT = oceQnet + TRELAX - oceQsw + [T_dilution_effect]*Cp surForcS = oceSflux + SRELAX + [S_dilution_effect] where T_dilution_effect = PmEpR*( temp_EvPrRn - SST ) and S_dilution_effect = PmEpR*( salt_EvPrRn - SSS ) oceFWflx = [PmEpR] TFLUX = surForcT + oceQsw + oceFreez + [PmEpR*SST]*Cp SFLUX = surForcS + [PmEpR*SSS] And with the default value: salt_EvPrRn=0. & temp_EvPrRn=UNSET_RL (=> no dilution effect on Temp.): surForcT = oceQnet + TRELAX - oceQsw surForcS = oceSflux + SRELAX - [PmEpR*SSS] Notes: 1) here PmEpR is assumed to be the fresh-water mass flux per surface area [units: kg/m^2/s] whereas the model variable EmPmR is still a volume flux per surface area [units: m/s]. 2) with Linear Free surface (nonlinFreeSurf=0), the term corresponding to w_surface*SST,SSS is missing in TFLUX,SFLUX [might be added later ?] to match exactly the Heat and Salt budget evolution. -------------------------------------------------------------------------- C pkg/diagnostics SIempmr, dumpfreq EmPmR, and tavefreq EmPmRtave C are identical but they differ from pkg/diagnostics EXFempmr, which C is EmPmR before impact of ice. C SIqnet, Qnet, and QNETtave are identical. C With #undef NONLIN_FRSURF SIqnet is identical to -(TFLUX-TRELAX). C Except over land and under sea ice, SIqnet is also identical to C EXFlwnet+EXFswnet-EXFhl-EXFhs. C SIqsw, Qsw, and QSWtave are identical. C Except under sea ice, SIqsw is also identical to EXFswnet. C pkg/diagnostics SIfu and TAUX, dumpfreq FU, and tavefreq FUtave C are identical but they differ from pkg/diagnostics EXFtaux, which C is stress before impact of ice. Also when using exf bulk C formulae, EXFtaux is defined on tracer rather than uvel points. C pkg/diagnostics SIfv and TAUY, dumpfreq FV, and tavefreq FVtave C are identical but they differ from pkg/diagnostics EXFtauy, which C is stress before impact of ice. Also when using exf bulk C formulae, EXFtauy is defined on tracer rather than vvel points.