Changeset 7158 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

Timestamp:

2016-10-29T01:21:05+02:00 (7 years ago)

Author:

clem

Message:

debug branch

File:

• branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90 (modified) (6 diffs)

branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim.F90

@@ -132 +132 @@
 
 # if defined key_agrif
-         IF( .NOT. Agrif_Root() ) THEN
-            lim_nbstep = MOD( lim_nbstep, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) + 1
-         ENDIF
+         IF( .NOT. Agrif_Root() )  lim_nbstep = MOD( lim_nbstep, Agrif_irhot() * Agrif_Parent(nn_fsbc) / nn_fsbc ) + 1
 # endif
 
@@ -180 +178 @@
                                       CALL lim_rst_opn( kt )   ! Open Ice restart file
          !
-#if defined key_agrif
-          IF( .NOT. Agrif_Root() )    CALL agrif_interp_lim3('T')
-#endif
          ! --- zap this if no ice dynamics --- !
          IF( .NOT. lk_c1d .AND. ln_limdyn ) THEN
@@ -201 +196 @@
          ! ---
 #if defined key_agrif
-          IF( .NOT. Agrif_Root() )    CALL agrif_interp_lim3('T')
+         IF( .NOT. Agrif_Root() )     CALL agrif_interp_lim3('T')
 #endif
 #if defined key_bdy
@@ -270 +265 @@
                                       CALL lim_var_agg( 2 )     ! necessary calls (at least for coupling)
                                       !
+# if defined key_agrif
+!!         IF( .NOT. Agrif_Root() )     CALL Agrif_ChildGrid_To_ParentGrid()  ! clem: should be called at the update frequency only (cf agrif_lim3_update)
+# endif
                                       CALL lim_sbc_flx( kt )    ! -- Update surface ocean mass, heat and salt fluxes
-                                      !
+# if defined key_agrif
+!!         IF( .NOT. Agrif_Root() )     CALL Agrif_ParentGrid_To_ChildGrid()  ! clem: should be called at the update frequency only (cf agrif_lim3_update)
+# endif
          IF( ln_limdiahsb )           CALL lim_diahsb( kt )     ! -- Diagnostics and outputs 
          !
@@ -438 +438 @@
       !
       ! sea-ice timestep and inverse
-      rdt_ice   = nn_fsbc * rdttra(1)  
+      rdt_ice   = REAL(nn_fsbc) * rdt  
       r1_rdtice = 1._wp / rdt_ice 
 
@@ -448 +448 @@
       IF( lwp .AND. ln_limdiachk )  CALL ctl_warn('online conservation check activated but it does not work with BDY')
 #endif
+      !
+      IF( lwp ) WRITE(numout,*) '   ice timestep rdt_ice  = ', rdt_ice
       !
    END SUBROUTINE ice_run