Changeset 5282 for branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
- Timestamp:
- 2015-05-18T17:19:50+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
r4990 r5282 82 82 WRITE(numout,*) 83 83 SELECT CASE ( jprstlib ) 84 CASE ( jprstdimg ) ; WRITE(numout,*) ' open ocean restart binary file: '//clname85 84 CASE DEFAULT ; WRITE(numout,*) ' open ocean restart NetCDF file: '//clname 86 85 END SELECT … … 110 109 !!---------------------------------------------------------------------- 111 110 112 CALL iom_rstput( kt, nitrst, numrow, 'rdt' , rdt ) ! dynamics time step 113 CALL iom_rstput( kt, nitrst, numrow, 'rdttra1', rdttra(1) ) ! surface tracer time step 111 CALL iom_rstput( kt, nitrst, numrow, 'rdt' , rdt ) ! dynamics and tracer time step 114 112 115 113 CALL iom_rstput( kt, nitrst, numrow, 'ub' , ub ) ! before fields … … 165 163 SELECT CASE ( jprstlib ) 166 164 CASE ( jpnf90 ) ; WRITE(numout,*) 'rst_read : read oce NetCDF restart file' 167 CASE ( jprstdimg ) ; WRITE(numout,*) 'rst_read : read oce binary restart file'168 165 END SELECT 169 166 IF ( snc4set%luse ) WRITE(numout,*) 'rst_read : configured with NetCDF4 support' … … 171 168 ENDIF 172 169 173 IF ( jprstlib == jprstdimg ) THEN174 ! eventually read netcdf file (monobloc) for restarting on different number of processors175 ! if {cn_ocerst_in}.nc exists, then set jlibalt to jpnf90176 INQUIRE( FILE = TRIM(cn_ocerst_in)//'.nc', EXIST = llok )177 IF ( llok ) THEN ; jlibalt = jpnf90 ; ELSE ; jlibalt = jprstlib ; ENDIF178 ENDIF179 170 CALL iom_open( cn_ocerst_in, numror, kiolib = jlibalt ) 180 171 ENDIF … … 189 180 !! ** Method : Read in restart.nc file fields which are necessary for restart 190 181 !!---------------------------------------------------------------------- 191 REAL(wp) :: zrdt , zrdttra1182 REAL(wp) :: zrdt 192 183 INTEGER :: jk 193 184 LOGICAL :: llok … … 201 192 IF( zrdt /= rdt ) neuler = 0 202 193 ENDIF 203 IF( iom_varid( numror, 'rdttra1', ldstop = .FALSE. ) > 0 ) THEN 204 CALL iom_get( numror, 'rdttra1', zrdttra1 ) 205 IF( zrdttra1 /= rdttra(1) ) neuler = 0 206 ENDIF 207 ! 194 ! 208 195 IF( iom_varid( numror, 'ub', ldstop = .FALSE. ) > 0 ) THEN 209 196 CALL iom_get( numror, jpdom_autoglo, 'ub' , ub ) ! before fields
Note: See TracChangeset
for help on using the changeset viewer.