Commit 580c08aca9af55392a331da2402211cae4a31c41

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

Correction d'un bog avec les termes lies a la mte (mte=.true.) dans les equation…

…s de reaction. Le probleme faisait en sorte que le modele n'etait pas conservatif.
Showing 1 changed file with 2 additions and 2 deletions   Show diff stats
src/extras/bio/bio_ismer.F90
@@ -673,9 +673,9 @@ @@ -673,9 +673,9 @@
673 dd(p2,l,ci) =amr2*gamma*ff2*(qn2+qa2)*cc(p2,ci) 673 dd(p2,l,ci) =amr2*gamma*ff2*(qn2+qa2)*cc(p2,ci)
674 674
675 dd(p1,d,ci) =amr1*mu11*(cc(p1,ci)+p0)/(k5+cc(p1,ci)+p0)*cc(p1,ci) & 675 dd(p1,d,ci) =amr1*mu11*(cc(p1,ci)+p0)/(k5+cc(p1,ci)+p0)*cc(p1,ci) &
676 - +(1.-beta)*cc(p1,ci)**2*(g1max*r11*fac1+g2max*r21*fac2) 676 + +(1.-beta)*cc(p1,ci)**2*(hmr1*g1max*r11*fac1+hmr2*g2max*r21*fac2)
677 dd(p2,d,ci) =amr2*mu12*(cc(p2,ci)+p0)/(k5+cc(p2,ci)+p0)*cc(p2,ci) & 677 dd(p2,d,ci) =amr2*mu12*(cc(p2,ci)+p0)/(k5+cc(p2,ci)+p0)*cc(p2,ci) &
678 - +(1.-beta)*g2max*r21*cc(p2,ci)**2*fac2 678 + +(1.-beta)*cc(p2,ci)**2*(hmr1*g1max*r11*fac1+hmr2*g2max*r21*fac2)
679 679
680 dd(b,d,ci) =(1.-beta)*g1max*r12*cc(b,ci)**2*fac1 680 dd(b,d,ci) =(1.-beta)*g1max*r12*cc(b,ci)**2*fac1
681 681