- 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/C14b/trcsms_c14b.F90
r1736 r2038 17 17 USE par_trc ! TOP parameters 18 18 USE trc ! TOP variables 19 USE trdm ld_trc_oce20 USE trdm ld_trc19 USE trdmod_oce 20 USE trdmod_trc 21 21 USE iom 22 22 … … 126 126 zpv , & !: piston velocity 127 127 zdemi, ztra 128 #if defined key_ trc_dia3d&& defined key_iomput128 #if defined key_diatrc && defined key_iomput 129 129 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d 130 130 #endif … … 259 259 qtr_c14(ji,jj) = -zpv * zsol * zpco2at & 260 260 & * ( trb(ji,jj,1,jpc14) - zatmbc14(ji,jj) ) & 261 #if defined key_ off_degrad261 #if defined key_degrad 262 262 & * facvol(ji,jj,1) & 263 263 #endif … … 270 270 qint_c14(ji,jj) = qint_c14(ji,jj) + qtr_c14(ji,jj) * rdt 271 271 272 # if defined key_ trc_diaadd&& ! defined key_iomput272 # if defined key_diatrc && ! defined key_iomput 273 273 ! Save 2D diagnostics 274 274 trc2d(ji,jj,jp_c14b0_2d ) = qtr_c14 (ji,jj) … … 282 282 DO jj = 1, jpj 283 283 DO ji = 1, jpi 284 #if ! defined key_ off_degrad284 #if ! defined key_degrad 285 285 ztra = trn(ji,jj,jk,jpc14) * xaccum 286 286 #else … … 288 288 #endif 289 289 tra(ji,jj,jk,jpc14) = tra(ji,jj,jk,jpc14) - ztra / rdt 290 #if defined key_ trc_dia3d290 #if defined key_diatrc 291 291 ! Save 3D diagnostics 292 292 # if ! defined key_iomput … … 300 300 END DO 301 301 302 #if defined key_ trc_diaadd&& defined key_iomput302 #if defined key_diatrc && defined key_iomput 303 303 CALL iom_put( "qtrC14b" , qtr_c14 ) 304 304 CALL iom_put( "qintC14b" , qint_c14 ) 305 305 #endif 306 #if defined key_ trc_dia3d&& defined key_iomput306 #if defined key_diatrc && defined key_iomput 307 307 CALL iom_put( "fdecay" , zw3d ) 308 308 #endif 309 309 IF( l_trdtrc ) THEN 310 CALL trd_mod_trc( tra(:,:,:,jpc14), jpc14, jptr c_trd_sms, kt ) ! save trends310 CALL trd_mod_trc( tra(:,:,:,jpc14), jpc14, jptra_trd_sms, kt ) ! save trends 311 311 END IF 312 312
Note: See TracChangeset
for help on using the changeset viewer.