Bugfix inactive top cell surface restoring

E3SM-Ocean-Discussion · May 2, 2022 · 2236e64 · 2236e64
1 parent 8b87c04
commit 2236e64
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 4 deletions.
diff --git a/components/mpas-ocean/src/shared/Makefile b/components/mpas-ocean/src/shared/Makefile
@@ -187,7 +187,7 @@ mpas_ocn_transport_tests.o: mpas_ocn_config.o
 
 mpas_ocn_effective_density_in_land_ice.o: mpas_ocn_constants.o mpas_ocn_config.o
 
-mpas_ocn_forcing_restoring.o: mpas_ocn_constants.o mpas_ocn_config.o
+mpas_ocn_forcing_restoring.o: mpas_ocn_constants.o mpas_ocn_config.o mpas_ocn_mesh.o
 
 mpas_ocn_tracer_surface_restoring.o: mpas_ocn_constants.o mpas_ocn_config.o mpas_ocn_framework_forcing.o mpas_ocn_diagnostics_variables.o
 

diff --git a/components/mpas-ocean/src/shared/mpas_ocn_tracer_surface_restoring.F b/components/mpas-ocean/src/shared/mpas_ocn_tracer_surface_restoring.F
@@ -30,6 +30,7 @@ module ocn_tracer_surface_restoring
    use ocn_config
    use ocn_framework_forcing
    use ocn_diagnostics_variables
+   use ocn_mesh
 
    implicit none
 
@@ -135,11 +136,11 @@ subroutine ocn_tracer_surface_restoring_compute(groupName, nTracers, nCells, tra
 
       err = 0
 
-      iLevel = 1  ! base surface flux restoring on tracer fields in the top layer
-
       !$omp parallel
-      !$omp do schedule(runtime) private(iTracer)
+      !$omp do schedule(runtime) private(iTracer, iLevel)
       do iCell=1,nCells
+        iLevel = minLevelCell(iCell)  ! base surface flux restoring on tracer fields in the top layer
+
         do iTracer=1,nTracers
 
           ! For monthly salinity restoring, tracersSurfaceRestoringValue contains the zero-mean deltaS