bug fix to KM equation

LMNS3d · Jan 22, 2020 · 97e6ab5 · 97e6ab5
1 parent 5844ab0
commit 97e6ab5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/simulation_code/m_bubbles.f90 b/src/simulation_code/m_bubbles.f90
@@ -401,7 +401,7 @@ FUNCTION f_rddot_KM( fpbdot, fgamma_gas, fCp, fCpbw, fRho, fR, fV, fR0, fC )
             cdot_star = tmp1 + ( 2.d0/Web + 4.d0*Re_inv*fV )*fV/(fR**2.d0)
             tmp1 = fV/fC
             tmp2 =  1.5D0*fV**2d0*( tmp1/3.D0-1.D0 ) +                      &
-                    (fCpbw - fCp)/fRho * (1.d0 - tmp1) +      &
+                    (fCpbw - fCp)/fRho * (1.d0 + tmp1) +      &
                     cdot_star * fR/fRho/fC
 
             f_rddot_KM = tmp2/( fR*(1.d0-tmp1) + 4.d0*Re_inv/(fRho*fC) )