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

EQUATIONS: Print warning in the case of the unphysical variables.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2015-02-22 19:14:40 -03:00
parent d9d73e2cf6
commit e077a985ea
1 changed files with 28 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

28
src/equations.F90
View File

@ -3821,6 +3821,7 @@ module equations
it = it - 1
end do
if (it <= 0) then
print *, ''
print *, 'no initial brackets found', wl, wu, fl, fu
info = .false.
return
@ -3886,6 +3887,19 @@ module equations
w = w - dw
vv = vv - dv
! check if the new enthalpy gives physical pressure and velocity
!
if (w < wl) then
print *, '[SRHD, 2D(W,v²)] Enthalpy too small: ', w, wl
info = .false.
return
end if
if (vv < 0.0d+00 .or. vv >= 1.0d+00) then
print *, '[SRHD, 2D(W,v²)] Unphysical velocity |v|²: ', vv
info = .false.
return
end if
! calculate the normalized error
!
err = max(abs(dw / w), abs(dv))
@ -4022,6 +4036,7 @@ module equations
it = it - 1
end do
if (it <= 0) then
print *, ''
print *, 'no initial brackets found', wl, wu, fl, fu
info = .false.
return
@ -4087,6 +4102,19 @@ module equations
w = w - dw
uu = uu - du
! check if the new enthalpy gives physical pressure
!
if (w < wl) then
print *, '[SRHD, 2D(W,u²)] Enthalpy too small: ', w, wl
info = .false.
return
end if
if (uu < 0.0d+00) then
print *, '[SRHD, 2D(W,u²)] Unphysical velocity |u|²: ', uu
info = .false.
return
end if
! calculate the normalized error
!
err = max(abs(dw / w), abs(du))