diff --git a/low/agora/config_agora_sigma00.dice b/low/agora/config_agora_sigma00.dice new file mode 100644 index 0000000..b771d82 --- /dev/null +++ b/low/agora/config_agora_sigma00.dice @@ -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 + diff --git a/low/agora/config_agora_sigma10.dice b/low/agora/config_agora_sigma10.dice new file mode 100644 index 0000000..2b46073 --- /dev/null +++ b/low/agora/config_agora_sigma10.dice @@ -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 + diff --git a/low/agora/config_agora_sigma20.dice b/low/agora/config_agora_sigma20.dice new file mode 100644 index 0000000..3f5e503 --- /dev/null +++ b/low/agora/config_agora_sigma20.dice @@ -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 + diff --git a/low/agora/config_agora_sigma30.dice b/low/agora/config_agora_sigma30.dice new file mode 100644 index 0000000..32fdd4f --- /dev/null +++ b/low/agora/config_agora_sigma30.dice @@ -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 + diff --git a/low/agora/nml/agora_sigma00.nml b/low/agora/nml/agora_sigma00.nml new file mode 100644 index 0000000..c851ac8 --- /dev/null +++ b/low/agora/nml/agora_sigma00.nml @@ -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 +/ diff --git a/low/agora/nml/agora_sigma10.nml b/low/agora/nml/agora_sigma10.nml new file mode 100644 index 0000000..a9bdf93 --- /dev/null +++ b/low/agora/nml/agora_sigma10.nml @@ -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 +/ diff --git a/low/agora/nml/agora_sigma20.nml b/low/agora/nml/agora_sigma20.nml new file mode 100644 index 0000000..e7df59c --- /dev/null +++ b/low/agora/nml/agora_sigma20.nml @@ -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 +/ diff --git a/low/agora/nml/agora_sigma30.nml b/low/agora/nml/agora_sigma30.nml new file mode 100644 index 0000000..c6d1f29 --- /dev/null +++ b/low/agora/nml/agora_sigma30.nml @@ -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 +/ diff --git a/low/agora/params/agora_sigma00.params b/low/agora/params/agora_sigma00.params new file mode 100644 index 0000000..55851e4 --- /dev/null +++ b/low/agora/params/agora_sigma00.params @@ -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 diff --git a/low/agora/params/agora_sigma10.params b/low/agora/params/agora_sigma10.params new file mode 100644 index 0000000..b65166c --- /dev/null +++ b/low/agora/params/agora_sigma10.params @@ -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 diff --git a/low/agora/params/agora_sigma20.params b/low/agora/params/agora_sigma20.params new file mode 100644 index 0000000..84f94a2 --- /dev/null +++ b/low/agora/params/agora_sigma20.params @@ -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 diff --git a/low/agora/params/agora_sigma30.params b/low/agora/params/agora_sigma30.params new file mode 100644 index 0000000..345fa92 --- /dev/null +++ b/low/agora/params/agora_sigma30.params @@ -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 diff --git a/low/pbox/config_pbox.dice b/low/pbox/config_pbox.dice new file mode 100644 index 0000000..8b47e12 --- /dev/null +++ b/low/pbox/config_pbox.dice @@ -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 + diff --git a/low/pbox/nml/agora_sigma00.nml b/low/pbox/nml/agora_sigma00.nml new file mode 100644 index 0000000..c851ac8 --- /dev/null +++ b/low/pbox/nml/agora_sigma00.nml @@ -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 +/ diff --git a/low/pbox/nml/agora_sigma10.nml b/low/pbox/nml/agora_sigma10.nml new file mode 100644 index 0000000..a9bdf93 --- /dev/null +++ b/low/pbox/nml/agora_sigma10.nml @@ -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 +/ diff --git a/low/pbox/nml/agora_sigma30.nml b/low/pbox/nml/agora_sigma30.nml new file mode 100644 index 0000000..c6d1f29 --- /dev/null +++ b/low/pbox/nml/agora_sigma30.nml @@ -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 +/ diff --git a/low/pbox/nml/base.nml b/low/pbox/nml/base.nml new file mode 100644 index 0000000..d938716 --- /dev/null +++ b/low/pbox/nml/base.nml @@ -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 +/ diff --git a/low/pbox/params/pbox.params b/low/pbox/params/pbox.params new file mode 100644 index 0000000..b5a61ea --- /dev/null +++ b/low/pbox/params/pbox.params @@ -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 diff --git a/med/agora/config_agora_sigma00.dice b/med/agora/config_agora_sigma00.dice new file mode 100644 index 0000000..b771d82 --- /dev/null +++ b/med/agora/config_agora_sigma00.dice @@ -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 + diff --git a/med/agora/config_agora_sigma10.dice b/med/agora/config_agora_sigma10.dice new file mode 100644 index 0000000..2b46073 --- /dev/null +++ b/med/agora/config_agora_sigma10.dice @@ -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 + diff --git a/med/agora/config_agora_sigma20.dice b/med/agora/config_agora_sigma20.dice new file mode 100644 index 0000000..3f5e503 --- /dev/null +++ b/med/agora/config_agora_sigma20.dice @@ -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 + diff --git a/med/agora/config_agora_sigma30.dice b/med/agora/config_agora_sigma30.dice new file mode 100644 index 0000000..32fdd4f --- /dev/null +++ b/med/agora/config_agora_sigma30.dice @@ -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 + diff --git a/med/agora/nml/agora_sigma00.nml b/med/agora/nml/agora_sigma00.nml new file mode 100644 index 0000000..058dcb4 --- /dev/null +++ b/med/agora/nml/agora_sigma00.nml @@ -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 +/ diff --git a/med/agora/nml/agora_sigma10.nml b/med/agora/nml/agora_sigma10.nml new file mode 100644 index 0000000..47b10d4 --- /dev/null +++ b/med/agora/nml/agora_sigma10.nml @@ -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 +/ diff --git a/med/agora/nml/agora_sigma20.nml b/med/agora/nml/agora_sigma20.nml new file mode 100644 index 0000000..c35062c --- /dev/null +++ b/med/agora/nml/agora_sigma20.nml @@ -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 +/ diff --git a/med/agora/nml/agora_sigma30.nml b/med/agora/nml/agora_sigma30.nml new file mode 100644 index 0000000..72834d6 --- /dev/null +++ b/med/agora/nml/agora_sigma30.nml @@ -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 +/ diff --git a/med/agora/params/agora_sigma00.params b/med/agora/params/agora_sigma00.params new file mode 100644 index 0000000..aa476d1 --- /dev/null +++ b/med/agora/params/agora_sigma00.params @@ -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 diff --git a/med/agora/params/agora_sigma10.params b/med/agora/params/agora_sigma10.params new file mode 100644 index 0000000..dadf2ef --- /dev/null +++ b/med/agora/params/agora_sigma10.params @@ -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 diff --git a/med/agora/params/agora_sigma20.params b/med/agora/params/agora_sigma20.params new file mode 100644 index 0000000..0e084b1 --- /dev/null +++ b/med/agora/params/agora_sigma20.params @@ -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 diff --git a/med/agora/params/agora_sigma30.params b/med/agora/params/agora_sigma30.params new file mode 100644 index 0000000..2f804a1 --- /dev/null +++ b/med/agora/params/agora_sigma30.params @@ -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