Changeset 3475 for branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90

Timestamp:

2012-09-13T19:00:23+02:00 (12 years ago)

Author:

cetlod

Message:

branch:2012/dev_r3438_LOCEAN15_PISLOB last updates from PISCES, hopefully... see ticket #972

File:

• branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90 (modified) (5 diffs)

branches/2012/dev_r3438_LOCEAN15_PISLOB/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zfechem.F90

@@ -34 +34 @@
    LOGICAL          ::  ln_ligvar  = .FALSE.    !: boolean for variable ligand concentration following Tagliabue and voelker
    REAL(wp), PUBLIC ::  xlam1      = 0.005_wp   !: scavenging rate of Iron 
-   REAL(wp), PUBLIC ::  xlamdust   = 100.0_wp   !: scavenging rate of Iron by dust 
+   REAL(wp), PUBLIC ::  xlamdust   = 150.0_wp   !: scavenging rate of Iron by dust 
    REAL(wp), PUBLIC ::  ligand     = 0.6E-9_wp  !: ligand concentration in the ocean 
 
@@ -236 +236 @@
                ! to later allocate scavenged iron to the different organic pools
                ! ---------------------------------------------------------
+               zdenom1 = xlam1 * trn(ji,jj,jk,jppoc) / zlam1b
 #if  defined key_kriest
-               zdenom1 = xlam1 * trn(ji,jj,jk,jppoc) / zlam1b 
-#else
-               zdenom1 = xlam1 * trn(ji,jj,jk,jppoc) / zlam1b
                zdenom2 = xlam1 * trn(ji,jj,jk,jpgoc) / zlam1b
 #endif
@@ -246 +244 @@
                !  due to increased lithogenic particles and let say it is unknown processes (precipitation, ...)
                !  -----------------------------------------------------------
+               zlamfac = MAX( 0.e0, ( gphit(ji,jj) + 55.) / 30. )
+               zlamfac = MIN( 1.  , zlamfac )
+               zdep    = MIN( 1., 1000. / fsdept(ji,jj,jk) )
                zlam1b  = xlam1 * MAX( 0.e0, ( trn(ji,jj,jk,jpfer) * 1.e9 - ztotlig(ji,jj,jk) ) )
-               zcoag   = zfeequi * zlam1b * zstep
+               zcoag   = zfeequi * zlam1b * zstep + 1E-4 * ( 1. - zlamfac ) * zdep * zstep *zfecoll
 
                !  Compute the coagulation of colloidal iron. This parameterization 
@@ -253 +254 @@
                !  It requires certainly some more work as it is very poorly constrained.
                !  ----------------------------------------------------------------
-               zlamfac = MAX( 0.e0, ( gphit(ji,jj) + 55.) / 30. )
-               zlamfac = MIN( 1.  , zlamfac )
-               zdep    = MIN( 1., 1000. / fsdept(ji,jj,jk) )
                zlam1a  = ( 0.369  * 0.3 * trn(ji,jj,jk,jpdoc) + 102.4  * trn(ji,jj,jk,jppoc) ) * xdiss(ji,jj,jk)    &
                    &   + ( 114.   * 0.3 * trn(ji,jj,jk,jpdoc) + 5.09E3 * trn(ji,jj,jk,jppoc) )
+               zaggdfea = zlam1a * zstep * zfecoll
 #if defined key_kriest
-               zlam1a   = zlam1a + 1E-4 * ( 1. - zlamfac ) * zdep
-               zaggdfea = zlam1a * zstep * zfecoll
                zaggdfeb = 0.
                !
@@ -267 +264 @@
 #else
                zlam1b = 3.53E3 *   trn(ji,jj,jk,jpgoc) * xdiss(ji,jj,jk) + 1E-4 * ( 1. - zlamfac ) * zdep 
-               zaggdfea = zlam1a * zstep * zfecoll
                zaggdfeb = zlam1b * zstep * zfecoll
                !