Changeset 13673 for NEMO/branches

Timestamp:

2020-10-26T08:51:22+01:00 (4 years ago)

Author:

jchanut

Message:

#2129, Correct use of Niglo0/Njglo0 in place of jpiglo/jpjglo.
Accomodate with 2DV domains with vertical North/South? walls.

File:

• NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_user.F90 (modified) (3 diffs)

NEMO/branches/2020/dev_r13312_AGRIF-03-04_jchanut_vinterp_tstep/src/NST/agrif_user.F90

@@ -69 +69 @@
       IF ( .NOT. lk_south ) THEN
          CALL Agrif_Set_NearCommonBorderY(.TRUE.)
+      ENDIF
+
+      IF ( .NOT. lk_north ) THEN
+         CALL Agrif_Set_DistantCommonBorderY(.TRUE.)
       ENDIF
 
@@ -799 +803 @@
 ! JC => side effects of lines below to be checked:
       lk_west  = .NOT. ( Agrif_Ix() == 1 )
-      lk_east  = .NOT. ( Agrif_Ix() + nbcellsx/AGRIF_Irhox() == Agrif_Parent(jpiglo) -1 )
+      lk_east  = .NOT. ( Agrif_Ix() + nbcellsx/AGRIF_Irhox() == Agrif_Parent(Ni0glo) -1 )
       lk_south = .NOT. ( Agrif_Iy() == 1 )
-      lk_north = .NOT. ( Agrif_Iy() + nbcellsy/AGRIF_Irhoy() == Agrif_Parent(jpjglo) -1 )
+      lk_north = .NOT. ( Agrif_Iy() + nbcellsy/AGRIF_Irhoy() == Agrif_Parent(Nj0glo) -1 )
       !
       ! Set the number of ghost cells according to periodicity
@@ -810 +814 @@
       IF(   jperio == 1  )   nbghostcells_x   = 0
       IF( .NOT. lk_south )   nbghostcells_y_s = 0
-      ! For 2DV domains:
-      IF (( nbcellsy <= 3 ).AND.(AGRIF_Irhoy()==1)) THEN
-         lk_north  = .FALSE. ; lk_south = .FALSE.
-         nbghostcells_y_s = nbghostcells
-      ENDIF
-      IF (( nbcellsx <= 3 ).AND.(AGRIF_Irhox()==1)) THEN
-         lk_east  = .FALSE. ; lk_north = .FALSE.
-      ENDIF
+      IF( .NOT. lk_north )   nbghostcells_y_n = 0
+      !
       ! Some checks
       IF( jpiglo /= nbcellsx + 2 + 2*nn_hls + nbghostcells_x   + nbghostcells_x   )   CALL ctl_stop( 'STOP',    &