Skip to content

Commit

Permalink
Reset IOF%seaice_melt to zero every timestep just like old net_melt
Browse files Browse the repository at this point in the history
This fixes a bug introduced in my previous two commits, in
which I was unintentionally adding IOF%seaice_melt rate on top of the
previous timesteps result, without realizing that it never gets set
back to zero!
  • Loading branch information
Henri Drake committed Dec 11, 2024
1 parent ec559ba commit 5ecae71
Showing 1 changed file with 5 additions and 0 deletions.
5 changes: 5 additions & 0 deletions src/SIS_slow_thermo.F90
Original file line number Diff line number Diff line change
Expand Up @@ -854,6 +854,11 @@ subroutine SIS2_thermodynamics(IST, dt_slow, CS, OSS, FIA, IOF, G, US, IG)
IOF%lprec_ocn_top(i,j) = IOF%lprec_ocn_top(i,j) + IST%part_size(i,j,k) * FIA%lprec_top(i,j,k)
enddo ; enddo ; enddo

! Reset sea ice melt rate to zero
do j=jsc,jec ; do i=isc,iec
IOF%seaice_melt(i,j) = 0.
enddo ; enddo

! Add fluxes of snow and other properties to the ocean due to recent ridging or drifting events.
if (allocated(IST%snow_to_ocn)) then
!$OMP do
Expand Down

0 comments on commit 5ecae71

Please sign in to comment.