set lamult to 1 below ice

NOAA-GFDL · Jul 9, 2021 · 7d85ab2 · 7d85ab2
1 parent 5f756ab
commit 7d85ab2
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/config_src/drivers/nuopc_cap/mom_surface_forcing_nuopc.F90 b/config_src/drivers/nuopc_cap/mom_surface_forcing_nuopc.F90
@@ -879,7 +879,11 @@ subroutine convert_IOB_to_forces(IOB, forces, index_bounds, Time, G, US, CS)
   ! wave to ocean coupling
   if ( associated(IOB%lamult)) then
     do j=js,je; do i=is,ie
-      forces%lamult(i,j) = IOB%lamult(i-i0,j-j0)
+      if (IOB%ice_fraction(i-i0,j-j0) <= 0.05 ) then
+        forces%lamult(i,j) = IOB%lamult(i-i0,j-j0)
+      else
+        forces%lamult(i,j) = 1.0
+      endif
     enddo ; enddo
     call pass_var(forces%lamult, G%domain, halo=1 )
   endif