Fix bug and make sure melt potential is always <= 0

NOAA-GFDL · Aug 6, 2018 · d0bc3bc · d0bc3bc
1 parent c6b2534
commit d0bc3bc
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/config_src/mct_driver/ocn_cap_methods.F90 b/config_src/mct_driver/ocn_cap_methods.F90
@@ -183,8 +183,8 @@ subroutine ocn_export(ind, ocn_public, grid, o2x, dt_int, ncouple_per_day)
         else
           ! Melt_potential: change from J/m^2 to W/m^2
           o2x(ind%o2x_Fioo_q, n) = -ocn_public%melt_potential(ig,jg) * grid%mask2dT(i,j) * I_time_int !* ncouple_per_day
-          ! reset melt_potential
-          ocn_public%melt_potential(ig,jg) = 0.0
+          ! make sure Melt_potential is always <= 0
+          if (o2x(ind%o2x_Fioo_q, n) > 0.0) o2x(ind%o2x_Fioo_q, n) = 0.0
         endif
         ! Make a copy of ssh in order to do a halo update. We use the usual MOM domain
         ! in order to update halos. i.e. does not use global indexing.