From c1c5580cd2141d0e15871a3970c84afa323b16bd Mon Sep 17 00:00:00 2001
From: Philippe WAUTELET <philippe.wautelet@aero.obs-mip.fr>
Date: Fri, 17 Dec 2021 13:47:25 +0100
Subject: [PATCH] Philippe 17/12/2021: bugfix: inversion of values was
incorrectly done for XMEAN_POVO
---
src/MNH/write_lfifm1_for_diag.f90 | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/src/MNH/write_lfifm1_for_diag.f90 b/src/MNH/write_lfifm1_for_diag.f90
index 60240e9ca..0a16f3078 100644
--- a/src/MNH/write_lfifm1_for_diag.f90
+++ b/src/MNH/write_lfifm1_for_diag.f90
@@ -2689,9 +2689,10 @@ IF ( LMEAN_POVO ) THEN
IWORK1(:,:)=0
ZWORK21(:,:)=0.
IF (XMEAN_POVO(1)>XMEAN_POVO(2)) THEN
- XMEAN_POVO(1) = ZX0D
- XMEAN_POVO(2) = XMEAN_POVO(1)
- ZX0D = XMEAN_POVO(2)
+ !Invert values (smallest must be first)
+ ZX0D = XMEAN_POVO(1)
+ XMEAN_POVO(1) = XMEAN_POVO(2)
+ XMEAN_POVO(2) = ZX0D
END IF
DO JK=IKB,IKE
WHERE((XPABST(:,:,JK)>XMEAN_POVO(1)).AND.(XPABST(:,:,JK)<XMEAN_POVO(2)))
--
GitLab