Changeset 6939 for branches/UKMO/CO6_KD490/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

Timestamp:

2016-09-19T16:36:50+02:00 (8 years ago)

Author:

deazer

Message:

Updated to include bug bix restart date label with respect to truncation error.

File:

• branches/UKMO/CO6_KD490/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90 (modified) (2 diffs)

branches/UKMO/CO6_KD490/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90

@@ -55 +55 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt     ! ocean time-step
-      INTEGER             ::   iyear, imonth, iday 
-      REAL (wp)           ::   zsec 
+      INTEGER             ::   iyear, imonth, iday
+      REAL (wp)           ::   zsec
+      REAL (wp)           ::   zfjulday
       !!
       CHARACTER(LEN=20)   ::   clkt     ! ocean time-step deine as a character
@@ -84 +85 @@
       IF( kt == nitrst - 1 .OR. nstock == 1 .OR. ( kt == nitend .AND. .NOT. lrst_oce ) ) THEN
          IF( nitrst <= nitend .AND. nitrst > 0 ) THEN 
-            IF ( ln_rstdate ) THEN 
-              CALL ju2ymds( fjulday + rdttra(1) / rday, iyear, imonth, iday, zsec )            
-              WRITE(clkt, '(i4.4,2i2.2)') iyear, imonth, iday 
-            ELSE 
-              ! beware of the format used to write kt (default is i8.8, that should be large enough...) 
-              IF( nitrst > 999999999 ) THEN   ;   WRITE(clkt, *       ) nitrst 
-              ELSE                            ;   WRITE(clkt, '(i8.8)') nitrst 
-              ENDIF 
+            IF ( ln_rstdate ) THEN
+               zfjulday = fjulday + rdttra(1) / rday
+               IF( ABS(zfjulday - REAL(NINT(zfjulday),wp)) < 0.1 / rday )   zfjulday = REAL(NINT(zfjulday),wp)   ! avoid truncation error
+               CALL ju2ymds( zfjulday, iyear, imonth, iday, zsec )           
+               WRITE(clkt, '(i4.4,2i2.2)') iyear, imonth, iday
+            ELSE
+               ! beware of the format used to write kt (default is i8.8, that should be large enough...)
+               IF( nitrst > 999999999 ) THEN   ;   WRITE(clkt, *       ) nitrst
+               ELSE                            ;   WRITE(clkt, '(i8.8)') nitrst
+               ENDIF
             ENDIF
             ! create the file