- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/TOP_SRC/trcrst.F90
r5513 r6808 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 … … 25 21 USE oce_trc 26 22 USE trc 27 USE trcnam_trp28 23 USE iom 29 24 USE daymod 25 30 26 IMPLICIT NONE 31 27 PRIVATE … … 36 32 PUBLIC trc_rst_cal 37 33 38 !! * Substitutions 39 # include "top_substitute.h90" 40 34 !!---------------------------------------------------------------------- 35 !! NEMO/TOP 3.7 , NEMO Consortium (2010) 36 !! $Id$ 37 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 38 !!---------------------------------------------------------------------- 41 39 CONTAINS 42 40 … … 132 130 !!---------------------------------------------------------------------- 133 131 ! 134 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 135 133 ! prognostic variables 136 134 ! -------------------- … … 200 198 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~~' 201 199 202 IF ( jprstlib == jprstdimg ) THEN203 ! eventually read netcdf file (monobloc) for restarting on different number of processors204 ! if {cn_trcrst_in}.nc exists, then set jlibalt to jpnf90205 INQUIRE( FILE = TRIM(cn_trcrst_indir)//'/'//TRIM(cn_trcrst_in)//'.nc', EXIST = llok )206 IF ( llok ) THEN ; jlibalt = jpnf90 ; ELSE ; jlibalt = jprstlib ; ENDIF207 ENDIF208 209 200 IF( ln_rsttr ) THEN 210 201 CALL iom_open( TRIM(cn_trcrst_indir)//'/'//cn_trcrst_in, numrtr, kiolib = jlibalt ) … … 236 227 ELSE 237 228 ndastp = ndate0 - 1 ! ndate0 read in the namelist in dom_nam 238 adatrj = ( REAL( nittrc000-1, wp ) * rdt tra(1)) / rday229 adatrj = ( REAL( nittrc000-1, wp ) * rdt ) / rday 239 230 ! note this is wrong if time step has changed during run 240 231 ENDIF … … 289 280 ! 290 281 DO jk = 1, jpk 291 zvol(:,:,jk) = e1e2t(:,:) * fse3t_a(:,:,jk) * tmask(:,:,jk)282 zvol(:,:,jk) = e1e2t(:,:) * e3t_a(:,:,jk) * tmask(:,:,jk) 292 283 END DO 293 284 !
Note: See TracChangeset
for help on using the changeset viewer.