Changeset 7180 for branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90

Timestamp:

2016-11-03T16:41:10+01:00 (7 years ago)

Author:

aumont

Message:

various bug fixes on iron chemistry

File:

• branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90 (modified) (3 diffs)

branches/CNRS/dev_r6270_PISCES_QUOTA/NEMOGCM/NEMO/TOP_SRC/PISCES/P4Z/p4zligand.F90

@@ -75 +75 @@
                ! ---------------------------------------------------------
                zstep   = xstep
-               zstep2  = zstep / 365 ! per year
 # if defined key_degrad
                zstep = zstep * facvol(ji,jj,jk)
-               zstep2 = zstep2 * facvol(ji,jj,jk)
 # endif
+               zstep2  = zstep / 365. ! per year
                ! production from remineralisation of organic matter
                zlgwp  = orem(ji,jj,jk) * rlig
@@ -88 +87 @@
                zlgwr = 1. / zlgwr * tgfunc(ji,jj,jk) * zstep2 * trn(ji,jj,jk,jplgw)
                ! photochem loss of weak ligand
-               zlablgw = MAX( 0., trn(ji,jj,jk, jpfer) * plig(ji,jj,jk) )
                zlgwpr = prlgw * zstep * etot(ji,jj,jk) * trn(ji,jj,jk,jplgw) * (1. - fr_i(ji,jj))
                tra(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zlgwp - zlgwr - zlgwpr
@@ -105 +103 @@
                ! lower in the aphotici zone
                ! ! 25 Wm-2 constant
-               zrfepa = rfep * (1- EXP(-1*etot(ji,jj,jk) / 25. ) ) * (1.- fr_i(ji,jj))
+               zrfepa = rfep * (1- EXP(-etot(ji,jj,jk) / 25. ) ) * (1.- fr_i(ji,jj))
                zrfepa = MAX( (zrfepa / 10.0), zrfepa ) ! min of 10 days lifetime
                zfepr = rfep * zstep * trn(ji,jj,jk,jpfep)