Juan 26/04/2017:ppm.f90 , in DIF2X/Y , initialize realistic value in all HALO pts
authorJuan Escobar <escj@aero.obs-mip.fr>
Wed, 26 Apr 2017 15:40:55 +0000 (17:40 +0200)
committerJuan Escobar <escj@aero.obs-mip.fr>
Wed, 26 Apr 2017 15:40:55 +0000 (17:40 +0200)
src/MNH/ppm.f90

index 7f52edc..22c8f6a 100644 (file)
@@ -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
 !
 !-------------------------------------------------------------------------------
 !