- Timestamp:
- 2019-11-22T15:29:17+01:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src
- Property svn:mergeinfo deleted
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/TOP/PISCES/P4Z/p5zlim.F90
r11536 r11949 99 99 CONTAINS 100 100 101 SUBROUTINE p5z_lim( kt, knt )101 SUBROUTINE p5z_lim( kt, knt, Kbb, Kmm ) 102 102 !!--------------------------------------------------------------------- 103 103 !! *** ROUTINE p5z_lim *** … … 110 110 ! 111 111 INTEGER, INTENT(in) :: kt, knt 112 INTEGER, INTENT(in) :: Kbb, Kmm ! time level indices 112 113 ! 113 114 INTEGER :: ji, jj, jk … … 134 135 ! Tuning of the iron concentration to a minimum level that is set to the detection limit 135 136 !------------------------------------- 136 zno3 = tr b(ji,jj,jk,jpno3) / 40.e-6137 zno3 = tr(ji,jj,jk,jpno3,Kbb) / 40.e-6 137 138 zferlim = MAX( 3e-11 * zno3 * zno3, 5e-12 ) 138 139 zferlim = MIN( zferlim, 7e-11 ) 139 tr b(ji,jj,jk,jpfer) = MAX( trb(ji,jj,jk,jpfer), zferlim )140 tr(ji,jj,jk,jpfer,Kbb) = MAX( tr(ji,jj,jk,jpfer,Kbb), zferlim ) 140 141 141 142 ! Computation of the mean relative size of each community 142 143 ! ------------------------------------------------------- 143 z1_trnphy = 1. / ( tr b(ji,jj,jk,jpphy) + rtrn )144 z1_trnpic = 1. / ( tr b(ji,jj,jk,jppic) + rtrn )145 z1_trndia = 1. / ( tr b(ji,jj,jk,jpdia) + rtrn )146 znanochl = tr b(ji,jj,jk,jpnch) * z1_trnphy147 zpicochl = tr b(ji,jj,jk,jppch) * z1_trnpic148 zdiatchl = tr b(ji,jj,jk,jpdch) * z1_trndia144 z1_trnphy = 1. / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn ) 145 z1_trnpic = 1. / ( tr(ji,jj,jk,jppic,Kbb) + rtrn ) 146 z1_trndia = 1. / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 147 znanochl = tr(ji,jj,jk,jpnch,Kbb) * z1_trnphy 148 zpicochl = tr(ji,jj,jk,jppch,Kbb) * z1_trnpic 149 zdiatchl = tr(ji,jj,jk,jpdch,Kbb) * z1_trndia 149 150 150 151 ! Computation of a variable Ks for iron on diatoms taking into account … … 182 183 ! Based on the different papers by Pahlow et al., and Smith et al. 183 184 ! ----------------------------------------------------------------- 184 znutlim = MAX( tr b(ji,jj,jk,jpnh4) / zconc0nnh4, &185 & tr b(ji,jj,jk,jpno3) / zconc0n)185 znutlim = MAX( tr(ji,jj,jk,jpnh4,Kbb) / zconc0nnh4, & 186 & tr(ji,jj,jk,jpno3,Kbb) / zconc0n) 186 187 fanano = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 187 znutlim = tr b(ji,jj,jk,jppo4) / zconc0npo4188 znutlim = tr(ji,jj,jk,jppo4,Kbb) / zconc0npo4 188 189 fananop = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 189 190 znutlim = biron(ji,jj,jk) / zconcnfe 190 191 fananof = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 191 znutlim = MAX( tr b(ji,jj,jk,jpnh4) / zconc0pnh4, &192 & tr b(ji,jj,jk,jpno3) / zconc0p)192 znutlim = MAX( tr(ji,jj,jk,jpnh4,Kbb) / zconc0pnh4, & 193 & tr(ji,jj,jk,jpno3,Kbb) / zconc0p) 193 194 fapico = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 194 znutlim = tr b(ji,jj,jk,jppo4) / zconc0ppo4195 znutlim = tr(ji,jj,jk,jppo4,Kbb) / zconc0ppo4 195 196 fapicop = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 196 197 znutlim = biron(ji,jj,jk) / zconcpfe 197 198 fapicof = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 198 znutlim = MAX( tr b(ji,jj,jk,jpnh4) / zconc1dnh4, &199 & tr b(ji,jj,jk,jpno3) / zconc1d )199 znutlim = MAX( tr(ji,jj,jk,jpnh4,Kbb) / zconc1dnh4, & 200 & tr(ji,jj,jk,jpno3,Kbb) / zconc1d ) 200 201 fadiat = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 201 znutlim = tr b(ji,jj,jk,jppo4) / zconc0dpo4202 znutlim = tr(ji,jj,jk,jppo4,Kbb) / zconc0dpo4 202 203 fadiatp = MAX(0.01, MIN(0.99, 1. / ( SQRT(znutlim) + 1.) ) ) 203 204 znutlim = biron(ji,jj,jk) / zconcdfe … … 206 207 ! Michaelis-Menten Limitation term for nutrients Small bacteria 207 208 ! ------------------------------------------------------------- 208 zbactnh4 = tr b(ji,jj,jk,jpnh4) / ( concbnh4 + trb(ji,jj,jk,jpnh4) )209 zbactno3 = tr b(ji,jj,jk,jpno3) / ( concbno3 + trb(ji,jj,jk,jpno3) ) * (1. - zbactnh4)209 zbactnh4 = tr(ji,jj,jk,jpnh4,Kbb) / ( concbnh4 + tr(ji,jj,jk,jpnh4,Kbb) ) 210 zbactno3 = tr(ji,jj,jk,jpno3,Kbb) / ( concbno3 + tr(ji,jj,jk,jpno3,Kbb) ) * (1. - zbactnh4) 210 211 ! 211 212 zlim1 = zbactno3 + zbactnh4 212 zlim2 = tr b(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + concbpo4)213 zlim2 = tr(ji,jj,jk,jppo4,Kbb) / ( tr(ji,jj,jk,jppo4,Kbb) + concbpo4) 213 214 zlim3 = biron(ji,jj,jk) / ( concbfe + biron(ji,jj,jk) ) 214 zlim4 = tr b(ji,jj,jk,jpdoc) / ( xkdoc + trb(ji,jj,jk,jpdoc) )215 zlim4 = tr(ji,jj,jk,jpdoc,Kbb) / ( xkdoc + tr(ji,jj,jk,jpdoc,Kbb) ) 215 216 xlimbacl(ji,jj,jk) = MIN( zlim1, zlim2, zlim3 ) 216 217 xlimbac (ji,jj,jk) = xlimbacl(ji,jj,jk) * zlim4 … … 219 220 ! ----------------------------------------------- 220 221 zfalim = (1.-fanano) / fanano 221 xnanonh4(ji,jj,jk) = (1. - fanano) * tr b(ji,jj,jk,jpnh4) / ( zfalim * zconc0nnh4 + trb(ji,jj,jk,jpnh4) )222 xnanono3(ji,jj,jk) = (1. - fanano) * tr b(ji,jj,jk,jpno3) / ( zfalim * zconc0n + trb(ji,jj,jk,jpno3) ) &222 xnanonh4(ji,jj,jk) = (1. - fanano) * tr(ji,jj,jk,jpnh4,Kbb) / ( zfalim * zconc0nnh4 + tr(ji,jj,jk,jpnh4,Kbb) ) 223 xnanono3(ji,jj,jk) = (1. - fanano) * tr(ji,jj,jk,jpno3,Kbb) / ( zfalim * zconc0n + tr(ji,jj,jk,jpno3,Kbb) ) & 223 224 & * (1. - xnanonh4(ji,jj,jk)) 224 225 ! 225 226 zfalim = (1.-fananop) / fananop 226 xnanopo4(ji,jj,jk) = (1. - fananop) * tr b(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + zfalim * zconc0npo4 )227 xnanodop(ji,jj,jk) = tr b(ji,jj,jk,jpdop) / ( trb(ji,jj,jk,jpdop) + xkdoc ) &227 xnanopo4(ji,jj,jk) = (1. - fananop) * tr(ji,jj,jk,jppo4,Kbb) / ( tr(ji,jj,jk,jppo4,Kbb) + zfalim * zconc0npo4 ) 228 xnanodop(ji,jj,jk) = tr(ji,jj,jk,jpdop,Kbb) / ( tr(ji,jj,jk,jpdop,Kbb) + xkdoc ) & 228 229 & * ( 1.0 - xnanopo4(ji,jj,jk) ) 229 230 xnanodop(ji,jj,jk) = 0. … … 232 233 xnanofer(ji,jj,jk) = (1. - fananof) * biron(ji,jj,jk) / ( biron(ji,jj,jk) + zfalim * zconcnfe ) 233 234 ! 234 zratiof = tr b(ji,jj,jk,jpnfe) * z1_trnphy235 zratiof = tr(ji,jj,jk,jpnfe,Kbb) * z1_trnphy 235 236 zqfemn = xcoef1 * znanochl + xcoef2 + xcoef3 * xnanono3(ji,jj,jk) 236 237 ! 237 zration = tr b(ji,jj,jk,jpnph) * z1_trnphy238 zration = tr(ji,jj,jk,jpnph,Kbb) * z1_trnphy 238 239 zration = MIN(xqnnmax(ji,jj,jk), MAX( 2. * xqnnmin(ji,jj,jk), zration )) 239 240 fvnuptk(ji,jj,jk) = 1. / zpsiuptk * rno3 * 2. * xqnnmin(ji,jj,jk) / (zration + rtrn) & … … 250 251 ! ---------------------------------------------------------------- 251 252 zfalim = (1.-fapico) / fapico 252 xpiconh4(ji,jj,jk) = (1. - fapico) * tr b(ji,jj,jk,jpnh4) / ( zfalim * zconc0pnh4 + trb(ji,jj,jk,jpnh4) )253 xpicono3(ji,jj,jk) = (1. - fapico) * tr b(ji,jj,jk,jpno3) / ( zfalim * zconc0p + trb(ji,jj,jk,jpno3) ) &253 xpiconh4(ji,jj,jk) = (1. - fapico) * tr(ji,jj,jk,jpnh4,Kbb) / ( zfalim * zconc0pnh4 + tr(ji,jj,jk,jpnh4,Kbb) ) 254 xpicono3(ji,jj,jk) = (1. - fapico) * tr(ji,jj,jk,jpno3,Kbb) / ( zfalim * zconc0p + tr(ji,jj,jk,jpno3,Kbb) ) & 254 255 & * (1. - xpiconh4(ji,jj,jk)) 255 256 ! 256 257 zfalim = (1.-fapicop) / fapicop 257 xpicopo4(ji,jj,jk) = (1. - fapicop) * tr b(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + zfalim * zconc0ppo4 )258 xpicodop(ji,jj,jk) = tr b(ji,jj,jk,jpdop) / ( trb(ji,jj,jk,jpdop) + xkdoc ) &258 xpicopo4(ji,jj,jk) = (1. - fapicop) * tr(ji,jj,jk,jppo4,Kbb) / ( tr(ji,jj,jk,jppo4,Kbb) + zfalim * zconc0ppo4 ) 259 xpicodop(ji,jj,jk) = tr(ji,jj,jk,jpdop,Kbb) / ( tr(ji,jj,jk,jpdop,Kbb) + xkdoc ) & 259 260 & * ( 1.0 - xpicopo4(ji,jj,jk) ) 260 261 xpicodop(ji,jj,jk) = 0. … … 263 264 xpicofer(ji,jj,jk) = (1. - fapicof) * biron(ji,jj,jk) / ( biron(ji,jj,jk) + zfalim * zconcpfe ) 264 265 ! 265 zratiof = tr b(ji,jj,jk,jppfe) * z1_trnpic266 zratiof = tr(ji,jj,jk,jppfe,Kbb) * z1_trnpic 266 267 zqfemp = xcoef1 * zpicochl + xcoef2 + xcoef3 * xpicono3(ji,jj,jk) 267 268 ! 268 zration = tr b(ji,jj,jk,jpnpi) * z1_trnpic269 zration = tr(ji,jj,jk,jpnpi,Kbb) * z1_trnpic 269 270 zration = MIN(xqnpmax(ji,jj,jk), MAX( 2. * xqnpmin(ji,jj,jk), zration )) 270 271 fvpuptk(ji,jj,jk) = 1. / zpsiuptk * rno3 * 2. * xqnpmin(ji,jj,jk) / (zration + rtrn) & … … 281 282 ! ------------------------------------------------------ 282 283 zfalim = (1.-fadiat) / fadiat 283 xdiatnh4(ji,jj,jk) = (1. - fadiat) * tr b(ji,jj,jk,jpnh4) / ( zfalim * zconc1dnh4 + trb(ji,jj,jk,jpnh4) )284 xdiatno3(ji,jj,jk) = (1. - fadiat) * tr b(ji,jj,jk,jpno3) / ( zfalim * zconc1d + trb(ji,jj,jk,jpno3) ) &284 xdiatnh4(ji,jj,jk) = (1. - fadiat) * tr(ji,jj,jk,jpnh4,Kbb) / ( zfalim * zconc1dnh4 + tr(ji,jj,jk,jpnh4,Kbb) ) 285 xdiatno3(ji,jj,jk) = (1. - fadiat) * tr(ji,jj,jk,jpno3,Kbb) / ( zfalim * zconc1d + tr(ji,jj,jk,jpno3,Kbb) ) & 285 286 & * (1. - xdiatnh4(ji,jj,jk)) 286 287 ! 287 288 zfalim = (1.-fadiatp) / fadiatp 288 xdiatpo4(ji,jj,jk) = (1. - fadiatp) * tr b(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + zfalim * zconc0dpo4 )289 xdiatdop(ji,jj,jk) = tr b(ji,jj,jk,jpdop) / ( trb(ji,jj,jk,jpdop) + xkdoc ) &289 xdiatpo4(ji,jj,jk) = (1. - fadiatp) * tr(ji,jj,jk,jppo4,Kbb) / ( tr(ji,jj,jk,jppo4,Kbb) + zfalim * zconc0dpo4 ) 290 xdiatdop(ji,jj,jk) = tr(ji,jj,jk,jpdop,Kbb) / ( tr(ji,jj,jk,jpdop,Kbb) + xkdoc ) & 290 291 & * ( 1.0 - xdiatpo4(ji,jj,jk) ) 291 292 xdiatdop(ji,jj,jk) = 0. … … 294 295 xdiatfer(ji,jj,jk) = (1. - fadiatf) * biron(ji,jj,jk) / ( biron(ji,jj,jk) + zfalim * zconcdfe ) 295 296 ! 296 zratiof = tr b(ji,jj,jk,jpdfe) * z1_trndia297 zratiof = tr(ji,jj,jk,jpdfe,Kbb) * z1_trndia 297 298 zqfemd = xcoef1 * zdiatchl + xcoef2 + xcoef3 * xdiatno3(ji,jj,jk) 298 299 ! 299 zration = tr b(ji,jj,jk,jpndi) * z1_trndia300 zration = tr(ji,jj,jk,jpndi,Kbb) * z1_trndia 300 301 zration = MIN(xqndmax(ji,jj,jk), MAX( 2. * xqndmin(ji,jj,jk), zration )) 301 302 fvduptk(ji,jj,jk) = 1. / zpsiuptk * rno3 * 2. * xqndmin(ji,jj,jk) / (zration + rtrn) & … … 305 306 & / (xqndmax(ji,jj,jk) - 2. * xqndmin(ji,jj,jk) ) ) & 306 307 & * xqndmax(ji,jj,jk) / (zration + rtrn) 307 zlim3 = tr b(ji,jj,jk,jpsil) / ( trb(ji,jj,jk,jpsil) + xksi(ji,jj) )308 zlim3 = tr(ji,jj,jk,jpsil,Kbb) / ( tr(ji,jj,jk,jpsil,Kbb) + xksi(ji,jj) ) 308 309 zlim4 = MAX( 0., ( zratiof - zqfemd ) / qfdopt ) 309 310 xlimdfe(ji,jj,jk) = MIN( 1., zlim4 ) … … 330 331 ! ------------------------------ 331 332 zfuptk = 0.23 * zfvn 332 zrpho = 2.24 * tr b(ji,jj,jk,jpnch) / ( trb(ji,jj,jk,jpnph) * rno3 * 15. + rtrn )333 zrpho = 2.24 * tr(ji,jj,jk,jpnch,Kbb) / ( tr(ji,jj,jk,jpnph,Kbb) * rno3 * 15. + rtrn ) 333 334 zrass = 1. - 0.2 - zrpho - zfuptk 334 335 xqpnmax(ji,jj,jk) = ( zfuptk + zrpho ) * 0.0128 * 16. + zrass * 1./ 7.2 * 16. 335 xqpnmax(ji,jj,jk) = xqpnmax(ji,jj,jk) * tr b(ji,jj,jk,jpnph) / ( trb(ji,jj,jk,jpphy) + rtrn ) + 0.13336 xqpnmax(ji,jj,jk) = xqpnmax(ji,jj,jk) * tr(ji,jj,jk,jpnph,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn ) + 0.13 336 337 xqpnmin(ji,jj,jk) = 0.13 + 0.23 * 0.0128 * 16. 337 338 … … 344 345 ! ------------------------------ 345 346 zfuptk = 0.35 * zfvn 346 zrpho = 2.24 * tr b(ji,jj,jk,jppch) / ( trb(ji,jj,jk,jpnpi) * rno3 * 15. + rtrn )347 zrpho = 2.24 * tr(ji,jj,jk,jppch,Kbb) / ( tr(ji,jj,jk,jpnpi,Kbb) * rno3 * 15. + rtrn ) 347 348 zrass = 1. - 0.4 - zrpho - zfuptk 348 349 xqppmax(ji,jj,jk) = (zrpho + zfuptk) * 0.0128 * 16. + zrass * 1./ 9. * 16. 349 xqppmax(ji,jj,jk) = xqppmax(ji,jj,jk) * tr b(ji,jj,jk,jpnpi) / ( trb(ji,jj,jk,jppic) + rtrn ) + 0.13350 xqppmax(ji,jj,jk) = xqppmax(ji,jj,jk) * tr(ji,jj,jk,jpnpi,Kbb) / ( tr(ji,jj,jk,jppic,Kbb) + rtrn ) + 0.13 350 351 xqppmin(ji,jj,jk) = 0.13 351 352 … … 354 355 zfvn = 2. * fvduptk(ji,jj,jk) 355 356 sized(ji,jj,jk) = MAX(1., MIN(xsizerd, 1.0 / ( MAX(rtrn, zfvn) ) ) ) 356 zcoef = tr b(ji,jj,jk,jpdia) - MIN(xsizedia, trb(ji,jj,jk,jpdia) )357 zcoef = tr(ji,jj,jk,jpdia,Kbb) - MIN(xsizedia, tr(ji,jj,jk,jpdia,Kbb) ) 357 358 sized(ji,jj,jk) = 1. + xsizerd * zcoef *1E6 / ( 1. + zcoef * 1E6 ) 358 359 … … 360 361 ! -------------------- 361 362 zfuptk = 0.2 * zfvn 362 zrpho = 2.24 * tr b(ji,jj,jk,jpdch) / ( trb(ji,jj,jk,jpndi) * rno3 * 15. + rtrn )363 zrpho = 2.24 * tr(ji,jj,jk,jpdch,Kbb) / ( tr(ji,jj,jk,jpndi,Kbb) * rno3 * 15. + rtrn ) 363 364 zrass = 1. - 0.2 - zrpho - zfuptk 364 365 xqpdmax(ji,jj,jk) = ( zfuptk + zrpho ) * 0.0128 * 16. + zrass * 1./ 7.2 * 16. 365 xqpdmax(ji,jj,jk) = xqpdmax(ji,jj,jk) * tr b(ji,jj,jk,jpndi) / ( trb(ji,jj,jk,jpdia) + rtrn ) + 0.13366 xqpdmax(ji,jj,jk) = xqpdmax(ji,jj,jk) * tr(ji,jj,jk,jpndi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) + 0.13 366 367 xqpdmin(ji,jj,jk) = 0.13 + 0.2 * 0.0128 * 16. 367 368 … … 375 376 DO jj = 1, jpj 376 377 DO ji = 1, jpi 377 zlim1 = tr b(ji,jj,jk,jpnh4) / ( trb(ji,jj,jk,jpnh4) + concnnh4 ) + trb(ji,jj,jk,jpno3) &378 & / ( tr b(ji,jj,jk,jpno3) + concnno3 ) * ( 1.0 - trb(ji,jj,jk,jpnh4) &379 & / ( tr b(ji,jj,jk,jpnh4) + concnnh4 ) )380 zlim2 = tr b(ji,jj,jk,jppo4) / ( trb(ji,jj,jk,jppo4) + concnpo4 )381 zlim3 = tr b(ji,jj,jk,jpfer) / ( trb(ji,jj,jk,jpfer) + 5.E-11 )382 ztem1 = MAX( 0., ts n(ji,jj,jk,jp_tem) )383 ztem2 = ts n(ji,jj,jk,jp_tem) - 10.378 zlim1 = tr(ji,jj,jk,jpnh4,Kbb) / ( tr(ji,jj,jk,jpnh4,Kbb) + concnnh4 ) + tr(ji,jj,jk,jpno3,Kbb) & 379 & / ( tr(ji,jj,jk,jpno3,Kbb) + concnno3 ) * ( 1.0 - tr(ji,jj,jk,jpnh4,Kbb) & 380 & / ( tr(ji,jj,jk,jpnh4,Kbb) + concnnh4 ) ) 381 zlim2 = tr(ji,jj,jk,jppo4,Kbb) / ( tr(ji,jj,jk,jppo4,Kbb) + concnpo4 ) 382 zlim3 = tr(ji,jj,jk,jpfer,Kbb) / ( tr(ji,jj,jk,jpfer,Kbb) + 5.E-11 ) 383 ztem1 = MAX( 0., ts(ji,jj,jk,jp_tem,Kmm) ) 384 ztem2 = ts(ji,jj,jk,jp_tem,Kmm) - 10. 384 385 zetot1 = MAX( 0., etot(ji,jj,jk) - 1.) / ( 4. + etot(ji,jj,jk) ) * 20. / ( 20. + etot(ji,jj,jk) ) 385 386 386 387 ! xfracal(ji,jj,jk) = caco3r * MIN( zlim1, zlim2, zlim3 ) & 387 388 xfracal(ji,jj,jk) = caco3r & 388 & * ztem1 / ( 1. + ztem1 ) * MAX( 1., tr b(ji,jj,jk,jpphy)*1E6 ) &389 & * ztem1 / ( 1. + ztem1 ) * MAX( 1., tr(ji,jj,jk,jpphy,Kbb)*1E6 ) & 389 390 & * ( 1. + EXP(-ztem2 * ztem2 / 25. ) ) & 390 391 & * zetot1 * MIN( 1., 50. / ( hmld(ji,jj) + rtrn ) ) … … 398 399 DO ji = 1, jpi 399 400 ! denitrification factor computed from O2 levels 400 nitrfac(ji,jj,jk) = MAX( 0.e0, 0.4 * ( 6.e-6 - tr b(ji,jj,jk,jpoxy) ) &401 & / ( oxymin + tr b(ji,jj,jk,jpoxy) ) )401 nitrfac(ji,jj,jk) = MAX( 0.e0, 0.4 * ( 6.e-6 - tr(ji,jj,jk,jpoxy,Kbb) ) & 402 & / ( oxymin + tr(ji,jj,jk,jpoxy,Kbb) ) ) 402 403 nitrfac(ji,jj,jk) = MIN( 1., nitrfac(ji,jj,jk) ) 403 404 END DO
Note: See TracChangeset
for help on using the changeset viewer.