From a9b417b069cdc1b42cc34cd26f94d0d824bc9f74 Mon Sep 17 00:00:00 2001 From: Tomas Torsvik Date: Thu, 11 Dec 2025 14:45:30 +0100 Subject: [PATCH 1/2] Namelist options for DMS tuning --- cime_config/namelist_definition_blom.xml | 30 ++++++++++++++++++++++++ hamocc/mo_param_bgc.F90 | 13 +++++----- 2 files changed, 37 insertions(+), 6 deletions(-) diff --git a/cime_config/namelist_definition_blom.xml b/cime_config/namelist_definition_blom.xml index ff94d855..f05b636b 100644 --- a/cime_config/namelist_definition_blom.xml +++ b/cime_config/namelist_definition_blom.xml @@ -6063,6 +6063,36 @@ Half-saturation constant for NO2 for nitrification on NO2 (kmol/m3) + + real + bgcparams + bgcparams + + None + + bacterial removal; constant for DMS scheme (Kloster et al. (2006), Table 1) + + + + real + bgcparams + bgcparams + + None + + production with delcar; constant for DMS scheme (Kloster et al. (2006), Table 1) + + + + real + bgcparams + bgcparams + + None + + production with delsil; constant for DMS scheme (Kloster et al. (2006), Table 1) + + real bgcparams diff --git a/hamocc/mo_param_bgc.F90 b/hamocc/mo_param_bgc.F90 index f4411554..fb364c65 100644 --- a/hamocc/mo_param_bgc.F90 +++ b/hamocc/mo_param_bgc.F90 @@ -496,9 +496,9 @@ module mo_param_bgc ! Set constants for dms scheme following Kloster et al. (2006), Table 1 real(rp), protected :: dmsp1 = 10._rp ! 2*5. production with temp real(rp), protected :: dmsp2 = 0.0011_rp - real(rp), protected :: dmsp3 = 0.0864_rp ! bacterial removal, but reduced 50% to increase dms emissions + real(rp), protected :: dmsp3 = 0.1296_rp ! bacterial removal, but reduced 50% to increase dms emissions real(rp), protected :: dmsp4 = 1.25_rp*0.10_rp ! production with delcar, but reduced by ~7% - real(rp), protected :: dmsp5 = 1.25_rp*0.02_rp ! production with delsil, but increased by a factor of ~2 + real(rp), protected :: dmsp5 = 1.36e-02_rp ! production with delsil real(rp), protected :: dmsp6 = 0.1e-07_rp ! half saturation microbial ! Scaling factor for pH dependency (used if with_dmsph=.true.) @@ -724,10 +724,11 @@ subroutine read_bgcnamelist() bkoxdnra_sed,bkdnra_sed,q10anh4nitr_sed, & bkoxamox_sed,bkanh4nitr_sed,q10ano2nitr_sed, & bkoxnitr_sed,bkano2nitr_sed,sed_alpha_poc,sed_qual_sc, & - sed_denit,sed_sulf, & - sed_O2thresh_hypoxic,sed_O2thresh_sulf,sed_NO3thresh_sulf,& - gammapsl,gammazsl,alphasl,alphasr,docl_remin,docsl_remin, & - docsr_remin,docr_remin,yield_n2o_inf,bkamoxn2o + sed_denit,sed_sulf,sed_O2thresh_hypoxic, & + sed_O2thresh_sulf,sed_NO3thresh_sulf,gammapsl, & + gammazsl,alphasl,alphasr,docl_remin,docsl_remin, & + docsr_remin,docr_remin,yield_n2o_inf,bkamoxn2o, & + dmsp3,dmsp4,dmsp5 if (mnproc.eq.1) then write(io_stdo_bgc,*) From 576caa77fc77db594a2e97faade0e9b13abdcffc Mon Sep 17 00:00:00 2001 From: Tomas Torsvik Date: Thu, 11 Dec 2025 15:22:06 +0100 Subject: [PATCH 2/2] Revert default values, to be updated later --- hamocc/mo_param_bgc.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/hamocc/mo_param_bgc.F90 b/hamocc/mo_param_bgc.F90 index fb364c65..b72e58fb 100644 --- a/hamocc/mo_param_bgc.F90 +++ b/hamocc/mo_param_bgc.F90 @@ -496,9 +496,9 @@ module mo_param_bgc ! Set constants for dms scheme following Kloster et al. (2006), Table 1 real(rp), protected :: dmsp1 = 10._rp ! 2*5. production with temp real(rp), protected :: dmsp2 = 0.0011_rp - real(rp), protected :: dmsp3 = 0.1296_rp ! bacterial removal, but reduced 50% to increase dms emissions + real(rp), protected :: dmsp3 = 0.0864_rp ! bacterial removal, but reduced 50% to increase dms emissions real(rp), protected :: dmsp4 = 1.25_rp*0.10_rp ! production with delcar, but reduced by ~7% - real(rp), protected :: dmsp5 = 1.36e-02_rp ! production with delsil + real(rp), protected :: dmsp5 = 1.25_rp*0.02_rp ! production with delsil, but increased by a factor of ~2 real(rp), protected :: dmsp6 = 0.1e-07_rp ! half saturation microbial ! Scaling factor for pH dependency (used if with_dmsph=.true.)