Changeset 7761 for trunk/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

Timestamp:

2017-03-06T18:58:35+01:00 (7 years ago)

Author:

clem

Message:

make AGRIF and LIM3 fully compatible

File:

• trunk/NEMOGCM/NEMO/NST_SRC/agrif_user.F90 (modified) (8 diffs)

trunk/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

@@ -37 +37 @@
       jpim1   = jpi-1 
       jpjm1   = jpj-1 
-      jpkm1   = jpk-1                                         
+      jpkm1   = MAX( 1, jpk-1 )                                         
       jpij    = jpi*jpj 
       nperio  = 0
@@ -83 +83 @@
 # if defined key_top
    CALL Agrif_InitValues_cont_top
+# endif
+# if defined key_lim3
+   CALL Agrif_InitValues_cont_lim3
 # endif
    !
@@ -290 +293 @@
          ! check if vmask agree with parent along northern and southern boundaries:
          CALL Agrif_Bc_variable(vmsk_id,calledweight=1.,procname=interpvmsk)
-    ! check if tmask and vertical scale factors agree with parent over first two coarse grid points:
+         ! check if tmask and vertical scale factors agree with parent over first two coarse grid points:
          CALL Agrif_Bc_variable(e3t_id,calledweight=1.,procname=interpe3t)
          !
          IF (lk_mpp) CALL mpp_sum( kindic_agr )
-         IF( kindic_agr /= 0 ) THEN                   
+         IF( kindic_agr /= 0 ) THEN
             CALL ctl_stop('Child Bathymetry is not correct near boundaries.')
          ELSE
@@ -457 +460 @@
 !   CALL Agrif_Set_Updatetype(vb2b_update_id,update1 = Agrif_Update_Full_Weighting, update2 = Agrif_Update_Average)
 !   CALL Agrif_Set_Updatetype(sshn_id, update = Agrif_Update_Full_Weighting)
- 
+
    !
 END SUBROUTINE agrif_declare_var
@@ -579 +582 @@
    CALL agrif_declare_var_lim3
 
-   ! Controls (clem)
+   ! Controls
+
+   ! clem: For some reason, nn_fsbc(child)/=1 does not work properly (signal is largely degraded by the agrif zoom)
+   !          the run must satisfy CFL=Uice/(dx/dt) < 0.6/nn_fsbc(child)
+   !          therefore, if nn_fsbc(child)>1 one must reduce the time-step in proportion to nn_fsbc(child), which is not acceptable
+   !       If a solution is found, the following stop could be removed
+   IF( nn_fsbc > 1 )  CALL ctl_stop('nn_fsbc(child) must be set to 1 otherwise agrif and lim3 do not work properly')
+
    ! stop if rhot * nn_fsbc(parent) /= N * nn_fsbc(child) with N being integer
    IF( MOD( Agrif_irhot() * Agrif_Parent(nn_fsbc), nn_fsbc ) /= 0 )  THEN
@@ -591 +601 @@
    ! First Interpolations (using "after" ice subtime step => lim_nbstep=1)
    !----------------------------------------------------------------------
-   lim_nbstep = 1
+!!   lim_nbstep = 1
+   lim_nbstep = ( Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) ! clem: to have calledweight=1 in interp (otherwise the western border of the zoom is wrong)
    CALL agrif_interp_lim3('U') ! interpolation of ice velocities
    CALL agrif_interp_lim3('V') ! interpolation of ice velocities
@@ -618 +629 @@
    ! 1. Declaration of the type of variable which have to be interpolated (parent=>child)
    !       agrif_declare_variable(position,1st point index,--,--,dimensions,name)
+   !           ex.:  position=> 1,1 = not-centered (in i and j)
+   !                            2,2 =     centered (    -     )
+   !                 index   => 1,1 = one ghost line
+   !                            2,2 = two ghost lines
    !-------------------------------------------------------------------------------------
    CALL agrif_declare_variable((/2,2,0/),(/3,3,0/),(/'x','y','N'/),(/1,1,1/),(/nlci,nlcj,jpl*(5+nlay_s+nlay_i)/),tra_ice_id )
@@ -637 +652 @@
    ! 4. Set update type in case 2 ways (child=>parent) (normal & tangent to the grid cell for velocities)
    !--------------------------------------------------
-   CALL Agrif_Set_Updatetype(tra_ice_id, update = AGRIF_Update_Average) ! clem je comprends pas average/copy
+   CALL Agrif_Set_Updatetype(tra_ice_id, update = AGRIF_Update_Average)
    CALL Agrif_Set_Updatetype(u_ice_id  ,update1 = Agrif_Update_Copy   , update2 = Agrif_Update_Average)
    CALL Agrif_Set_Updatetype(v_ice_id  ,update1 = Agrif_Update_Average, update2 = Agrif_Update_Copy   )