Changeset 15265 for NEMO/trunk/tests/DOME

Timestamp:

2021-09-16T13:13:13+02:00 (3 years ago)

Author:

jchanut

Message:

Improve ghost cell initialization with AGRIF + minor changes such as missing _wp, tests namelists updates, etc... can be assigned to #2638

Location:

NEMO/trunk/tests/DOME

Files:

• EXPREF/1_namelist_cfg (modified) (1 diff)
• MY_SRC/usrdef_nam.F90 (modified) (1 diff)

NEMO/trunk/tests/DOME/EXPREF/1_namelist_cfg

@@ -7 +7 @@
 &namagrif      !  AGRIF zoom                                            ("key_agrif")
 !-----------------------------------------------------------------------
-   ln_agrif_2way   = .true. !  activate two way nesting
-   ln_init_chfrpar = .false.!  initialize child grids from parent
    ln_vert_remap   = .true. !  use vertical remapping
-   ln_chk_bathy    = .true. !  =T  check the parent bathymetry
-   ln_spc_dyn      = .false.
    rn_sponge_tra   = 0.0    !  coefficient for tracer   sponge layer []
-   rn_sponge_dyn   = 0.002  !  coefficient for dynamics sponge layer []
-   rn_trelax_tra   = 0.01   !  inverse of relaxation time (in steps) for tracers []
-   rn_trelax_dyn   = 0.01   !  inverse of relaxation time (in steps) for dynamics []
 /
 !-----------------------------------------------------------------------

NEMO/trunk/tests/DOME/MY_SRC/usrdef_nam.F90

@@ -96 +96 @@
          ighost_s = nbghostcells
          ighost_n = nbghostcells
-
-         IF  ( Agrif_Ix() == 1 ) ighost_w = 1 
-         IF  ( Agrif_Ix() + nbcellsx/AGRIF_Irhox() == Agrif_Parent(Ni0glo)-1 ) ighost_e = 1 
-         IF  ( Agrif_Iy() == 1 ) ighost_s = 1 
-         IF  ( Agrif_Iy() + nbcellsy/AGRIF_Irhoy() == Agrif_Parent(Nj0glo)-1 ) ighost_n = 1 
+         ! In case one sets zoom boundaries over domain edges: 
+         IF  ( Agrif_Ix() == 2 - Agrif_Parent(nbghostcells_x_w) ) ighost_w = 1 
+         IF  ( Agrif_Ix() + nbcellsx/AGRIF_Irhox() == Agrif_Parent(Ni0glo)-Agrif_Parent(nbghostcells_x_w) ) ighost_e = 1 
+         IF  ( Agrif_Iy() == 2 - Agrif_Parent(nbghostcells_y_s) ) ighost_s = 1 
+         IF  ( Agrif_Iy() + nbcellsy/AGRIF_Irhoy() == Agrif_Parent(Nj0glo)-Agrif_Parent(nbghostcells_y_s) ) ighost_n = 1 
          kpi  = nbcellsx + ighost_w + ighost_e
          kpj  = nbcellsy + ighost_s + ighost_n