Fix Redi resolution taper

E3SM-Project · Apr 21, 2022 · 492d8dc · 492d8dc
1 parent 88fd9c0
commit 492d8dc
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 5 deletions.
diff --git a/components/mpas-ocean/src/shared/mpas_ocn_gm.F b/components/mpas-ocean/src/shared/mpas_ocn_gm.F
@@ -1096,9 +1096,9 @@ subroutine ocn_GM_init(domain, err)!{{{
                else if (dcEdge(iEdge) >= config_GM_horizontal_ramp_max) then
                   gmHorizontalTaper(iEdge) = 1.0_RKIND
                else
-                  gmHorizontalTaper(iEdge) = 1.0_RKIND &
-                     /(config_GM_horizontal_ramp_max - config_GM_horizontal_ramp_min) &
-                     *(dcEdge(iEdge) - config_GM_horizontal_ramp_min)
+                  gmHorizontalTaper(iEdge) = &
+                      (dcEdge(iEdge) - config_GM_horizontal_ramp_min) &
+                     /(config_GM_horizontal_ramp_max - config_GM_horizontal_ramp_min)
                end if
             end do
             !$omp end do

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
@@ -554,8 +554,8 @@ subroutine ocn_tracer_hmix_Redi_init(domain, err)!{{{
                else if (dcEdge(iEdge) >= config_Redi_horizontal_ramp_max) then
                   RediHorizontalTaper(iEdge) = 1.0_RKIND
                else
-                  RediHorizontalTaper(iEdge) = RediHorizontalTaper(iEdge) &
-                     *(dcEdge(iEdge) - config_Redi_horizontal_ramp_min) &
+                  RediHorizontalTaper(iEdge) = &
+                      (dcEdge(iEdge) - config_Redi_horizontal_ramp_min) &
                      /(config_Redi_horizontal_ramp_max - config_Redi_horizontal_ramp_min)
                end if
             end do