Changeset 3798

Timestamp:

05/02/18 11:50:27 (6 years ago)

Author:

mafoipsl

Message:

Consolidation of Reproducibility_after_restart ie it's now working for all PeriodLength?
and test included to verify start on XXXX0101.

File:

• CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/DRIVER/opa9.driver (modified) (2 diffs)

CONFIG/UNIFORM/v6/IPSLCM6/GENERAL/DRIVER/opa9.driver

@@ -121 +121 @@
       IGCM_debug_Print 1 'Reproducibility_after_Restart forced ie kt forced to kt red in restart file : ' ${opa9_UserChoices_Reproducibility_after_restart}
 
-      ## This option is uncompatible with PeriodLength=1M. See later patch overwrite ndastp and kt in ocean restart file
-
-      if [ X${config_UserChoices_PeriodLength} != X"1Y" ] ; then
-        IGCM_debug_Exit " PeriodLength = ${config_UserChoices_PeriodLength} is uncompatible with Reproducibility_after_Restart=y, please use 1Y"
-      fi
-
     elif ( [ "${CumulPeriod}" -eq 1 ] && [ "${config_OCE_Restart}" = "n" ] ) ; then
       OPA_LRSTAR=.FALSE.
@@ -170 +164 @@
       ##-- Getting from restart file the last time step of the run at which we want to re-start
       NEMO_KT_FROM_RESTART=$(ncdump -v kt ${NEMO_RESTART_FILE}|grep 'kt ='|awk '{print $3}' )
-      ## Force the date in the restart file to be the last day before the begiining of the new run
-      (( yearrst = year - 1 ))
-      IGCM_debug_Print 1 "cmd done to modify NEMO restartfile:   ncap2 -s ndastp=${yearrst}1231 ${NEMO_RESTART_FILE}  tmp.nc"
-      ## Take care : 1231 is forced here. Still to be done for monthly run or 
-      ## make nrstdt=1 working for all NEMO component. Not done for sea ice yet.
-      ncap2 -s "ndastp=${yearrst}1231" ${NEMO_RESTART_FILE}  tmp.nc
-      mv tmp.nc ${NEMO_RESTART_FILE}
+
+      if ( [ "${CumulPeriod}" -eq 1 ] ) ; then
+        ## Force the date in the restart file to be the last day before the begining of the new run
+        (( yearrst = year - 1 ))
+        IGCM_debug_Print 1 "cmd done to modify NEMO restartfile:   ncap2 -s ndastp=${yearrst}1231 ${NEMO_RESTART_FILE}  tmp.nc"
+        ## Take care : 1231 is forced here. Still to be done for monthly run or 
+        ## make nrstdt=1 working for all NEMO component. Not done for sea ice yet.
+        ncap2 -s "ndastp=${yearrst}1231" ${NEMO_RESTART_FILE}  tmp.nc
+        mv tmp.nc ${NEMO_RESTART_FILE}
+        ##-- Getting from restart file the "new" restart date of run at which we want to re-start
+        NEMO_NDASTP_FROM_RESTART=$(ncdump -v ndastp ${NEMO_RESTART_FILE}|grep 'ndastp ='|awk '{print $3}' )
+        IGCM_debug_Print 1 "NEMO_NDASTP (from Restart) for 1st Period       : ${NEMO_NDASTP_FROM_RESTART}"
+
+        IGCM_date_GetYearMonthDay ${PeriodDateBegin} Verif_repro_year Verif_repro_month Verif_repro_day
+
+        if [ X${Verif_repro_month}${Verif_repro_day} != X"0101" ] ; then
+          IGCM_debug_Exit " PeriodDateBegin = ${PeriodDateBegin}is uncompatible with Reproducibility_after_Restart=y, please start on 0101"
+        fi
+
+      fi
+
       (( NEMO_NIT000 = NEMO_KT_FROM_RESTART + 1 ))
       (( NEMO_NITEND = NEMO_KT_FROM_RESTART + OPA_NITEND - OPA_NIT000 + 1 ))
       IGCM_debug_Print 1 "NEMO_NIT000 (from Restart)        : ${NEMO_NIT000}"
       IGCM_debug_Print 1 "NEMO_NITEND (from Restart + ... ) : ${NEMO_NITEND}"
-      ##-- Getting from restart file the "new" restart date of run at which we want to re-start
-      NEMO_NDASTP_FROM_RESTART=$(ncdump -v ndastp ${NEMO_RESTART_FILE}|grep 'ndastp ='|awk '{print $3}' )
-      IGCM_debug_Print 1 "NEMO_NDASTP (from Restart)       : ${NEMO_NDASTP_FROM_RESTART}"
       ##-- We force one restart file at the end of the trunk
       OPA_NSTOCK="${NEMO_NITEND}"