From 74f2451675ef715bba59cfea95933723450a1868 Mon Sep 17 00:00:00 2001
From: Juan Escobar <escj@aero.obs-mip.fr>
Date: Wed, 26 Apr 2017 17:40:55 +0200
Subject: [PATCH] Juan 26/04/2017:ppm.f90 , in DIF2X/Y , initialize realistic
 value in all HALO pts

---
 src/MNH/ppm.f90 | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/src/MNH/ppm.f90 b/src/MNH/ppm.f90
index 7f52edc77..22c8f6a78 100644
--- a/src/MNH/ppm.f90
+++ b/src/MNH/ppm.f90
@@ -582,6 +582,7 @@ CONTAINS
 !!    
 !!    18.3.2006.  T. Maric - original version
 !!    07/2010     J.Escobar : Correction for reproducility
+!!    04/2017     J.Escobar : initialize realistic value in all HALO pts
 !-------------------------------------------------------------------------------
 !
 !
@@ -604,7 +605,9 @@ INTEGER :: IIE,IJE        ! End useful area in x,y directions
 !*       1.0.     COMPUTE THE DOMAIN DIMENSIONS
 !                 -----------------------------
 !
-CALL GET_INDICE_ll(IIB,IJB,IIE,IJE)
+!!$CALL GET_INDICE_ll(IIB,IJB,IIE,IJE)
+IIB=2 ; IIE = SIZE(PQ,1) -1
+IJB=2 ; IJE = SIZE(PQ,2) -1
 !
 !-------------------------------------------------------------------------------
 !
@@ -1008,6 +1011,7 @@ CONTAINS
 !!    
 !!    18.3.2006.  T. Maric - original version, works only for periodic boundary
 !!                           conditions and on one domain
+!!    04/2017     J.Escobar : initialize realistic value in all HALO pts
 !!
 !-------------------------------------------------------------------------------
 !
@@ -1031,7 +1035,9 @@ INTEGER :: IIE,IJE        ! End useful area in x,y directions
 !*       1.0.     COMPUTE THE DOMAIN DIMENSIONS
 !                 -----------------------------
 !
-CALL GET_INDICE_ll(IIB,IJB,IIE,IJE)
+!!$CALL GET_INDICE_ll(IIB,IJB,IIE,IJE)
+IIB=2 ; IIE = SIZE(PQ,1) -1
+IJB=2 ; IJE = SIZE(PQ,2) -1
 !
 !-------------------------------------------------------------------------------
 !
-- 
GitLab