Changeset 6336

Timestamp:

2016-02-19T14:53:05+01:00 (9 years ago)

Author:

mathiot

Message:

correction of truncation error for rst date stamp

Location:

branches/UKMO/restart_datestamp/NEMOGCM/NEMO/OPA_SRC

Files:

• ICB/icbrst.F90 (modified) (2 diffs)
• IOM/restart.F90 (modified) (2 diffs)

branches/UKMO/restart_datestamp/NEMOGCM/NEMO/OPA_SRC/ICB/icbrst.F90

@@ -235 +235 @@
       INTEGER             ::   iyear, imonth, iday
       REAL (wp)           ::   zsec
+      REAL (wp)           ::   zfjulday
       CHARACTER(len=256)  :: cl_path
       CHARACTER(len=256)  :: cl_filename
@@ -246 +247 @@
       IF( cl_path(LEN_TRIM(cl_path):) /= '/' ) cl_path = TRIM(cl_path) // '/'
       IF ( ln_rstdate ) THEN
-         CALL ju2ymds( fjulday + rdttra(1) / rday, iyear, imonth, iday, zsec )           
+         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

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

@@ -57 +57 @@
       INTEGER             ::   iyear, imonth, iday
       REAL (wp)           ::   zsec
+      REAL (wp)           ::   zfjulday
       !!
       CHARACTER(LEN=20)   ::   clkt     ! ocean time-step deine as a character
@@ -85 +86 @@
          IF( nitrst <= nitend .AND. nitrst > 0 ) THEN 
             IF ( ln_rstdate ) THEN
-               CALL ju2ymds( fjulday + ((kt-nit000)*rdttra(1)) / rday, iyear, imonth, iday, zsec )           
+               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