reduce wind inpit and wave dissipation in presence of sea ice

umwm · Feb 8, 2019 · ed40d41 · ed40d41
1 parent 76ce5af
commit ed40d41
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/umwm_source_functions.F90 b/src/umwm_source_functions.F90
@@ -68,6 +68,7 @@ subroutine sin_d12()
 
  do concurrent (o = 1:om, p = 1:pm, i = istart:iend)
  ssin(o,p,i) = twopi * rhorat(i) * ssin(o,p,i) * fkovg(o,i)
+ ssin(o,p,i) = (1 - fice(i)) * ssin(o,p,i)
  end do
 
  ! prevent negative sin for diagnostic tail
@@ -100,6 +101,7 @@ subroutine sds_d12
 
  do concurrent(o = 1:om, p = 1:pm, i = istart:iend)
  sds(o,p,i) = twopisds_fac * f(o) * dummy(o,p,i) * (e(o,p,i) * k4(o,i))**sds_power
+ sds(o,p,i) = (1 - fice(i)) * sds(o,p,i) 
  end do
 
  end subroutine sds_d12