From b4f557ab3a0b68c9303a38afe3759bb3061af15f Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 24 Nov 2023 18:02:47 +0100
Subject: [PATCH] Didier R. 24/11/2023: fix computation of diag theta_es
---
src/MNH/write_lfifm1_for_diag.f90 | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/src/MNH/write_lfifm1_for_diag.f90 b/src/MNH/write_lfifm1_for_diag.f90
index 84ff78bda..96a3954b2 100644
--- a/src/MNH/write_lfifm1_for_diag.f90
+++ b/src/MNH/write_lfifm1_for_diag.f90
@@ -2454,7 +2454,7 @@ IF (LMOIST_ES .AND. (NRR>0)) THEN
(3.5*ALOG(XTHT(:,:,:)*( XPABST(:,:,:)/XP00 )**(XRD/XCPD) ) &
- ALOG( XPABST(:,:,:)*0.01*ZWORK31(:,:,:) / ( 0.622+ZWORK31(:,:,:) ) ) &
-4.805 ) ) + 55.
- ZTHETAES(:,:,:)= XTHT(:,:,:) * EXP( (3376. / ZTHETAE(:,:,:) - 2.54) &
+ ZTHETAES(:,:,:)= XTHT(:,:,:) * EXP( (3376. / ZTHETAES(:,:,:) - 2.54) &
*ZWORK31(:,:,:) *(1. +0.81 *ZWORK31(:,:,:)) )
TZFIELD = TFIELDMETADATA( &
CMNHNAME = 'THETAES', &
--
GitLab