- 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/PISCES/p4zsink.F90
r1970 r2038 97 97 REAL(wp) :: zdiv , zdiv1, zdiv2, zdiv3, zdiv4, zdiv5 98 98 REAL(wp) :: zval1, zval2, zval3, zval4 99 #if defined key_ trc_diaadd99 #if defined key_diatrc 100 100 REAL(wp) :: zrfact2 101 101 INTEGER :: ik1 … … 207 207 & * (zfm*xkr_mass_max**2-xkr_mass_min**2) & 208 208 & * (zeps-1.)**2/(zdiv2*zdiv3)) & 209 # if defined key_ off_degrad209 # if defined key_degrad 210 210 & *facvol(ji,jj,jk) & 211 211 # endif … … 219 219 & -zfm*xkr_mass_max**3*(1.+3.*((zeps-1.)/ & 220 220 & (zeps-2.)+(zeps-1.)/zdiv3)+(zeps-1.)/zdiv1)) & 221 # if defined key_ off_degrad221 # if defined key_degrad 222 222 & *facvol(ji,jj,jk) & 223 223 # endif … … 225 225 226 226 zagg3 = ( 0.163*trn(ji,jj,jk,jpnum)**2*zfm**2*8. * xkr_mass_max**3 & 227 # if defined key_ off_degrad227 # if defined key_degrad 228 228 & *facvol(ji,jj,jk) & 229 229 # endif … … 242 242 & ((zfm*zfm*xkr_mass_max**2*zsm-xkr_mass_min**2) & 243 243 & *xkr_eta)/(zdiv*zdiv3*zdiv5) ) & 244 # if defined key_ off_degrad244 # if defined key_degrad 245 245 & *facvol(ji,jj,jk) & 246 246 # endif … … 252 252 & /zdiv3-(xkr_mass_min**2-zfm*zsm*xkr_mass_max**2) & 253 253 & /zdiv) & 254 # if defined key_ off_degrad254 # if defined key_degrad 255 255 & *facvol(ji,jj,jk) & 256 256 # endif … … 266 266 zaggdoc = ( 0.4 * trn(ji,jj,jk,jpdoc) & 267 267 & + 1018. * trn(ji,jj,jk,jppoc) ) * xstep & 268 # if defined key_ off_degrad268 # if defined key_degrad 269 269 & * facvol(ji,jj,jk) & 270 270 # endif … … 281 281 END DO 282 282 283 #if defined key_ trc_diaadd283 #if defined key_diatrc 284 284 zrfact2 = 1.e3 * rfact2r 285 285 ik1 = iksed + 1 … … 474 474 REAL(wp) :: zagg , zaggfe, zaggdoc, zaggdoc2 475 475 REAL(wp) :: zfact, zwsmax 476 #if defined key_ trc_dia3d476 #if defined key_diatrc 477 477 REAL(wp) :: zrfact2 478 478 INTEGER :: ik1 … … 552 552 zfact = xstep * xdiss(ji,jj,jk) 553 553 ! Part I : Coagulation dependent on turbulence 554 # if defined key_ off_degrad554 # if defined key_degrad 555 555 zagg1 = 940.* zfact * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jppoc) * facvol(ji,jj,jk) 556 556 zagg2 = 1.054e4 * zfact * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpgoc) * facvol(ji,jj,jk) … … 563 563 564 564 ! Aggregation of small into large particles 565 # if defined key_ off_degrad565 # if defined key_degrad 566 566 zagg3 = 0.66 * xstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jpgoc) * facvol(ji,jj,jk) 567 567 zagg4 = 0.e0 * xstep * trn(ji,jj,jk,jppoc) * trn(ji,jj,jk,jppoc) * facvol(ji,jj,jk) … … 574 574 575 575 ! Aggregation of DOC to small particles 576 #if defined key_ off_degrad576 #if defined key_degrad 577 577 zaggdoc = ( 80.* trn(ji,jj,jk,jpdoc) + 698. * trn(ji,jj,jk,jppoc) ) & 578 578 & * facvol(ji,jj,jk) * zfact * trn(ji,jj,jk,jpdoc) … … 595 595 END DO 596 596 597 #if defined key_ trc_diaadd597 #if defined key_diatrc 598 598 zrfact2 = 1.e3 * rfact2r 599 599 ik1 = iksed + 1 … … 651 651 652 652 DO jk = 1, jpkm1 653 # if defined key_ off_degrad653 # if defined key_degrad 654 654 zwsink2(:,:,jk+1) = -pwsink(:,:,jk) / rday * tmask(:,:,jk+1) * facvol(:,:,jk) 655 655 # else
Note: See TracChangeset
for help on using the changeset viewer.