From 47a5deb345aa0e93f272982235a00ccb0c9a45fa Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Fri, 24 Feb 2023 13:16:13 +0100
Subject: [PATCH] Quentin 24/02/2023: bugfix, no apply remove_indices on
gradient_m
---
src/MNH/gradient_m.f90 | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/src/MNH/gradient_m.f90 b/src/MNH/gradient_m.f90
index 1c4660151..feb4c58f3 100644
--- a/src/MNH/gradient_m.f90
+++ b/src/MNH/gradient_m.f90
@@ -742,8 +742,8 @@ IKT=SIZE(PY,3)
IKTB=1+JPVEXT_TURB
IKTE=IKT-JPVEXT_TURB
-PGZ_M_W(:,:,:) = (PY(:,:,:)-PY(:,:,IKTB-KL:IKTE-KL)) &
- / PDZZ(:,:,:)
+PGZ_M_W(:,:,IKTB:IKTE) = (PY(:,:,IKTB:IKTE)-PY(:,:,IKTB-KL:IKTE-KL)) &
+ / PDZZ(:,:,IKTB:IKTE)
PGZ_M_W(:,:,KKU)= (PY(:,:,KKU)-PY(:,:,KKU-KL)) &
/ PDZZ(:,:,KKU)
PGZ_M_W(:,:,KKA)= PGZ_M_W(:,:,KKU) ! -999.
--
GitLab