Changeset 2804 for branches/2011/dev_r2802_LOCEAN10_agrif_lim/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

Timestamp:

2011-07-18T17:05:52+02:00 (13 years ago)

Author:

rblod

Message:

dev_r2802_LOCEAN10_agrif_lim: first implementation see ticket #848

File:

• branches/2011/dev_r2802_LOCEAN10_agrif_lim/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90 (modified) (5 diffs)

branches/2011/dev_r2802_LOCEAN10_agrif_lim/NEMOGCM/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -48 +48 @@
    USE prtctl           ! Print control
 
+#if defined key_agrif
+   USE agrif_ice
+   USE agrif_lim2_update
+#endif
+
    IMPLICIT NONE
    PRIVATE
@@ -110 +115 @@
          !
          CALL ice_init_2
+#if defined key_agrif
+         IF( .NOT. Agrif_Root() ) CALL Agrif_InitValues_cont_lim2  ! AGRIF: set the meshes
+# endif
       ENDIF
 
@@ -115 +123 @@
       IF( MOD( kt-1, nn_fsbc ) == 0 ) THEN     !  Ice time-step only  !
          !                                     !----------------------!
+#if defined key_agrif
+         IF( .NOT. Agrif_Root() ) childfreq = MOD(childfreq,Agrif_rhot()) + 1
+#endif
+         !  
          !  Bulk Formulea !
          !----------------!
@@ -197 +209 @@
          IF( .NOT. lk_c1d ) THEN                       ! Ice dynamics & transport (except in 1D case)
                            CALL lim_dyn_2      ( kt )      ! Ice dynamics    ( rheology/dynamics )
+
                            CALL lim_trp_2      ( kt )      ! Ice transport   ( Advection/diffusion )
            IF( ln_limdmp ) CALL lim_dmp_2      ( kt )      ! Ice damping 
@@ -223 +236 @@
       !                                        !  at all ocean time step  !
       !                                        !--------------------------!
+#if defined key_agrif && defined key_lim2
+      IF( .NOT. Agrif_Root() )   CALL agrif_update_lim2( kt )
+#endif
       !                                               
       !                                              ! Update surface ocean stresses (only in ice-dynamic case)