From acbab971f3b220499d69dc682882320203b18692 Mon Sep 17 00:00:00 2001 From: Grzegorz Kowal Date: Wed, 15 Mar 2023 16:57:18 -0300 Subject: [PATCH] USER_PROBLEM: Print the Prandtl number as well. Signed-off-by: Grzegorz Kowal --- sources/user_problem.F90 | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/sources/user_problem.F90 b/sources/user_problem.F90 index 7263773..08ab662 100644 --- a/sources/user_problem.F90 +++ b/sources/user_problem.F90 @@ -35,6 +35,7 @@ module user_problem real(kind=8), save :: beta = 1.00d+00 real(kind=8), save :: zeta = 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 :: bamp = 1.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 :: valf = 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 :: tdec = 1.00d+00 @@ -165,6 +167,7 @@ module user_problem status = 1 return end if + call get_parameter("viscosity" , nu ) call get_parameter("resistivity", eta) if (eta < 0.0d+00) then if (verbose) & @@ -213,6 +216,7 @@ module user_problem csnd = sqrt(csnd2) valf = sqrt(bamp**2 / dens) lund = valf / max(tiny(eta), eta) + prdl = nu / max(tiny(eta), eta) dlta = lund**(-alpha) ! fill up the boundary values @@ -411,6 +415,7 @@ module user_problem if (resistive) then call print_parameter(verbose, '( ) S (Lundquist number)', lund) end if + call print_parameter(verbose, '( ) P (Prandtl number)', prdl) call print_parameter(verbose, '(*) delta (thickness)', dlta) call print_parameter(verbose, '(*) perturbation', perturbation) call print_parameter(verbose, '(*) bper' , bper)