Juan 29/03/2019: check very small domain size(=0 ) & force N/S/E/W check on getting LB bounds

src/LIB/SURCOUCHE/src/mode_distriblb.f90

@@ -103,12 +103,14 @@ END SUBROUTINE GET_DISTRIB_LB
 !!    MODIFICATIONS
 !!    -------------
 !!      Original    23/09/98
+!!      Modif 
+!!     J.Escobar 28/03/2019: for very small domain , force N/S/E/W check on getting LB bounds
 !-------------------------------------------------------------------------------
 !
 USE MODD_PARAMETERS_ll,ONLY : JPHEXT
 USE MODD_VAR_ll,       ONLY : TCRRT_PROCONF
 USE MODD_STRUCTURE_ll,  ONLY : MODELSPLITTING_ll
-USE MODE_TOOLS_ll,     ONLY : GET_INTERSECTION_ll,GET_GLOBALDIMS_ll,LWEST_ll
+USE MODE_TOOLS_ll,     ONLY : GET_INTERSECTION_ll,GET_GLOBALDIMS_ll,LWEST_ll,LEAST_ll
 
 !*       0.    DECLARATIONS
 !              ------------ 
@@ -191,7 +193,7 @@ ELSE
   IYEND=IJMAX_ll+ 2 * JPHEXT
 ENDIF
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,YMODE,IINFO,KIP)
-IF (IINFO /= 1) THEN  ! no empty intersection
+IF (IINFO/=1 .AND. LWEST_ll(KIP) ) THEN  ! no empty intersection
   IF (HCOORD == 'LOC') THEN
     KIB=IXORI
     KIE=IXENDI
@@ -242,7 +244,7 @@ ELSE
   IYEND=IJMAX_ll+ 2 * JPHEXT
 ENDIF
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,YMODE,IINFO,KIP)
-IF (IINFO /=1) THEN
+IF (IINFO/=1  .AND. LEAST_ll(KIP) ) THEN
   IF (HCOORD == 'LOC') THEN
     IF (KIB == 0) KIB=1+KIE
     KIE=KIE+1+IXENDI-IXORI
@@ -317,7 +319,7 @@ END SUBROUTINE GET_DISTRIBX_LB
 USE MODD_PARAMETERS_ll,ONLY : JPHEXT
 USE MODD_VAR_ll,       ONLY : TCRRT_PROCONF
 USE MODD_STRUCTURE_ll,  ONLY : MODELSPLITTING_ll
-USE MODE_TOOLS_ll,     ONLY : GET_INTERSECTION_ll,GET_GLOBALDIMS_ll,LSOUTH_ll
+USE MODE_TOOLS_ll,     ONLY : GET_INTERSECTION_ll,GET_GLOBALDIMS_ll,LSOUTH_ll,LNORTH_ll
 !*       0.    DECLARATIONS
 !              ------------ 
 !*       0.1   declarations of arguments
@@ -396,7 +398,7 @@ ELSE
   IYEND=JPHEXT !1
 ENDIF
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,YMODE,IINFO,KIP)
-IF (IINFO /= 1) THEN  ! no empty intersection
+IF (IINFO/=1 .AND. LSOUTH_ll(KIP) ) THEN  ! no empty intersection
   IF (HCOORD == 'LOC') THEN
     KIB=IXORI
     KIE=IXENDI
@@ -442,7 +444,7 @@ ELSE
   IYEND=IJMAX_ll + 2 * JPHEXT - JPHEXT + JPHEXT  ! + 2 * JPHEXT 
 ENDIF
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,YMODE,IINFO,KIP)
-IF (IINFO /=1) THEN
+IF (IINFO/=1 .AND. LNORTH_ll(KIP) ) THEN
   IF (HCOORD == 'LOC') THEN
     KIB=IXORI
     KIE=IXENDI

src/LIB/SURCOUCHE/src/mode_toolsz_ll.f90

@@ -202,6 +202,7 @@
     !!    -------------
     !     Original 01/05/98
     !     R.Guivarch 29/11/99 : x and y splitting : HSPLITTING
+    !     J.Escobar 28/03/2019: check very small domain(0 size) 
     !
     !-------------------------------------------------------------------------------
     !
@@ -332,6 +333,14 @@
        CALL CARTESIANZ(TPROC,NB_PROC,X_DIM,Y_DIM,Z_DIM,X_DOMAINS,Y_DOMAINS,Z_DOMAINS,321)
        !
     END IF
+    IF (      ( (1+TPROC(IP)%NXEND-TPROC(IP)%NXOR) == 0 )  &
+         .OR. ( (1+TPROC(IP)%NYEND-TPROC(IP)%NYOR) == 0 )  ) THEN
+       PRINT*, "/!\ SPLITZ: some proc have 0 size local domaine , to much processors used for domaine size ", &
+            " IP="  ,IP , &
+            " DIMX=",1+TPROC(IP)%NXEND-TPROC(IP)%NXOR, &
+            " DIMY=",1+TPROC(IP)%NYEND-TPROC(IP)%NYOR
+       CALL ABORT()
+    END IF
     !
     !*       3.     shift from physical to extended domain
     !
@@ -352,6 +361,7 @@
                 PRINT*,"NYOR=",TPROC(IK)%NYOR," NYEND=",TPROC(IK)%NYEND," TAILLE=",1+TPROC(IK)%NYEND-TPROC(IK)%NYOR
                 PRINT*,"NZOR=",TPROC(IK)%NZOR," NZEND=",TPROC(IK)%NZEND," TAILLE=",1+TPROC(IK)%NZEND-TPROC(IK)%NZOR
              END DO
+
           ENDIF
        END IF
     END IF

src/MNH/get_sizex_lb.f90

@@ -83,6 +83,7 @@ END MODULE MODI_GET_SIZEX_LB
 !!      Original    23/09/98
 !!   J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1 
 !!  Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
+!!   J.Escobar 28/03/2019: for very small domain , force N/S/E/W check on getting LB bounds
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -139,7 +140,7 @@ IF (KRIMX /=0) THEN
   IYOR=1 
   IYEND=KJMAX_ll+ 2 * JPHEXT
   CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-  IF (IINFO/=1) THEN  ! no empty intersection
+  IF (IINFO/=1 .AND. LWEST_ll() ) THEN  ! no empty intersection
     KISIZEXF=KISIZEXF + (IXENDI - IXORI +1)
     KJSIZEXF= IYENDI - IYORI +1 
   ENDIF
@@ -150,7 +151,7 @@ IF (KRIMX /=0) THEN
   IYOR=1
   IYEND=KJMAX_ll+ 2 * JPHEXT
   CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-  IF (IINFO/=1) THEN ! no empty intersection
+  IF (IINFO/=1 .AND. LEAST_ll() ) THEN ! no empty intersection
     KISIZEXF=KISIZEXF + (IXENDI - IXORI +1) ! added to the western side 
     KJSIZEXF= IYENDI - IYORI +1 
     KISIZEXFU=KISIZEXFU + ( IXENDI - IXORI +1)
@@ -162,7 +163,7 @@ IF (KRIMX /=0) THEN
   IYOR=1 
   IYEND=KJMAX_ll+ 2 * JPHEXT
   CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-  IF (IINFO /= 1) THEN ! no empty intersection
+  IF (IINFO/=1 .AND. LWEST_ll() ) THEN ! no empty intersection
     KISIZEXFU=KISIZEXFU + (IXENDI - IXORI +1)
     KJSIZEXFU= IYENDI - IYORI +1 
   ENDIF
@@ -176,7 +177,7 @@ IXEND=JPHEXT+2 ! 3
 IYOR=1
 IYEND=KJMAX_ll+ 2 * JPHEXT
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO /=1 ) THEN ! no empty intersection
+IF (IINFO /=1 .AND. LWEST_ll() ) THEN ! no empty intersection
   KISIZEX4=KISIZEX4 + ( IXENDI - IXORI +1)
   KJSIZEX4= IYENDI - IYORI +1 
 ENDIF
@@ -186,7 +187,7 @@ IXEND=KIMAX_ll+ 2 * JPHEXT - JPHEXT + JPHEXT  ! + 2*JPHEXT
 IYOR=1
 IYEND=KJMAX_ll+ 2 * JPHEXT
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO/=1) THEN ! no empty intersection
+IF (IINFO/=1 .AND. LEAST_ll() ) THEN ! no empty intersection
   KISIZEX4=KISIZEX4 + (IXENDI - IXORI +1)
   KJSIZEX4= IYENDI - IYORI +1 
 ENDIF
@@ -199,7 +200,7 @@ IXEND=JPHEXT ! 1
 IYOR=1
 IYEND=KJMAX_ll+ 2 * JPHEXT
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO /=1) THEN   ! no empty intersection
+IF (IINFO/=1 .AND. LWEST_ll() ) THEN   ! no empty intersection
   KISIZEX2=KISIZEX2 + ( IXENDI - IXORI +1)
   KJSIZEX2= IYENDI - IYORI +1 
 ENDIF
@@ -209,7 +210,7 @@ IXEND=KIMAX_ll + 2 * JPHEXT - JPHEXT + JPHEXT !  + 2 * JPHEXT
 IYOR=1
 IYEND=KJMAX_ll+ 2 * JPHEXT
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO /=1) THEN    ! no empty intersection
+IF (IINFO/=1 .AND. LEAST_ll() ) THEN    ! no empty intersection
   KISIZEX2=KISIZEX2 + ( IXENDI - IXORI +1)
   KJSIZEX2= IYENDI - IYORI +1 
 ENDIF

src/MNH/get_sizey_lb.f90

@@ -83,6 +83,7 @@ END MODULE MODI_GET_SIZEY_LB
 !!      Original    23/09/98
 !!   J.Escobar : 15/09/2015 : WENO5 & JPHEXT <> 1 
 !!  Philippe Wautelet: 05/2016-04/2018: new data structures and calls for I/O
+!!   J.Escobar 28/03/2019: for very small domain , force N/S/E/W check on getting LB bounds
 !-------------------------------------------------------------------------------
 !
 !*       0.    DECLARATIONS
@@ -139,7 +140,7 @@ IF (KRIMY /=0) THEN
   IYOR=1 
   IYEND=KRIMY+JPHEXT ! +1
   CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-  IF (IINFO/=1) THEN  ! no empty intersection
+  IF (IINFO/=1 .AND. LSOUTH_ll() ) THEN  ! no empty intersection
     KISIZEYF= IXENDI - IXORI +1
     KJSIZEYF= KJSIZEYF + (IYENDI - IYORI +1)
   ENDIF
@@ -150,7 +151,7 @@ IF (KRIMY /=0) THEN
   IYOR=KJMAX_ll + 2 * JPHEXT-KRIMY-JPHEXT+1  ! -KRIMY
   IYEND=KJMAX_ll+ 2 * JPHEXT
   CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-  IF (IINFO/=1) THEN ! no empty intersection
+  IF (IINFO/=1 .AND. LNORTH_ll() ) THEN ! no empty intersection
     KISIZEYF=IXENDI - IXORI +1 
     KJSIZEYF= KJSIZEYF + (IYENDI - IYORI +1 )! added to the southern  side 
     KISIZEYFV= IXENDI - IXORI +1
@@ -162,7 +163,7 @@ IF (KRIMY /=0) THEN
   IYOR=2
   IYEND=KRIMY+JPHEXT+1   !+2
   CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-  IF (IINFO /= 1) THEN ! no empty intersection
+  IF (IINFO/=1 .AND. LSOUTH_ll() ) THEN ! no empty intersection
     KISIZEYFV= IXENDI - IXORI +1
     KJSIZEYFV= KJSIZEYFV + (IYENDI - IYORI +1 )
   ENDIF
@@ -176,7 +177,7 @@ IXEND=KIMAX_ll+ 2 * JPHEXT
 IYOR=2  !2
 IYEND=JPHEXT+2 !3
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO /=1 ) THEN ! no empty intersection
+IF (IINFO /=1  .AND. LSOUTH_ll() ) THEN ! no empty intersection
   KISIZEY4= IXENDI - IXORI +1
   KJSIZEY4= KJSIZEY4 + (IYENDI - IYORI +1)
 ENDIF
@@ -186,7 +187,7 @@ IXEND=KIMAX_ll+ 2 * JPHEXT
 IYOR=KJMAX_ll + 2 * JPHEXT - JPHEXT      ! + JPHEXT
 IYEND=KJMAX_ll+ 2 * JPHEXT - JPHEXT + JPHEXT  ! + 2*JPHEXT
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO/=1) THEN ! no empty intersection
+IF (IINFO/=1 .AND. LNORTH_ll() ) THEN ! no empty intersection
   KISIZEY4=IXENDI - IXORI +1
   KJSIZEY4=KJSIZEY4 + ( IYENDI - IYORI +1 )
 ENDIF
@@ -199,7 +200,7 @@ IXEND=KIMAX_ll+ 2 * JPHEXT
 IYOR=1   ! 1
 IYEND=JPHEXT  ! 1
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO /=1) THEN   ! no empty intersection
+IF (IINFO /=1 .AND. LSOUTH_ll() ) THEN   ! no empty intersection
   KISIZEY2= IXENDI - IXORI +1
   KJSIZEY2=KJSIZEY2 + (IYENDI - IYORI +1 )
 ENDIF
@@ -209,7 +210,7 @@ IXEND=KIMAX_ll + 2 * JPHEXT
 IYOR=KJMAX_ll + 2 * JPHEXT - JPHEXT + 1 !  + 2 * JPHEXT
 IYEND=KJMAX_ll+ 2 * JPHEXT - JPHEXT + JPHEXT !  + 2 * JPHEXT
 CALL GET_INTERSECTION_ll(IXOR,IYOR,IXEND,IYEND,IXORI,IYORI,IXENDI,IYENDI,"EXTE",IINFO)
-IF (IINFO /=1) THEN    ! no empty intersection
+IF (IINFO /=1 .AND. LNORTH_ll() ) THEN    ! no empty intersection
   KISIZEY2=  IXENDI - IXORI +1
   KJSIZEY2= KJSIZEY2 + (IYENDI - IYORI +1 )
 ENDIF