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

USER_PROBLEM: Guide field does not determine the Alfven speed.

Browse Source 
The Alfven speed is calculated from the upstream reconnecting magnetic
field component only.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2023-03-15 16:52:50 -03:00
parent 7c65c0faa4
commit 3deaa1e9fb
1 changed files with 1 additions and 1 deletions

2
sources/user_problem.F90

@ -211,7 +211,7 @@ module user_problem
ptot = pres + pmag
csnd2 = pres / dens
csnd = sqrt(csnd2)
valf = sqrt(2.0d+00 * pmag / dens)
valf = sqrt(bamp**2 / dens)
lund = valf / max(tiny(eta), eta)
dlta = lund**(-alpha)