Mercurial > hgweb > burst-coupling

changeset 76:c8ed926f0017

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Tuned biogeochemistry parameters from Joerg's _17 case (used for CMIP simulations).
author Marco van Hulten <marco@hulten.org>
date Fri, 26 Apr 2019 16:37:56 +0200
parents a345c6f2748f
children 79f08d3ffa55
files beleg_bgc.F90
diffstat 1 files changed, 12 insertions(+), 13 deletions(-) [+]
line wrap: on
line diff
--- a/beleg_bgc.F90	Wed Apr 10 17:43:26 2019 +0200
+++ b/beleg_bgc.F90	Fri Apr 26 16:37:56 2019 +0200
@@ -151,10 +151,10 @@
 
 !ik addded parameter definition; taken from OCPROD.F
       remido=0.004*dtb      !1/d -remineralization rate (of DOM)
-      dyphy=0.008*dtb       !1/d -mortality rate of phytoplankton 
-      grazra=1.0*dtb        !1/d -grazing rate
+      dyphy=0.004*dtb       !1/d -mortality rate of phytoplankton 
+      grazra=1.2*dtb        !1/d -grazing rate
       spemor=3.*1.e6*dtb    !1/d -mortality rate
-      gammap=0.03*dtb       !1/d -exudation rate
+      gammap=0.04*dtb       !1/d -exudation rate
       gammaz=0.06*dtb       !1/d -excretion rate
       ecan=0.95             ! fraction of mortality as PO_4
       pi_alpha=0.02*0.4     ! initial slope of production vs irradiance curve (alpha) (0.002 for 10 steps per day)
@@ -162,8 +162,8 @@
       zinges = 0.5          !dimensionless fraction -assimilation efficiency
       epsher = 0.9          !dimensionless fraction -fraction of grazing egested
 #elif defined(WLIN)
-      zinges = 0.55         !dimensionless fraction -assimilation efficiency
-      epsher = 0.85         !dimensionless fraction -fraction of grazing egest      
+      zinges = 0.7          !dimensionless fraction -assimilation efficiency
+      epsher = 0.9          !dimensionless fraction -fraction of grazing egested
 #else
       zinges = 0.6          !dimensionless fraction -assimilation efficiency
       epsher = 0.8          !dimensionless fraction -fraction of grazing egest      
@@ -182,8 +182,8 @@
 
 ! half sat. constants, note that the units are kmol/m3 !
       bkphy  = 4.e-8    !i.e. 0.04 mmol P/m3
-      bkzoo  = 4.e-8    !i.e. 0.04 mmol P/m3
-      bkopal = 1.e-6    !i.e. 1.0 mmol Si/m3
+      bkzoo  = 8.e-8    !i.e. 0.08 mmol P/m3
+      bkopal = 5.e-6    !i.e. 4.0 mmol Si/m3
 
 !sinking speed
       wpoc  =  5.*dtb       !m/d  iris : 5.
@@ -197,12 +197,11 @@
 
       
 ! deep see remineralisation constants
-
       drempoc  = 0.025*dtb    !1/d
       dremdoc  = 0.004*dtb    !1/d
-      dphymor  = 0.008 *dtb   !1/d
+      dphymor  = 0.004*dtb    !1/d
       dzoomor  = 3.*1.e6*dtb  !1/d -mortality rate
-      dremopal = 0.005*dtb    !1/d
+      dremopal = 0.003*dtb    !1/d
       dremn2o  = 0.01*dtb     !1/d
       dremsul  = 0.005*dtb    ! remineralization rate for sulphate reduction 
       
@@ -248,8 +247,8 @@
       ropal = 10.5 ! opal to organic phosphorous production ratio      
       calmax= 0.20
 #elif defined(WLIN)
-      rcalc = 40.  ! calcium carbonate to organic phosphorous production ratio
-      ropal = 30.  ! opal to organic phosphorous production ratio      
+      rcalc = 38.  ! calcium carbonate to organic phosphorous production ratio
+      ropal = 45.  ! opal to organic phosphorous production ratio      
 #else
       rcalc = 40.  ! iris 40 !calcium carbonate to organic phosphorous production ratio
       ropal = 30.  ! iris 25 !opal to organic phosphorous production ratio      
@@ -271,7 +270,7 @@
 
 !ik weight percent iron in dust deposition times Fe solubility
 ! the latter three values come from Johnson et al., 1997
-      fetune=0.5                  ! factor introduced to tune deposistion/solubility
+      fetune=0.6                  ! factor introduced to tune deposistion/solubility
       perc_diron = fetune * 0.035 * 0.01 / 55.85
       riron= 5.*rcar*1.e-6        ! fe to P ratio in organic matter
       fesoly=0.5*1.e-9            ! max. diss. iron concentration in deep water