Changeset 14276 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedini.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/SED/sedini.F90 (modified) (2 diffs)

NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/SED/sedini.F90

@@ -371 +371 @@
       denssol = 2.6
 
-      ! Initialization of diffusion coefficient as function of porosity [cm**2/s]
-      !--------------------------------------------------------------------
-!      DO jn = 1, jpsol
-!         DO jk = 1, jpksed
-!            DO ji = 1, jpoce
-!               diff(ji,jk,jn) = dcoef / ( 1.0 - 2.0 * log(por(jk)) )
-!            END DO
-!         END DO
-!      END DO
-
       ! Accumulation rate from Burwicz et al. (2011). This is used to
       ! compute the flux of clays and minerals
@@ -387 +377 @@
           ztmp1 = 0.117 / ( 1.0 + ( zkbot(ji) / 200.)**3 )
           ztmp2 = 0.006 / ( 1.0 + ( zkbot(ji) / 4000.)**10 )
-          wacc(ji) = ztmp1 + ztmp2
+          wacc(ji) = (ztmp1 + ztmp2)/10.0
+          wacc(ji) = ztmp2 / 10.0
       END DO