Opened 10 years ago

Closed 10 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: release-3.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)

ChangesetAuthorTimeChangeLog
1892acc2010-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 10 years ago by acc

  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.