added low and med
This commit is contained in:
hlasco
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma00.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma00
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 73.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.6911
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.3089
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma10.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma10
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 71.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.70
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.30
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma20.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma20
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 71.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.70
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.30
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma30.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma30
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 71.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.70
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.30
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=20.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens','temp'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
ycentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
zcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma00.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma10.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma20.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma30.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 100000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 200000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.6
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 100000
|
||||||
|
npart_pot2 200000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 0.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 1.0
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.01
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 12500
|
||||||
|
npart_pot4 25000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 100000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 200000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.6
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 100000
|
||||||
|
npart_pot2 200000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 10.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 0.34
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.01
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 12500
|
||||||
|
npart_pot4 25000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 100000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 200000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.6
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 100000
|
||||||
|
npart_pot2 200000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 20.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 0.34
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.01
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 12500
|
||||||
|
npart_pot4 25000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 100000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 200000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.6
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 100000
|
||||||
|
npart_pot2 200000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 30.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 0.34
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.01
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 12500
|
||||||
|
npart_pot4 25000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/pbox.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/pbox
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 0.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 73.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.6911
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.3089
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=20.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens','temp'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
ycentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
zcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma00.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma10.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.5
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.false.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-5
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma30.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.05
|
||||||
|
tend=10.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=10.
|
||||||
|
movie=.true.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=12
|
||||||
|
movie_vars_txt='stars','dens'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=0.05,0.0,0.0,0.0 0.05,0.0,0.0,0.0
|
||||||
|
ycentre_frame=0.05,0.0,0.0,0.0 0.05,0.0,0.0,0.0
|
||||||
|
zcentre_frame=0.05,0.0,0.0,0.0 0.05,0.0,0.0,0.0
|
||||||
|
deltax_frame=0.05,0.0 0.05,0.0
|
||||||
|
deltay_frame=0.05,0.0 0.05,0.0
|
||||||
|
deltaz_frame=0.01,0.0 0.01,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/pbox/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=12
|
||||||
|
ngridmax=5000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=0.1
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
!&BOUNDARY_PARAMS
|
||||||
|
!nboundary=6
|
||||||
|
!bound_type= 0, 0, 0, 0, 0, 0
|
||||||
|
!ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
!ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
!jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
!jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
!kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
!kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
!!no_inflow=.true.
|
||||||
|
!/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-8
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='pbox.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,85 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 0.0012
|
||||||
|
# Halo spin parameter
|
||||||
|
# lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 0.1
|
||||||
|
#boxsize2 0.05
|
||||||
|
boxsize2 0.025
|
||||||
|
#boxsize4 0.0125
|
||||||
|
boxsize3 0.00625
|
||||||
|
#boxsize6 0.003125
|
||||||
|
#boxsize4 0.0015625
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 1
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 0.1
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.05
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 1: Gas
|
||||||
|
###########################
|
||||||
|
mass_frac1 1.0
|
||||||
|
npart1 5000000
|
||||||
|
npart_pot1 10000000
|
||||||
|
#accept_min1 0.90
|
||||||
|
#accept_max1 0.98
|
||||||
|
model1 15
|
||||||
|
scale_length1 0.1
|
||||||
|
#cut1 10
|
||||||
|
flatx_cut1 0.1
|
||||||
|
flaty_cut1 0.1
|
||||||
|
flatz_cut1 0.1
|
||||||
|
sigma_cut1 0.000000001
|
||||||
|
mcmc_step1 0.15
|
||||||
|
type1 0
|
||||||
|
stream_fraction1 0.00
|
||||||
|
metal1 0.02041
|
||||||
|
dens_fluct1 1
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient1 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init1 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma1 20.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj1 0.1
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss1 0.005
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed1 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq1 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro1 0.3
|
||||||
|
compute_vel1 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly1 1.0
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma00.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma00
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 73.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.6911
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.3089
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma10.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma10
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 71.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.70
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.30
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma20.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma20
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 71.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.70
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.30
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,39 @@
|
|||||||
|
# List of galaxy parameters file
|
||||||
|
# Putting identical galaxies one beside the other allows DICE
|
||||||
|
# to re-use the previous computation,
|
||||||
|
# and thus reduces the global execution time
|
||||||
|
# x[kpc] y[kpc] z[kpc] vx[km/s] vy[km/s] vz[km/s] spin[deg] incl[deg]
|
||||||
|
Galaxy params/agora_sigma30.params 0. 0. 0. 0. 0. 0. 0. 0.
|
||||||
|
# Number of OpenMP threads
|
||||||
|
Nthreads 16
|
||||||
|
# Compute the mean inter-particle distance
|
||||||
|
# Useful to define the softening length in SPH codes
|
||||||
|
MeanPartDist 0
|
||||||
|
# Maximum number of component for a galaxy model
|
||||||
|
MaxCompNumber 64
|
||||||
|
# Output quantities in the rz plane
|
||||||
|
OutputRz 0
|
||||||
|
# Name of the final IC file
|
||||||
|
# The file extension wil be added by DICE
|
||||||
|
Filename out/agora_sigma30
|
||||||
|
# File format of the IC file [Gadget1,Gadget2]
|
||||||
|
ICformat Gadget2
|
||||||
|
# Cosmological parameters for Virial quantities
|
||||||
|
# Default values are taken from Planck cosmology
|
||||||
|
# Redshift
|
||||||
|
Redshift 1.0
|
||||||
|
# Hubble constant
|
||||||
|
H0 71.0
|
||||||
|
# Dark energy density parameter
|
||||||
|
OmegaL 0.70
|
||||||
|
# Matter density parameter
|
||||||
|
OmegaM 0.30
|
||||||
|
# Spatial curvature density paramter
|
||||||
|
OmegaK 0.00
|
||||||
|
# Normalize the sum of mass fractions to one
|
||||||
|
NormMassFact 0
|
||||||
|
# Size of the GSL integration space
|
||||||
|
GslWorkspaceSize 10000
|
||||||
|
# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
|
||||||
|
GslIntegrationScheme 3
|
||||||
|
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.25
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.true.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=14
|
||||||
|
movie_vars_txt='stars','dens','temp'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
ycentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
zcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=14
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-6
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma00.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.25
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.true.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=14
|
||||||
|
movie_vars_txt='stars','dens','temp'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=14
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-6
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma10.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.25
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.true.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=14
|
||||||
|
movie_vars_txt='stars','dens','temp'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=14
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-6
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma20.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,145 @@
|
|||||||
|
&RUN_PARAMS
|
||||||
|
cosmo=.false.
|
||||||
|
pic=.true.
|
||||||
|
poisson=.true.
|
||||||
|
hydro=.true.
|
||||||
|
nrestart=0
|
||||||
|
nremap=25
|
||||||
|
nsubcycle=1,1,2
|
||||||
|
ncontrol=1
|
||||||
|
rt=.false.
|
||||||
|
/
|
||||||
|
|
||||||
|
&OUTPUT_PARAMS
|
||||||
|
delta_tout=0.25
|
||||||
|
tend=50.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&MOVIE_PARAMS
|
||||||
|
imovout=2000
|
||||||
|
imov=0
|
||||||
|
tendmov=50.
|
||||||
|
movie=.true.
|
||||||
|
nw_frame=1080
|
||||||
|
nh_frame=1080
|
||||||
|
levelmax_frame=14
|
||||||
|
movie_vars_txt='stars','dens','temp'
|
||||||
|
proj_axis='zx'
|
||||||
|
xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
|
||||||
|
deltax_frame=30,0.0 30,0.0
|
||||||
|
deltay_frame=30,0.0 30,0.0
|
||||||
|
deltaz_frame=30,0.0 30,0.0
|
||||||
|
! zoom_only_frame=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&INIT_PARAMS
|
||||||
|
filetype='dice'
|
||||||
|
initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
|
||||||
|
/
|
||||||
|
|
||||||
|
&UNITS_PARAMS
|
||||||
|
units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
|
||||||
|
units_time = 0.470430312423675E+15 ! G=1
|
||||||
|
units_length = 0.308567758128200E+22 ! kpc
|
||||||
|
/
|
||||||
|
|
||||||
|
&AMR_PARAMS
|
||||||
|
levelmin=7
|
||||||
|
levelmax=14
|
||||||
|
ngridmax=3000000
|
||||||
|
npartmax=5000000
|
||||||
|
boxlen=320.0
|
||||||
|
/
|
||||||
|
|
||||||
|
&POISSON_PARAMS
|
||||||
|
epsilon=1e-4
|
||||||
|
/
|
||||||
|
|
||||||
|
&BOUNDARY_PARAMS
|
||||||
|
nboundary=6
|
||||||
|
bound_type= 2, 2, 2, 2, 2, 2
|
||||||
|
ibound_min=-1, 1, -1, -1, -1, -1
|
||||||
|
ibound_max=-1, 1, 1, 1, 1, 1
|
||||||
|
jbound_min= 0, 0, -1, 1, -1, -1
|
||||||
|
jbound_max= 0, 0, -1, 1, 1, 1
|
||||||
|
kbound_min= 0, 0, 0, 0, -1, 1
|
||||||
|
kbound_max= 0, 0, 0, 0, -1, 1
|
||||||
|
no_inflow=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&COOLING_PARAMS
|
||||||
|
cooling=.true.
|
||||||
|
metal=.true.
|
||||||
|
z_ave=1.0 ! To explore explore [0.001, 1]
|
||||||
|
self_shielding=.true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&SF_PARAMS
|
||||||
|
sf_virial=.true. ! Enable turbulent SF
|
||||||
|
sf_model=5 ! See mom2/star_formation.f90
|
||||||
|
n_star=0.01 ! SF density threshold. [0.01, 10]
|
||||||
|
eps_star=1.0 ! SF efficency. Should be 1 with sf_virial=.true.
|
||||||
|
T2_star=10. ! ISM polytropic temperature
|
||||||
|
g_star=1.0 ! ISM polytropic index
|
||||||
|
m_star=10.0 ! Star particle mass
|
||||||
|
sf_log_properties=.false. ! do not set this to .true.
|
||||||
|
/
|
||||||
|
|
||||||
|
&FEEDBACK_PARAMS
|
||||||
|
momentum_feedback=1 ! SN feedback model (from patch/mom2)
|
||||||
|
delayed_cooling=.false. ! Cooling starts from t=0
|
||||||
|
yield=0.177703 ! SN metal yield
|
||||||
|
eta_sn=0.162637 ! SN efficiency
|
||||||
|
t_sne=3.0 ! time for a star to turn into SN
|
||||||
|
f_w=0.
|
||||||
|
/
|
||||||
|
|
||||||
|
&HYDRO_PARAMS
|
||||||
|
gamma=1.6666
|
||||||
|
courant_factor=0.5
|
||||||
|
slope_type=1
|
||||||
|
pressure_fix=.true. ! Scheme to prevent negative T
|
||||||
|
beta_fix=0.5 ! Controls threshold at which energy is
|
||||||
|
! truncated
|
||||||
|
riemann='hllc'
|
||||||
|
/
|
||||||
|
|
||||||
|
&REFINE_PARAMS
|
||||||
|
mass_sph=1e-6
|
||||||
|
m_refine=10*8.
|
||||||
|
interpol_type=0
|
||||||
|
interpol_var=1
|
||||||
|
/
|
||||||
|
|
||||||
|
&DICE_PARAMS
|
||||||
|
ic_file='agora_sigma30.g2' ! Name of the initial conditions file
|
||||||
|
ic_nfile=1 ! If greater than one, look for files with name matching ic_file//'.n'
|
||||||
|
ic_ifout=1 ! Change ramses output index for restarts
|
||||||
|
ic_format='Gadget2' ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
|
||||||
|
ic_center=0.0,0.0,0.0 ! Shift center parameter. ICs are automatically shifted with boxlen/2
|
||||||
|
ic_scale_pos=1.0 ! Scaling factor for the position vector
|
||||||
|
ic_scale_vel=1.0 ! Scaling factor for the velocity vector
|
||||||
|
ic_scale_mass=1.0 ! Scaling factor for the mass
|
||||||
|
ic_scale_u=1.0 ! Scaling factor for the internal energy
|
||||||
|
ic_scale_age=1.0 ! Scaling factor for the particles age
|
||||||
|
ic_scale_metal=1.0 ! Scaling factor for the metallicity
|
||||||
|
ic_head_name='HEAD' ! Name of the Header datablock (Gadget2 format only)
|
||||||
|
ic_pos_name='POS ' ! Name of the position vector datablock (Gadget2 format only)
|
||||||
|
ic_vel_name='VEL ' ! Name of the velocity vector datablock (Gadget2 format only)
|
||||||
|
ic_mass_name='MASS' ! Name of the mass datablock (Gadget2 format only)
|
||||||
|
ic_id_name='ID ' ! Name of the particle identifier datablock (Gadget2 format only)
|
||||||
|
ic_u_name='U ' ! Name of the internal energy datablock (Gadget2 format only)
|
||||||
|
ic_metal_name='Z ' ! Name of the metallicity datablock (Gadget2 format only)
|
||||||
|
ic_age_name='AGE ' ! Name of the particle age datablock (Gadget2 format only)
|
||||||
|
IG_rho=1.0D-6 ! Intergalactic gas density
|
||||||
|
IG_T2=1.0D7 ! Intergalactic gas temperature
|
||||||
|
IG_metal=0.0 ! Intergalactic gas metallicity
|
||||||
|
amr_struct=.false. ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
|
||||||
|
gadget_scale_l=3.085677581282D21 ! Gadget file length unit
|
||||||
|
gadget_scale_v=1.0D5 ! Gadget file velocity unit
|
||||||
|
gadget_scale_m=1.9891D43 ! Gadget file mass unit
|
||||||
|
gadget_scale_t=3.15360e+13 ! Gadget file time unit
|
||||||
|
ic_skip_type = -1 ! Skip specific particle type
|
||||||
|
/
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 1000000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 2000000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.3
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 1000000
|
||||||
|
npart_pot2 2000000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 0.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 1.0
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.01
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 125000
|
||||||
|
npart_pot4 250000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 1000000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 2000000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.3
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 1000000
|
||||||
|
npart_pot2 2000000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 10.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 0.34
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.005
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 125000
|
||||||
|
npart_pot4 250000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 1000000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 2000000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.3
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 1000000
|
||||||
|
npart_pot2 2000000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 20.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 0.34
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.005
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 125000
|
||||||
|
npart_pot4 250000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
@@ -0,0 +1,218 @@
|
|||||||
|
###################
|
||||||
|
# Global parameters
|
||||||
|
###################
|
||||||
|
# Virial velocity of the galaxy [km/s]
|
||||||
|
v200 0.0
|
||||||
|
# Virial mass of the galaxy [1e10 Msol]
|
||||||
|
# Overrides the v200 parameter
|
||||||
|
m200 112.127
|
||||||
|
# Halo spin parameter
|
||||||
|
lambda 0.04
|
||||||
|
# Refinement level of the potential grid
|
||||||
|
level_coarse 7
|
||||||
|
# Refinement level of the plane plane density grid
|
||||||
|
level_grid_mid_dens 7
|
||||||
|
# Refinement level of the turbulence grid
|
||||||
|
level_grid_turb 7
|
||||||
|
# Refinement level of the gaussian field grid
|
||||||
|
level_grid_dens_fluct 7
|
||||||
|
# Size of the potential grid [kpc]
|
||||||
|
boxsize1 320.0
|
||||||
|
boxsize2 160.0
|
||||||
|
boxsize3 80.0
|
||||||
|
boxsize4 40.0
|
||||||
|
boxsize5 20.0
|
||||||
|
boxsize6 10.0
|
||||||
|
boxsize7 5.0
|
||||||
|
# Dispersion for the Gaussian field fluctuations
|
||||||
|
dens_fluct_sigma 0.50
|
||||||
|
# Physical injection scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_inj 2.00
|
||||||
|
# Physical dissipation scale of the random field fluctuations [kpc]
|
||||||
|
dens_fluct_scale_diss 0.25
|
||||||
|
dens_fluct_seed 1212
|
||||||
|
# Seed for the random number generator
|
||||||
|
seed 1246
|
||||||
|
# Switch to MCMC ntry algorithm to position particles for a value > 1
|
||||||
|
mcmc_ntry 1
|
||||||
|
# Number of iterations to reach hydrostatic equilibrium (zero to deactivate)
|
||||||
|
hydro_eq_niter 3
|
||||||
|
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Components parameters
|
||||||
|
#######################
|
||||||
|
|
||||||
|
#######################
|
||||||
|
# Component 1: Halo
|
||||||
|
#######################
|
||||||
|
# Fraction of the virial mass in the component 1
|
||||||
|
mass_frac1 0.9578424465115449
|
||||||
|
# Number of particles for the component 1
|
||||||
|
npart1 1000000
|
||||||
|
# Target mass of individual particles [Msol]
|
||||||
|
#part_mass1 2e6
|
||||||
|
# Number of particles for the potential computation
|
||||||
|
npart_pot1 2000000
|
||||||
|
# Target mass of individual particles for potential computation [Msol]
|
||||||
|
#part_mass_pot1 1e6
|
||||||
|
# Component 1 density model
|
||||||
|
# Available models:
|
||||||
|
# 1 = Exponential disk + sech-z profile
|
||||||
|
# 2 = Myamoto-Nagai profile
|
||||||
|
# 3 = Exponential disk + exponential-z profile
|
||||||
|
# 4 = Hernquist profile
|
||||||
|
# 5 = Plummer profile
|
||||||
|
# 6 = Jaffe profile
|
||||||
|
# 7 = Isothermal profile
|
||||||
|
# 8 = NFW profile
|
||||||
|
# 9 = Burkert
|
||||||
|
# 10 = Einasto profile
|
||||||
|
# 11 = Mestel profile
|
||||||
|
# 12 = Kalnajs profile
|
||||||
|
# 13 = Sersic profile
|
||||||
|
# 14 = Toomre-Kuzmin profile
|
||||||
|
# 15 = Uniform profile
|
||||||
|
# 16 = Pseudo-isothermal profile
|
||||||
|
model1 8
|
||||||
|
# Scale length of the density profile [kpc]
|
||||||
|
scale_length1 1.0
|
||||||
|
# Core radius in the density profile [kpc]
|
||||||
|
rcore1 0.
|
||||||
|
# Concentration parameter
|
||||||
|
# if positive, the scale length is recomputed to match the concentration
|
||||||
|
# the concentration parameter of an NFW halo with a mass M200*mass_frac
|
||||||
|
concentration1 10
|
||||||
|
# Gravitational softening for the poisson solver [kpc]
|
||||||
|
softening1 0.3
|
||||||
|
# Component 1 radial density cut [kpc]
|
||||||
|
cut1 0.0
|
||||||
|
# Component 1 thickness parameters
|
||||||
|
flatx1 1.00
|
||||||
|
flaty1 1.00
|
||||||
|
flatz1 1.00
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm,
|
||||||
|
# expressed in units of the component scale length. Default value is 0.5.
|
||||||
|
mcmc_step1 0.50
|
||||||
|
# Maximum velocity for the component 1 particles in expressed in units of escape velocity
|
||||||
|
vmax_esc1 5.0
|
||||||
|
# Particles type (GADGET format -- 0=Gas,1=Halo,2=Disk,3=Bulge,4=Stars)
|
||||||
|
type1 1
|
||||||
|
# Radius at which the density profile should have the exact same
|
||||||
|
# value as the NFW profile with the previously defined concentration
|
||||||
|
# and with a mass M200*mass_frac [kpc]
|
||||||
|
radius_nfw1 -1.0
|
||||||
|
# Cut the density function in the Jeans equation integration
|
||||||
|
jeans_mass_cut1 0
|
||||||
|
# Compute component velocity (debug option)
|
||||||
|
compute_vel1 1
|
||||||
|
# Number of integral of motion for the Jeans equations
|
||||||
|
# 0 = Do not use Jeans equation for sigma_r and sigma_z (model specidifed by sigmar_model and sigmaz_model)
|
||||||
|
# 1 = Spherically symmetric Jeans equation
|
||||||
|
# 2 = Jeans equations with 2 integrals of motion
|
||||||
|
# 3 = Jeans equations with 3 integrals of motion (solved on a 2D grid in the r-z plane)
|
||||||
|
jeans_dim1 1
|
||||||
|
# Method for streaming velocity computation
|
||||||
|
# 0 = - User defined fixed fraction of the circular velocity profile (stream_fraction keyword)
|
||||||
|
# 1 = Bullock 2001 - Streaming velocity profile following the cumulative mass profile
|
||||||
|
# 2 = Springel 1999 - Streaming velocity profile following is a fixed fraction of the rotation curve
|
||||||
|
# 3 = - Solid Body rotation
|
||||||
|
stream_method1 1
|
||||||
|
# Streaming fraction of the component 1 (only for steam_method=4)
|
||||||
|
stream_fraction1 0.00
|
||||||
|
# Alpha coefficient for the generalized normal distribution [Gaussian=2.0]
|
||||||
|
# for the random generation of velocities
|
||||||
|
ggd_beta1 2.0
|
||||||
|
# Minimum acceptance for the MCMC chain
|
||||||
|
accept_min1 0.95
|
||||||
|
# Maximum acceptance for the MCMC chain
|
||||||
|
accept_max1 0.99
|
||||||
|
|
||||||
|
################################
|
||||||
|
# Component 2: Thin stellar disk
|
||||||
|
################################
|
||||||
|
mass_frac2 0.03066166043860979
|
||||||
|
npart2 1000000
|
||||||
|
npart_pot2 2000000
|
||||||
|
model2 3
|
||||||
|
# If the value is zero and the particle type is not 1 the size is determined using the spin conservation
|
||||||
|
# (Fitting formula from Mo, Mao & White 1998)
|
||||||
|
scale_length2 3.43218
|
||||||
|
cut2 30.
|
||||||
|
flatz2 0.1
|
||||||
|
mcmc_step2 0.3
|
||||||
|
type2 1
|
||||||
|
stream_fraction2 1.00
|
||||||
|
# Epicycle approximation in the Jeans equations
|
||||||
|
epicycle2 1
|
||||||
|
# Minimal value for the Toomre parameter
|
||||||
|
Q_lim2 1.5
|
||||||
|
# Fixed value for the Toomre parameter
|
||||||
|
Q_fixed2 0.0
|
||||||
|
# Additional term for the Toomre parameter
|
||||||
|
Q_boost2 0.0
|
||||||
|
# Past constant Star Formation Rate [Msol/yr]
|
||||||
|
# Negative value recomptes the SFR according to Bouché et al. 2010
|
||||||
|
SFR2 -1.
|
||||||
|
# Minimum age of the stars [Myr]
|
||||||
|
min_age2 50.
|
||||||
|
# Metallicity of the stars [Zsolar]
|
||||||
|
metal2 0.02041
|
||||||
|
# Gaussian fluctuations in the density field
|
||||||
|
dens_fluct2 0
|
||||||
|
# Use the density cut during the Jeans equation integration
|
||||||
|
jeans_mass_cut2 1
|
||||||
|
compute_vel2 1
|
||||||
|
|
||||||
|
###########################
|
||||||
|
# Component 3: Gaseous disk
|
||||||
|
###########################
|
||||||
|
mass_frac3 0.007663631417945721
|
||||||
|
npart3 1000000
|
||||||
|
npart_pot3 2000000
|
||||||
|
model3 3
|
||||||
|
scale_length3 3.43218
|
||||||
|
cut3 30.
|
||||||
|
flatz3 0.1
|
||||||
|
mcmc_step3 0.15
|
||||||
|
type3 0
|
||||||
|
stream_fraction3 1.00
|
||||||
|
metal3 0.02041
|
||||||
|
# Metallicity follows density gradient
|
||||||
|
metal_gradient3 0
|
||||||
|
# Temperature of the gas particles [K]
|
||||||
|
t_init3 1e4
|
||||||
|
# Turbulent velocity dispersion [km/s]
|
||||||
|
turb_sigma3 30.0
|
||||||
|
# Turbulence injection scale [kpc]
|
||||||
|
turb_scale_inj3 0.34
|
||||||
|
# Turbulence dissipation scale [kpc]
|
||||||
|
turb_scale_diss3 0.005
|
||||||
|
# Seed for the turbulent gaussian field
|
||||||
|
turb_seed3 1234
|
||||||
|
# Compute hydrostatic equilibrium
|
||||||
|
hydro_eq3 0
|
||||||
|
# Gaussian step for the MCMC Metropolis-Hasting particle postionning algorithm
|
||||||
|
# within the hydro equilibrium algorithm
|
||||||
|
mcmc_step_hydro3 0.3
|
||||||
|
compute_vel3 1
|
||||||
|
# Polytropic index for the gas (1 = isothermal)
|
||||||
|
gamma_poly3 1.0
|
||||||
|
|
||||||
|
############################
|
||||||
|
# Component 4: Stellar bulge
|
||||||
|
############################
|
||||||
|
mass_frac4 0.0038322616318995422
|
||||||
|
npart4 125000
|
||||||
|
npart_pot4 250000
|
||||||
|
model4 4
|
||||||
|
scale_length4 0.343218
|
||||||
|
flatz4 0.80
|
||||||
|
mcmc_step4 0.3
|
||||||
|
type4 1
|
||||||
|
stream_fraction4 1.
|
||||||
|
min_age4 50.
|
||||||
|
metal4 0.02041
|
||||||
|
Q_lim4 1.5
|
||||||
|
epicycle4 1
|
||||||
|
compute_vel4 1
|
||||||
Reference in New Issue
Block a user