From 61159ba6aaad20ee39b889832219781142c2bc29 Mon Sep 17 00:00:00 2001 From: apcraig Date: Thu, 8 Feb 2018 07:02:23 +0000 Subject: [PATCH] comment out brine fix for bfb with one test failing --- columnphysics/icepack_brine.F90 | 9 +++++---- 1 file changed, 5 insertions(+), 4 deletions(-) diff --git a/columnphysics/icepack_brine.F90 b/columnphysics/icepack_brine.F90 index 91c646bd1..657b9965a 100644 --- a/columnphysics/icepack_brine.F90 +++ b/columnphysics/icepack_brine.F90 @@ -545,11 +545,12 @@ subroutine update_hbrine (meltb, meltt, & dhrunoff = -dhS_top*aice0 hbrocn = max(c0,hbrocn - dhrunoff) exp_arg = darcy_coeff/bphi_min*dt - if (exp_arg > exp_argmax) then - hbrocn_new = c0 - else +! tcx tcraig avoids underflows but is not bit-for-bit +! if (exp_arg > exp_argmax) then +! hbrocn_new = c0 +! else hbrocn_new = hbrocn*exp(-exp_arg) - endif +! endif hbr = max(hbrmin, h_ocn + hbrocn_new) hbrocn_new = hbr-h_ocn darcy_V = -SIGN((hbrocn-hbrocn_new)/dt*bphi_min, hbrocn)