For equalGM, change taper to GM too

E3SM-Project · May 25, 2022 · 142ee46 · 142ee46
1 parent 4cb9c1b
commit 142ee46
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/components/mpas-ocean/src/shared/mpas_ocn_gm.F b/components/mpas-ocean/src/shared/mpas_ocn_gm.F
@@ -770,7 +770,7 @@ subroutine ocn_GM_compute_Bolus_velocity(statePool, &
             !$omp parallel
             !$omp do schedule(runtime)
             do iEdge=1, nEdges
-               RediKappa(iEdge) = RediHorizontalTaper(iEdge) * gmBolusKappa(iEdge)
+               RediKappa(iEdge) = gmHorizontalTaper(iEdge) * gmBolusKappa(iEdge)
             end do
             !$omp end do
             !$omp end parallel

diff --git a/components/mpas-ocean/src/shared/mpas_ocn_tracer_hmix_redi.F b/components/mpas-ocean/src/shared/mpas_ocn_tracer_hmix_redi.F
@@ -538,6 +538,8 @@ subroutine ocn_tracer_hmix_Redi_init(domain, err)!{{{
          ! initialize Redi kappa array
          if (config_Redi_closure == 'constant'.or. &
              config_Redi_closure == 'equalGM') then
+             ! For 'equalGM' RediKappa is updated every step. Just
+             ! initialize here.
             !$omp parallel
             !$omp do schedule(runtime)
             do iEdge = 1, nEdges