- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/TOP_SRC/PISCES/P2Z/p2zbio.F90
r5215 r6808 59 59 REAL(wp) :: fdbod ! zooplankton mortality fraction that goes to detritus 60 60 61 !! * Substitution62 # include " top_substitute.h90"61 !! * Substitutions 62 # include "vectopt_loop_substitute.h90" 63 63 !!---------------------------------------------------------------------- 64 64 !! NEMO/TOP 3.3 , NEMO Consortium (2010) … … 66 66 !! Software governed by the CeCILL licence (NEMOGCM/NEMO_CeCILL.txt) 67 67 !!---------------------------------------------------------------------- 68 69 68 CONTAINS 70 69 … … 186 185 ! closure : flux grazing is redistributed below level jpkbio 187 186 zzoobod = tmminz * zzoo * zzoo 188 xksi(ji,jj) = xksi(ji,jj) + (1-fdbod) * zzoobod * fse3t(ji,jj,jk)187 xksi(ji,jj) = xksi(ji,jj) + (1-fdbod) * zzoobod * e3t_n(ji,jj,jk) 189 188 zboddet = fdbod * zzoobod 190 189 … … 241 240 IF( ln_diatrc .OR. lk_iomput ) THEN 242 241 ! convert fluxes in per day 243 ze3t = fse3t(ji,jj,jk) * 86400.242 ze3t = e3t_n(ji,jj,jk) * 86400._wp 244 243 zw2d(ji,jj,1) = zw2d(ji,jj,1) + zno3phy * ze3t 245 244 zw2d(ji,jj,2) = zw2d(ji,jj,2) + znh4phy * ze3t … … 362 361 IF( ln_diatrc .OR. lk_iomput ) THEN 363 362 ! convert fluxes in per day 364 ze3t = fse3t(ji,jj,jk) * 86400.363 ze3t = e3t_n(ji,jj,jk) * 86400._wp 365 364 zw2d(ji,jj,1) = zw2d(ji,jj,1) + zno3phy * ze3t 366 365 zw2d(ji,jj,2) = zw2d(ji,jj,2) + znh4phy * ze3t … … 381 380 zw2d(ji,jj,17) = zw2d(ji,jj,17) + zdetdom * ze3t 382 381 ! 383 zw3d(ji,jj,jk,1) = zno3phy * 86400 384 zw3d(ji,jj,jk,2) = znh4phy * 86400 385 zw3d(ji,jj,jk,3) = znh4no3 * 86400 382 zw3d(ji,jj,jk,1) = zno3phy * 86400._wp 383 zw3d(ji,jj,jk,2) = znh4phy * 86400._wp 384 zw3d(ji,jj,jk,3) = znh4no3 * 86400._wp 386 385 ! 387 386 ENDIF … … 599 598 600 599 !!====================================================================== 601 END MODULE 600 END MODULE p2zbio
Note: See TracChangeset
for help on using the changeset viewer.