gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

EQUATIONS: Avoid division errors in nr_iterate_srhd_adi_2d().

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2015-02-19 10:01:05 -02:00
parent 020142137d
commit f46981fb03
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/equations.F90
View File

@ -3646,7 +3646,7 @@ module equations
!
logical :: keep
integer :: it, cn
real(kind=8) :: ww, vm, gm
real(kind=8) :: ww, vm, vs
real(kind=8) :: pr, dpw, dpv
real(kind=8) :: f, dfw, dfv, df
real(kind=8) :: g, dgw, dgv, dg
@ -3686,7 +3686,7 @@ module equations
!
ww = w * w
vm = 1.0d+00 - vv
gm = 1.0d+00 / sqrt(vm)
vs = sqrt(vm)
! calculate the thermal pressure and its derivatives
!
@ -3694,13 +3694,13 @@ module equations
! dP/dW = (γ - 1)/γ (1 - |V|²)
! dP/d|V|² = (γ - 1)/γ (- W - 1/2 D Γ)
!
pr = gammaxi * (w - dn * gm) * vm
pr = gammaxi * (w * vm - dn * vs)
dpw = gammaxi * vm
dpv = gammaxi * (- w + 0.5d+00 * dn * gm)
dpv = gammaxi * (- w + 0.5d+00 * dn / vs)
! calculate F(W,|V|²) and G(W,|V|²)
!
f = w - pr - en
f = w - en - pr
g = vv * ww - mm
! calculate dF(W,|V|²)/dW and dF(W,|V|²)/d|V|²