Changeset 10356 for NEMO/trunk/src/TOP/PISCES/SED/sedchem.F90

Timestamp:

2018-11-22T12:39:19+01:00 (5 years ago)

Author:

cetlod

Message:

Bugfix in sediment model, see ticket #2171

File:

• NEMO/trunk/src/TOP/PISCES/SED/sedchem.F90 (modified) (2 diffs)

NEMO/trunk/src/TOP/PISCES/SED/sedchem.F90

@@ -678 +678 @@
          &         - 0.07711*zsal + 0.0041249*zsal15
 
+         ! CONVERT FROM DIFFERENT PH SCALES
+         total2free  = 1.0/(1.0 + zst/zcks)
+         free2SWS    = 1. + zst/zcks + zft/(zckf*total2free)
+         total2SWS   = total2free * free2SWS
+         SWS2total   = 1.0 / total2SWS
+
+
          ! K1, K2 OF CARBONIC ACID, KB OF BORIC ACID, KW (H2O) (LIT.?)
          zak1    = 10**(zck1) * total2SWS
@@ -743 +750 @@
          aksis(ji) = zaksi * EXP( zbuf1 * zcpexp + zbuf2 * zcpexp2 )
 
-         ! CONVERT FROM DIFFERENT PH SCALES
-         total2free  = 1.0/(1.0 + zst/aks3s(ji))
-         free2SWS    = 1. + zst/aks3s(ji) + zft/akf3s(ji)
-         total2SWS   = total2free * free2SWS
-         SWS2total   = 1.0 / total2SWS
-
          ! Convert to total scale
          ak1s(ji)  = ak1s(ji)  * SWS2total