Changeset 617 for trunk/NEMO/TOP_SRC/SMS/p4zlys.F

Timestamp:

2007-02-21T14:03:10+01:00 (17 years ago)

Author:

opalod

Message:

* empty log message *

File:

• trunk/NEMO/TOP_SRC/SMS/p4zlys.F (modified) (3 diffs)

trunk/NEMO/TOP_SRC/SMS/p4zlys.F

@@ -44 +44 @@
       INTEGER ji, jj, jk, jn
       REAL zbot, zalk, zdic, zph, remco3, ah2
-      REAL delco3, excess, dispot, zfact, zalka
+      REAL dispot, zfact, zalka
+      REAL omegaca, excess, excess0
+      REAL co3(jpi,jpj,jpk)
 C
 C
@@ -114 +116 @@
 C ------------------------------------------------
 C
-            delco3 = co3(ji,jj,jk)-aksp(ji,jj,jk)/calcon
+            omegaca = ( calcon * co3(ji,jj,jk) )/aksp(ji,jj,jk)
+
 C
 C* 2.2  SET DEGREE OF UNDER-/SUPERSATURATION
 C ------------------------------------------
 C
-            excess = max(0.,-delco3)
+            excess0 = max(0.,(1.-omegaca))
+            excess = excess0**nca
+
 C
 C* 2.3  AMOUNT CACO3 (12C) THAT RE-ENTERS SOLUTION
@@ -126 +131 @@
 C --------------------------------------------------------------
 C
-            dispot = trn(ji,jj,jk,jpcal)*
-     &        excess/(dispo0+excess)
+            dispot = kdca * excess * trn(ji,jj,jk,jpcal)
 #    if defined key_off_degrad
      &        *facvol(ji,jj,jk)
 #    endif
+
 C
 C* 2.4  CHANGE OF [CO3--] , [ALK], PARTICULATE [CACO3],