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

USER_PROBLEM: Print the Prandtl number as well.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2023-03-15 16:57:18 -03:00
parent 3deaa1e9fb
commit acbab971f3
1 changed files with 5 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
sources/user_problem.F90
View File

@ -35,6 +35,7 @@ module user_problem
real(kind=8), save :: beta = 1.00d+00 real(kind=8), save :: beta = 1.00d+00
real(kind=8), save :: zeta = 0.00d+00 real(kind=8), save :: zeta = 0.00d+00
real(kind=8), save :: eta = 0.00d+00 real(kind=8), save :: eta = 0.00d+00
real(kind=8), save :: nu = 0.00d+00
real(kind=8), save :: dens = 1.00d+00 real(kind=8), save :: dens = 1.00d+00
real(kind=8), save :: bamp = 1.00d+00 real(kind=8), save :: bamp = 1.00d+00
real(kind=8), save :: bgui = 0.00d+00 real(kind=8), save :: bgui = 0.00d+00
@ -44,6 +45,7 @@ module user_problem
real(kind=8), save :: ptot = 1.00d+00 real(kind=8), save :: ptot = 1.00d+00
real(kind=8), save :: valf = 1.00d+00 real(kind=8), save :: valf = 1.00d+00
real(kind=8), save :: lund = 1.00d+00 real(kind=8), save :: lund = 1.00d+00
real(kind=8), save :: prdl = 0.00d+00
real(kind=8), save :: dlta = 1.00d-16 real(kind=8), save :: dlta = 1.00d-16
real(kind=8), save :: tdec = 1.00d+00 real(kind=8), save :: tdec = 1.00d+00
@ -165,6 +167,7 @@ module user_problem
status = 1 status = 1
return return
end if end if
call get_parameter("viscosity" , nu )
call get_parameter("resistivity", eta) call get_parameter("resistivity", eta)
if (eta < 0.0d+00) then if (eta < 0.0d+00) then
if (verbose) & if (verbose) &
@ -213,6 +216,7 @@ module user_problem
csnd = sqrt(csnd2) csnd = sqrt(csnd2)
valf = sqrt(bamp**2 / dens) valf = sqrt(bamp**2 / dens)
lund = valf / max(tiny(eta), eta) lund = valf / max(tiny(eta), eta)
prdl = nu / max(tiny(eta), eta)
dlta = lund**(-alpha) dlta = lund**(-alpha)
! fill up the boundary values ! fill up the boundary values
@ -411,6 +415,7 @@ module user_problem
if (resistive) then if (resistive) then
call print_parameter(verbose, '( ) S (Lundquist number)', lund) call print_parameter(verbose, '( ) S (Lundquist number)', lund)
end if end if
call print_parameter(verbose, '( ) P (Prandtl number)', prdl)
call print_parameter(verbose, '(*) delta (thickness)', dlta) call print_parameter(verbose, '(*) delta (thickness)', dlta)
call print_parameter(verbose, '(*) perturbation', perturbation) call print_parameter(verbose, '(*) perturbation', perturbation)
call print_parameter(verbose, '(*) bper' , bper) call print_parameter(verbose, '(*) bper' , bper)