Changeset 2517 for branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

Timestamp:

2010-12-23T17:34:46+01:00 (14 years ago)

Author:

cetlod

Message:

v3.3beta:Ensure restartability of ORCA2_OFF_PISCES & define lk_offline flag for OFFLINE mode

File:

• branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/trcrst.F90 (modified) (9 diffs)

branches/nemo_v3_3_beta/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

@@ -34 +34 @@
    USE trcrst_c14b     ! C14 bomb restart
    USE trcrst_my_trc   ! MY_TRC   restart
-#if defined key_offline
-    USE daymod
-#endif
+   USE daymod
    IMPLICIT NONE
    PRIVATE
@@ -63 +61 @@
       !!----------------------------------------------------------------------
       !
-# if ! defined key_offline
-      IF( kt == nit000 ) lrst_trc = .FALSE. 
-# else
-      IF( kt == nit000 ) THEN
-        lrst_trc = .FALSE. 
-        nitrst = nitend  
-      ENDIF
-
-      IF( MOD( kt - 1, nstock ) == 0 ) THEN
-         ! we use kt - 1 and not kt - nit000 to keep the same periodicity from the beginning of the experiment
-         nitrst = kt + nstock - 1                  ! define the next value of nitrst for restart writing
-         IF( nitrst > nitend )   nitrst = nitend   ! make sure we write a restart at the end of the run
-      ENDIF
-# endif
-     ! to get better performances with NetCDF format:
-     ! we open and define the tracer restart file one tracer time step before writing the data (-> at nitrst - 2*nn_dttrc + 1)
-     ! except if we write tracer restart files every tracer time step or if a tracer restart file was writen at nitend - 2*nn_dttrc + 1
-     IF( kt == nitrst - 2*nn_dttrc + 1 .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc + 1 .AND. .NOT. lrst_trc ) ) THEN
+      IF( lk_offline ) THEN
+         IF( kt == nit000 ) THEN
+            lrst_trc = .FALSE.
+            nitrst = nitend
+         ENDIF
+
+         IF( MOD( kt - 1, nstock ) == 0 ) THEN
+            ! we use kt - 1 and not kt - nit000 to keep the same periodicity from the beginning of the experiment
+            nitrst = kt + nstock - 1                  ! define the next value of nitrst for restart writing
+            IF( nitrst > nitend )   nitrst = nitend   ! make sure we write a restart at the end of the run
+         ENDIF
+      ELSE
+         IF( kt == nit000 ) lrst_trc = .FALSE.
+      ENDIF
+
+      ! to get better performances with NetCDF format:
+      ! we open and define the tracer restart file one tracer time step before writing the data (-> at nitrst - 2*nn_dttrc + 1)
+      ! except if we write tracer restart files every tracer time step or if a tracer restart file was writen at nitend - 2*nn_dttrc + 1
+      IF( kt == nitrst - 2*nn_dttrc + 1 .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc + 1 .AND. .NOT. lrst_trc ) ) THEN
          ! beware of the format used to write kt (default is i8.8, that should be large enough)
          IF( nitrst > 1.0e9 ) THEN   ;   WRITE(clkt,*       ) nitrst
@@ -102 +101 @@
       !!----------------------------------------------------------------------
       INTEGER  ::  jn     
-      INTEGER  ::  iarak0 
-      REAL(wp) ::  zarak0
       INTEGER  ::  jlibalt = jprstlib
       LOGICAL  ::  llok
@@ -126 +123 @@
       CALL trc_rst_cal( nit000, 'READ' )   ! calendar
 
-      IF( ln_trcadv_cen2 .OR. ln_trcadv_tvd ) THEN   ;   iarak0 = 1
-      ELSE                                           ;   iarak0 = 0
-      ENDIF
-      CALL iom_get( numrtr, 'arak0', zarak0 )
-
-      IF( iarak0 /= NINT( zarak0 ) )   &                           ! Control of the scheme
-         & CALL ctl_stop( ' ===>>>> : problem with advection scheme', &
-         & ' it must be the same type for both restart and previous run', &
-         & ' centered or euler '  )
-      IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '    arakawa option      : ', NINT( zarak0 )
-
       ! READ prognostic variables and computes diagnostic variable
       DO jn = 1, jptra
@@ -171 +156 @@
 
       CALL trc_rst_cal( kt, 'WRITE' )   ! calendar
-
-      IF( ln_trcadv_cen2 .OR. ln_trcadv_tvd ) THEN   ;   zarak0 = 1.
-      ELSE                                           ;   zarak0 = 0.
-      ENDIF
-      CALL iom_rstput( kt, nitrst, numrtw, 'arak0', zarak0 )
-
+      CALL iom_rstput( kt, nitrst, numrtw, 'rdttrc1', rdttrc(1) )   ! surface passive tracer time step
       ! prognostic variables 
       ! -------------------- 
@@ -232 +212 @@
       CHARACTER(len=*), INTENT(in) ::   cdrw       ! "READ"/"WRITE" flag
       !
-      REAL(wp) ::  zkt
-#if defined key_offline
+      REAL(wp) ::  zkt, zrdttrc1
       REAL(wp) ::  zndastp
-#endif
 
       ! Time domain : restart
@@ -257 +235 @@
             &   CALL ctl_stop( ' ===>>>> : problem with nit000 for the restart',                 &
             &                  ' verify the restart file or rerun with nn_rsttr = 0 (namelist)' )
-#if defined key_offline
-         ! define ndastp and adatrj
-         IF ( nn_rsttr == 2 ) THEN
-            CALL iom_get( numrtr, 'ndastp', zndastp ) 
-            ndastp = NINT( zndastp )
-            CALL iom_get( numrtr, 'adatrj', adatrj  )
-         ELSE
-            ndastp = ndate0 - 1     ! ndate0 read in the namelist in dom_nam
-            adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday
-            ! note this is wrong if time step has changed during run
+         IF( lk_offline ) THEN      ! set the date in offline mode
+            ! Check dynamics and tracer time-step consistency and force Euler restart if changed
+            IF( iom_varid( numrtr, 'rdttrc1', ldstop = .FALSE. ) > 0 )   THEN
+               CALL iom_get( numrtr, 'rdttrc1', zrdttrc1 )
+               IF( zrdttrc1 /= rdttrc(1) )   neuler = 0
+            ENDIF
+            !                          ! define ndastp and adatrj
+            IF ( nn_rsttr == 2 ) THEN
+               CALL iom_get( numrtr, 'ndastp', zndastp ) 
+               ndastp = NINT( zndastp )
+               CALL iom_get( numrtr, 'adatrj', adatrj  )
+            ELSE
+               ndastp = ndate0 - 1     ! ndate0 read in the namelist in dom_nam
+               adatrj = ( REAL( nit000-1, wp ) * rdttra(1) ) / rday
+               ! note this is wrong if time step has changed during run
+            ENDIF
+            !
+            IF(lwp) THEN
+              WRITE(numout,*) ' *** Info used values : '
+              WRITE(numout,*) '   date ndastp                                      : ', ndastp
+              WRITE(numout,*) '   number of elapsed days since the begining of run : ', adatrj
+              WRITE(numout,*)
+            ENDIF
+            !
+            CALL day_init          ! compute calendar
+            !
          ENDIF
          !
-         IF(lwp) THEN
-           WRITE(numout,*) ' *** Info used values : '
-           WRITE(numout,*) '   date ndastp                                      : ', ndastp
-           WRITE(numout,*) '   number of elapsed days since the begining of run : ', adatrj
-           WRITE(numout,*)
-         ENDIF
-         !
-         CALL day_init          ! compute calendar
-         !
-#endif
-
       ELSEIF( TRIM(cdrw) == 'WRITE' ) THEN
          !
@@ -287 +270 @@
             IF(lwp) WRITE(numout,*) '~~~~~~~'
          ENDIF
-         ! calendar control
          CALL iom_rstput( kt, nitrst, numrtw, 'kt'     , REAL( kt    , wp) )   ! time-step
          CALL iom_rstput( kt, nitrst, numrtw, 'ndastp' , REAL( ndastp, wp) )   ! date
@@ -304 +286 @@
       !!----------------------------------------------------------------------
 
-      INTEGER  :: ji, jj, jk, jn
+      INTEGER  :: jn
       REAL(wp) :: zdiag_var, zdiag_varmin, zdiag_varmax, zdiag_tot
-      REAL(wp) :: zder, zvol
+      REAL(wp) :: zder
       !!----------------------------------------------------------------------