Changeset 5038 for branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90
- Timestamp:
- 2015-01-20T15:26:13+01:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2014/dev_r4621_NOC4_BDY_VERT_INTERP/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90
r4147 r5038 21 21 USE lbclnk ! 22 22 USE prtctl_trc ! Print control for debbuging 23 USE trd mod_oce24 USE trd mod_trc23 USE trd_oce 24 USE trdtrc 25 25 USE iom 26 26 … … 110 110 IF( nn_timing == 1 ) CALL timing_start('p2z_bio') 111 111 ! 112 IF( ln_diatrc ) THEN112 IF( ln_diatrc .OR. lk_iomput ) THEN 113 113 CALL wrk_alloc( jpi, jpj, 17, zw2d ) 114 114 CALL wrk_alloc( jpi, jpj, jpk, 3, zw3d ) … … 122 122 123 123 xksi(:,:) = 0.e0 ! zooplakton closure ( fbod) 124 IF( ln_diatrc ) THEN124 IF( ln_diatrc .OR. lk_iomput ) THEN 125 125 zw2d (:,:,:) = 0.e0 126 126 zw3d(:,:,:,:) = 0.e0 … … 239 239 ! trend number 17 in p2zexp 240 240 ENDIF 241 IF( ln_diatrc ) THEN241 IF( ln_diatrc .OR. lk_iomput ) THEN 242 242 ! convert fluxes in per day 243 243 ze3t = fse3t(ji,jj,jk) * 86400. … … 360 360 ! trend number 17 in p2zexp 361 361 ENDIF 362 IF( ln_diatrc ) THEN362 IF( ln_diatrc .OR. lk_iomput ) THEN 363 363 ! convert fluxes in per day 364 364 ze3t = fse3t(ji,jj,jk) * 86400. … … 390 390 END DO 391 391 392 IF( ln_diatrc ) THEN 393 ! 392 IF( ln_diatrc .OR. lk_iomput ) THEN 394 393 DO jl = 1, 17 395 394 CALL lbc_lnk( zw2d(:,:,jl),'T', 1. ) … … 398 397 CALL lbc_lnk( zw3d(:,:,:,jl),'T', 1. ) 399 398 END DO 400 IF( lk_iomput ) THEN 399 ENDIF 400 IF( lk_iomput ) THEN 401 401 ! Save diagnostics 402 CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 403 CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 404 CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 405 CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 406 CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 407 CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 408 CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 409 CALL iom_put( "TZOODET", zw2d(:,:,8) ) 410 CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 411 CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 412 CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 413 CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 414 CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 415 CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 416 CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 417 CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 418 ! 419 CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 420 CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 421 CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 422 ! 423 ELSE 402 CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 403 CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 404 CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 405 CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 406 CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 407 CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 408 CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 409 CALL iom_put( "TZOODET", zw2d(:,:,8) ) 410 CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 411 CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 412 CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 413 CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 414 CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 415 CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 416 CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 417 CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 418 ! 419 CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 420 CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 421 CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 422 ! 423 ELSE 424 IF( ln_diatrc ) THEN 424 425 ! 425 426 trc2d(:,:,jp_pcs0_2d ) = zw2d(:,:,1) … … 457 458 IF( l_trdtrc ) THEN 458 459 DO jl = jp_pcs0_trd, jp_pcs1_trd 459 CALL trd_ mod_trc( trbio(:,:,:,jl), jl, kt ) ! handle the trend460 CALL trd_trc( trbio(:,:,:,jl), jl, kt ) ! handle the trend 460 461 END DO 461 462 ENDIF … … 467 468 ENDIF 468 469 ! 469 IF( ln_diatrc ) THEN470 IF( ln_diatrc .OR. lk_iomput ) THEN 470 471 CALL wrk_dealloc( jpi, jpj, 17, zw2d ) 471 472 CALL wrk_dealloc( jpi, jpj, jpk, 3, zw3d ) … … 500 501 READ ( numnatp_cfg, namlobphy, IOSTAT = ios, ERR = 902 ) 501 502 902 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobphy in configuration namelist', lwp ) 502 WRITE ( numonp, namlobphy )503 IF(lwm) WRITE ( numonp, namlobphy ) 503 504 504 505 IF(lwp) THEN … … 519 520 READ ( numnatp_cfg, namlobnut, IOSTAT = ios, ERR = 904 ) 520 521 904 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobnut in configuration namelist', lwp ) 521 WRITE ( numonp, namlobnut )522 IF(lwm) WRITE ( numonp, namlobnut ) 522 523 523 524 IF(lwp) THEN … … 537 538 READ ( numnatp_cfg, namlobzoo, IOSTAT = ios, ERR = 906 ) 538 539 906 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobzoo in configuration namelist', lwp ) 539 WRITE ( numonp, namlobzoo )540 IF(lwm) WRITE ( numonp, namlobzoo ) 540 541 541 542 IF(lwp) THEN … … 560 561 READ ( numnatp_cfg, namlobdet, IOSTAT = ios, ERR = 908 ) 561 562 908 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdet in configuration namelist', lwp ) 562 WRITE ( numonp, namlobdet )563 IF(lwm) WRITE ( numonp, namlobdet ) 563 564 564 565 IF(lwp) THEN … … 576 577 READ ( numnatp_cfg, namlobdom, IOSTAT = ios, ERR = 910 ) 577 578 910 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namlobdom in configuration namelist', lwp ) 578 WRITE ( numonp, namlobdom )579 IF(lwm) WRITE ( numonp, namlobdom ) 579 580 580 581 IF(lwp) THEN
Note: See TracChangeset
for help on using the changeset viewer.