add WOMBATlite generic tracer module

zero det_sediment and caco3_sediment tracers below k=1
COSIMA · May 17, 2024 · a5e8780 · a5e8780
1 parent 857fcda
commit a5e8780
Show file tree

Hide file tree

Showing 2 changed files with 16 additions and 0 deletions.
diff --git a/MOM6/CMakeLists.txt b/MOM6/CMakeLists.txt
@@ -317,6 +317,7 @@ target_sources(OM3_mom6 PRIVATE
   GFDL_generic_tracers/generic_tracers/generic_miniBLING.F90
   GFDL_generic_tracers/generic_tracers/generic_SF6.F90
   GFDL_generic_tracers/generic_tracers/generic_TOPAZ.F90
+  GFDL_generic_tracers/generic_tracers/generic_WOMBATlite.F90
 
   GFDL_generic_tracers/mocsy/src/mocsy_buffesm.F90
   GFDL_generic_tracers/mocsy/src/mocsy_constants.F90

diff --git a/MOM6/patches/MOM_generic_tracer.F90.patch b/MOM6/patches/MOM_generic_tracer.F90.patch
@@ -0,0 +1,15 @@
+diff --git a/MOM6/src/tracer/MOM_generic_tracer.F90 b/MOM6/src/tracer/MOM_generic_tracer.F90.new
+index 131110e..7334a56 100644
+--- a/MOM6/src/tracer/MOM_generic_tracer.F90
++++ b/MOM6/src/tracer/MOM_generic_tracer.F90.new
+@@ -352,7 +352,9 @@ contains
+           enddo ; enddo ; enddo
+
+           !jgj: Reset CASED to 0 below K=1
+-          if ( (trim(g_tracer_name) == 'cased') .or. (trim(g_tracer_name) == 'ca13csed') ) then
++          ! dts: also WOMBAT sediment tracers
++          if ( (trim(g_tracer_name) == 'cased') .or. (trim(g_tracer_name) == 'ca13csed') .or. &
++            (trim(g_tracer_name) == 'det_sediment') .or. (trim(g_tracer_name) == 'caco3_sediment')) then
+             do k=2,nk ; do j=jsc,jec ; do i=isc,iec
+               if (tr_ptr(i,j,k) /= CS%tracer_land_val) then
+                 tr_ptr(i,j,k) = 0.0