Changeset 7378 for branches/2016

Timestamp:

2016-11-30T11:15:30+01:00 (7 years ago)

Author:

clem

Message:

debug agrif

File:

• branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/NST_SRC/agrif_user.F90 (modified) (5 diffs)

branches/2016/dev_CNRS_2016/NEMOGCM/NEMO/NST_SRC/agrif_user.F90

@@ -39 +39 @@
       jpkm1   = jpk-1                                         
       jpij    = jpi*jpj 
- !SF      jpidta  = jpiglo
- !SF     jpjdta  = jpjglo
       nperio  = 0
       jperio  = 0
@@ -239 +237 @@
          WRITE(cl_check2,*)  NINT(rdt)
          WRITE(cl_check3,*)  NINT(Agrif_Parent(rdt)/Agrif_Rhot())
-         CALL ctl_warn( 'incompatible time step between grids',   &
+         CALL ctl_stop( 'incompatible time step between ocean grids',   &
                &               'parent grid value : '//cl_check1    ,   & 
                &               'child  grid value : '//cl_check2    ,   & 
-               &               'value on child grid will be changed to : '//cl_check3 )
-         rdt=Agrif_Parent(rdt)/Agrif_Rhot()
+               &               'value on child grid should be changed to : '//cl_check3 )
       ENDIF
 
@@ -566 +563 @@
 #  endif
 
+#if defined key_lim3
+SUBROUTINE Agrif_InitValues_cont_lim3
+   !!----------------------------------------------------------------------
+   !!                 *** ROUTINE Agrif_InitValues_cont_lim3 ***
+   !!
+   !! ** Purpose :: Initialisation of variables to be interpolated for LIM3
+   !!----------------------------------------------------------------------
+   USE Agrif_Util
+   USE sbc_oce, ONLY : nn_fsbc  ! clem: necessary otherwise Agrif_Parent(nn_fsbc) = nn_fsbc
+   USE ice
+   USE agrif_ice
+   USE in_out_manager
+   USE agrif_lim3_update
+   USE agrif_lim3_interp
+   USE lib_mpp
+   !
+   IMPLICIT NONE
+   !!----------------------------------------------------------------------
+   !
+   ! Declaration of the type of variable which have to be interpolated (parent=>child)
+   !----------------------------------------------------------------------------------
+   CALL agrif_declare_var_lim3
+
+   ! Controls (clem)
+   ! 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
+      CALL ctl_stop('rhot * nn_fsbc(parent) /= N * nn_fsbc(child), therefore nn_fsbc(child) should be set to 1 or nn_fsbc(parent)')
+   ENDIF
+
+   ! stop if update frequency is different from nn_fsbc
+   IF( nbclineupdate > nn_fsbc )  CALL ctl_stop('With ice model on child grid, nn_cln_update should be set to 1 or nn_fsbc')
+
+
+   ! First Interpolations (using "after" ice subtime step => lim_nbstep=1)
+   !----------------------------------------------------------------------
+   lim_nbstep = 1
+   CALL agrif_interp_lim3('U') ! interpolation of ice velocities
+   CALL agrif_interp_lim3('V') ! interpolation of ice velocities
+   CALL agrif_interp_lim3('T') ! interpolation of ice tracers
+   lim_nbstep = 0
+   
+   ! Update in case 2 ways
+   !----------------------
+   CALL agrif_update_lim3(0)
+
+   !
+END SUBROUTINE Agrif_InitValues_cont_lim3
+
+SUBROUTINE agrif_declare_var_lim3
+   !!----------------------------------------------------------------------
+   !!                 *** ROUTINE agrif_declare_var_lim3 ***
+   !!
+   !! ** Purpose :: Declaration of variables to be interpolated for LIM3
+   !!----------------------------------------------------------------------
+   USE Agrif_Util
+   USE ice
+
+   IMPLICIT NONE
+   !!----------------------------------------------------------------------
+   !
+   ! 1. Declaration of the type of variable which have to be interpolated (parent=>child)
+   !       agrif_declare_variable(position,1st point index,--,--,dimensions,name)
+   !-------------------------------------------------------------------------------------
+   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 )
+   CALL agrif_declare_variable((/1,2/)    ,(/2,3/),(/'x','y'/)    ,(/1,1/)  ,(/nlci,nlcj/)                      ,u_ice_id   )
+   CALL agrif_declare_variable((/2,1/)    ,(/3,2/),(/'x','y'/)    ,(/1,1/)  ,(/nlci,nlcj/)                      ,v_ice_id   )
+
+   ! 2. Set interpolations (normal & tangent to the grid cell for velocities)
+   !-----------------------------------
+   CALL Agrif_Set_bcinterp(tra_ice_id,  interp = AGRIF_linear)
+   CALL Agrif_Set_bcinterp(u_ice_id  , interp1 = Agrif_linear,interp2 = AGRIF_ppm   )
+   CALL Agrif_Set_bcinterp(v_ice_id  , interp1 = AGRIF_ppm   ,interp2 = Agrif_linear)
+
+   ! 3. Set location of interpolations
+   !----------------------------------
+   CALL Agrif_Set_bc(tra_ice_id,(/0,1/))
+   CALL Agrif_Set_bc(u_ice_id  ,(/0,1/))
+   CALL Agrif_Set_bc(v_ice_id  ,(/0,1/))
+
+   ! 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(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   )
+
+END SUBROUTINE agrif_declare_var_lim3
+#endif
+
 
 # if defined key_top
@@ -622 +707 @@
          WRITE(cl_check2,*)  rdt
          WRITE(cl_check3,*)  rdt*Agrif_Rhot()
-         CALL ctl_warn( 'incompatible time step between grids',   &
+         CALL ctl_stop( 'incompatible time step between grids',   &
                &               'parent grid value : '//cl_check1    ,   & 
                &               'child  grid value : '//cl_check2    ,   & 
-               &               'value on child grid will be changed to  &
+               &               'value on child grid should be changed to  &
                &               :'//cl_check3  )
-         rdt=rdt*Agrif_Rhot()
       ENDIF
 
@@ -779 +863 @@
    IF( agrif_oce_alloc()  > 0 )   CALL ctl_warn('agrif agrif_oce_alloc: allocation of arrays failed')
 # if defined key_lim2
-   IF( agrif_ice_alloc()  > 0 )   CALL ctl_stop('agrif agrif_ice_alloc: allocation of arrays failed')
+   IF( agrif_ice_alloc()  > 0 )   CALL ctl_stop('agrif agrif_ice_alloc: allocation of arrays failed') ! only for LIM2 (not LIM3)
 # endif
    !