From af5a47534d339075f2b15659abb1757001e01831 Mon Sep 17 00:00:00 2001
From: Pierre Tulet <pierre.tulet@aero.obs-mip.fr>
Date: Mon, 22 May 2023 17:20:00 +0200
Subject: [PATCH] =?UTF-8?q?P.Tulet=2022/02/23:=20correction=20du=20calcul?=
=?UTF-8?q?=20de=20la=20viscosit=C3=A9=20pour=20SST=20<283.15=20et=20SST?=
=?UTF-8?q?=20>=20303.15.=20=20correction=20of=20the=20viscosity=20calcula?=
=?UTF-8?q?tion=20for=20SST=20<283.15=20and=20SST=20>=20303.15.?=
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit
---
src/SURFEX/coupling_sltn.F90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/SURFEX/coupling_sltn.F90 b/src/SURFEX/coupling_sltn.F90
index 3db0c7b1a..cdb6c10ee 100644
--- a/src/SURFEX/coupling_sltn.F90
+++ b/src/SURFEX/coupling_sltn.F90
@@ -192,7 +192,7 @@ IF ((CEMISPARAM_SLT .eq. "Ova14").OR.(CEMISPARAM_SLT .eq. "OvB21a").OR.(CEMISPAR
! Tableau d'interpolation pour calculer ZNUWATER en fonction de la SST
! Cas ou 0 < SST < 10 C
- WHERE ((PSST(:) >= 273.15).AND.(PSST(:) < 283.15))
+ WHERE (PSST(:) < 283.15)
ZVISCO(:) = ZNUWATER(1) + (PSST(:) - ZWT(1)) * (ZNUWATER(2)-ZNUWATER(1)) / &
(ZWT(2) - ZWT(1))
ENDWHERE
@@ -210,7 +210,7 @@ IF ((CEMISPARAM_SLT .eq. "Ova14").OR.(CEMISPARAM_SLT .eq. "OvB21a").OR.(CEMISPAR
ENDWHERE
! Cas ou 30 < SST < 40 C
- WHERE ((PSST(:) >= 303.15).AND.(PSST(:) < 313.15))
+ WHERE (PSST(:) >= 303.15)
ZVISCO(:) = ZNUWATER(4) + (PSST(:) - ZWT(4)) * (ZNUWATER(5)-ZNUWATER(4)) / &
(ZWT(5) - ZWT(4))
ENDWHERE
--
GitLab