- Timestamp:
- 2009-11-20T13:21:07+01:00 (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/TOP_SRC/PISCES/p4zrem.F90
r1736 r1744 41 41 & denitr !: denitrification array 42 42 43 REAL(wp) :: & 44 xstep !: Time step duration for biology 43 45 44 46 !!* Substitution … … 70 72 REAL(wp) :: zofer2, zdenom, zdenom2 71 73 #endif 72 REAL(wp) :: zlamfac, z step, zonitr74 REAL(wp) :: zlamfac, zonitr 73 75 REAL(wp), DIMENSION(jpi,jpj) :: ztempbac 74 76 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zdepbac, zfesatur, zolimi … … 78 80 79 81 80 IF( ( kt * jnt ) == nittrc000 ) CALL p4z_rem_init ! Initialization (first time-step only) 81 82 zstep = rfact2 / rday ! Time step duration for the biology 83 82 IF( ( kt * jnt ) == nittrc000 ) THEN 83 CALL p4z_rem_init ! Initialization (first time-step only) 84 xstep = rfact2 / rday ! Time step duration for the biology 85 nitrfac(:,:,:) = 0.0 86 denitr (:,:,:) = 0.0 87 ENDIF 88 89 90 ! Initialisation of temprary arrys 91 zdepbac (:,:,:) = 0.0 92 zfesatur(:,:,:) = 0.0 93 zolimi (:,:,:) = 0.0 94 ztempbac(:,:) = 0.0 84 95 85 96 ! Computation of the mean phytoplankton concentration as … … 124 135 ! of the bacterial activity. 125 136 ! ---------------------------------------------------------------- 126 zremik = xremik * zstep / 1.e-6 * xlimbac(ji,jj,jk) &137 zremik = xremik * xstep / 1.e-6 * xlimbac(ji,jj,jk) & 127 138 # if defined key_off_degrad 128 139 & * facvol(ji,jj,jk) & 129 140 # endif 130 141 & * zdepbac(ji,jj,jk) 131 zremik = MAX( zremik, 5.5e-4 * zstep )142 zremik = MAX( zremik, 5.5e-4 * xstep ) 132 143 133 144 ! Ammonification in oxic waters with oxygen consumption … … 144 155 END DO 145 156 146 zolimi (:,:,:) = MAX( 0.e0, zolimi (:,:,:) ) 147 denitr (:,:,:) = MAX( 0.e0, denitr (:,:,:) ) 157 DO jk = 1, jpkm1 158 DO jj = 1, jpj 159 DO ji = 1, jpi 160 zolimi (ji,jj,jk) = MAX( 0.e0, zolimi (ji,jj,jk) ) 161 denitr (ji,jj,jk) = MAX( 0.e0, denitr (ji,jj,jk) ) 162 END DO 163 END DO 164 END DO 148 165 149 166 DO jk = 1, jpkm1 … … 154 171 ! below 2 umol/L. Inhibited at strong light 155 172 ! ---------------------------------------------------------- 156 zonitr = nitrif * zstep * trn(ji,jj,jk,jpnh4) / ( 1.+ emoy(ji,jj,jk) ) &173 zonitr = nitrif * xstep * trn(ji,jj,jk,jpnh4) / ( 1.+ emoy(ji,jj,jk) ) & 157 174 # if defined key_off_degrad 158 175 & * facvol(ji,jj,jk) & … … 216 233 ! POC disaggregation by turbulence and bacterial activity. 217 234 ! ------------------------------------------------------------- 218 zremip = xremip * zstep * tgfunc(ji,jj,jk) &235 zremip = xremip * xstep * tgfunc(ji,jj,jk) & 219 236 # if defined key_off_degrad 220 237 & * facvol(ji,jj,jk) & … … 273 290 znusil = 0.225 * ( 1. + tn(ji,jj,jk) / 15.) * zsatur + 0.775 * zsatur2**9 274 291 # if defined key_off_degrad 275 zsiremin = xsirem * zstep * znusil * facvol(ji,jj,jk)292 zsiremin = xsirem * xstep * znusil * facvol(ji,jj,jk) 276 293 # else 277 zsiremin = xsirem * zstep * znusil294 zsiremin = xsirem * xstep * znusil 278 295 # endif 279 296 zosil = zsiremin * trn(ji,jj,jk,jpdsi) … … 334 351 335 352 # if defined key_off_degrad 336 zscave = zfeequi * zlam1b * zstep * facvol(ji,jj,jk)353 zscave = zfeequi * zlam1b * xstep * facvol(ji,jj,jk) 337 354 # else 338 zscave = zfeequi * zlam1b * zstep355 zscave = zfeequi * zlam1b * xstep 339 356 # endif 340 357 … … 358 375 359 376 # if defined key_off_degrad 360 zaggdfe = zlam1b * zstep * 0.5 * ( trn(ji,jj,jk,jpfer) - zfeequi ) * facvol(ji,jj,jk)377 zaggdfe = zlam1b * xstep * 0.5 * ( trn(ji,jj,jk,jpfer) - zfeequi ) * facvol(ji,jj,jk) 361 378 # else 362 zaggdfe = zlam1b * zstep * 0.5 * ( trn(ji,jj,jk,jpfer) - zfeequi )379 zaggdfe = zlam1b * xstep * 0.5 * ( trn(ji,jj,jk,jpfer) - zfeequi ) 363 380 # endif 364 381
Note: See TracChangeset
for help on using the changeset viewer.