DISPATCH
hydro_parameters.f90
1 #define MAXBOUND 8
2 #define MAXREGION 8
3 #define MAXLEVEL 30
4 
5 module hydro_parameters
6  use amr_parameters
7 
8  ! Number of independant variables
9 #ifndef NVAR
10  integer,parameter::nvar=ndim+2
11 #else
12  integer,parameter::nvar=nvar
13 #endif
14  ! Size of hydro kernel
15  integer,parameter::iu1=-1
16  integer,parameter::iu2=+4
17  integer,parameter::ju1=(1-ndim/2)-1*(ndim/2)
18  integer,parameter::ju2=(1-ndim/2)+4*(ndim/2)
19  integer,parameter::ku1=(1-ndim/3)-1*(ndim/3)
20  integer,parameter::ku2=(1-ndim/3)+4*(ndim/3)
21  integer,parameter::if1=1
22  integer,parameter::if2=3
23  integer,parameter::jf1=1
24  integer,parameter::jf2=(1-ndim/2)+3*(ndim/2)
25  integer,parameter::kf1=1
26  integer,parameter::kf2=(1-ndim/3)+3*(ndim/3)
27 
28  ! Imposed boundary condition variables
29  real(dp),dimension(1:MAXBOUND,1:nvar)::boundary_var
30  real(dp),dimension(1:MAXBOUND)::d_bound=0.0d0
31  real(dp),dimension(1:MAXBOUND)::p_bound=0.0d0
32  real(dp),dimension(1:MAXBOUND)::u_bound=0.0d0
33  real(dp),dimension(1:MAXBOUND)::v_bound=0.0d0
34  real(dp),dimension(1:MAXBOUND)::w_bound=0.0d0
35 
36  ! Refinement parameters for hydro
37  real(dp)::err_grad_d=-1.0 ! Density gradient
38  real(dp)::err_grad_u= 0.0 ! Velocity gradient
39  real(dp)::err_grad_p=-1.0 ! Pressure gradient
40  real(dp)::floor_d=1.d-10 ! Density floor
41  real(dp)::floor_u=1.d-10 ! Velocity floor
42  real(dp)::floor_p=1.d-10 ! Pressure floor
43  real(dp)::mass_sph=0.0d0 ! mass_sph
44  real(dp),dimension(1:MAXLEVEL)::jeans_refine=-1.0
45 
46  ! Initial conditions hydro variables
47  real(dp),dimension(1:MAXREGION)::d_region=0.
48  real(dp),dimension(1:MAXREGION)::u_region=0.
49  real(dp),dimension(1:MAXREGION)::v_region=0.
50  real(dp),dimension(1:MAXREGION)::w_region=0.
51  real(dp),dimension(1:MAXREGION)::p_region=0.
52 
53  ! Hydro solver parameters
54  integer ::niter_riemann=10
55  logical:: isothermal=.false.
56  real(dp)::slope_type=1
57  real(dp)::gamma=1.4d0
58  real(dp)::courant_factor=0.5d0
59  real(dp)::difmag=0.0d0
60  real(dp)::smallc=1.d-10
61  real(dp)::smallr=1.d-10
62  character(LEN=10)::scheme='muscl'
63  character(LEN=10)::riemann='llf'
64 
65  ! Interpolation parameters
66  integer ::interpol_var=0
67  integer ::interpol_type=1
68 
69  ! Passive variables index
70  integer::imetal=6
71  integer::idelay=6
72  integer::ixion=6
73  integer::ichem=6
74 
75 end module hydro_parameters