Changeset 15287 for NEMO/trunk/src/TOP

Timestamp:

2021-09-24T13:11:02+02:00 (3 years ago)

Author:

cetlod

Message:

trunk : bugfixes in PISCES

Location:

NEMO/trunk/src/TOP/PISCES/P4Z

Files:

• p4zlys.F90 (modified) (2 diffs)
• p4zsed.F90 (modified) (2 diffs)

NEMO/trunk/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
@@ -98 +98 @@
          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.2**(nca - 0.11)
+            zexcess = zexcess0**0.11
+            zdispot = zkd * zexcess * tr(ji,jj,jk,jpcal,Kbb)
+        ENDIF
+
         !  CHANGE OF [CO3--] , [ALK], PARTICULATE [CACO3],
         !       AND [SUM(CO2)] DUE TO CACO3 DISSOLUTION/PRECIPITATION

NEMO/trunk/src/TOP/PISCES/P4Z/p4zsed.F90

@@ -265 +265 @@
             tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0
             tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zfact / 3.0
+            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zfact * 2.0 / 3.0            
             tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) - zfact * 2.0 / 3.0
             tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + zfact * 1.0 / 3.0
@@ -282 +283 @@
             tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) + zfact / 3.0
             tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * zfact / 3.0
+            tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zfact * 2.0 / 3.0            
             tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) - 16.0 / 46.0 * zfact * ( 1.0 - 1.0 / 3.0 ) &
             &                     * ztrpo4(ji,jj,jk) / (ztrpo4(ji,jj,jk) + ztrdop(ji,jj,jk) + rtrn)