153 |
surForcS = oceSflux + SRELAX - [PmEpR*SSS] |
surForcS = oceSflux + SRELAX - [PmEpR*SSS] |
154 |
|
|
155 |
Notes: |
Notes: |
156 |
1) here PmEpR is assumed to be the fresh-water mass flux |
1) Here PmEpR is assumed to be the fresh-water mass flux |
157 |
per surface area [units: kg/m^2/s] whereas the model variable |
per surface area [units: kg/m^2/s] whereas the model variable |
158 |
EmPmR is still a volume flux per surface area [units: m/s]. |
EmPmR is still a volume flux per surface area [units: m/s]. |
159 |
2) with Linear Free surface (nonlinFreeSurf=0), the term |
2) With Linear FreeSurf, there is a small loss of heat/salt/tracer at the |
160 |
corresponding to w_surface*SST,SSS is missing in TFLUX,SFLUX |
surface (d/dt(Eta) inconsistent with fixed surf grid cell thickness: |
161 |
[might be added later ?] to match exactly the Heat and Salt |
a) if linFSConserveTr is set to true (file "data"), then the model |
162 |
budget evolution. |
corrects for this (see S/R calc_wsurf_tr.F). |
163 |
-------------------------------------------------------------------------- |
b) if linFSConserveTr=False (=the default), the term corresponding |
164 |
|
to w_surface*SST,SSS is missing in TFLUX,SFLUX to match exactly the |
165 |
|
Heat and Salt budget evolution so need to account for this term in |
166 |
|
the budget, which can be obtained from the 1rst level |
167 |
|
of WTHMASS / WSLTMASS diagnostics (diagnostics ADVr_TH / ADVr_SLT |
168 |
|
is not the right one for this purpose, it's just zero at k=1). |
169 |
|
3) There is just one more undocumented/hidden case |
170 |
|
(nonlinFreeSurf=0 & select_rStar=-1) which is not covered. |