jprod_lithdet update for COBALTv3 (#41)

* jprod_lithdet update for COBALTv3 * change unit of lithogenic detrital flux and production integral in upper 100m * Update ocean.stats --------- Co-authored-by: Fei Da <[email protected]> Co-authored-by: Fei Da <[email protected]> Co-authored-by: Andrew Ross <[email protected]> Co-authored-by: Yi-Cheng Teng - NOAA GFDL <[email protected]>
NOAA-CEFI-Regional-Ocean-Modeling · Apr 12, 2024 · f7ac0b7 · f7ac0b7
1 parent de7eb09
commit f7ac0b7
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 10 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.2897044161471751E-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.2185516137658286E-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.4391883490880411E-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.9283640532311142E-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.6652617822038532E-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.2253748503107894E-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.7781350394219996E-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.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
diff --git a/generic_tracers/generic_COBALT.F90 b/generic_tracers/generic_COBALT.F90
@@ -2688,7 +2688,7 @@ subroutine generic_COBALT_register_diag(diag_list)
     bact(1)%id_juptake_ldon_100 = register_diag_field(package_name, vardesc_temp%name, axes(1:2),&
          init_time, vardesc_temp%longname,vardesc_temp%units, missing_value = missing_value1)
 
-    vardesc_temp = vardesc("jprod_lithdet_100","Lithogenic detritus production integral in upper 100m",'h','1','s','mol m-2 s-1','f')
+    vardesc_temp = vardesc("jprod_lithdet_100","Lithogenic detritus production integral in upper 100m",'h','1','s','g m-2 s-1','f')
     cobalt%id_jprod_lithdet_100 = register_diag_field(package_name, vardesc_temp%name, axes(1:2),&
          init_time, vardesc_temp%longname,vardesc_temp%units, missing_value = missing_value1)
 
@@ -2920,7 +2920,7 @@ subroutine generic_COBALT_register_diag(diag_list)
     cobalt%id_fcadet_arag_100 = register_diag_field(package_name, vardesc_temp%name, axes(1:2),&
          init_time, vardesc_temp%longname,vardesc_temp%units, missing_value = missing_value1)
 
-    vardesc_temp = vardesc("flithdet_100","Lithogenic detritus sinking flux @ 100m",'h','1','s','mol m-2 s-1','f')
+    vardesc_temp = vardesc("flithdet_100","Lithogenic detritus sinking flux @ 100m",'h','1','s','g m-2 s-1','f')
     cobalt%id_flithdet_100 = register_diag_field(package_name, vardesc_temp%name, axes(1:2),&
          init_time, vardesc_temp%longname,vardesc_temp%units, missing_value = missing_value1)
 
@@ -8074,7 +8074,8 @@ subroutine generic_COBALT_update_from_source(tracer_list,Temp,Salt,rho_dzt,dzt,h
 
     do k = 1, nk ; do j = jsc, jec ; do i = isc, iec   !{
        cobalt%jprod_lithdet(i,j,k)=( cobalt%total_filter_feeding(i,j,k)/ &
-                                   ( phyto(LARGE)%f_n(i,j,k) + phyto(DIAZO)%f_n(i,j,k) + epsln ) * &
+                                   ( phyto(LARGE)%f_n(i,j,k) + phyto(DIAZO)%f_n(i,j,k) + &
+                                     0.8*phyto(MEDIUM)%f_n(i,j,k) + 0.3*phyto(SMALL)%f_n(i,j,k) + epsln) * &
                                     cobalt%phi_lith + cobalt%k_lith ) * cobalt%f_lith(i,j,k)
     enddo; enddo ; enddo !} i,j,k