Changeset 15548 for NEMO/branches/2021/ticket2632_r14588_theta_sbcblk/src/TOP/PISCES/P4Z/p4zligand.F90

Timestamp:

2021-11-28T18:59:49+01:00 (3 years ago)

Author:

gsamson

Message:

update branch to the head of the trunk (r15547); ticket #2632

Location:

NEMO/branches/2021/ticket2632_r14588_theta_sbcblk

Files:

• . (modified) (1 prop)
• src/TOP/PISCES/P4Z/p4zligand.F90 (modified) (7 diffs)

NEMO/branches/2021/ticket2632_r14588_theta_sbcblk

@@ -9 +9 @@
 
 # SETTE
-^/utils/CI/sette@14244        sette
+^/utils/CI/sette@HEAD        sette
+

@@ -9 +9 @@
 
 # SETTE
-^/utils/CI/sette@14244        sette
+^/utils/CI/sette@HEAD        sette
+

NEMO/branches/2021/ticket2632_r14588_theta_sbcblk/src/TOP/PISCES/P4Z/p4zligand.F90

@@ -25 +25 @@
    REAL(wp), PUBLIC ::  rlig     !: Remin ligand production
    REAL(wp), PUBLIC ::  prlgw    !: Photochemical of weak ligand
+   REAL(wp), PUBLIC ::  xklig    !: 1/2 saturation constant of photolysis
 
    !! * Substitutions
@@ -45 +46 @@
       !
       INTEGER  ::   ji, jj, jk
-      REAL(wp) ::   zlgwp, zlgwpr, zlgwr, zlablgw
-      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zligrem, zligpr, zligprod
+      REAL(wp) ::   zlgwp, zlgwpr, zlgwr, zlablgw 
+      REAL(wp) ::   zlam1a, zlam1b, zaggliga, zligco
+      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zligrem, zligpr, zligprod, zlcoll3d
       CHARACTER (len=25) ::   charout
       !!---------------------------------------------------------------------
@@ -52 +54 @@
       IF( ln_timing )   CALL timing_start('p4z_ligand')
       !
-      DO_3D( 1, 1, 1, 1, 1, jpkm1 )
+      DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1)
          !
          ! ------------------------------------------------------------------
@@ -62 +64 @@
          ! This is based on the idea that as LGW is lower
          ! there is a larger fraction of refractory OM
-         zlgwr = max( rlgs , rlgw * exp( -2 * (tr(ji,jj,jk,jplgw,Kbb)*1e9) ) ) ! years
-         zlgwr = 1. / zlgwr * tgfunc(ji,jj,jk) * ( xstep / nyear_len(1) ) * blim(ji,jj,jk) * tr(ji,jj,jk,jplgw,Kbb)
+         zlgwr = ( 1.0 / rlgs * MAX(0., tr(ji,jj,jk,jplgw,Kbb) - xfecolagg(ji,jj,jk) * 1.0E-9 )    &
+         &       + 1.0 / rlgw * xfecolagg(ji,jj,jk) * 1.0E-9 ) / ( rtrn + tr(ji,jj,jk,jplgw,Kbb) )
+         zlgwr = zlgwr * tgfunc(ji,jj,jk) * ( xstep / nyear_len(1) ) * blim(ji,jj,jk) * tr(ji,jj,jk,jplgw,Kbb)
          ! photochem loss of weak ligand
-         zlgwpr = prlgw * xstep * etot(ji,jj,jk) * tr(ji,jj,jk,jplgw,Kbb) * (1. - fr_i(ji,jj))
-         tr(ji,jj,jk,jplgw,Krhs) = tr(ji,jj,jk,jplgw,Krhs) + zlgwp - zlgwr - zlgwpr
+         zlgwpr = prlgw * xstep * etot(ji,jj,jk) * tr(ji,jj,jk,jplgw,Kbb)**3 * (1. - fr_i(ji,jj))   &
+         &        / ( tr(ji,jj,jk,jplgw,Kbb)**2 + (xklig)**2)
+         ! Coagulation of ligands due to various processes (Brownian, shear, diff. sedimentation
+         ! xcoagfe is computed in p4zfechem
+         ! -------------------------------------------------------------------------------------
+         ! 50% of the ligands are supposed to be in the colloidal size fraction
+         ! as for FeL
+         zligco   = 0.5 * MAX(0., tr(ji,jj,jk,jplgw,Kbb) - xfecolagg(ji,jj,jk) * 1.0E-9 )
+         zaggliga = xcoagfe(ji,jj,jk) * xstep * zligco
+
+         tr(ji,jj,jk,jplgw,Krhs) = tr(ji,jj,jk,jplgw,Krhs) + zlgwp - zlgwr - zlgwpr - zaggliga
+         !
          zligrem(ji,jj,jk)   = zlgwr
          zligpr(ji,jj,jk)    = zlgwpr
-         zligprod(ji,jj,jk) = zlgwp
-         !
+         zligprod(ji,jj,jk)  = zlgwp
+         zlcoll3d(ji,jj,jk)  = zaggliga
       END_3D
       !
@@ -84 +97 @@
          IF( iom_use( "LPRODR" ) ) THEN
            zligprod(:,:,jpk) = 0. ; CALL iom_put( "LPRODR", zligprod(:,:,:) * 1e9 * 1.e+3 * rfact2r * tmask(:,:,:) )
+         ENDIF
+         IF( iom_use( "LGWCOLL" ) ) THEN
+            zlcoll3d(:,:,jpk) = 0. ; CALL iom_put( "LGWCOLL", zlcoll3d(:,:,:) * 1.e9 *  1.e+3 * rfact2r * tmask(:,:,:) )
          ENDIF
       ENDIF
@@ -110 +126 @@
       INTEGER ::   ios   ! Local integer 
       !
-      NAMELIST/nampislig/ rlgw, prlgw, rlgs, rlig
+      NAMELIST/nampislig/ rlgw, prlgw, rlgs, rlig, xklig
       !!----------------------------------------------------------------------
       !
@@ -130 +146 @@
          WRITE(numout,*) '      Photolysis of weak ligand                    prlgw =', prlgw
          WRITE(numout,*) '      Lifetime (years) of strong ligands           rlgs  =', rlgs
+         WRITE(numout,*) '      1/2 saturation for photolysis                xklig =', xklig
       ENDIF
       !

@@ -9 +9 @@
 
 # SETTE
-^/utils/CI/sette@14244        sette
+^/utils/CI/sette@HEAD        sette
+