Commit b0a803da18416bfc5f5c66a9f4bf9d94bb532d5a

Authored by dumoda01
1 parent 8473db4a
Exists in master and in 1 other branch snow

Complement de modification a la revision precedente sur l'harmonisation de la pr…

…ise en charge du nitrate et de l'ammonium par le module bio.
Showing 1 changed file with 22 additions and 7 deletions   Show diff stats
src/extras/bio/bio.F90
... ... @@ -565,16 +565,31 @@
565 565 if (bio_eulerian) then
566 566 do j=1,numcc
567 567  
568   - if(j==1) then !CHG3
569   - do k=1,nlev
570   - cc(j,k) = nit(k)
  568 + if (j==1 .and. bio_model==6) then
  569 + do i=1,nlev
  570 + cc(j,i) = nit(i)
  571 + end do
  572 + else if (j==9 .and. bio_model==6) then
  573 + do i=1,nlev
  574 + cc(j,i) = amm(i)
  575 + end do
  576 + else if (j==5 .and. bio_model==4) then
  577 + do i=1,nlev
  578 + cc(j,i) = nit(i)
571 579 end do
572   - else if (j==9) then !CHG5
573   - do k=1,nlev
574   - cc(j,k) = amm(k)
  580 + else if (j==6 .and. bio_model==4) then
  581 + do i=1,nlev
  582 + cc(j,i) = amm(i)
  583 + end do
  584 + else if (j==7 .and. bio_model==2) then
  585 + do i=1,nlev
  586 + cc(7,i) = nit(i)
  587 + end do
  588 + else if (j==6 .and. bio_model==2) then
  589 + do i=1,nlev
  590 + cc(6,i) = amm(i)
575 591 end do
576 592 else
577   -
578 593 ! do advection step due to settling or rising
579 594 call adv_center(nlev,dt,h,h,ws(j,:),flux, &
580 595 flux,_ZERO_,_ZERO_,w_adv_discr,adv_mode_1,cc(j,:))
... ...