Opened 14 years ago
Closed 14 years ago
#675 closed Bug (fixed)
Missing salt flux into the ocean when using vvl and sea-ice
| Reported by: | acc | Owned by: | acc |
|---|---|---|---|
| Priority: | low | Milestone: | |
| Component: | OCE | Version: | v3.2 |
| Severity: | Keywords: | ||
| Cc: |
Description
In tra_sbc (trasbc.F90), when running with the vvl option, the salinity trend due to the surface boundary condition is set to zero. The dilution/concentration effect due to water exchange at the surface is supposedly accounted for by using the "true layer thickness" in the non-linear free surface case. However, it is emp that is used for the layer thickness computation not emps. Therefore the contribution due to the formation/melting of sea-ice and any SSS restoration that is not converted to a freshwater flux (nn_sssr=1) will be lost. This omission will eventually be corrected by the reformulation of the heat and mass fluxes underway in ticket #665. Until that work is complete, the missing trend can be added by replacing line 136:
zsa = 0.0e0
with:
zsa = (emps(ji,jj) - emp(ji,jj)) * zsrau * sn(ji,jj,1) * zse3t
Commit History (1)
| Changeset | Author | Time | ChangeLog |
|---|---|---|---|
| 1892 | acc | 2010-06-01T16:00:31+02:00 | Added missing salt flux in trasbc.F90 (vvl case only), see ticket #675 |
Change History (1)
comment:1 Changed 14 years ago by acc
- Resolution set to fixed
- Status changed from new to closed
