WOMBAT: Made the minimum of sediment dissolution configurable at run time with new parameter "omegamax_sed"

Author: pearseb

generic_tracers/generic_WOMBATlite.F90

@@ -185,6 +185,7 @@ module generic_WOMBATlite
         wcaco3, &
         caco3lrem, &
         caco3lrem_sed, &
+        omegamax_sed, &
         f_inorg, &
         disscal, &
         dissara, &
@@ -1472,6 +1473,12 @@ subroutine user_add_params
     !-----------------------------------------------------------------------
     call g_tracer_add_param('caco3lrem_sed', wombat%caco3lrem_sed, 0.01/86400.0)
 
+    ! Ceiling of omega in the sediments (controls rate of CaCO3 dissolution) [0-1]
+    ! - if == 1.0, then there may be at minimum no dissolution of CaCO3
+    ! - if < 1.0, then there is always some dissolution of CaCO3 when when supersaturated
+    !-----------------------------------------------------------------------
+    call g_tracer_add_param('omegamax_sed', wombat%omegamax_sed, 0.7)
+
     ! CaCO3 inorganic fraction [1]
     !-----------------------------------------------------------------------
     call g_tracer_add_param('f_inorg', wombat%f_inorg, 0.04)
@@ -3300,11 +3307,11 @@ subroutine generic_WOMBATlite_update_from_source(tracer_list, Temp, Salt,  &
       wombat%det_sed_remin(i,j) = wombat%detlrem_sed * fbc * wombat%p_det_sediment(i,j,1) ! [mol/m2/s]
       wombat%detfe_sed_remin(i,j) = wombat%detlrem_sed * fbc * wombat%p_detfe_sediment(i,j,1) ! [mol/m2/s]
       if (do_caco3_dynamics) then
-        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) * &
-                                            max(0.1, (1.0 - wombat%sedomega_cal(i,j)))**(4.5)
+        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) &
+                                      * max((1.0-wombat%omegamax_sed), (1.0-wombat%sedomega_cal(i,j)))**(4.5)
       else
-        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) * &
-                                            max(0.1, (1.0 - 0.2081))**(4.5)
+        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) &
+                                      * (1.0 - 0.2081)**(4.5)
       endif
 
       ! Remineralisation of sediments to supply nutrient fields.

generic_tracers/generic_WOMBATmid.F90

@@ -222,6 +222,7 @@ module generic_WOMBATmid
         wcaco3, &
         caco3lrem, &
         caco3lrem_sed, &
+        omegamax_sed, &
         f_inorg, &
         disscal, &
         dissara, &
@@ -2085,6 +2086,12 @@ subroutine user_add_params
     !-----------------------------------------------------------------------
     call g_tracer_add_param('caco3lrem_sed', wombat%caco3lrem_sed, 0.01/86400.0)
 
+    ! Ceiling of omega in the sediments (controls rate of CaCO3 dissolution) [0-1]
+    ! - if == 1.0, then there may be at minimum no dissolution of CaCO3
+    ! - if < 1.0, then there is always some dissolution of CaCO3 when when supersaturated
+    !-----------------------------------------------------------------------
+    call g_tracer_add_param('omegamax_sed', wombat%omegamax_sed, 0.7)
+
     ! CaCO3 inorganic fraction [1]
     !-----------------------------------------------------------------------
     call g_tracer_add_param('f_inorg', wombat%f_inorg, 0.04)
@@ -3453,7 +3460,6 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt,  &
       ! 1. Light limitation of chlorophyll production
       ! 2. minimum and optimal rates of chlorophyll growth
       ! 3. Calculate mg Chl m-3 s-1
-    
       !!!~~~ Phytoplankton ~~~!!!
       pchl_pisl = phy_pisl / ( wombat%phy_mumax(i,j,k) * 86400.0 * & 
                   (1. - min(wombat%phy_lnit(i,j,k), wombat%phy_lfer(i,j,k))) + epsi )
@@ -4524,11 +4530,11 @@ subroutine generic_WOMBATmid_update_from_source(tracer_list, Temp, Salt,  &
       wombat%det_sed_remin(i,j) = wombat%detlrem_sed * fbc * wombat%p_det_sediment(i,j,1) ! [mol/m2/s]
       wombat%detfe_sed_remin(i,j) = wombat%detlrem_sed * fbc * wombat%p_detfe_sediment(i,j,1) ! [mol/m2/s]
       if (do_caco3_dynamics) then
-        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) * &
-                                            max(0.1, (1.0 - wombat%sedomega_cal(i,j)))**(4.5)
+        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) &
+                                      * max((1.0 - wombat%omegamax_sed), (1.0 - wombat%sedomega_cal(i,j)))**(4.5)
       else
-        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) * &
-                                            max(0.1, (1.0 - 0.2081))**(4.5)
+        wombat%caco3_sed_remin(i,j) = wombat%caco3lrem_sed * fbc * wombat%p_caco3_sediment(i,j,1) &
+                                      * (1.0 - 0.2081)**(4.5)
       endif
       if (do_open_n_cycle) then
         ! sedimentary denitrification (Bohlen et al., 2012 Global Biogeochemical Cycles)