Changeset 1655

Timestamp:

2009-10-14T16:26:01+02:00 (15 years ago)

Author:

cetlod

Message:

Call to write restart in TOP not properly done, see ticket:43

File:

• trunk/NEMO/TOP_SRC/trcrst.F90 (modified) (2 diffs)

trunk/NEMO/TOP_SRC/trcrst.F90

@@ -57 +57 @@
       !!----------------------------------------------------------------------
       !
+# if ! defined key_off_tra
       IF( kt == nit000 ) lrst_trc = .FALSE. 
-# if defined key_off_tra
+# 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
@@ -64 +70 @@
          IF( nitrst > nitend )   nitrst = nitend   ! make sure we write a restart at the end of the run
       ENDIF
-#endif
-      
-      IF( kt == nitrst - ndttrc .OR. nitend - nit000 + 1 < 2 * ndttrc ) THEN
-         ! beware if model runs less than 2*ndttrc time step
+# 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*ndttrc + 1)
+     ! except if we write tracer restart files every tracer time step or if a tracer restart file was writen at nitend - 2*ndttrc + 1
+     IF( kt == nitrst - 2*ndttrc + 1 .OR. nstock == ndttrc .OR. ( kt == nitend - ndttrc + 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