Changeset 2038 for branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/trcrst.F90
- Timestamp:
- 2010-08-02T12:57:40+02:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/DEV_r2006_merge_TRA_TRC/NEMO/TOP_SRC/trcrst.F90
r1970 r2038 25 25 USE oce_trc 26 26 USE trc 27 USE trc trp_lec27 USE trcnam_trp 28 28 USE lib_mpp 29 29 USE iom … … 33 33 USE trcrst_c14b ! C14 bomb restart 34 34 USE trcrst_my_trc ! MY_TRC restart 35 #if defined key_off _tra35 #if defined key_offline 36 36 USE daymod 37 37 #endif … … 67 67 !!---------------------------------------------------------------------- 68 68 ! 69 # if ! defined key_off _tra69 # if ! defined key_offline 70 70 IF( kt == nit000 ) lrst_trc = .FALSE. 71 71 # else … … 82 82 # endif 83 83 ! to get better performances with NetCDF format: 84 ! we open and define the tracer restart file one tracer time step before writing the data (-> at nitrst - 2*n dttrc + 1)85 ! except if we write tracer restart files every tracer time step or if a tracer restart file was writen at nitend - 2*n dttrc + 186 IF( kt == nitrst - 2*n dttrc + 1 .OR. nstock == ndttrc .OR. ( kt == nitend - ndttrc + 1 .AND. .NOT. lrst_trc ) ) THEN84 ! we open and define the tracer restart file one tracer time step before writing the data (-> at nitrst - 2*nn_dttrc + 1) 85 ! 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 86 IF( kt == nitrst - 2*nn_dttrc + 1 .OR. nstock == nn_dttrc .OR. ( kt == nitend - nn_dttrc + 1 .AND. .NOT. lrst_trc ) ) THEN 87 87 ! beware of the format used to write kt (default is i8.8, that should be large enough) 88 88 IF( nitrst > 1.0e9 ) THEN ; WRITE(clkt,* ) nitrst … … 222 222 !! 223 223 !! According to namelist parameter nrstdt, 224 !! n rsttr = 0 no control on the date (nittrc000 is arbitrary).225 !! n rsttr = 1 we verify that nit000 is equal to the last224 !! nn_rsttr = 0 no control on the date (nittrc000 is arbitrary). 225 !! nn_rsttr = 1 we verify that nit000 is equal to the last 226 226 !! time step of previous run + 1. 227 227 !! In both those options, the exact duration of the experiment … … 230 230 !! This is valid is the time step has remained constant. 231 231 !! 232 !! n rsttr = 2 the duration of the experiment in days (adatrj)232 !! nn_rsttr = 2 the duration of the experiment in days (adatrj) 233 233 !! has been stored in the restart file. 234 234 !!---------------------------------------------------------------------- … … 237 237 ! 238 238 REAL(wp) :: zkt 239 #if defined key_off _tra239 #if defined key_offline 240 240 REAL(wp) :: zndastp 241 241 #endif … … 250 250 WRITE(numout,*) ' previous time-step : ', NINT( zkt ) 251 251 WRITE(numout,*) ' *** restart option' 252 SELECT CASE ( n rsttr )253 CASE ( 0 ) ; WRITE(numout,*) ' n rsttr = 0 : no control of nittrc000'254 CASE ( 1 ) ; WRITE(numout,*) ' n rsttr = 1 : no control the date at nit000 (use ndate0 read in the namelist)'255 CASE ( 2 ) ; WRITE(numout,*) ' n rsttr = 2 : calendar parameters read in restart'252 SELECT CASE ( nn_rsttr ) 253 CASE ( 0 ) ; WRITE(numout,*) ' nn_rsttr = 0 : no control of nittrc000' 254 CASE ( 1 ) ; WRITE(numout,*) ' nn_rsttr = 1 : no control the date at nit000 (use ndate0 read in the namelist)' 255 CASE ( 2 ) ; WRITE(numout,*) ' nn_rsttr = 2 : calendar parameters read in restart' 256 256 END SELECT 257 257 WRITE(numout,*) 258 258 ENDIF 259 259 ! Control of date 260 IF( nittrc000 - NINT( zkt ) /= 1 .AND. n rsttr /= 0 ) &260 IF( nittrc000 - NINT( zkt ) /= 1 .AND. nn_rsttr /= 0 ) & 261 261 & CALL ctl_stop( ' ===>>>> : problem with nit000 for the restart', & 262 & ' verify the restart file or rerun with n rsttr = 0 (namelist)' )263 #if defined key_off _tra262 & ' verify the restart file or rerun with nn_rsttr = 0 (namelist)' ) 263 #if defined key_offline 264 264 ! define ndastp and adatrj 265 IF ( n rsttr == 2 ) THEN265 IF ( nn_rsttr == 2 ) THEN 266 266 CALL iom_get( numrtr, 'ndastp', zndastp ) 267 267 ndastp = NINT( zndastp ) … … 329 329 DO ji = 1, jpi 330 330 zvol = cvol(ji,jj,jk) 331 # if defined key_ off_degrad331 # if defined key_degrad 332 332 zvol = zvol * facvol(ji,jj,jk) 333 333 # endif
Note: See TracChangeset
for help on using the changeset viewer.