Refactor of vertical reconstruction adding six new schemes

.testing/Makefile

@@ -269,16 +269,16 @@ $(BUILD)/timing/Makefile: MOM_ACFLAGS += --with-driver=timing_tests
 # Build executables
 .NOTPARALLEL:$(foreach e,$(UNIT_EXECS),$(BUILD)/unit/$(e))
 $(BUILD)/unit/test_%: $(BUILD)/unit/Makefile FORCE
-	cd $(@D) && $(TIME) $(MAKE) $(@F) -j
+	cd $(@D) && $(TIME) $(MAKE) $(@F)
 $(BUILD)/unit/Makefile: $(foreach e,$(UNIT_EXECS),../config_src/drivers/unit_tests/$(e).F90)
 
 .NOTPARALLEL:$(foreach e,$(TIMING_EXECS),$(BUILD)/timing/$(e))
 $(BUILD)/timing/time_%: $(BUILD)/timing/Makefile FORCE
-	cd $(@D) && $(TIME) $(MAKE) $(@F) -j
+	cd $(@D) && $(TIME) $(MAKE) $(@F)
 $(BUILD)/timing/Makefile: $(foreach e,$(TIMING_EXECS),../config_src/drivers/timing_tests/$(e).F90)
 
 $(BUILD)/%/MOM6: $(BUILD)/%/Makefile FORCE
-	cd $(@D) && $(TIME) $(MAKE) $(@F) -j
+	cd $(@D) && $(TIME) $(MAKE) $(@F)
 
 # Target codebase should use its own build system
 $(BUILD)/target/MOM6: $(BUILD)/target FORCE | $(TARGET_CODEBASE)

config_src/drivers/timing_tests/time_MOM_remapping.F90

@@ -9,8 +9,30 @@ program time_MOM_remapping
 implicit none
 
 type(remapping_CS) :: CS
-integer, parameter :: nk=75, nij=20*20, nits=10, nsamp=100, nschemes = 2
-character(len=10) :: scheme_labels(nschemes)
+integer, parameter :: nk=75, nij=20*20, nits=10, nsamp=100, nschemes = 22
+character(len=16) :: scheme_labels(nschemes) = [ character(len=16) :: &
+      'PCM', &
+      'C_PCM', &
+      'PLM', &
+      'C_MPLM_WA', &
+      'C_EMPLM_WA', &
+      'C_PLM_HYBGEN', &
+      'C_PLM_CW', &
+      'C_PLM_CWK', &
+      'C_MPLM_WA_POLY', &
+      'C_EMPLM_WA_POLY', &
+      'C_MPLM_CWK', &
+      'PPM_H4', &
+      'PPM_IH4', &
+      'PQM_IH4IH3', &
+      'PPM_CW', &
+      'PPM_HYBGEN', &
+      'C_PPM_H4_2018', &
+      'C_PPM_H4_2019', &
+      'C_PPM_HYBGEN', &
+      'C_PPM_CW', &
+      'C_PPM_CWK', &
+      'C_EPPM_CWK' ]
 real, dimension(nschemes) :: timings ! Time for nits of nij calls for each scheme [s]
 real, dimension(nschemes) :: tmean ! Mean time for a call [s]
 real, dimension(nschemes) :: tstd ! Standard deviation of time for a call [s]
@@ -30,9 +52,6 @@ program time_MOM_remapping
 seed(:) = 102030405
 call random_seed(put=seed)
 
-scheme_labels(1) = 'PCM'
-scheme_labels(2) = 'PLM'
-
 ! Set up some test data (note: using k,i indexing rather than i,k)
 allocate( u0(nk,nij), h0(nk,nij), u1(nk,nij), h1(nk,nij) )
 call random_number(u0) ! In range 0-1
@@ -59,7 +78,8 @@ program time_MOM_remapping
 do isamp = 1, nsamp
   ! Time reconstruction + remapping
   do ischeme = 1, nschemes
-    call initialize_remapping(CS, remapping_scheme=trim(scheme_labels(ischeme)))
+    call initialize_remapping(CS, remapping_scheme=trim(scheme_labels(ischeme)), nk=nk, &
+                              h_neglect=1.e-10)
     call cpu_time(start)
     do iter = 1, nits ! Make many passes to reduce sampling error
       do ij = 1, nij ! Calling nij times to make similar to cost in MOM_ALE()

config_src/drivers/unit_tests/test_MOM_remapping.F90

@@ -2,6 +2,18 @@ program test_MOM_remapping
 
 use MOM_remapping, only : remapping_unit_tests
 
-if (remapping_unit_tests(.true.)) stop 1
+integer :: n !< Number of arguments, or tests
+character(len=12) :: cmd_ln_arg !< Command line argument (if any)
+
+n = command_argument_count()
+
+if (n==1) then
+  call get_command_argument(1, cmd_ln_arg)
+  read(cmd_ln_arg,*) n
+else
+  n = 3000 ! Fallback value if no argument provided
+endif
+
+if (remapping_unit_tests(.true., num_comp_samp=n)) stop 1
 
 end program test_MOM_remapping

docs/discrete_space.rst

@@ -17,4 +17,5 @@ algorithm.
     api/generated/pages/Discrete_Coriolis
     api/generated/pages/Discrete_PG
     api/generated/pages/Energetic_Consistency
+    api/generated/pages/Vertical_Reconstruction
     api/generated/pages/Discrete_OBC

docs/zotero.bib

@@ -2946,3 +2946,17 @@ @article{Young1994
 	pages={1812--1826},
 	year={1994}
 }
+
+@article{van_leer_1977,
+	title = {Towards the ultimate conservative difference scheme. {IV}. {A} new approach to numerical convection},
+	volume = {23},
+	issn = {0021-9991},
+	doi = {10.1016/0021-9991(77)90095-X},
+	number = {3},
+	journal = {Journal of Computational Physics},
+	author = {Van Leer, Bram},
+	month = mar,
+	year = {1977},
+	pages = {276--299},
+}
+

src/ALE/MOM_ALE.F90

@@ -248,14 +248,14 @@ subroutine ALE_init( param_file, GV, US, max_depth, CS)
                  default=default_answer_date, do_not_log=.not.GV%Boussinesq)
   if (.not.GV%Boussinesq) CS%answer_date = max(CS%answer_date, 20230701)
 
-  call initialize_remapping( CS%remapCS, string, &
+  call initialize_remapping( CS%remapCS, string, nk=GV%ke, &
                              boundary_extrapolation=init_boundary_extrap, &
                              check_reconstruction=check_reconstruction, &
                              check_remapping=check_remapping, &
                              force_bounds_in_subcell=force_bounds_in_subcell, &
                              om4_remap_via_sub_cells=om4_remap_via_sub_cells, &
                              answer_date=CS%answer_date)
-  call initialize_remapping( CS%vel_remapCS, vel_string, &
+  call initialize_remapping( CS%vel_remapCS, vel_string, nk=GV%ke, &
                              boundary_extrapolation=init_boundary_extrap, &
                              check_reconstruction=check_reconstruction, &
                              check_remapping=check_remapping, &