Skip to content

Commit

Permalink
Fixes flux at BLD base
Browse files Browse the repository at this point in the history
vanroekel committed Nov 23, 2024
1 parent 5f15586 commit be5c95a
Showing 1 changed file with 3 additions and 2 deletions.
Original file line number Diff line number Diff line change
@@ -114,7 +114,8 @@ subroutine ocn_brine_rejection_tendency(tend, layerThickness, seaIceSalinityFlux

do iCell = 1, nCells
if (seaIceSalinityFlux(iCell) > 0) then !only compute for salinity into ocean
sis_flux(:) = 0.0_RKIND
!this is set to 1 such that across the BLD we don't have a huge salt flux divergence.
sis_flux(:) = 1.0_RKIND
D = boundaryLayerDepth(iCell)
z = 0.0_RKIND
A = (config_brine_param_n + 1) / D**(config_brine_param_n+1)
@@ -127,7 +128,7 @@ subroutine ocn_brine_rejection_tendency(tend, layerThickness, seaIceSalinityFlux
end do
! need to pick up the last bit that goes over
sis_flux(k) = 1.0_RKIND ! this should finish

do k=minLevelCell(iCell),maxLevelCell(iCell)
! need to flip over so we take k+1 - k
tend(index_salinity, k, iCell) = tend(index_salinity, k, iCell) + &

0 comments on commit be5c95a

Please sign in to comment.