Changeset 12477 for NEMO/branches/UKMO/r12083_restart_datestamp/src/TOP

Timestamp:

2020-02-27T13:20:58+01:00 (4 years ago)

Author:

jcastill

Message:

Changes as in the original branch, plus changes for bgc restart (in branch AMM15_v3_6_STABLE_package_collate)

File:

• NEMO/branches/UKMO/r12083_restart_datestamp/src/TOP/trcrst.F90 (modified) (3 diffs)

NEMO/branches/UKMO/r12083_restart_datestamp/src/TOP/trcrst.F90

@@ -20 +20 @@
    !!----------------------------------------------------------------------
    USE oce_trc
+   USE dom_oce, ONLY: fjulday 
+   USE phycst, ONLY: rday 
+   USE in_out_manager, ONLY: ln_rstdate
    USE trc
    USE iom
+   USE ioipsl, ONLY: ju2ymds
    USE daymod
    USE lib_mpp
@@ -47 +51 @@
       !!----------------------------------------------------------------------
       INTEGER, INTENT(in) ::   kt       ! number of iteration
+      INTEGER             ::   iyear, imonth, iday 
+      REAL (wp)           ::   zsec 
+      REAL (wp)           ::   zfjulday
       !
       CHARACTER(LEN=20)   ::   clkt     ! ocean time-step define as a character
@@ -78 +85 @@
       ! 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 .OR. nn_stock == nn_dttrc .OR. ( kt == nitend - nn_dttrc .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
-         ELSE                        ;   WRITE(clkt,'(i8.8)') nitrst
+      IF( kt == nitrst - 2*nn_dttrc + 1 .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc .AND. .NOT. lrst_trc ) ) THEN 
+         if ( ln_rstdate ) THEN 
+            zfjulday = fjulday + rdttrc(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 > 1.0e9 ) THEN   ;   WRITE(clkt,*       ) nitrst 
+            ELSE                        ;   WRITE(clkt,'(i8.8)') nitrst 
+            ENDIF
          ENDIF
          ! create the file