1 |
Jun 21, 2009 |
2 |
Some redundancies and relations for model diagnostics: |
3 |
|
4 |
1) Qnet = QNETtave = - oceQnet = SIqnet |
5 |
|
6 |
2) For open-ocean, i.e., when there is no sea ice: |
7 |
Qnet = QNETtave = - oceQnet = SIqnet = EXFqnet = SIqneto = SIatmQnt |
8 |
|
9 |
3) EXFqnet = EXFlwnet + EXFswnet - EXFhl - EXFhs |
10 |
|
11 |
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 |
|
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 + [PmEpR*SST]*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 |
-------------------------------------------------------------------------- |