From 4e72eaa44304f0b5f9379dc64519a9235393f6a8 Mon Sep 17 00:00:00 2001 From: Alexander Hager Date: Wed, 19 Jul 2023 14:56:46 -0600 Subject: [PATCH] Zero diffusivity at boundary edges Explicitly sets diffusivity to zero at boundary edges, instead of just waterFluxDiff --- .../src/mode_forward/mpas_li_subglacial_hydro.F | 1 + 1 file changed, 1 insertion(+) diff --git a/components/mpas-albany-landice/src/mode_forward/mpas_li_subglacial_hydro.F b/components/mpas-albany-landice/src/mode_forward/mpas_li_subglacial_hydro.F index f73177ab3907..5a4049f146e3 100644 --- a/components/mpas-albany-landice/src/mode_forward/mpas_li_subglacial_hydro.F +++ b/components/mpas-albany-landice/src/mode_forward/mpas_li_subglacial_hydro.F @@ -1018,6 +1018,7 @@ subroutine calc_edge_quantities(block, err) endif end do where (waterFluxMask == 2) + diffusivity = 0.0_RKIND waterFluxAdvec = 0.0_RKIND waterFluxDiffu = 0.0_RKIND waterVelocity = 0.0_RKIND