Changeset 14963 for NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zfechem.F90

Timestamp:

2021-06-09T18:14:48+02:00 (3 years ago)

Author:

aumont

Message:

Code updates + debugging + configuration files updates

File:

• NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zfechem.F90 (modified) (8 diffs)

NEMO/branches/2021/dev_r14383_PISCES_NEWDEV_PISCO/src/TOP/PISCES/P4Z/p4zfechem.F90

@@ -56 +56 @@
       INTEGER  ::   ji, jj, jk, jic, jn
       REAL(wp) ::   zlam1a, zlam1b
-      REAL(wp) ::   zkeq, zfesatur, zfecoll, fe3sol, zligco
+      REAL(wp) ::   zkeq, zfesatur, fe3sol, zligco
       REAL(wp) ::   zscave, zaggdfea, zaggdfeb, ztrc, zdust, zklight
       REAL(wp) ::   ztfe, zhplus, zxlam, zaggliga, zaggligb
@@ -62 +62 @@
       REAL(wp) ::   zrfact2
       CHARACTER (len=25) :: charout
-      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zTL1, zFe3, ztotlig, zfeprecip, zFeL1
+      REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zTL1, zFe3, ztotlig, zfeprecip, zFeL1, zfecoll
       REAL(wp), DIMENSION(jpi,jpj,jpk) ::   zcoll3d, zscav3d, zlcoll3d
       !!---------------------------------------------------------------------
@@ -75 +75 @@
       zscav3d  (:,:,jpk) = 0.
       zlcoll3d (:,:,jpk) = 0.
-
+      zfecoll  (:,:,jpk) = 0.
+      xfecolagg(:,:,jpk) = 0.
+      xcoagfe  (:,:,jpk) = 0.
+      !
       ! Total ligand concentration : Ligands can be chosen to be constant or variable
       ! Parameterization from Pham and Ito (2018)
       ! -------------------------------------------------
+      xfecolagg(:,:,:) = ligand * 1E9 + MAX(0., chemo2(:,:,:) - tr(:,:,:,jpoxy,Kbb) ) / 400.E-6
       IF( ln_ligvar ) THEN
-         ztotlig(:,:,:) =  0.09 * 0.667 * tr(:,:,:,jpdoc,Kbb) * 1E6 &
-           &             + ligand * 1E9 &
-           &             + MAX(0., chemo2(:,:,:) - tr(:,:,:,jpoxy,Kbb) ) / 400.E-6
+         ztotlig(:,:,:) =  0.09 * 0.667 * tr(:,:,:,jpdoc,Kbb) * 1E6 + xfecolagg(:,:,:)
          ztotlig(:,:,:) =  MIN( ztotlig(:,:,:), 10. )
       ELSE
@@ -111 +113 @@
       !
       zdust = 0.         ! if no dust available
+
+      ! Computation of the colloidal fraction that is subjecto to coagulation
+      ! The assumption is that 50% of complexed iron is colloidal. Furthermore
+      ! The refractory part is supposed to be non sticky. The refractory
+      ! fraction is supposed to equal to the background concentration + 
+      ! the fraction that accumulates in the deep ocean. AOU is taken as a 
+      ! proxy of that accumulation following numerous studies showing 
+      ! some relationship between weak ligands and AOU.
+      ! An issue with that parameterization is that when ligands are not
+      ! prognostic or non variable, all the colloidal fraction is supposed
+      ! to coagulate
+      ! ----------------------------------------------------------------------
+      IF (ln_ligand) THEN
+         zfecoll(:,:,:) = 0.5 * zFeL1(:,:,:) * MAX(0., tr(:,:,:,jplgw,Kbb) - xfecolagg(:,:,:) * 1.0E-9 ) / ( tr(:,:,:,jplgw,Kbb) + rtrn ) 
+      ELSE
+         IF (ln_ligvar) THEN
+            zfecoll(:,:,:) = 0.5 * zFeL1(:,:,:) * MAX(0., tr(:,:,:,jplgw,Kbb) - xfecolagg(:,:,:) * 1.0E-9 ) / ( tr(:,:,:,jplgw,Kbb) + rtrn )   
+         ELSE
+            zfecoll(:,:,:) = 0.5 * zFeL1(:,:,:)
+         ENDIF
+      ENDIF
+
       DO_3D( 1, 1, 1, 1, 1, jpkm1 )
          ! Scavenging rate of iron. This scavenging rate depends on the load of particles of sea water. 
@@ -121 +145 @@
          &         + fesol(ji,jj,jk,5) / zhplus )
          !
-         zfecoll = 0.5 * zFeL1(ji,jj,jk)
          ! precipitation of Fe3+, creation of nanoparticles
          zprecip = MAX( 0., ( zFe3(ji,jj,jk) - fe3sol ) ) * kfep * xstep * ( 1.0 - nitrfac(ji,jj,jk) ) 
@@ -144 +167 @@
              &    + ( 2.49  * tr(ji,jj,jk,jppoc,Kbb) )     &
              &    + ( 127.8 * 0.3 * tr(ji,jj,jk,jpdoc,Kbb) + 725.7 * tr(ji,jj,jk,jppoc,Kbb) )
-         zaggdfea = zlam1a * xstep * zfecoll
+         zaggdfea = zlam1a * xstep * zfecoll(ji,jj,jk)
                !
          zlam1b   = ( 1.94 * xdiss(ji,jj,jk) + 1.37 ) * tr(ji,jj,jk,jpgoc,Kbb)
-         zaggdfeb = zlam1b * xstep * zfecoll
-
+         zaggdfeb = zlam1b * xstep * zfecoll(ji,jj,jk)
+         xcoagfe(ji,jj,jk) = zlam1a + zlam1b
          !
          tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) - zscave - zaggdfea - zaggdfeb &
@@ -179 +202 @@
       biron(:,:,:) = tr(:,:,:,jpfer,Kbb) 
       !
-      IF( ln_ligand ) THEN
-         !
-         DO_3D( 1, 1, 1, 1, 1, jpkm1 )
-             ! Coagulation of ligands due to various processes (Brownian, shear, diff. sedimentation
-             ! Coefficients are taken from p4zagg
-             ! -------------------------------------------------------------------------------------
-             zlam1a   = ( 12.0  * 0.3 * tr(ji,jj,jk,jpdoc,Kbb) + 9.05  * tr(ji,jj,jk,jppoc,Kbb) ) * xdiss(ji,jj,jk)    &
-                 &    + ( 2.49  * tr(ji,jj,jk,jppoc,Kbb) )     &
-                 &    + ( 127.8 * 0.3 * tr(ji,jj,jk,jpdoc,Kbb) + 725.7 * tr(ji,jj,jk,jppoc,Kbb) )
-             !
-             zlam1b   = ( 1.94 * xdiss(ji,jj,jk) + 1.37 ) * tr(ji,jj,jk,jpgoc,Kbb)
-             ! 50% of the ligands are supposed to be in the colloidal size fraction
-             ! as for FeL
-             zligco   = 0.5 * tr(ji,jj,jk,jplgw,Kbb)
-             zaggliga = zlam1a * xstep * zligco 
-             zaggligb = zlam1b * xstep * zligco
-             !
-             tr(ji,jj,jk,jplgw,Krhs) = tr(ji,jj,jk,jplgw,Krhs)  - zaggliga - zaggligb
-             zlcoll3d(ji,jj,jk)  = zaggliga + zaggligb
-         END_3D
-      ENDIF
-
       !  Output of some diagnostics variables
       !     ---------------------------------
@@ -213 +214 @@
          IF( iom_use("FECOLL") )  CALL iom_put("FECOLL" , zcoll3d(:,:,:)  * 1e9 * tmask(:,:,:) * zrfact2 )
          IF( iom_use("FEPREC") )  CALL iom_put("FEPREC" , zfeprecip(:,:,:) *1e9*tmask(:,:,:)*zrfact2 )
-         IF( iom_use("LGWCOLL"))  CALL iom_put("LGWCOLL", zlcoll3d(:,:,:) * 1e9 * tmask(:,:,:) * zrfact2 )
       ENDIF