Changeset 6160 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM
- Timestamp:
- 2015-12-22T15:32:58+01:00 (9 years ago)
- Location:
- branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
r6154 r6160 76 76 IF( nitrst > nitend ) nitrst = nitend ! make sure we write a restart at the end of the run 77 77 ENDIF 78 !! JPALM test -79 IF( kt == nitrst - 1 .OR. nstock == 1 .OR. ( kt == nitend .AND. .NOT. lrst_oce ) ) THEN80 IF( nitrst <= nitend .AND. nitrst > 0 ) THEN81 IF ( ln_rstdate ) THEN82 WRITE(numout,*) 'JPALM - we want all informations for restart file name -- '83 WRITE(numout,*) 'JPALM - 1 - conditions for openning rst file :'84 WRITE(numout,*) 'JPALM - oce restart kt = ',kt,'==', nitrst - 185 WRITE(numout,*) 'JPALM - oce restart nstock = ',nstock,'== 1 '86 WRITE(numout,*) 'JPALM - oce restart kt = ',kt,'==',nitend87 WRITE(numout,*) 'JPALM - oce restart nitrst = ',nitrst,'<=',nitend,' and > 0'88 WRITE(numout,*) '-----------------------------'89 CALL ju2ymds( fjulday + rdttra(1) / rday, iyear, imonth, iday, zsec )90 WRITE(clkt, '(i4.4,2i2.2)') iyear, imonth, iday91 WRITE(numout,*) 'JPALM - 2 - date calculated by ju2ymds: '92 WRITE(numout,*) 'JPALM - oce restart fjulday = ',fjulday93 WRITE(numout,*) 'JPALM - oce restart rdttra/rday = ',rdttra(1)/rday94 WRITE(numout,*) 'JPALM - oce restart rdttra = ',rdttra(1)95 WRITE(numout,*) 'JPALM - oce restart rday = ',rday96 WRITE(numout,*) 'JPALM - oce restart iyear = ',iyear97 WRITE(numout,*) 'JPALM - oce restart imonth = ',imonth98 WRITE(numout,*) 'JPALM - oce restart iday = ',iday99 100 ENDIF101 ENDIF102 ENDIF103 !!104 78 ! to get better performances with NetCDF format: 105 79 ! we open and define the ocean restart file one time step before writing the data (-> at nitrst - 1) -
branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r6155 r6160 86 86 IF( kt == nitrst - 2*nn_dttrc .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc .AND. .NOT. lrst_trc ) ) THEN 87 87 IF ( ln_rstdate ) THEN 88 WRITE(numout,*) '-----------------------------' 89 WRITE(numout,*) 'JPALM - we want all informations for restart file name -- ' 90 WRITE(numout,*) 'JPALM - 1 - conditions for openning rst file :' 91 WRITE(numout,*) 'JPALM - trc restart kt = ',kt,'==', nitrst - 2*nn_dttrc 92 WRITE(numout,*) 'JPALM - trc restart nstock = ',nstock,'==',nn_dttrc 93 WRITE(numout,*) 'JPALM - trc restart kt = ',kt,'==',nitend - nn_dttrc 94 WRITE(numout,*) '-----------------------------' 88 !! JPALM -- 22-12-2015 -- modif to get the good date on restart trc file name 89 !! -- the condition to open the rst file is not the same than for the dynamic rst. 90 !! -- here it - for an obscure reason - is open 2 time-step before the restart writing process 91 !! instead of 1. 92 !! -- i am not sure if someone forgot +1 in the if loop condition as 93 !! it is writen in all comments nitrst - 2*nn_dttrc + 1 and the condition is nitrst - 2*nn_dttrc 94 !! -- nevertheless we didn't wanted to broke something already working 95 !! and just adapted the part we added. 96 !! -- So instead of calling ju2ymds( fjulday + (rdttra(1)) .... 97 !! we call ju2ymds( fjulday + (2*rdttra(1)) ..... 98 !!-------------------------------------------------------------------- 95 99 CALL ju2ymds( fjulday + (2*rdttra(1)) / rday, iyear, imonth, iday, zsec ) 96 100 WRITE(clkt, '(i4.4,2i2.2)') iyear, imonth, iday 97 WRITE(numout,*) 'JPALM - 2 - date calculated by ju2ymds: '98 WRITE(numout,*) 'JPALM - trc restart fjulday = ',fjulday99 WRITE(numout,*) 'JPALM - trc restart rdttra/rday = ',rdttra(1)/rday100 WRITE(numout,*) 'JPALM - trc restart rdttra = ',rdttra(1)101 WRITE(numout,*) 'JPALM - trc restart rday = ',rday102 WRITE(numout,*) 'JPALM - trc restart iyear = ',iyear103 WRITE(numout,*) 'JPALM - trc restart imonth = ',imonth104 WRITE(numout,*) 'JPALM - trc restart iday = ',iday105 101 ELSE 106 102 ! beware of the format used to write kt (default is i8.8, that should be large enough)
Note: See TracChangeset
for help on using the changeset viewer.