USER_PROBLEM: Calculate Alfvén speed using bamp only.

This is the Alfvén speed along the X direction defining the velocity
scale, time scale, and the Lundquist number.

Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2022-11-27 18:16:24 -03:00
parent 20492048aa
commit 5793fe7810

View File

@ -145,7 +145,7 @@ module user_problem
csnd2 = pres / dens
end if
csnd = sqrt(csnd2)
valf = sqrt(2.0d+00 * pmag / dens)
valf = bamp / sqrt(dens)
lund = valf / max(tiny(eta), eta)
! determine if to append or create another file