Changeset 14786 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zlys.F90

Timestamp:

2021-05-05T09:09:10+02:00 (3 years ago)

Author:

aumont

Message:

Various bug fixes in PISCES

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zlys.F90 (modified) (3 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zlys.F90

@@ -60 +60 @@
       INTEGER  ::   ji, jj, jk, jn
       REAL(wp) ::   zdispot, zfact, zcalcon
-      REAL(wp) ::   zomegaca, zexcess, zexcess0
+      REAL(wp) ::   zomegaca, zexcess, zexcess0, zkd
       CHARACTER (len=25) ::   charout
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zco3, zcaldiss, zhinit, zhi, zco3sat
@@ -68 +68 @@
       !
       zhinit  (:,:,:) = hi(:,:,:) * 1000. / ( rhop(:,:,:) + rtrn )
+      zco3    (:,:,:) = 0._wp     ;   zco3sat(:,:,:) = 0._wp
+      zcaldiss(:,:,:) = 0._wp
       !
       !     -------------------------------------------
@@ -99 +101 @@
          excess(ji,jj,jk) = 1._wp - zomegaca
          zexcess0 = MAX( 0., excess(ji,jj,jk) )
-         zexcess  = zexcess0**nca
 
-         ! AMOUNT CACO3 (12C) THAT RE-ENTERS SOLUTION
-         !       (ACCORDING TO THIS FORMULATION ALSO SOME PARTICULATE
-         !       CACO3 GETS DISSOLVED EVEN IN THE CASE OF OVERSATURATION)
-         zdispot = kdca * zexcess * tr(ji,jj,jk,jpcal,Kbb)
+         IF (zomegaca < 0.8) THEN
+            zexcess  = zexcess0**nca
+            ! AMOUNT CACO3 THAT RE-ENTERS SOLUTION
+            zdispot = kdca * zexcess * tr(ji,jj,jk,jpcal,Kbb)
+         ELSE
+
+            zkd = kdca * 0.1**(nca - 0.4)
+            zexcess  = zexcess0**0.4
+            zdispot = zkd * zexcess * tr(ji,jj,jk,jpcal,Kbb)
+        ENDIF
+
         !  CHANGE OF [CO3--] , [ALK], PARTICULATE [CACO3],
         !       AND [SUM(CO2)] DUE TO CACO3 DISSOLUTION/PRECIPITATION