- Timestamp:
- 2017-10-30T12:44:29+01:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8600_xios_read/NEMOGCM/NEMO/OPA_SRC/IOM/restart.F90
r8612 r8668 202 202 IF(.NOT.lxios_set) lxios_read = lxios_read.AND.lxios_sini 203 203 IF( lxios_read) THEN 204 rxios_context = 'nemo_rst'204 crxios_context = 'nemo_rst' 205 205 if(.NOT.lxios_set) then 206 206 IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS' 207 CALL iom_init( rxios_context )207 CALL iom_init( crxios_context ) 208 208 lxios_set = .TRUE. 209 209 endif 210 210 ENDIF 211 211 IF( TRIM(Agrif_CFixed()) /= '0' .AND. lxios_read) THEN 212 CALL iom_init( rxios_context )212 CALL iom_init( crxios_context ) 213 213 IF(lwp) WRITE(numout,*) 'Enable restart reading by XIOS for AGRIF' 214 214 lxios_set = .TRUE. … … 237 237 ! Check dynamics and tracer time-step consistency and force Euler restart if changed 238 238 IF( iom_varid( numror, 'rdt', ldstop = .FALSE. ) > 0 ) THEN 239 CALL iom_get( numror, 'rdt', zrdt, l rxios = lxios_read )239 CALL iom_get( numror, 'rdt', zrdt, ldxios = lxios_read ) 240 240 IF( zrdt /= rdt ) neuler = 0 241 241 ENDIF 242 242 243 243 ! Diurnal DSST 244 IF( ln_diurnal ) CALL iom_get( numror, jpdom_autoglo, 'Dsst' , x_dsst, l rxios = lxios_read )244 IF( ln_diurnal ) CALL iom_get( numror, jpdom_autoglo, 'Dsst' , x_dsst, ldxios = lxios_read ) 245 245 IF ( ln_diurnal_only ) THEN 246 246 IF(lwp) WRITE( numout, * ) & 247 247 & "rst_read:- ln_diurnal_only set, setting rhop=rau0" 248 248 rhop = rau0 249 CALL iom_get( numror, jpdom_autoglo, 'tn' , tsn(:,:,1,jp_tem), l rxios = lxios_read )249 CALL iom_get( numror, jpdom_autoglo, 'tn' , tsn(:,:,1,jp_tem), ldxios = lxios_read ) 250 250 RETURN 251 251 ENDIF 252 252 253 253 IF( iom_varid( numror, 'ub', ldstop = .FALSE. ) > 0 ) THEN 254 CALL iom_get( numror, jpdom_autoglo, 'ub' , ub, l rxios = lxios_read ) ! before fields255 CALL iom_get( numror, jpdom_autoglo, 'vb' , vb, l rxios = lxios_read )256 CALL iom_get( numror, jpdom_autoglo, 'tb' , tsb(:,:,:,jp_tem), l rxios = lxios_read )257 CALL iom_get( numror, jpdom_autoglo, 'sb' , tsb(:,:,:,jp_sal), l rxios = lxios_read )258 CALL iom_get( numror, jpdom_autoglo, 'sshb' , sshb, l rxios = lxios_read )254 CALL iom_get( numror, jpdom_autoglo, 'ub' , ub, ldxios = lxios_read ) ! before fields 255 CALL iom_get( numror, jpdom_autoglo, 'vb' , vb, ldxios = lxios_read ) 256 CALL iom_get( numror, jpdom_autoglo, 'tb' , tsb(:,:,:,jp_tem), ldxios = lxios_read ) 257 CALL iom_get( numror, jpdom_autoglo, 'sb' , tsb(:,:,:,jp_sal), ldxios = lxios_read ) 258 CALL iom_get( numror, jpdom_autoglo, 'sshb' , sshb, ldxios = lxios_read ) 259 259 ELSE 260 260 neuler = 0 261 261 ENDIF 262 262 ! 263 CALL iom_get( numror, jpdom_autoglo, 'un' , un, l rxios = lxios_read ) ! now fields264 CALL iom_get( numror, jpdom_autoglo, 'vn' , vn, l rxios = lxios_read )265 CALL iom_get( numror, jpdom_autoglo, 'tn' , tsn(:,:,:,jp_tem), l rxios = lxios_read )266 CALL iom_get( numror, jpdom_autoglo, 'sn' , tsn(:,:,:,jp_sal), l rxios = lxios_read )267 CALL iom_get( numror, jpdom_autoglo, 'sshn' , sshn, l rxios = lxios_read )263 CALL iom_get( numror, jpdom_autoglo, 'un' , un, ldxios = lxios_read ) ! now fields 264 CALL iom_get( numror, jpdom_autoglo, 'vn' , vn, ldxios = lxios_read ) 265 CALL iom_get( numror, jpdom_autoglo, 'tn' , tsn(:,:,:,jp_tem), ldxios = lxios_read ) 266 CALL iom_get( numror, jpdom_autoglo, 'sn' , tsn(:,:,:,jp_sal), ldxios = lxios_read ) 267 CALL iom_get( numror, jpdom_autoglo, 'sshn' , sshn, ldxios = lxios_read ) 268 268 IF( iom_varid( numror, 'rhop', ldstop = .FALSE. ) > 0 ) THEN 269 CALL iom_get( numror, jpdom_autoglo, 'rhop' , rhop, l rxios = lxios_read ) ! now potential density269 CALL iom_get( numror, jpdom_autoglo, 'rhop' , rhop, ldxios = lxios_read ) ! now potential density 270 270 ELSE 271 271 CALL eos( tsn, rhd, rhop, gdept_n(:,:,:) )
Note: See TracChangeset
for help on using the changeset viewer.