- Timestamp:
- 2016-11-28T17:04:10+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_INGV_UKMO_2016/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r5836 r7351 14 14 !!---------------------------------------------------------------------- 15 15 !!---------------------------------------------------------------------- 16 !! trc_rst : Restart for passive tracer 17 !!---------------------------------------------------------------------- 18 !!---------------------------------------------------------------------- 19 !! 'key_top' TOP models 20 !!---------------------------------------------------------------------- 16 !! trc_rst : Restart for passive tracer 21 17 !! trc_rst_opn : open restart file 22 18 !! trc_rst_read : read restart file … … 27 23 USE iom 28 24 USE daymod 25 29 26 IMPLICIT NONE 30 27 PRIVATE … … 35 32 PUBLIC trc_rst_cal 36 33 37 !! * Substitutions 38 # include "domzgr_substitute.h90" 39 34 !!---------------------------------------------------------------------- 35 !! NEMO/TOP 3.7 , NEMO Consortium (2010) 36 !! $Id$ 37 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 38 !!---------------------------------------------------------------------- 40 39 CONTAINS 41 40 … … 131 130 !!---------------------------------------------------------------------- 132 131 ! 133 CALL iom_rstput( kt, nitrst, numrtw, 'rdttrc1', rdttrc (1) ) ! surfacepassive tracer time step132 CALL iom_rstput( kt, nitrst, numrtw, 'rdttrc1', rdttrc ) ! passive tracer time step 134 133 ! prognostic variables 135 134 ! -------------------- … … 199 198 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 200 199 201 IF ( jprstlib == jprstdimg ) THEN202 ! eventually read netcdf file (monobloc) for restarting on different number of processors203 ! if {cn_trcrst_in}.nc exists, then set jlibalt to jpnf90204 INQUIRE( FILE = TRIM(cn_trcrst_indir)//'/'//TRIM(cn_trcrst_in)//'.nc', EXIST = llok )205 IF ( llok ) THEN ; jlibalt = jpnf90 ; ELSE ; jlibalt = jprstlib ; ENDIF206 ENDIF207 208 200 IF( ln_rsttr ) THEN 209 201 CALL iom_open( TRIM(cn_trcrst_indir)//'/'//cn_trcrst_in, numrtr, kiolib = jlibalt ) … … 235 227 ELSE 236 228 ndastp = ndate0 - 1 ! ndate0 read in the namelist in dom_nam 237 adatrj = ( REAL( nittrc000-1, wp ) * rdt tra(1)) / rday229 adatrj = ( REAL( nittrc000-1, wp ) * rdt ) / rday 238 230 ! note this is wrong if time step has changed during run 239 231 ENDIF … … 288 280 ! 289 281 DO jk = 1, jpk 290 zvol(:,:,jk) = e1e2t(:,:) * fse3t_a(:,:,jk) * tmask(:,:,jk)282 zvol(:,:,jk) = e1e2t(:,:) * e3t_a(:,:,jk) * tmask(:,:,jk) 291 283 END DO 292 284 ! … … 303 295 IF(lwp) WRITE(numout,9000) jn, TRIM( ctrcnm(jn) ), zmean, zmin, zmax, zdrift 304 296 END DO 305 WRITE(numout,*)297 IF(lwp) WRITE(numout,*) 306 298 9000 FORMAT(' tracer nb :',i2,' name :',a10,' mean :',e18.10,' min :',e18.10, & 307 299 & ' max :',e18.10,' drift :',e18.10, ' %')
Note: See TracChangeset
for help on using the changeset viewer.