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 + [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 |
|
|
-------------------------------------------------------------------------- |