MITgcm/doc/diags_changes.txt

Jun 21, 2009
Some redundancies and relations for model diagnostics:

1) Qnet = QNETtave = - oceQnet = SIqnet

2) For open-ocean, i.e., when there is no sea ice:
   Qnet = QNETtave = - oceQnet = SIqnet = EXFqnet = SIqneto = SIatmQnt

3) EXFqnet = EXFlwnet + EXFswnet - EXFhl - EXFhs

4) Qsw = QSWtave = - oceQsw = SIqsw

5) For open-ocean, i.e., when there is no sea ice:
   Qsw = QSWtave = - oceQsw = SIqsw = EXFswnet

6) EmPmR = EmPmRtave = -oceFWflx = SIempmr

================================================================

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.
--------------------------------------------------------------------------
1	dimitri	1.1	Jun 21, 2009
2			Some redundancies and relations for model diagnostics:
3
4			1) Qnet = QNETtave = - oceQnet = SIqnet
5
6	dimitri	1.2	2) For open-ocean, i.e., when there is no sea ice:
7			Qnet = QNETtave = - oceQnet = SIqnet = EXFqnet = SIqneto = SIatmQnt
8	dimitri	1.1
9	dimitri	1.2	3) EXFqnet = EXFlwnet + EXFswnet - EXFhl - EXFhs
10	dimitri	1.1
11	dimitri	1.2	4) Qsw = QSWtave = - oceQsw = SIqsw
12
13			5) For open-ocean, i.e., when there is no sea ice:
14			Qsw = QSWtave = - oceQsw = SIqsw = EXFswnet
15
16			6) EmPmR = EmPmRtave = -oceFWflx = SIempmr
17	dimitri	1.1
18			================================================================
19
20			Nov 19, 2006 (after tag checkpoint58r_post)
21
22			I) Some diagnostics have been renamed (essentially, to better
23			match the content):
24
25			1) name: description (oceanic set-up):
26			PRESSURE Cell-Center Height
27			replaced by:
28			RCENTER Cell-Center Height
29
30			2) name: description (oceanic set-up):
31			TICE heat from melt/freeze of sea-ice, >0 increases theta
32			replaced by:
33			oceFreez heating from freezing of sea-water (allowFreezing=T)
34
35			3) name: description (oceanic set-up):
36			TAUX zonal surface wind stress, >0 increases uVel
37			TAUY meridional surf. wind stress, >0 increases vVel
38			replaced by:
39			oceTAUX zonal surface wind stress, >0 increases uVel
40			oceTAUY meridional surf. wind stress, >0 increases vVel
41
42			4) name: description (oceanic set-up):
43			SWFLUX net upward SW radiation, >0 increases theta
44			replaced by:
45			oceQsw net Short-Wave radiation (+=down), >0 increases theta
46
47			5) name: description (oceanic set-up):
48			DIFx_TH Zonal Diffusive Flux of Pot.Temperature
49			DIFy_TH Meridional Diffusive Flux of Pot.Temperature
50			replaced by:
51			DFxE_TH Zonal Diffusive Flux of Pot.Temperature
52			DFyE_TH Meridional Diffusive Flux of Pot.Temperature
53
54			6) name: description (oceanic set-up):
55			DIFx_SLT Zonal Diffusive Flux of Salinity
56			DIFy_SLT Meridional Diffusive Flux of Salinity
57			replaced by:
58			DFxE_SLT Zonal Diffusive Flux of Salinity
59			DFyE_SLT Meridional Diffusive Flux of Salinity
60
61			-----------------------------------------------
62			II) Change description of existing diagnostics:
63
64			name: old description (oceanic set-up):
65			TFLUX net surface heat flux, >0 increases theta
66			SFLUX net surface salt flux, >0 increases salt
67			name: new description (oceanic set-up):
68			TFLUX total heat flux (match heat-content variations), >0 increases theta
69			SFLUX total salt flux (match salt-content variations), >0 increases salt
70
71			-----------------------------------------------
72			III) New diagnostics have been added:
73
74			name: description (oceanic set-up):
75			atmPload Atmospheric pressure loading
76			sIceLoad sea-ice loading (in Mass of ice+snow / area unit)
77			oceQnet net surface heat flux into the ocean (+=down), >0 increases theta
78			oceFWflx net surface Fresh-Water flux into the ocean (+=down), >0 decreases salinity
79			oceSflux net surface Salt flux into the ocean (+=down), >0 increases salinity
80			surForcT model surface forcing for Temperature, >0 increases theta
81			surForcS model surface forcing for Salinity, >0 increases salinity
82
83			--------------------------------------------------------------------------
84			Relation between surForcT,surForcS and others surface forcing diagnostics:
85			[x] = average of model variable "x" over the diagnostic time interval
86
87			a) if useRealFreshWaterFlux=F or (nonlinFreeSurf=0 & usingZCoords=T)
88			surForcT = oceQnet + TRELAX - oceQsw
89			surForcS = oceSflux + SRELAX - [PmEpR*So]
90			(with So = local Sea-Surface Salinity (SSS) if convertFW2Salt=-1
91			and So = convertFW2Salt otherwise)
92			oceFWflx = [PmEpR]
93			TFLUX = surForcT + oceQsw + oceFreez
94			SFLUX = surForcS
95
96			b) if useRealFreshWaterFlux=T & (nonlinFreeSurf>0 or usingPCoords=T),
97			In general:
98			surForcT = oceQnet + TRELAX - oceQsw + [T_dilution_effect]*Cp
99			surForcS = oceSflux + SRELAX + [S_dilution_effect]
100			where T_dilution_effect = PmEpR*( temp_EvPrRn - SST )
101			and S_dilution_effect = PmEpR*( salt_EvPrRn - SSS )
102			oceFWflx = [PmEpR]
103			TFLUX = surForcT + oceQsw + oceFreez + [PmEpRSST]Cp
104			SFLUX = surForcS + [PmEpR*SSS]
105
106			And with the default value: salt_EvPrRn=0. & temp_EvPrRn=UNSET_RL
107			(=> no dilution effect on Temp.):
108			surForcT = oceQnet + TRELAX - oceQsw
109			surForcS = oceSflux + SRELAX - [PmEpR*SSS]
110
111			Notes:
112			1) here PmEpR is assumed to be the fresh-water mass flux
113			per surface area [units: kg/m^2/s] whereas the model variable
114			EmPmR is still a volume flux per surface area [units: m/s].
115			2) with Linear Free surface (nonlinFreeSurf=0), the term
116			corresponding to w_surface*SST,SSS is missing in TFLUX,SFLUX
117			[might be added later ?] to match exactly the Heat and Salt
118			budget evolution.
119			--------------------------------------------------------------------------