- 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/LOBSTER/trcbio.F90
r1953 r2038 20 20 USE lbclnk ! 21 21 USE prtctl_trc ! Print control for debbuging 22 USE trdmld_trc 23 USE trdmld_trc_oce 22 USE trdmod_trc 24 23 USE iom 25 24 … … 57 56 !! source sink 58 57 !! 59 !! IF 'key_ trc_diabio' defined , the biogeochemical trends58 !! IF 'key_diabio' defined , the biogeochemical trends 60 59 !! for passive tracers are saved for futher diagnostics. 61 60 !!--------------------------------------------------------------------- … … 71 70 REAL(wp) :: zfilpz, zfildz, zphya, zzooa, zno3a 72 71 REAL(wp) :: znh4a, zdeta, zdoma, zzoobod, zboddet, zdomaju 73 #if defined key_ trc_diaadd72 #if defined key_diatrc 74 73 REAL(wp) :: ze3t 75 74 #endif 76 #if defined key_ trc_diaadd&& defined key_iomput75 #if defined key_diatrc && defined key_iomput 77 76 REAL(wp), DIMENSION(jpi,jpj,17) :: zw2d 78 # if defined key_trc_dia3d79 77 REAL(wp), DIMENSION(jpi,jpj,jpk,3) :: zw3d 80 # endif81 78 #endif 82 79 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: ztrbio … … 91 88 92 89 fbod(:,:) = 0.e0 93 #if defined key_ trc_diaadd&& ! defined key_iomput90 #if defined key_diatrc && ! defined key_iomput 94 91 DO jl = jp_lob0_2d, jp_lob1_2d 95 92 trc2d(:,:,jl) = 0.e0 96 93 END DO 97 94 #endif 98 #if defined key_ trc_diaadd&& defined key_iomput95 #if defined key_diatrc && defined key_iomput 99 96 zw2d(:,:,:) = 0.e0 100 # if defined key_trc_dia3d101 97 zw3d(:,:,:,:) = 0.e0 102 # endif103 98 #endif 104 99 … … 117 112 118 113 ! negative trophic variables DO not contribute to the fluxes 119 zdet = MAX( 0.e0, trn(ji,jj,jk,jp det) )120 zzoo = MAX( 0.e0, trn(ji,jj,jk,jp zoo) )121 zphy = MAX( 0.e0, trn(ji,jj,jk,jp phy) )122 zno3 = MAX( 0.e0, trn(ji,jj,jk,jp no3) )123 znh4 = MAX( 0.e0, trn(ji,jj,jk,jp nh4) )124 zdom = MAX( 0.e0, trn(ji,jj,jk,jp dom) )114 zdet = MAX( 0.e0, trn(ji,jj,jk,jp_lob_det) ) 115 zzoo = MAX( 0.e0, trn(ji,jj,jk,jp_lob_zoo) ) 116 zphy = MAX( 0.e0, trn(ji,jj,jk,jp_lob_phy) ) 117 zno3 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_no3) ) 118 znh4 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_nh4) ) 119 zdom = MAX( 0.e0, trn(ji,jj,jk,jp_lob_dom) ) 125 120 126 121 ! Limitations … … 194 189 195 190 ! tracer flux at totox-point added to the general trend 196 tra(ji,jj,jk,jp det) = tra(ji,jj,jk,jpdet) + zdeta197 tra(ji,jj,jk,jp zoo) = tra(ji,jj,jk,jpzoo) + zzooa198 tra(ji,jj,jk,jp phy) = tra(ji,jj,jk,jpphy) + zphya199 tra(ji,jj,jk,jp no3) = tra(ji,jj,jk,jpno3) + zno3a200 tra(ji,jj,jk,jp nh4) = tra(ji,jj,jk,jpnh4) + znh4a201 tra(ji,jj,jk,jp dom) = tra(ji,jj,jk,jpdom) + zdoma202 203 #if defined key_ trc_diabio191 tra(ji,jj,jk,jp_lob_det) = tra(ji,jj,jk,jp_lob_det) + zdeta 192 tra(ji,jj,jk,jp_lob_zoo) = tra(ji,jj,jk,jp_lob_zoo) + zzooa 193 tra(ji,jj,jk,jp_lob_phy) = tra(ji,jj,jk,jp_lob_phy) + zphya 194 tra(ji,jj,jk,jp_lob_no3) = tra(ji,jj,jk,jp_lob_no3) + zno3a 195 tra(ji,jj,jk,jp_lob_nh4) = tra(ji,jj,jk,jp_lob_nh4) + znh4a 196 tra(ji,jj,jk,jp_lob_dom) = tra(ji,jj,jk,jp_lob_dom) + zdoma 197 198 #if defined key_diabio 204 199 trbio(ji,jj,jk,jp_lob0_trd ) = zno3phy 205 200 trbio(ji,jj,jk,jp_lob0_trd + 1) = znh4phy … … 238 233 ENDIF 239 234 240 #if defined key_ trc_diaadd235 #if defined key_diatrc 241 236 ! convert fluxes in per day 242 237 ze3t = fse3t(ji,jj,jk) * 86400. … … 282 277 zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t 283 278 #endif 284 #if defined key_ trc_dia3d279 #if defined key_diatrc 285 280 # if ! defined key_iomput 286 281 trc3d(ji,jj,jk,jp_lob0_3d ) = zno3phy * 86400 … … 307 302 ! trophic variables( det, zoo, phy, no3, nh4, dom) 308 303 ! negative trophic variables DO not contribute to the fluxes 309 zdet = MAX( 0.e0, trn(ji,jj,jk,jp det) )310 zzoo = MAX( 0.e0, trn(ji,jj,jk,jp zoo) )311 zphy = MAX( 0.e0, trn(ji,jj,jk,jp phy) )312 zno3 = MAX( 0.e0, trn(ji,jj,jk,jp no3) )313 znh4 = MAX( 0.e0, trn(ji,jj,jk,jp nh4) )314 zdom = MAX( 0.e0, trn(ji,jj,jk,jp dom) )304 zdet = MAX( 0.e0, trn(ji,jj,jk,jp_lob_det) ) 305 zzoo = MAX( 0.e0, trn(ji,jj,jk,jp_lob_zoo) ) 306 zphy = MAX( 0.e0, trn(ji,jj,jk,jp_lob_phy) ) 307 zno3 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_no3) ) 308 znh4 = MAX( 0.e0, trn(ji,jj,jk,jp_lob_nh4) ) 309 zdom = MAX( 0.e0, trn(ji,jj,jk,jp_lob_dom) ) 315 310 316 311 ! Limitations … … 363 358 364 359 ! tracer flux at totox-point added to the general trend 365 tra(ji,jj,jk,jp det) = tra(ji,jj,jk,jpdet) + zdeta366 tra(ji,jj,jk,jp zoo) = tra(ji,jj,jk,jpzoo) + zzooa367 tra(ji,jj,jk,jp phy) = tra(ji,jj,jk,jpphy) + zphya368 tra(ji,jj,jk,jp no3) = tra(ji,jj,jk,jpno3) + zno3a369 tra(ji,jj,jk,jp nh4) = tra(ji,jj,jk,jpnh4) + znh4a370 tra(ji,jj,jk,jp dom) = tra(ji,jj,jk,jpdom) + zdoma360 tra(ji,jj,jk,jp_lob_det) = tra(ji,jj,jk,jp_lob_det) + zdeta 361 tra(ji,jj,jk,jp_lob_zoo) = tra(ji,jj,jk,jp_lob_zoo) + zzooa 362 tra(ji,jj,jk,jp_lob_phy) = tra(ji,jj,jk,jp_lob_phy) + zphya 363 tra(ji,jj,jk,jp_lob_no3) = tra(ji,jj,jk,jp_lob_no3) + zno3a 364 tra(ji,jj,jk,jp_lob_nh4) = tra(ji,jj,jk,jp_lob_nh4) + znh4a 365 tra(ji,jj,jk,jp_lob_dom) = tra(ji,jj,jk,jp_lob_dom) + zdoma 371 366 ! 372 #if defined key_ trc_diabio367 #if defined key_diabio 373 368 trbio(ji,jj,jk,jp_lob0_trd ) = zno3phy 374 369 trbio(ji,jj,jk,jp_lob0_trd + 1) = znh4phy … … 406 401 ! trend number 17 in trcexp 407 402 ENDIF 408 #if defined key_ trc_diaadd && defined key_trc_dia3d403 #if defined key_diatrc 409 404 # if ! defined key_iomput 410 405 trc3d(ji,jj,jk,jp_lob0_3d ) = zno3phy * 86400 … … 421 416 END DO 422 417 423 #if defined key_ trc_diaadd418 #if defined key_diatrc 424 419 ! Lateral boundary conditions 425 420 # if ! defined key_iomput … … 452 447 #endif 453 448 454 #if defined key_ trc_diaadd && defined key_trc_dia3d449 #if defined key_diatrc 455 450 ! Lateral boundary conditions 456 451 # if ! defined key_iomput … … 469 464 #endif 470 465 471 #if defined key_ trc_diabio466 #if defined key_diabio 472 467 ! Lateral boundary conditions on trcbio 473 468 DO jl = jp_lob0_trd, jp_lob1_trd
Note: See TracChangeset
for help on using the changeset viewer.