From 16dd9128b3d3a07447498a1d536e165d597a6318 Mon Sep 17 00:00:00 2001
From: Daley Calvert <dcalve@users.noreply.github.com>
Date: Wed, 15 Mar 2023 10:47:06 +0000
Subject: [PATCH] #175: Correct isotherm depth calculation
---
src/OCE/DIA/diahth.F90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/OCE/DIA/diahth.F90 b/src/OCE/DIA/diahth.F90
index dc71bfa7..f66bf3c9 100644
--- a/src/OCE/DIA/diahth.F90
+++ b/src/OCE/DIA/diahth.F90
@@ -306,7 +306,7 @@ CONTAINS
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) ) &
+ & * ( 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