Changeset 13233 for NEMO/branches/2019/dev_r11708_aumont_PISCES_QUOTA/src/TOP/PISCES/P4Z/p4zrem.F90

Timestamp:

2020-07-02T20:34:16+02:00 (8 months ago)

Author:

aumont

Message:

update of the PISCES comments

File:

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

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

@@ -92 +92 @@
       ! Computation of the mean bacterial concentration
       ! this parameterization has been deduced from a model version
-      ! that was modeling explicitely bacteria
-      ! -------------------------------------------------------
+      ! that was modeling explicitely bacteria. This is a very old param 
+      ! that will be very soon updated based on results from a much more
+      ! recent version of PISCES with bacteria.
+      ! ----------------------------------------------------------------
       DO jk = 1, jpkm1
          DO jj = 1, jpj
@@ -111 +113 @@
       END DO
 
-      IF( ln_p4z ) THEN
+      IF( ln_p4z ) THEN ! Standard PISCES code
          DO jk = 1, jpkm1
             DO jj = 1, jpj
@@ -152 +154 @@
             END DO
          END DO
-      ELSE
+      ELSE ! PISCES-QUOTA
          DO jk = 1, jpkm1
             DO jj = 1, jpj
@@ -231 +233 @@
 
                ! Bacterial uptake of iron. No iron is available in DOC. So
-               ! Bacteries are obliged to take up iron from the water. Some
+               ! Bacteria are obliged to take up iron from the water. Some
                ! studies (especially at Papa) have shown this uptake to be significant
                ! ---------------------------------------------------------------------
@@ -262 +264 @@
          DO jj = 1, jpj
             DO ji = 1, jpi
-               ! Remineralization rate of BSi depedant on T and saturation
+
+               ! Remineralization rate of BSi dependent on T and saturation
                ! The parameterization is taken from Ridgwell et al. (2002) 
                ! ---------------------------------------------------------
@@ -269 +272 @@
                zsatur2  = ( 1. + tsn(ji,jj,jk,jp_tem) / 400.)**37
                znusil   = 0.225  * ( 1. + tsn(ji,jj,jk,jp_tem) / 15.) * zsatur + 0.775 * zsatur2 * zsatur**9.25
+
+               ! Two fractions of bSi are considered : a labile one and a more
+               ! refractory one based on the commonly observed two step 
+               ! dissolution of bSi (initial rapid dissolution followed by 
+               ! more slowly dissolution).
                ! Computation of the vertical evolution of the labile fraction
                ! of bSi. This is computed assuming steady state.
+               ! --------------------------------------------------------------
                IF ( gdept_n(ji,jj,jk) > zdep ) THEN
                   zfacsib(ji,jj,jk) = zfacsib(ji,jj,jk-1) * EXP( -0.5 * ( xsiremlab - xsirem )  &
@@ -280 +289 @@
                zsiremin = ( xsiremlab * zfacsi(ji,jj,jk) + xsirem * ( 1. - zfacsi(ji,jj,jk) ) ) * xstep * znusil
                zosil    = zsiremin * trb(ji,jj,jk,jpgsi)
+
                ! Update of the TRA arrays
                tra(ji,jj,jk,jpgsi) = tra(ji,jj,jk,jpgsi) - zosil