Changeset 1457 for trunk/NEMO/TOP_SRC/LOBSTER/trcbio.F90
- Timestamp:
- 2009-05-23T10:16:38+02:00 (15 years ago)
- File:
-
- 1 edited
-
trunk/NEMO/TOP_SRC/LOBSTER/trcbio.F90 (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/LOBSTER/trcbio.F90
r1264 r1457 22 22 USE trdmld_trc 23 23 USE trdmld_trc_oce 24 USE iom 24 25 25 26 IMPLICIT NONE … … 29 30 30 31 !!* Substitution 31 # include " domzgr_substitute.h90"32 # include "top_substitute.h90" 32 33 !!---------------------------------------------------------------------- 33 34 !! NEMO/TOP 2.0 , LOCEAN-IPSL (2007) … … 71 72 REAL(wp) :: znh4a, zdeta, zdoma, zzoobod, zboddet, zdomaju 72 73 #if defined key_trc_diaadd 73 REAL(wp), DIMENSION(jpi,jpj,jpk) :: ze3t 74 REAL(wp) :: ze3t 75 #endif 76 #if defined key_trc_diaadd && defined key_iomput 77 REAL(wp), DIMENSION(jpi,jpj,17) :: zw2d 78 # if defined key_trc_dia3d 79 REAL(wp), DIMENSION(jpi,jpj,jpk,3) :: zw3d 80 # endif 74 81 #endif 75 82 REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE :: ztrbio … … 83 90 ENDIF 84 91 85 #if defined key_trc_diaadd86 ! convert fluxes in per day87 ze3t(:,:,:) = 0.e088 DO jk = 1, jpkbm189 ze3t(:,:,jk) = fse3t(:,:,jk) * 86400.90 END DO91 #endif92 93 92 fbod(:,:) = 0.e0 94 #if defined key_trc_diaadd 93 #if defined key_trc_diaadd && ! defined key_iomput 95 94 DO jl = jp_lob0_2d, jp_lob1_2d 96 95 trc2d(:,:,jl) = 0.e0 97 96 END DO 97 #endif 98 #if defined key_trc_diaadd && defined key_iomput 99 zw2d(:,:,:) = 0.e0 100 # if defined key_trc_dia3d 101 zw3d(:,:,:,:) = 0.e0 102 # endif 98 103 #endif 99 104 … … 107 112 ! ! -------------------------- ! 108 113 DO jj = 2, jpjm1 109 DO ji = 2, jpim1 !!gm use of fs_2 fs_jpm1 required here 110 114 DO ji = fs_2, fs_jpim1 111 115 ! trophic variables( det, zoo, phy, no3, nh4, dom) 112 116 ! ------------------------------------------------ … … 126 130 zlno3 = zno3 * EXP( -psinut * znh4 ) / ( akno3 + zno3 ) 127 131 zlnh4 = znh4 / (znh4+aknh4) 128 129 132 130 133 ! sinks and sources … … 189 192 zdeta = zphydet + zzoodet - zdetzoo - zdetnh4 - zdetdom + zboddet 190 193 zdoma = zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju 194 195 ! tracer flux at totox-point added to the general trend 196 tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta 197 tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa 198 tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya 199 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a 200 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a 201 tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma 191 202 192 203 #if defined key_trc_diabio … … 206 217 trbio(ji,jj,jk,jp_lob0_trd + 14) = zdetnh4 207 218 trbio(ji,jj,jk,jp_lob0_trd + 15) = zdetdom 208 #endif209 #if defined key_trc_diaadd210 trc2d(ji,jj,jp_lob0_2d ) = trc2d(ji,jj, jp_lob0_2d ) + zno3phy * ze3t(ji,jj,jk)211 trc2d(ji,jj,jp_lob0_2d + 1) = trc2d(ji,jj, jp_lob0_2d + 1) + znh4phy * ze3t(ji,jj,jk)212 trc2d(ji,jj,jp_lob0_2d + 2) = trc2d(ji,jj, jp_lob0_2d + 2) + zphydom * ze3t(ji,jj,jk)213 trc2d(ji,jj,jp_lob0_2d + 3) = trc2d(ji,jj, jp_lob0_2d + 3) + zphynh4 * ze3t(ji,jj,jk)214 trc2d(ji,jj,jp_lob0_2d + 4) = trc2d(ji,jj, jp_lob0_2d + 4) + zphyzoo * ze3t(ji,jj,jk)215 trc2d(ji,jj,jp_lob0_2d + 5) = trc2d(ji,jj, jp_lob0_2d + 5) + zphydet * ze3t(ji,jj,jk)216 trc2d(ji,jj,jp_lob0_2d + 6) = trc2d(ji,jj, jp_lob0_2d + 6) + zdetzoo * ze3t(ji,jj,jk)217 ! trend number 8 is in trcsed.F218 trc2d(ji,jj,jp_lob0_2d + 8) = trc2d(ji,jj,jp_lob0_2d + 8) + zzoodet * ze3t(ji,jj,jk)219 trc2d(ji,jj,jp_lob0_2d + 9) = trc2d(ji,jj,jp_lob0_2d + 9) + zzoobod * ze3t(ji,jj,jk)220 trc2d(ji,jj,jp_lob0_2d + 10) = trc2d(ji,jj,jp_lob0_2d + 10) + zzoonh4 * ze3t(ji,jj,jk)221 trc2d(ji,jj,jp_lob0_2d + 11) = trc2d(ji,jj,jp_lob0_2d + 11) + zzoodom * ze3t(ji,jj,jk)222 trc2d(ji,jj,jp_lob0_2d + 12) = trc2d(ji,jj,jp_lob0_2d + 12) + znh4no3 * ze3t(ji,jj,jk)223 trc2d(ji,jj,jp_lob0_2d + 13) = trc2d(ji,jj,jp_lob0_2d + 13) + zdomnh4 * ze3t(ji,jj,jk)224 trc2d(ji,jj,jp_lob0_2d + 14) = trc2d(ji,jj,jp_lob0_2d + 14) + zdetnh4 * ze3t(ji,jj,jk)225 trc2d(ji,jj,jp_lob0_2d + 15) = trc2d(ji,jj,jp_lob0_2d + 15) + ( zno3phy + znh4phy - zphynh4 &226 & - zphydom - zphyzoo - zphydet ) * ze3t(ji,jj,jk)227 trc2d(ji,jj,jp_lob0_2d + 16) = trc2d(ji,jj,jp_lob0_2d + 16) + ( zphyzoo + zdetzoo - zzoodet &228 & - zzoobod - zzoonh4 - zzoodom ) * ze3t(ji,jj,jk)229 trc2d(ji,jj,jp_lob0_2d + 17) = trc2d(ji,jj,jp_lob0_2d + 17) + zdetdom * ze3t(ji,jj,jk)230 ! trend number 19 is in trcexp.F231 trc3d(ji,jj,jk,jp_lob0_3d ) = zno3phy * 86400232 trc3d(ji,jj,jk,jp_lob0_3d + 1) = znh4phy * 86400233 trc3d(ji,jj,jk,jp_lob0_3d + 2) = znh4no3 * 86400234 219 #endif 235 220 IF( l_trdtrc ) THEN … … 253 238 ENDIF 254 239 255 256 ! tracer flux at totox-point added to the general trend 257 tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta 258 tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa 259 tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya 260 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a 261 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a 262 tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma 263 240 #if defined key_trc_diaadd 241 ! convert fluxes in per day 242 ze3t = fse3t(ji,jj,jk) * 86400. 243 #if ! defined key_iomput 244 trc2d(ji,jj,jp_lob0_2d ) = trc2d(ji,jj, jp_lob0_2d ) + zno3phy * ze3t 245 trc2d(ji,jj,jp_lob0_2d + 1) = trc2d(ji,jj, jp_lob0_2d + 1) + znh4phy * ze3t 246 trc2d(ji,jj,jp_lob0_2d + 2) = trc2d(ji,jj, jp_lob0_2d + 2) + zphydom * ze3t 247 trc2d(ji,jj,jp_lob0_2d + 3) = trc2d(ji,jj, jp_lob0_2d + 3) + zphynh4 * ze3t 248 trc2d(ji,jj,jp_lob0_2d + 4) = trc2d(ji,jj, jp_lob0_2d + 4) + zphyzoo * ze3t 249 trc2d(ji,jj,jp_lob0_2d + 5) = trc2d(ji,jj, jp_lob0_2d + 5) + zphydet * ze3t 250 trc2d(ji,jj,jp_lob0_2d + 6) = trc2d(ji,jj, jp_lob0_2d + 6) + zdetzoo * ze3t 251 ! trend number 8 is in trcsed.F 252 trc2d(ji,jj,jp_lob0_2d + 8) = trc2d(ji,jj,jp_lob0_2d + 8) + zzoodet * ze3t 253 trc2d(ji,jj,jp_lob0_2d + 9) = trc2d(ji,jj,jp_lob0_2d + 9) + zzoobod * ze3t 254 trc2d(ji,jj,jp_lob0_2d + 10) = trc2d(ji,jj,jp_lob0_2d + 10) + zzoonh4 * ze3t 255 trc2d(ji,jj,jp_lob0_2d + 11) = trc2d(ji,jj,jp_lob0_2d + 11) + zzoodom * ze3t 256 trc2d(ji,jj,jp_lob0_2d + 12) = trc2d(ji,jj,jp_lob0_2d + 12) + znh4no3 * ze3t 257 trc2d(ji,jj,jp_lob0_2d + 13) = trc2d(ji,jj,jp_lob0_2d + 13) + zdomnh4 * ze3t 258 trc2d(ji,jj,jp_lob0_2d + 14) = trc2d(ji,jj,jp_lob0_2d + 14) + zdetnh4 * ze3t 259 trc2d(ji,jj,jp_lob0_2d + 15) = trc2d(ji,jj,jp_lob0_2d + 15) + ( zno3phy + znh4phy - zphynh4 & 260 & - zphydom - zphyzoo - zphydet ) * ze3t 261 trc2d(ji,jj,jp_lob0_2d + 16) = trc2d(ji,jj,jp_lob0_2d + 16) + ( zphyzoo + zdetzoo - zzoodet & 262 & - zzoobod - zzoonh4 - zzoodom ) * ze3t 263 trc2d(ji,jj,jp_lob0_2d + 17) = trc2d(ji,jj,jp_lob0_2d + 17) + zdetdom * ze3t 264 ! trend number 19 is in trcexp.F 265 #else 266 zw2d(ji,jj,1) = zw2d(ji,jj,1) + zno3phy * ze3t 267 zw2d(ji,jj,2) = zw2d(ji,jj,2) + znh4phy * ze3t 268 zw2d(ji,jj,3) = zw2d(ji,jj,3) + zphydom * ze3t 269 zw2d(ji,jj,4) = zw2d(ji,jj,4) + zphynh4 * ze3t 270 zw2d(ji,jj,5) = zw2d(ji,jj,5) + zphyzoo * ze3t 271 zw2d(ji,jj,6) = zw2d(ji,jj,6) + zphydet * ze3t 272 zw2d(ji,jj,7) = zw2d(ji,jj,7) + zdetzoo * ze3t 273 zw2d(ji,jj,8) = zw2d(ji,jj,8) + zzoodet * ze3t 274 zw2d(ji,jj,9) = zw2d(ji,jj,9) + zzoobod * ze3t 275 zw2d(ji,jj,10) = zw2d(ji,jj,10) + zzoonh4 * ze3t 276 zw2d(ji,jj,11) = zw2d(ji,jj,11) + zzoodom * ze3t 277 zw2d(ji,jj,12) = zw2d(ji,jj,12) + znh4no3 * ze3t 278 zw2d(ji,jj,13) = zw2d(ji,jj,13) + zdomnh4 * ze3t 279 zw2d(ji,jj,14) = zw2d(ji,jj,14) + zdetnh4 * ze3t 280 zw2d(ji,jj,15) = zw2d(ji,jj,15) + ( zno3phy + znh4phy - zphynh4 - zphydom - zphyzoo - zphydet ) * ze3t 281 zw2d(ji,jj,16) = zw2d(ji,jj,16) + ( zphyzoo + zdetzoo - zzoodet - zzoobod - zzoonh4 - zzoodom ) * ze3t 282 zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t 283 #endif 284 #if defined key_trc_dia3d 285 # if ! defined key_iomput 286 trc3d(ji,jj,jk,jp_lob0_3d ) = zno3phy * 86400 287 trc3d(ji,jj,jk,jp_lob0_3d + 1) = znh4phy * 86400 288 trc3d(ji,jj,jk,jp_lob0_3d + 2) = znh4no3 * 86400 289 # else 290 zw3d(ji,jj,jk,1) = zno3phy * 86400 291 zw3d(ji,jj,jk,2) = znh4phy * 86400 292 zw3d(ji,jj,jk,3) = znh4no3 * 86400 293 # endif 294 #endif 295 #endif 264 296 END DO 265 297 END DO 266 298 END DO 267 299 268 !269 270 !!gm do loop until jpkm1 only!271 300 ! ! -------------------------- ! 272 DO jk = jpkb, jpk 301 DO jk = jpkb, jpkm1 ! Upper ocean (bio-layers) ! 273 302 ! ! -------------------------- ! 274 275 303 DO jj = 2, jpjm1 276 DO ji = 2,jpim1 !!gm use of fs_2 & fs_jpim1 required 277 304 DO ji = fs_2, fs_jpim1 278 305 ! remineralisation of all quantities towards nitrate 279 306 … … 335 362 zdoma = zphydom + zzoodom + zdetdom - zdomnh4 - zdomaju 336 363 364 ! tracer flux at totox-point added to the general trend 365 tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta 366 tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa 367 tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya 368 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a 369 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a 370 tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma 371 ! 337 372 #if defined key_trc_diabio 338 373 trbio(ji,jj,jk,jp_lob0_trd ) = zno3phy … … 351 386 trbio(ji,jj,jk,jp_lob0_trd + 14) = zdetnh4 352 387 trbio(ji,jj,jk,jp_lob0_trd + 15) = zdetdom 353 #endif354 #if defined key_trc_diaadd355 trc2d(ji,jj, jp_lob0_2d ) = trc2d(ji,jj, jp_lob0_2d ) + zno3phy * ze3t(ji,jj,jk)356 trc2d(ji,jj, jp_lob0_2d + 1) = trc2d(ji,jj, jp_lob0_2d + 1) + znh4phy * ze3t(ji,jj,jk)357 trc2d(ji,jj, jp_lob0_2d + 2) = trc2d(ji,jj, jp_lob0_2d + 2) + zphydom * ze3t(ji,jj,jk)358 trc2d(ji,jj, jp_lob0_2d + 3) = trc2d(ji,jj, jp_lob0_2d + 3) + zphynh4 * ze3t(ji,jj,jk)359 trc2d(ji,jj, jp_lob0_2d + 4) = trc2d(ji,jj, jp_lob0_2d + 4) + zphyzoo * ze3t(ji,jj,jk)360 trc2d(ji,jj, jp_lob0_2d + 5) = trc2d(ji,jj, jp_lob0_2d + 5) + zphydet * ze3t(ji,jj,jk)361 trc2d(ji,jj, jp_lob0_2d + 6) = trc2d(ji,jj, jp_lob0_2d + 6) + zdetzoo * ze3t(ji,jj,jk)362 ! trend number 8 is in trcsed.F363 trc2d(ji,jj,jp_lob0_2d + 8) = trc2d(ji,jj,jp_lob0_2d + 8) + zzoodet * ze3t(ji,jj,jk)364 trc2d(ji,jj,jp_lob0_2d + 9) = trc2d(ji,jj,jp_lob0_2d + 9) + zzoobod * ze3t(ji,jj,jk)365 trc2d(ji,jj,jp_lob0_2d + 10) = trc2d(ji,jj,jp_lob0_2d + 10) + zzoonh4 * ze3t(ji,jj,jk)366 trc2d(ji,jj,jp_lob0_2d + 11) = trc2d(ji,jj,jp_lob0_2d + 11) + zzoodom * ze3t(ji,jj,jk)367 trc2d(ji,jj,jp_lob0_2d + 12) = trc2d(ji,jj,jp_lob0_2d + 12) + znh4no3 * ze3t(ji,jj,jk)368 trc2d(ji,jj,jp_lob0_2d + 13) = trc2d(ji,jj,jp_lob0_2d + 13) + zdomnh4 * ze3t(ji,jj,jk)369 trc2d(ji,jj,jp_lob0_2d + 14) = trc2d(ji,jj,jp_lob0_2d + 14) + zdetnh4 * ze3t(ji,jj,jk)370 371 trc2d(ji,jj,jp_lob0_2d + 15) = trc2d(ji,jj,jp_lob0_2d + 15) + ( zno3phy + znh4phy - zphynh4 &372 & - zphydom - zphyzoo - zphydet ) * ze3t(ji,jj,jk)373 trc2d(ji,jj,jp_lob0_2d + 16) = trc2d(ji,jj,jp_lob0_2d + 16) + ( zphyzoo + zdetzoo - zzoodet &374 & - zzoobod - zzoonh4 - zzoodom ) * ze3t(ji,jj,jk)375 trc2d(ji,jj,jp_lob0_2d + 17) = trc2d(ji,jj,jp_lob0_2d + 17) + zdetdom * ze3t(ji,jj,jk)376 377 trc3d(ji,jj,jk,jp_lob0_3d ) = zno3phy * 86400378 trc3d(ji,jj,jk,jp_lob0_3d + 1) = znh4phy * 86400379 trc3d(ji,jj,jk,jp_lob0_3d + 2) = znh4no3 * 86400380 388 #endif 381 389 IF( l_trdtrc ) THEN … … 398 406 ! trend number 17 in trcexp 399 407 ENDIF 400 401 402 ! tracer flux at totox-point added to the general trend 403 tra(ji,jj,jk,jpdet) = tra(ji,jj,jk,jpdet) + zdeta 404 tra(ji,jj,jk,jpzoo) = tra(ji,jj,jk,jpzoo) + zzooa 405 tra(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zphya 406 tra(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) + zno3a 407 tra(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) + znh4a 408 tra(ji,jj,jk,jpdom) = tra(ji,jj,jk,jpdom) + zdoma 409 ! 408 #if defined key_trc_diaadd && defined key_trc_dia3d 409 # if ! defined key_iomput 410 trc3d(ji,jj,jk,jp_lob0_3d ) = zno3phy * 86400 411 trc3d(ji,jj,jk,jp_lob0_3d + 1) = znh4phy * 86400 412 trc3d(ji,jj,jk,jp_lob0_3d + 2) = znh4no3 * 86400 413 # else 414 zw3d(ji,jj,jk,1) = zno3phy * 86400 415 zw3d(ji,jj,jk,2) = znh4phy * 86400 416 zw3d(ji,jj,jk,3) = znh4no3 * 86400 417 # endif 418 #endif 410 419 END DO 411 420 END DO … … 413 422 414 423 #if defined key_trc_diaadd 415 ! Lateral boundary conditions on trc2d and trc3d 424 ! Lateral boundary conditions 425 # if ! defined key_iomput 416 426 DO jl = jp_lob0_2d, jp_lob1_2d 417 427 CALL lbc_lnk( trc2d(:,:,jl),'T', 1. ) 418 428 END DO 429 # else 430 DO jl = 1, 17 431 CALL lbc_lnk( zw2d(:,:,jl),'T', 1. ) 432 END DO 433 ! Save diagnostics 434 CALL iom_put( "TNO3PHY", zw2d(:,:,1) ) 435 CALL iom_put( "TNH4PHY", zw2d(:,:,2) ) 436 CALL iom_put( "TPHYDOM", zw2d(:,:,3) ) 437 CALL iom_put( "TPHYNH4", zw2d(:,:,4) ) 438 CALL iom_put( "TPHYZOO", zw2d(:,:,5) ) 439 CALL iom_put( "TPHYDET", zw2d(:,:,6) ) 440 CALL iom_put( "TDETZOO", zw2d(:,:,7) ) 441 CALL iom_put( "TZOODET", zw2d(:,:,8) ) 442 CALL iom_put( "TZOOBOD", zw2d(:,:,9) ) 443 CALL iom_put( "TZOONH4", zw2d(:,:,10) ) 444 CALL iom_put( "TZOODOM", zw2d(:,:,11) ) 445 CALL iom_put( "TNH4NO3", zw2d(:,:,12) ) 446 CALL iom_put( "TDOMNH4", zw2d(:,:,13) ) 447 CALL iom_put( "TDETNH4", zw2d(:,:,14) ) 448 CALL iom_put( "TPHYTOT", zw2d(:,:,15) ) 449 CALL iom_put( "TZOOTOT", zw2d(:,:,16) ) 450 CALL iom_put( "TDETDOM", zw2d(:,:,17) ) 451 # endif 452 #endif 453 454 #if defined key_trc_diaadd && defined key_trc_dia3d 455 ! Lateral boundary conditions 456 # if ! defined key_iomput 419 457 DO jl = jp_lob0_3d, jp_lob1_3d 420 458 CALL lbc_lnk( trc3d(:,:,1,jl),'T', 1. ) 421 459 END DO 460 # else 461 DO jl = 1, 3 462 CALL lbc_lnk( zw3d(:,:,:,jl),'T', 1. ) 463 END DO 464 ! save diagnostics 465 CALL iom_put( "FNO3PHY", zw3d(:,:,:,1) ) 466 CALL iom_put( "FNH4PHY", zw3d(:,:,:,2) ) 467 CALL iom_put( "FNH4NO3", zw3d(:,:,:,3) ) 468 # endif 422 469 #endif 423 470
Note: See TracChangeset
for help on using the changeset viewer.
