Recalibrate calcite and aragonite detritus flux (#47)

* recalibrate calcite and aragonite detritus flux * Update ocean.stats
NOAA-CEFI-Regional-Ocean-Modeling · Apr 15, 2024 · 864b0db · 864b0db
1 parent f7ac0b7
commit 864b0db
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 9 deletions.
diff --git a/.github/ref/OM4.single_column.COBALT.p4/docker-linux-gnu/ocean.stats b/.github/ref/OM4.single_column.COBALT.p4/docker-linux-gnu/ocean.stats
@@ -2,10 +2,10 @@
             [days]                 [m2 s-2]           [Nondim]       [m]             [kg]         [PSU]      [degC]       [Nondim]        [PSU]        [degC]
      0,  731595.000,     0, En 4.0531850015098383E-27, CFL  0.00000, SL -0.0000E+00, M 6.97073E+15, S 35.2427, T  6.2462, Me  0.00E+00, Se  0.00E+00, Te  0.00E+00
      6,  731595.250,     0, En 2.0223578954766472E-06, CFL  0.00281, SL -3.4291E-06, M 6.97073E+15, S 35.2427, T  6.2460, Me -8.57E-10, Se -2.85E-07, Te -2.66E-06
-    12,  731595.500,     0, En 4.2897044161823591E-06, CFL  0.00294, SL -6.4127E-06, M 6.97073E+15, S 35.2427, T  6.2457, Me -7.46E-10, Se -6.81E-07, Te -6.81E-06
-    18,  731595.750,     0, En 1.2185516288182668E-05, CFL  0.00794, SL -9.2141E-06, M 6.97073E+15, S 35.2427, T  6.2459, Me -7.00E-10, Se -9.89E-07, Te -1.00E-05
-    24,  731596.000,     0, En 2.4391883629687393E-05, CFL  0.00870, SL -1.3774E-05, M 6.97073E+15, S 35.2427, T  6.2457, Me -1.14E-09, Se -1.40E-06, Te -1.42E-05
-    30,  731596.250,     0, En 2.9283640524189851E-05, CFL  0.00724, SL -1.8458E-05, M 6.97073E+15, S 35.2427, T  6.2454, Me -1.17E-09, Se -1.93E-06, Te -1.96E-05
-    36,  731596.500,     0, En 3.6652617806009578E-05, CFL  0.00746, SL -2.3211E-05, M 6.97073E+15, S 35.2427, T  6.2451, Me -1.19E-09, Se -2.37E-06, Te -2.42E-05
-    42,  731596.750,     0, En 4.2253748566466643E-05, CFL  0.00812, SL -2.7653E-05, M 6.97073E+15, S 35.2427, T  6.2454, Me -1.11E-09, Se -2.79E-06, Te -2.86E-05
-    48,  731597.000,     0, En 4.7781350583834568E-05, CFL  0.00941, SL -3.1143E-05, M 6.97073E+15, S 35.2427, T  6.2453, Me -8.73E-10, Se -3.05E-06, Te -3.14E-05
+    12,  731595.500,     0, En 4.2897044161040916E-06, CFL  0.00294, SL -6.4127E-06, M 6.97073E+15, S 35.2427, T  6.2457, Me -7.46E-10, Se -6.81E-07, Te -6.81E-06
+    18,  731595.750,     0, En 1.2185516095946735E-05, CFL  0.00794, SL -9.2141E-06, M 6.97073E+15, S 35.2427, T  6.2459, Me -7.00E-10, Se -9.89E-07, Te -1.00E-05
+    24,  731596.000,     0, En 2.4391883408641453E-05, CFL  0.00870, SL -1.3774E-05, M 6.97073E+15, S 35.2427, T  6.2457, Me -1.14E-09, Se -1.40E-06, Te -1.42E-05
+    30,  731596.250,     0, En 2.9283640487163710E-05, CFL  0.00724, SL -1.8458E-05, M 6.97073E+15, S 35.2427, T  6.2454, Me -1.17E-09, Se -1.93E-06, Te -1.96E-05
+    36,  731596.500,     0, En 3.6652617800868474E-05, CFL  0.00746, SL -2.3211E-05, M 6.97073E+15, S 35.2427, T  6.2451, Me -1.19E-09, Se -2.37E-06, Te -2.42E-05
+    42,  731596.750,     0, En 4.2253748371831278E-05, CFL  0.00812, SL -2.7653E-05, M 6.97073E+15, S 35.2427, T  6.2454, Me -1.11E-09, Se -2.79E-06, Te -2.86E-05
+    48,  731597.000,     0, En 4.7781350107303387E-05, CFL  0.00941, SL -3.1143E-05, M 6.97073E+15, S 35.2427, T  6.2453, Me -8.73E-10, Se -3.05E-06, Te -3.14E-05
diff --git a/generic_tracers/generic_COBALT.F90 b/generic_tracers/generic_COBALT.F90
@@ -4818,8 +4818,8 @@ subroutine user_add_params
     call g_tracer_add_param('si_2_n_static_Md', phyto(MEDIUM)%si_2_n_static, 2.0)           ! mol Si mol N-1
     call g_tracer_add_param('si_2_n_max_Lg', phyto(LARGE)%si_2_n_max, 3.0)                  ! mol Si mol N-1
     call g_tracer_add_param('si_2_n_max_Lg', phyto(MEDIUM)%si_2_n_max, 1.0)                 ! mol Si mol N-1
-    call g_tracer_add_param('ca_2_n_arag', cobalt%ca_2_n_arag, 0.050 * 106.0 / 16.0)        ! mol Ca mol N-1
-    call g_tracer_add_param('ca_2_n_calc', cobalt%ca_2_n_calc, 0.015 * 106.0 / 16.0)        ! mol Ca mol N-1
+    call g_tracer_add_param('ca_2_n_arag', cobalt%ca_2_n_arag, 0.055 * 106.0 / 16.0)        ! mol Ca mol N-1
+    call g_tracer_add_param('ca_2_n_calc', cobalt%ca_2_n_calc, 0.050 * 106.0 / 16.0)        ! mol Ca mol N-1
     call g_tracer_add_param('caco3_sat_max', cobalt%caco3_sat_max,10.0)                     ! dimensionless
     !
     !-----------------------------------------------------------------------