Changeset 1459 for trunk/NEMO/TOP_SRC/C14b
- Timestamp:
- 2009-05-23T10:53:26+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/C14b/trcsms_c14b.F90
r1329 r1459 19 19 USE trdmld_trc_oce 20 20 USE trdmld_trc 21 USE iom 21 22 22 23 IMPLICIT NONE … … 125 126 zpv , & !: piston velocity 126 127 zdemi, ztra 128 #if defined key_trc_dia3d && defined key_iomput 129 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 130 #endif 127 131 !!---------------------------------------------------------------------- 128 132 … … 266 270 qint_c14(ji,jj) = qint_c14(ji,jj) + qtr_c14(ji,jj) * rdt 267 271 268 # if defined key_trc_diaadd 272 # if defined key_trc_diaadd && ! defined key_iomput 269 273 ! Save 2D diagnostics 270 274 trc2d(ji,jj,jp_c14b0_2d ) = qtr_c14 (ji,jj) … … 284 288 #endif 285 289 tra(ji,jj,jk,jpc14) = tra(ji,jj,jk,jpc14) - ztra / rdt 286 287 # if defined key_trc_dia3d 290 #if defined key_trc_dia3d 288 291 ! Save 3D diagnostics 292 # if ! defined key_iomput 289 293 trc3d(ji,jj,jk,jp_c14b0_3d ) = ztra ! radioactive decay 294 # else 295 zw3d(ji,jj,jk) = ztra ! radioactive decay 290 296 # endif 297 #endif 291 298 END DO 292 299 END DO 293 300 END DO 294 301 302 #if defined key_trc_diaadd && defined key_iomput 303 CALL iom_put( "qtrC14b" , qtr_c14 ) 304 CALL iom_put( "qintC14b" , qint_c14 ) 305 #endif 306 #if defined key_trc_dia3d && defined key_iomput 307 CALL iom_put( "fdecay" , zw3d ) 308 #endif 295 309 IF( l_trdtrc ) THEN 296 310 CALL trd_mod_trc( tra(:,:,:,jpc14), jpc14, jptrc_trd_sms, kt ) ! save trends
Note: See TracChangeset
for help on using the changeset viewer.