From 93f6eda1e000b9abc86ca48cd427e1ae1684e9aa Mon Sep 17 00:00:00 2001
From: Isabella Ascione <isabella.ascione@metoffice.gov.uk>
Date: Mon, 2 Sep 2024 14:00:26 +0000
Subject: [PATCH] Resolve "bug fixes in GOSI10 4.2.2 icethd_pnd"
---
src/ICE/icethd_pnd.F90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/ICE/icethd_pnd.F90 b/src/ICE/icethd_pnd.F90
index 825c4345..f268859a 100644
--- a/src/ICE/icethd_pnd.F90
+++ b/src/ICE/icethd_pnd.F90
@@ -686,8 +686,8 @@ CONTAINS
! differential growth of base of surface floating ice layer
zdTice = MAX( - ( t_su(ji,jj,jl) - zTd ) , 0._wp ) ! > 0
zomega = rcnd_i * zdTice / zrhoi_L
- zdHui = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_i(ji,jj,jl) )**2 ) &
- - v_il(ji,jj,jl) / a_i(ji,jj,jl)
+ zdHui = SQRT( 2._wp * zomega * rDt_ice + ( v_il(ji,jj,jl) / a_ip(ji,jj,jl) )**2 ) &
+ - v_il(ji,jj,jl) / a_ip(ji,jj,jl)
zdvice = min( zdHui*a_ip(ji,jj,jl) , v_ip(ji,jj,jl) )
IF ( zdvice > epsi10 ) THEN
--
GitLab