From 2babb173cc37d0d0fa7f4044758c3e416cc07f0e Mon Sep 17 00:00:00 2001
From: Quentin Rodier <quentin.rodier@meteo.fr>
Date: Thu, 3 Jun 2021 18:03:02 +0200
Subject: [PATCH] Quentin 03/06/2021: bugfix add the pressure gradient to the
sources at the border if L2D (removed by IBM contrib, commit fde83373)
---
src/MNH/pressurez.f90 | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/src/MNH/pressurez.f90 b/src/MNH/pressurez.f90
index b16ee01a5..7cf850831 100644
--- a/src/MNH/pressurez.f90
+++ b/src/MNH/pressurez.f90
@@ -563,12 +563,10 @@ END IF
!
IF ( HLBCX(1) /= 'CYCL' ) THEN
IF(LWEST_ll()) THEN
-! ZPHIT(IIB-1,:, :) = ZPHIT(IIB,:, :)
ZPHIT(IIB-1,:,IKB-1) = ZPHIT(IIB,:,IKB)
ZPHIT(IIB-1,:,IKE+1) = ZPHIT(IIB,:,IKE)
ENDIF
IF(LEAST_ll()) THEN
-! ZPHIT(IIE+1,:, :) = ZPHIT(IIE,:, :)
ZPHIT(IIE+1,:,IKB-1) = ZPHIT(IIE,:,IKB)
ZPHIT(IIE+1,:,IKE+1) = ZPHIT(IIE,:,IKE)
ENDIF
@@ -576,17 +574,20 @@ ENDIF
!
IF ( HLBCY(1) /= 'CYCL' ) THEN
IF (LSOUTH_ll()) THEN
-! ZPHIT(:,IJB-1, :) = ZPHIT(:,IJB, :)
ZPHIT(:,IJB-1,IKB-1) = ZPHIT(:,IJB,IKB)
ZPHIT(:,IJB-1,IKE+1) = ZPHIT(:,IJB,IKE)
ENDIF
IF (LNORTH_ll()) THEN
-! ZPHIT(:,IJE+1, :) = ZPHIT(:,IJE, :)
ZPHIT(:,IJE+1,IKB-1) = ZPHIT(:,IJE,IKB)
ZPHIT(:,IJE+1,IKE+1) = ZPHIT(:,IJE,IKE)
ENDIF
ENDIF
!
+IF ( L2D ) THEN
+ IF (LSOUTH_ll()) ZPHIT(:,IJB-1,:) = ZPHIT(:,IJB,:)
+ IF (LNORTH_ll()) ZPHIT(:,IJE+1,:) = ZPHIT(:,IJB,:)
+ENDIF
+!
IF (LIBM) THEN
!
IF ( HLBCX(1) == 'CYCL' ) THEN
--
GitLab