Changeset 14276 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zlys.F90

Timestamp:

2021-01-07T23:09:56+01:00 (3 years ago)

Author:

aumont

Message:

numerous updates to PISCES, PISCES-QUOTA and the sediment module

File:

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

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zlys.F90

@@ -57 +57 @@
       !
       INTEGER  ::   ji, jj, jk, jn
-      REAL(wp) ::   zdispot, zfact, zcalcon
+      REAL(wp) ::   zdispot, zfact, zcalcon, zkd
       REAL(wp) ::   zomegaca, zexcess, zexcess0
       CHARACTER (len=25) ::   charout
@@ -65 +65 @@
       IF( ln_timing )  CALL timing_start('p4z_lys')
       !
-      zco3    (:,:,:) = 0.
-      zcaldiss(:,:,:) = 0.
+      zco3    (:,:,:) = 0._wp
+      zcaldiss(:,:,:) = 0._wp     ;     zco3sat(:,:,:) = 0._wp
       zhinit  (:,:,:) = hi(:,:,:) * 1000. / ( rhop(:,:,:) + rtrn )
       !
@@ -106 +106 @@
                excess(ji,jj,jk) = 1._wp - zomegaca
                zexcess0 = MAX( 0., excess(ji,jj,jk) )
-               zexcess  = zexcess0**nca
+               IF (zomegaca < 0.8) THEN
+                  zexcess  = zexcess0**nca
+                 ! AMOUNT CACO3 THAT RE-ENTERS SOLUTION
+                  zdispot = kdca * zexcess * trb(ji,jj,jk,jpcal)
+               ELSE
+                  zkd = kdca * 0.1**(nca - 0.4)
+                  zexcess  = zexcess0**0.4
+                  zdispot = zkd * zexcess * trb(ji,jj,jk,jpcal)
+               ENDIF
 
-               ! AMOUNT CACO3 THAT RE-ENTERS SOLUTION
-               zdispot = kdca * zexcess * trb(ji,jj,jk,jpcal)
                ! CHANGE OF [CO3--] , [ALK], PARTICULATE [CACO3],
                ! AND [SUM(CO2)] DUE TO CACO3 DISSOLUTION/PRECIPITATION