From 5932d0686cb746b61ad108bd567b9f800f16b49d Mon Sep 17 00:00:00 2001
From: Simon Mueller <11-smueller@users.noreply.forge.nemo-ocean.eu>
Date: Fri, 5 May 2023 10:49:08 +0000
Subject: [PATCH] Resolve "Error in `dia_hth_dep` when calculating 26/28 deg
isotherm depth" ('branch_4.2')
---
src/OCE/DIA/diahth.F90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/OCE/DIA/diahth.F90 b/src/OCE/DIA/diahth.F90
index 26c97b7d..cdb657b2 100644
--- a/src/OCE/DIA/diahth.F90
+++ b/src/OCE/DIA/diahth.F90
@@ -326,8 +326,8 @@ CONTAINS
iid = iktem(ji,jj)
IF( iid /= 1 ) THEN
zztmp = gdept(ji,jj,iid ,Kmm) & ! linear interpolation
- & + ( gdept(ji,jj,iid+1,Kmm) - gdept(ji,jj,iid,Kmm) ) &
- & * ( 20.*tmask(ji,jj,iid+1) - ts(ji,jj,iid,jp_tem,Kmm) ) &
+ & + ( gdept(ji,jj,iid+1,Kmm) - gdept(ji,jj,iid,Kmm) ) &
+ & * ( ptem * tmask(ji,jj,iid+1) - ts(ji,jj,iid,jp_tem,Kmm) ) &
& / ( ts(ji,jj,iid+1,jp_tem,Kmm) - ts(ji,jj,iid,jp_tem,Kmm) + (1.-tmask(ji,jj,1)) )
pdept(ji,jj) = MIN( zztmp , zzdep) * tmask(ji,jj,1) ! bound by the ocean depth
ELSE
--
GitLab