Skip to content
Snippets Groups Projects
Commit 9f98c5ba authored by ESCOBAR MUNOZ Juan's avatar ESCOBAR MUNOZ Juan Committed by WAUTELET Philippe
Browse files

Juan 2/05/2016 : add STOP in case of problem with decomposition

parent d1425cd5
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
src/MNH/mode_extrapol.f90
+ 10
− 0
View file @ 9f98c5ba
...@@ -152,6 +152,7 @@ CONTAINS ...@@ -152,6 +152,7 @@ CONTAINS
!! ------------- !! -------------
!! !!
!! Original 18/02/2015 !! Original 18/02/2015
!! J.Escobar 2/05/2016 : add STOP in case of problem with decomposition
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
USE MODD_LBC_n USE MODD_LBC_n
USE MODE_MODELN_HANDLER USE MODE_MODELN_HANDLER
...@@ -224,6 +225,7 @@ CONTAINS ...@@ -224,6 +225,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, case not supported : & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, case not supported : &
& the child grid has to be one point larger or one point smaller in X dim" & the child grid has to be one point larger or one point smaller in X dim"
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO3D'
! IF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler ! IF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler
! PTAB(1,:,:) = 2. * PTAB(2,:,:) - PTAB(3,:,:) ! PTAB(1,:,:) = 2. * PTAB(2,:,:) - PTAB(3,:,:)
! ELSEIF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)=='CYCL' ) THEN ! ELSEIF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)=='CYCL' ) THEN
...@@ -238,6 +240,7 @@ CONTAINS ...@@ -238,6 +240,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, IDIMX_C = ", & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, IDIMX_C = ", &
IDIMX_C, ", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT IDIMX_C, ", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO3D'
ENDIF ENDIF
ENDIF ENDIF
IF ( IDIMY_C > IJE - IJB + 1 + 2*JPHEXT ) THEN IF ( IDIMY_C > IJE - IJB + 1 + 2*JPHEXT ) THEN
...@@ -252,6 +255,7 @@ CONTAINS ...@@ -252,6 +255,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, case not supported : & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, case not supported : &
& the child grid has to be one point larger or one point smaller in Y dim" & the child grid has to be one point larger or one point smaller in Y dim"
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO3D'
! IF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler ! IF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler
! PTAB(:,1,:) = 2. * PTAB(:,2,:) - PTAB(:,3,:) ! PTAB(:,1,:) = 2. * PTAB(:,2,:) - PTAB(:,3,:)
! ELSEIF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)=='CYCL' ) THEN ! ELSEIF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)=='CYCL' ) THEN
...@@ -266,6 +270,7 @@ CONTAINS ...@@ -266,6 +270,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, IDIMY_C = ", & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO3D, IDIMY_C = ", &
IDIMY_C, ", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT IDIMY_C, ", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO3D'
ENDIF ENDIF
ENDIF ENDIF
! !
...@@ -289,6 +294,7 @@ CONTAINS ...@@ -289,6 +294,7 @@ CONTAINS
!! ------------- !! -------------
!! !!
!! Original 18/02/2015 !! Original 18/02/2015
!! J.Escobar 2/05/2016 : add STOP in case of problem with decomposition
!------------------------------------------------------------------------------- !-------------------------------------------------------------------------------
USE MODD_LBC_n USE MODD_LBC_n
USE MODE_MODELN_HANDLER USE MODE_MODELN_HANDLER
...@@ -358,6 +364,7 @@ CONTAINS ...@@ -358,6 +364,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, case not supported : & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, case not supported : &
& the child grid has to be one point larger or one point smaller in X dim" & the child grid has to be one point larger or one point smaller in X dim"
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO2D'
! IF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler ! IF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler
! PTAB(1,:) = 2. * PTAB(2,:) - PTAB(3,:) ! PTAB(1,:) = 2. * PTAB(2,:) - PTAB(3,:)
! ELSEIF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)=='CYCL' ) THEN ! ELSEIF ( IIB>1 .AND. LWEST_ll() .AND. CLBCX(1)=='CYCL' ) THEN
...@@ -372,6 +379,7 @@ CONTAINS ...@@ -372,6 +379,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, IDIMX_C = ", IDIMX_C, & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, IDIMX_C = ", IDIMX_C, &
", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT ", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO2D'
ENDIF ENDIF
ENDIF ENDIF
IF ( IDIMY_C > IJE - IJB + 1 + 2*JPHEXT ) THEN IF ( IDIMY_C > IJE - IJB + 1 + 2*JPHEXT ) THEN
...@@ -390,6 +398,7 @@ CONTAINS ...@@ -390,6 +398,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, case not supported : & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, case not supported : &
& the child grid has to be one point larger or one point smaller in Y dim" & the child grid has to be one point larger or one point smaller in Y dim"
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO2D'
! IF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler ! IF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)/='CYCL' ) THEN !du cote ouest, on a un point dans le 'pseudo halo' a extrapoler
! PTAB(:,1) = 2. * PTAB(:,2) - PTAB(:,3) ! PTAB(:,1) = 2. * PTAB(:,2) - PTAB(:,3)
! ELSEIF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)=='CYCL' ) THEN ! ELSEIF ( IJB>1 .AND. LNORTH_ll() .AND. CLBCY(1)=='CYCL' ) THEN
...@@ -404,6 +413,7 @@ CONTAINS ...@@ -404,6 +413,7 @@ CONTAINS
WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, IDIMY_C = ", IDIMY_C, & WRITE(*,*) "ERROR in EXTRAPOL_ON_PSEUDO_HALO2D, IDIMY_C = ", IDIMY_C, &
", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT ", IIE - IIB + 1 + 2*JPHEXT = ", IIE - IIB + 1 + 2*JPHEXT
CALL ABORT CALL ABORT
STOP 'ERROR in EXTRAPOL_ON_PSEUDO_HALO2D'
ENDIF ENDIF
ENDIF ENDIF
! !
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Documentation | Mentions légales | CGU | Accessibilité : non conforme