- Timestamp:
- 2019-11-22T15:29:17+01:00 (5 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/p5zprod.F90
r11536 r11949 57 57 CONTAINS 58 58 59 SUBROUTINE p5z_prod( kt , knt )59 SUBROUTINE p5z_prod( kt , knt, Kbb, Kmm, Krhs ) 60 60 !!--------------------------------------------------------------------- 61 61 !! *** ROUTINE p5z_prod *** … … 68 68 ! 69 69 INTEGER, INTENT(in) :: kt, knt 70 INTEGER, INTENT(in) :: Kbb, Kmm, Krhs ! time level indices 70 71 ! 71 72 INTEGER :: ji, jj, jk … … 136 137 IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 137 138 zval = MAX( 1., zstrn(ji,jj) ) 138 IF( gdepw _n(ji,jj,jk+1) <= hmld(ji,jj) ) THEN139 IF( gdepw(ji,jj,jk+1,Kmm) <= hmld(ji,jj) ) THEN 139 140 zval = zval * MIN(1., heup_01(ji,jj) / ( hmld(ji,jj) + rtrn )) 140 141 ENDIF … … 160 161 IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 161 162 ! Computation of the P-I slope for nanos and diatoms 162 ztn = MAX( 0., ts n(ji,jj,jk,jp_tem) - 15. )163 ztn = MAX( 0., ts(ji,jj,jk,jp_tem,Kmm) - 15. ) 163 164 zadap = xadap * ztn / ( 2.+ ztn ) 164 165 ! 165 zpislopeadn(ji,jj,jk) = pislopen * tr b(ji,jj,jk,jpnch) &166 & /( tr b(ji,jj,jk,jpphy) * 12. + rtrn)166 zpislopeadn(ji,jj,jk) = pislopen * tr(ji,jj,jk,jpnch,Kbb) & 167 & /( tr(ji,jj,jk,jpphy,Kbb) * 12. + rtrn) 167 168 zpislopeadp(ji,jj,jk) = pislopep * ( 1. + zadap * EXP( -0.25 * epico(ji,jj,jk) ) ) & 168 & * tr b(ji,jj,jk,jppch) /( trb(ji,jj,jk,jppic) * 12. + rtrn)169 zpislopeadd(ji,jj,jk) = pisloped * tr b(ji,jj,jk,jpdch) &170 & /( tr b(ji,jj,jk,jpdia) * 12. + rtrn)169 & * tr(ji,jj,jk,jppch,Kbb) /( tr(ji,jj,jk,jppic,Kbb) * 12. + rtrn) 170 zpislopeadd(ji,jj,jk) = pisloped * tr(ji,jj,jk,jpdch,Kbb) & 171 & /( tr(ji,jj,jk,jpdia,Kbb) * 12. + rtrn) 171 172 ! 172 173 zpislopen = zpislopeadn(ji,jj,jk) / ( zprbio(ji,jj,jk) * rday * xlimphy(ji,jj,jk) + rtrn ) … … 203 204 ! Si/C is arbitrariliy increased for very high Si concentrations 204 205 ! to mimic the very high ratios observed in the Southern Ocean (silpot2) 205 zlim = tr b(ji,jj,jk,jpsil) / ( trb(ji,jj,jk,jpsil) + xksi1 )206 zlim = tr(ji,jj,jk,jpsil,Kbb) / ( tr(ji,jj,jk,jpsil,Kbb) + xksi1 ) 206 207 zsilim = MIN( zprdia(ji,jj,jk) / ( zprmaxd(ji,jj,jk) + rtrn ), xlimsi(ji,jj,jk) ) 207 208 zsilfac = 3.4 * EXP( -4.23 * zsilim ) * MAX( 0.e0, MIN( 1., 2.2 * ( zlim - 0.5 ) ) ) + 1.e0 208 zsiborn = tr b(ji,jj,jk,jpsil) * trb(ji,jj,jk,jpsil) * trb(ji,jj,jk,jpsil)209 zsiborn = tr(ji,jj,jk,jpsil,Kbb) * tr(ji,jj,jk,jpsil,Kbb) * tr(ji,jj,jk,jpsil,Kbb) 209 210 IF (gphit(ji,jj) < -30 ) THEN 210 211 zsilfac2 = 1. + 2. * zsiborn / ( zsiborn + xksi2**3 ) … … 236 237 IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 237 238 ! production terms for nanophyto. 238 zprorcan(ji,jj,jk) = zprbio(ji,jj,jk) * xlimphy(ji,jj,jk) * tr b(ji,jj,jk,jpphy) * rfact2239 zprorcan(ji,jj,jk) = zprbio(ji,jj,jk) * xlimphy(ji,jj,jk) * tr(ji,jj,jk,jpphy,Kbb) * rfact2 239 240 ! 240 zration = tr b(ji,jj,jk,jpnph) / ( trb(ji,jj,jk,jpphy) + rtrn )241 zratiop = tr b(ji,jj,jk,jppph) / ( trb(ji,jj,jk,jpphy) + rtrn )242 zratiof = tr b(ji,jj,jk,jpnfe) / ( trb(ji,jj,jk,jpphy) + rtrn )243 zprnutmax = zprnut(ji,jj,jk) * fvnuptk(ji,jj,jk) / rno3 * tr b(ji,jj,jk,jpphy) * rfact2241 zration = tr(ji,jj,jk,jpnph,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn ) 242 zratiop = tr(ji,jj,jk,jppph,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn ) 243 zratiof = tr(ji,jj,jk,jpnfe,Kbb) / ( tr(ji,jj,jk,jpphy,Kbb) + rtrn ) 244 zprnutmax = zprnut(ji,jj,jk) * fvnuptk(ji,jj,jk) / rno3 * tr(ji,jj,jk,jpphy,Kbb) * rfact2 244 245 ! Uptake of nitrogen 245 246 zrat = MIN( 1., zration / (xqnnmax(ji,jj,jk) + rtrn) ) … … 273 274 IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 274 275 ! production terms for picophyto. 275 zprorcap(ji,jj,jk) = zprpic(ji,jj,jk) * xlimpic(ji,jj,jk) * tr b(ji,jj,jk,jppic) * rfact2276 zprorcap(ji,jj,jk) = zprpic(ji,jj,jk) * xlimpic(ji,jj,jk) * tr(ji,jj,jk,jppic,Kbb) * rfact2 276 277 ! 277 zration = tr b(ji,jj,jk,jpnpi) / ( trb(ji,jj,jk,jppic) + rtrn )278 zratiop = tr b(ji,jj,jk,jpppi) / ( trb(ji,jj,jk,jppic) + rtrn )279 zratiof = tr b(ji,jj,jk,jppfe) / ( trb(ji,jj,jk,jppic) + rtrn )280 zprnutmax = zprnut(ji,jj,jk) * fvpuptk(ji,jj,jk) / rno3 * tr b(ji,jj,jk,jppic) * rfact2278 zration = tr(ji,jj,jk,jpnpi,Kbb) / ( tr(ji,jj,jk,jppic,Kbb) + rtrn ) 279 zratiop = tr(ji,jj,jk,jpppi,Kbb) / ( tr(ji,jj,jk,jppic,Kbb) + rtrn ) 280 zratiof = tr(ji,jj,jk,jppfe,Kbb) / ( tr(ji,jj,jk,jppic,Kbb) + rtrn ) 281 zprnutmax = zprnut(ji,jj,jk) * fvpuptk(ji,jj,jk) / rno3 * tr(ji,jj,jk,jppic,Kbb) * rfact2 281 282 ! Uptake of nitrogen 282 283 zrat = MIN( 1., zration / (xqnpmax(ji,jj,jk) + rtrn) ) … … 310 311 IF( etot_ndcy(ji,jj,jk) > 1.E-3 ) THEN 311 312 ! production terms for diatomees 312 zprorcad(ji,jj,jk) = zprdia(ji,jj,jk) * xlimdia(ji,jj,jk) * tr b(ji,jj,jk,jpdia) * rfact2313 zprorcad(ji,jj,jk) = zprdia(ji,jj,jk) * xlimdia(ji,jj,jk) * tr(ji,jj,jk,jpdia,Kbb) * rfact2 313 314 ! Computation of the respiration term according to pahlow 314 315 ! & oschlies (2013) 315 316 ! 316 zration = tr b(ji,jj,jk,jpndi) / ( trb(ji,jj,jk,jpdia) + rtrn )317 zratiop = tr b(ji,jj,jk,jppdi) / ( trb(ji,jj,jk,jpdia) + rtrn )318 zratiof = tr b(ji,jj,jk,jpdfe) / ( trb(ji,jj,jk,jpdia) + rtrn )319 zprnutmax = zprnut(ji,jj,jk) * fvduptk(ji,jj,jk) / rno3 * tr b(ji,jj,jk,jpdia) * rfact2317 zration = tr(ji,jj,jk,jpndi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 318 zratiop = tr(ji,jj,jk,jppdi,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 319 zratiof = tr(ji,jj,jk,jpdfe,Kbb) / ( tr(ji,jj,jk,jpdia,Kbb) + rtrn ) 320 zprnutmax = zprnut(ji,jj,jk) * fvduptk(ji,jj,jk) / rno3 * tr(ji,jj,jk,jpdia,Kbb) * rfact2 320 321 ! Uptake of nitrogen 321 322 zrat = MIN( 1., zration / (xqndmax(ji,jj,jk) + rtrn) ) … … 350 351 znanotot = enanom(ji,jj,jk) / ( zmxl_chl(ji,jj,jk) + rtrn ) 351 352 zprod = rday * (zpronewn(ji,jj,jk) + zproregn(ji,jj,jk)) * zprchln(ji,jj,jk) * xlimphy(ji,jj,jk) 352 thetannm_n = MIN ( thetannm, ( thetannm / (1. - 1.14 / 43.4 *ts n(ji,jj,jk,jp_tem))) &353 thetannm_n = MIN ( thetannm, ( thetannm / (1. - 1.14 / 43.4 *ts(ji,jj,jk,jp_tem,Kmm))) & 353 354 & * (1. - 1.14 / 43.4 * 20.)) 354 355 zprochln = thetannm_n * zprod / ( zpislopeadn(ji,jj,jk) * znanotot + rtrn ) … … 357 358 zpicotot = epicom(ji,jj,jk) / ( zmxl_chl(ji,jj,jk) + rtrn ) 358 359 zprod = rday * (zpronewp(ji,jj,jk) + zproregp(ji,jj,jk)) * zprchlp(ji,jj,jk) * xlimpic(ji,jj,jk) 359 thetanpm_n = MIN ( thetanpm, ( thetanpm / (1. - 1.14 / 43.4 *ts n(ji,jj,jk,jp_tem))) &360 thetanpm_n = MIN ( thetanpm, ( thetanpm / (1. - 1.14 / 43.4 *ts(ji,jj,jk,jp_tem,Kmm))) & 360 361 & * (1. - 1.14 / 43.4 * 20.)) 361 362 zprochlp = thetanpm_n * zprod / ( zpislopeadp(ji,jj,jk) * zpicotot + rtrn ) … … 364 365 zdiattot = ediatm(ji,jj,jk) / ( zmxl_chl(ji,jj,jk) + rtrn ) 365 366 zprod = rday * (zpronewd(ji,jj,jk) + zproregd(ji,jj,jk)) * zprchld(ji,jj,jk) * xlimdia(ji,jj,jk) 366 thetandm_n = MIN ( thetandm, ( thetandm / (1. - 1.14 / 43.4 *ts n(ji,jj,jk,jp_tem))) &367 thetandm_n = MIN ( thetandm, ( thetandm / (1. - 1.14 / 43.4 *ts(ji,jj,jk,jp_tem,Kmm))) & 367 368 & * (1. - 1.14 / 43.4 * 20.)) 368 369 zprochld = thetandm_n * zprod / ( zpislopeadd(ji,jj,jk) * zdiattot + rtrn ) 369 370 zprochld = MAX(zprochld, chlcmin * 12. * zprorcad(ji,jj,jk) ) 370 371 ! Update the arrays TRA which contain the Chla sources and sinks 371 tr a(ji,jj,jk,jpnch) = tra(ji,jj,jk,jpnch) + zprochln * texcretn372 tr a(ji,jj,jk,jpdch) = tra(ji,jj,jk,jpdch) + zprochld * texcretd373 tr a(ji,jj,jk,jppch) = tra(ji,jj,jk,jppch) + zprochlp * texcretp372 tr(ji,jj,jk,jpnch,Krhs) = tr(ji,jj,jk,jpnch,Krhs) + zprochln * texcretn 373 tr(ji,jj,jk,jpdch,Krhs) = tr(ji,jj,jk,jpdch,Krhs) + zprochld * texcretd 374 tr(ji,jj,jk,jppch,Krhs) = tr(ji,jj,jk,jppch,Krhs) + zprochlp * texcretp 374 375 ENDIF 375 376 END DO … … 386 387 zdocprod = excretd * zprorcad(ji,jj,jk) + excretn * zprorcan(ji,jj,jk) & 387 388 & + excretp * zprorcap(ji,jj,jk) 388 tr a(ji,jj,jk,jppo4) = tra(ji,jj,jk,jppo4) - zpropo4n(ji,jj,jk) - zpropo4d(ji,jj,jk) &389 tr(ji,jj,jk,jppo4,Krhs) = tr(ji,jj,jk,jppo4,Krhs) - zpropo4n(ji,jj,jk) - zpropo4d(ji,jj,jk) & 389 390 & - zpropo4p(ji,jj,jk) 390 tr a(ji,jj,jk,jpno3) = tra(ji,jj,jk,jpno3) - zpronewn(ji,jj,jk) - zpronewd(ji,jj,jk) &391 tr(ji,jj,jk,jpno3,Krhs) = tr(ji,jj,jk,jpno3,Krhs) - zpronewn(ji,jj,jk) - zpronewd(ji,jj,jk) & 391 392 & - zpronewp(ji,jj,jk) 392 tr a(ji,jj,jk,jpnh4) = tra(ji,jj,jk,jpnh4) - zproregn(ji,jj,jk) - zproregd(ji,jj,jk) &393 tr(ji,jj,jk,jpnh4,Krhs) = tr(ji,jj,jk,jpnh4,Krhs) - zproregn(ji,jj,jk) - zproregd(ji,jj,jk) & 393 394 & - zproregp(ji,jj,jk) 394 tr a(ji,jj,jk,jpphy) = tra(ji,jj,jk,jpphy) + zprorcan(ji,jj,jk) * texcretn &395 tr(ji,jj,jk,jpphy,Krhs) = tr(ji,jj,jk,jpphy,Krhs) + zprorcan(ji,jj,jk) * texcretn & 395 396 & - zpsino3 * zpronewn(ji,jj,jk) - zpsinh4 * zproregn(ji,jj,jk) & 396 397 & - zrespn(ji,jj,jk) 397 zcroissn(ji,jj,jk) = tr a(ji,jj,jk,jpphy) / rfact2/ (trb(ji,jj,jk,jpphy) + rtrn)398 tr a(ji,jj,jk,jpnph) = tra(ji,jj,jk,jpnph) + zprontot * texcretn399 tr a(ji,jj,jk,jppph) = tra(ji,jj,jk,jppph) + zpropo4n(ji,jj,jk) * texcretn &398 zcroissn(ji,jj,jk) = tr(ji,jj,jk,jpphy,Krhs) / rfact2/ (tr(ji,jj,jk,jpphy,Kbb) + rtrn) 399 tr(ji,jj,jk,jpnph,Krhs) = tr(ji,jj,jk,jpnph,Krhs) + zprontot * texcretn 400 tr(ji,jj,jk,jppph,Krhs) = tr(ji,jj,jk,jppph,Krhs) + zpropo4n(ji,jj,jk) * texcretn & 400 401 & + zprodopn(ji,jj,jk) * texcretn 401 tr a(ji,jj,jk,jpnfe) = tra(ji,jj,jk,jpnfe) + zprofen(ji,jj,jk) * texcretn402 tr a(ji,jj,jk,jppic) = tra(ji,jj,jk,jppic) + zprorcap(ji,jj,jk) * texcretp &402 tr(ji,jj,jk,jpnfe,Krhs) = tr(ji,jj,jk,jpnfe,Krhs) + zprofen(ji,jj,jk) * texcretn 403 tr(ji,jj,jk,jppic,Krhs) = tr(ji,jj,jk,jppic,Krhs) + zprorcap(ji,jj,jk) * texcretp & 403 404 & - zpsino3 * zpronewp(ji,jj,jk) - zpsinh4 * zproregp(ji,jj,jk) & 404 405 & - zrespp(ji,jj,jk) 405 zcroissp(ji,jj,jk) = tr a(ji,jj,jk,jppic) / rfact2/ (trb(ji,jj,jk,jppic) + rtrn)406 tr a(ji,jj,jk,jpnpi) = tra(ji,jj,jk,jpnpi) + zproptot * texcretp407 tr a(ji,jj,jk,jpppi) = tra(ji,jj,jk,jpppi) + zpropo4p(ji,jj,jk) * texcretp &406 zcroissp(ji,jj,jk) = tr(ji,jj,jk,jppic,Krhs) / rfact2/ (tr(ji,jj,jk,jppic,Kbb) + rtrn) 407 tr(ji,jj,jk,jpnpi,Krhs) = tr(ji,jj,jk,jpnpi,Krhs) + zproptot * texcretp 408 tr(ji,jj,jk,jpppi,Krhs) = tr(ji,jj,jk,jpppi,Krhs) + zpropo4p(ji,jj,jk) * texcretp & 408 409 & + zprodopp(ji,jj,jk) * texcretp 409 tr a(ji,jj,jk,jppfe) = tra(ji,jj,jk,jppfe) + zprofep(ji,jj,jk) * texcretp410 tr a(ji,jj,jk,jpdia) = tra(ji,jj,jk,jpdia) + zprorcad(ji,jj,jk) * texcretd &410 tr(ji,jj,jk,jppfe,Krhs) = tr(ji,jj,jk,jppfe,Krhs) + zprofep(ji,jj,jk) * texcretp 411 tr(ji,jj,jk,jpdia,Krhs) = tr(ji,jj,jk,jpdia,Krhs) + zprorcad(ji,jj,jk) * texcretd & 411 412 & - zpsino3 * zpronewd(ji,jj,jk) - zpsinh4 * zproregd(ji,jj,jk) & 412 413 & - zrespd(ji,jj,jk) 413 zcroissd(ji,jj,jk) = tr a(ji,jj,jk,jpdia) / rfact2 / (trb(ji,jj,jk,jpdia) + rtrn)414 tr a(ji,jj,jk,jpndi) = tra(ji,jj,jk,jpndi) + zprodtot * texcretd415 tr a(ji,jj,jk,jppdi) = tra(ji,jj,jk,jppdi) + zpropo4d(ji,jj,jk) * texcretd &414 zcroissd(ji,jj,jk) = tr(ji,jj,jk,jpdia,Krhs) / rfact2 / (tr(ji,jj,jk,jpdia,Kbb) + rtrn) 415 tr(ji,jj,jk,jpndi,Krhs) = tr(ji,jj,jk,jpndi,Krhs) + zprodtot * texcretd 416 tr(ji,jj,jk,jppdi,Krhs) = tr(ji,jj,jk,jppdi,Krhs) + zpropo4d(ji,jj,jk) * texcretd & 416 417 & + zprodopd(ji,jj,jk) * texcretd 417 tr a(ji,jj,jk,jpdfe) = tra(ji,jj,jk,jpdfe) + zprofed(ji,jj,jk) * texcretd418 tr a(ji,jj,jk,jpdsi) = tra(ji,jj,jk,jpdsi) + zprorcad(ji,jj,jk) * zysopt(ji,jj,jk) * texcretd419 tr a(ji,jj,jk,jpdoc) = tra(ji,jj,jk,jpdoc) + excretd * zprorcad(ji,jj,jk) + excretn * zprorcan(ji,jj,jk) &418 tr(ji,jj,jk,jpdfe,Krhs) = tr(ji,jj,jk,jpdfe,Krhs) + zprofed(ji,jj,jk) * texcretd 419 tr(ji,jj,jk,jpdsi,Krhs) = tr(ji,jj,jk,jpdsi,Krhs) + zprorcad(ji,jj,jk) * zysopt(ji,jj,jk) * texcretd 420 tr(ji,jj,jk,jpdoc,Krhs) = tr(ji,jj,jk,jpdoc,Krhs) + excretd * zprorcad(ji,jj,jk) + excretn * zprorcan(ji,jj,jk) & 420 421 & + excretp * zprorcap(ji,jj,jk) 421 tr a(ji,jj,jk,jpdon) = tra(ji,jj,jk,jpdon) + excretd * zprodtot + excretn * zprontot &422 tr(ji,jj,jk,jpdon,Krhs) = tr(ji,jj,jk,jpdon,Krhs) + excretd * zprodtot + excretn * zprontot & 422 423 & + excretp * zproptot 423 tr a(ji,jj,jk,jpdop) = tra(ji,jj,jk,jpdop) + excretd * zpropo4d(ji,jj,jk) + excretn * zpropo4n(ji,jj,jk) &424 tr(ji,jj,jk,jpdop,Krhs) = tr(ji,jj,jk,jpdop,Krhs) + excretd * zpropo4d(ji,jj,jk) + excretn * zpropo4n(ji,jj,jk) & 424 425 & - texcretn * zprodopn(ji,jj,jk) - texcretd * zprodopd(ji,jj,jk) + excretp * zpropo4p(ji,jj,jk) & 425 426 & - texcretp * zprodopp(ji,jj,jk) 426 tr a(ji,jj,jk,jpoxy) = tra(ji,jj,jk,jpoxy) + o2ut * ( zproregn(ji,jj,jk) + zproregd(ji,jj,jk) &427 tr(ji,jj,jk,jpoxy,Krhs) = tr(ji,jj,jk,jpoxy,Krhs) + o2ut * ( zproregn(ji,jj,jk) + zproregd(ji,jj,jk) & 427 428 & + zproregp(ji,jj,jk) ) + ( o2ut + o2nit ) * ( zpronewn(ji,jj,jk) & 428 429 & + zpronewd(ji,jj,jk) + zpronewp(ji,jj,jk) ) & 429 430 & - o2ut * ( zrespn(ji,jj,jk) + zrespp(ji,jj,jk) + zrespd(ji,jj,jk) ) 430 431 zfeup = texcretn * zprofen(ji,jj,jk) + texcretd * zprofed(ji,jj,jk) + texcretp * zprofep(ji,jj,jk) 431 tr a(ji,jj,jk,jpfer) = tra(ji,jj,jk,jpfer) - zfeup432 tr a(ji,jj,jk,jpsil) = tra(ji,jj,jk,jpsil) - texcretd * zprorcad(ji,jj,jk) * zysopt(ji,jj,jk)433 tr a(ji,jj,jk,jpdic) = tra(ji,jj,jk,jpdic) - zprorcan(ji,jj,jk) - zprorcad(ji,jj,jk) - zprorcap(ji,jj,jk) &432 tr(ji,jj,jk,jpfer,Krhs) = tr(ji,jj,jk,jpfer,Krhs) - zfeup 433 tr(ji,jj,jk,jpsil,Krhs) = tr(ji,jj,jk,jpsil,Krhs) - texcretd * zprorcad(ji,jj,jk) * zysopt(ji,jj,jk) 434 tr(ji,jj,jk,jpdic,Krhs) = tr(ji,jj,jk,jpdic,Krhs) - zprorcan(ji,jj,jk) - zprorcad(ji,jj,jk) - zprorcap(ji,jj,jk) & 434 435 & + zpsino3 * zpronewn(ji,jj,jk) + zpsinh4 * zproregn(ji,jj,jk) & 435 436 & + zpsino3 * zpronewp(ji,jj,jk) + zpsinh4 * zproregp(ji,jj,jk) & 436 437 & + zpsino3 * zpronewd(ji,jj,jk) + zpsinh4 * zproregd(ji,jj,jk) & 437 438 & + zrespn(ji,jj,jk) + zrespd(ji,jj,jk) + zrespp(ji,jj,jk) 438 tr a(ji,jj,jk,jptal) = tra(ji,jj,jk,jptal) + rno3 * ( zpronewn(ji,jj,jk) + zpronewd(ji,jj,jk) &439 tr(ji,jj,jk,jptal,Krhs) = tr(ji,jj,jk,jptal,Krhs) + rno3 * ( zpronewn(ji,jj,jk) + zpronewd(ji,jj,jk) & 439 440 & + zpronewp(ji,jj,jk) ) - rno3 * ( zproregn(ji,jj,jk) + zproregd(ji,jj,jk) & 440 441 & + zproregp(ji,jj,jk) ) … … 450 451 zdocprod = excretd * zprorcad(ji,jj,jk) + excretn * zprorcan(ji,jj,jk) + excretp * zprorcap(ji,jj,jk) 451 452 zfeup = texcretn * zprofen(ji,jj,jk) + texcretd * zprofed(ji,jj,jk) + texcretp * zprofep(ji,jj,jk) 452 tr a(ji,jj,jk,jplgw) = tra(ji,jj,jk,jplgw) + zdocprod * ldocp - zfeup * plig(ji,jj,jk) * lthet453 tr(ji,jj,jk,jplgw,Krhs) = tr(ji,jj,jk,jplgw,Krhs) + zdocprod * ldocp - zfeup * plig(ji,jj,jk) * lthet 453 454 zpligprod1(ji,jj,jk) = zdocprod * ldocp 454 455 zpligprod2(ji,jj,jk) = zfeup * plig(ji,jj,jk) * lthet … … 557 558 WRITE(charout, FMT="('prod')") 558 559 CALL prt_ctl_trc_info(charout) 559 CALL prt_ctl_trc(tab4d=tr a, mask=tmask, clinfo=ctrcnm)560 CALL prt_ctl_trc(tab4d=tr(:,:,:,:,Krhs), mask=tmask, clinfo=ctrcnm) 560 561 ENDIF 561 562 !
Note: See TracChangeset
for help on using the changeset viewer.