Changeset 14655

Timestamp:

2021-03-30T16:41:30+02:00 (3 years ago)

Author:

jchanut

Message:

AGRIF: Set ghost cell number to 0 over root grid to properly define analytical grids in zooms, #2638

File:

• NEMO/branches/2021/dev_r14608_AGRIF_domcfg/src/NST/agrif_user.F90 (modified) (1 diff)

NEMO/branches/2021/dev_r14608_AGRIF_domcfg/src/NST/agrif_user.F90

@@ -850 +850 @@
       ENDIF
 
-      ! Set the number of ghost cells according to periodicity
-      nbghostcells_x   = nbghostcells
-      nbghostcells_y_s = nbghostcells
-      nbghostcells_y_n = nbghostcells
-
-      !IF ((jperio == 1).OR.(jperio == 4) ) THEN   
-      IF ( l_Iperio ) nbghostcells_x = 0
-      lk_south = .NOT. ( Agrif_Iy() == 1 )
-      IF( .NOT. lk_south )   nbghostcells_y_s = 0
-
 ! JC => side effects of lines below to be checked:
       IF (.not.agrif_root()) THEN
-        lk_west  = .NOT. ( Agrif_Ix() == 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(Nj0glo) -1 )
-      !
-
-      !
-        IF(    l_Iperio    )   nbghostcells_x   = 0
-        IF( .NOT. lk_south )   nbghostcells_y_s = 0
-        IF( .NOT. lk_north )   nbghostcells_y_n = 0
-      !
-      ! Some checks
-        IF( (.NOT.ln_vert_remap).AND.(jpkglo>Agrif_Parent(jpkglo)) )                    CALL ctl_stop( 'STOP',    &
+         nbghostcells_x   = nbghostcells
+         nbghostcells_y_s = nbghostcells
+         nbghostcells_y_n = nbghostcells
+
+         lk_west  = .NOT. ( Agrif_Ix() == 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(Nj0glo) -1 )
+         !
+         ! Correct number of ghost cells according to periodicity
+         !
+         IF(    l_Iperio    )   nbghostcells_x   = 0
+         IF( .NOT. lk_south )   nbghostcells_y_s = 0
+         IF( .NOT. lk_north )   nbghostcells_y_n = 0
+         !
+         ! Some checks
+         IF( (.NOT.ln_vert_remap).AND.(jpkglo>Agrif_Parent(jpkglo)) )                    CALL ctl_stop( 'STOP',    &
            &   'agrif_nemo_init: Agrif children must have less or equal number of vertical levels without ln_vert_remap defined' ) 
-        IF( jpiglo /= nbcellsx + 2 + 2*nn_hls + nbghostcells_x   + nbghostcells_x   )   CALL ctl_stop( 'STOP',    &
+         IF( jpiglo /= nbcellsx + 2 + 2*nn_hls + nbghostcells_x   + nbghostcells_x   )   CALL ctl_stop( 'STOP',    &
            &   'agrif_nemo_init: Agrif children requires jpiglo == nbcellsx + 2 + 2*nn_hls + 2*nbghostcells_x' )
-        IF( jpjglo /= nbcellsy + 2 + 2*nn_hls + nbghostcells_y_s + nbghostcells_y_n )   CALL ctl_stop( 'STOP',    &
+         IF( jpjglo /= nbcellsy + 2 + 2*nn_hls + nbghostcells_y_s + nbghostcells_y_n )   CALL ctl_stop( 'STOP',    &
            &   'agrif_nemo_init: Agrif children requires jpjglo == nbcellsy + 2 + 2*nn_hls + nbghostcells_y_s + nbghostcells_y_n' )
-        IF( ln_use_jattr )   CALL ctl_stop( 'STOP', 'agrif_nemo_init:Agrif children requires ln_use_jattr = .false. ' )
-      ENDIF
-
+         IF( ln_use_jattr )   CALL ctl_stop( 'STOP', 'agrif_nemo_init:Agrif children requires ln_use_jattr = .false. ' )
+      ELSE
+         ! Root grid
+         nbghostcells_x   = 0
+         nbghostcells_y_s = 0
+         nbghostcells_y_n = 0
+      ENDIF
       !
       !