Changeset 1239 for trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

Timestamp:

2008-12-31T10:35:35+01:00 (15 years ago)

Author:

ctlod

Message:

this version allows to stabilize the computation of En production term, see ticket: #281

File:

• trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/SBC/sbcice_lim_2.F90

@@ -183 +183 @@
          !  Ice model step  !
          ! ---------------- !
+
+#if defined key_zdftke2
+         IF ( kt .NE. nitend+1 ) THEN
                                         CALL lim_rst_opn_2  ( kt )      ! Open Ice restart file
+         ENDIF
+#else
+                                        CALL lim_rst_opn_2  ( kt )      ! Open Ice restart file
+#endif
          IF( .NOT. lk_c1d ) THEN                                        ! Ice dynamics & transport (not in 1D case)
                                         CALL lim_dyn_2      ( kt )           ! Ice dynamics    ( rheology/dynamics )
@@ -197 +204 @@
                                         CALL lim_thd_2      ( kt )      ! Ice thermodynamics 
                                         CALL lim_sbc_2      ( kt )      ! Ice/Ocean Mass & Heat fluxes 
+#if defined key_zdftke2
+         IF( kt .NE. nitend+1 ) THEN
+            IF( ( MOD( kt+nn_fsbc-1, ninfo ) == 0 .OR. ntmoy == 1 ) .AND. .NOT. lk_mpp )   &
+               &                        CALL lim_dia_2      ( kt )      ! Ice Diagnostics 
+                                        CALL lim_wri_2      ( kt )      ! Ice outputs 
+            IF( lrst_ice )              CALL lim_rst_write_2( kt )      ! Ice restart file 
+            !
+         ENDIF
+#else
          IF( ( MOD( kt+nn_fsbc-1, ninfo ) == 0 .OR. ntmoy == 1 ) .AND. .NOT. lk_mpp )   &
-            &                           CALL lim_dia_2      ( kt )      ! Ice Diagnostics 
-                                        CALL lim_wri_2      ( kt )      ! Ice outputs 
-         IF( lrst_ice )                 CALL lim_rst_write_2( kt )      ! Ice restart file 
+            &                           CALL lim_dia_2      ( kt )      ! Ice Diagnostics
+                                        CALL lim_wri_2      ( kt )      ! Ice outputs
+         IF( lrst_ice )                 CALL lim_rst_write_2( kt )      ! Ice restart file
          !
+#endif
       ENDIF
       !