diff --git a/.gitignore b/.gitignore index cf324c2..3f31d5c 100644 --- a/.gitignore +++ b/.gitignore @@ -1,2 +1,3 @@ *.g2 *.log +*.*Store diff --git a/low/agora/config_agora_sigma30.dice b/high/agora/config_agora_sigma00.dice similarity index 87% rename from low/agora/config_agora_sigma30.dice rename to high/agora/config_agora_sigma00.dice index 32fdd4f..b771d82 100644 --- a/low/agora/config_agora_sigma30.dice +++ b/high/agora/config_agora_sigma00.dice @@ -3,7 +3,7 @@ # 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. +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 @@ -15,7 +15,7 @@ MaxCompNumber 64 OutputRz 0 # Name of the final IC file # The file extension wil be added by DICE -Filename out/agora_sigma30 +Filename out/agora_sigma00 # File format of the IC file [Gadget1,Gadget2] ICformat Gadget2 # Cosmological parameters for Virial quantities @@ -23,11 +23,11 @@ ICformat Gadget2 # Redshift Redshift 1.0 # Hubble constant -H0 71.0 +H0 73.0 # Dark energy density parameter -OmegaL 0.70 +OmegaL 0.6911 # Matter density parameter -OmegaM 0.30 +OmegaM 0.3089 # Spatial curvature density paramter OmegaK 0.00 # Normalize the sum of mass fractions to one diff --git a/low/agora/config_agora_sigma10.dice b/high/agora/config_agora_sigma20.dice similarity index 92% rename from low/agora/config_agora_sigma10.dice rename to high/agora/config_agora_sigma20.dice index 2b46073..3f5e503 100644 --- a/low/agora/config_agora_sigma10.dice +++ b/high/agora/config_agora_sigma20.dice @@ -3,7 +3,7 @@ # 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. +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 @@ -15,7 +15,7 @@ MaxCompNumber 64 OutputRz 0 # Name of the final IC file # The file extension wil be added by DICE -Filename out/agora_sigma10 +Filename out/agora_sigma20 # File format of the IC file [Gadget1,Gadget2] ICformat Gadget2 # Cosmological parameters for Virial quantities diff --git a/med/agora/nml/agora_sigma20.nml b/high/agora/nml/base.nml similarity index 92% rename from med/agora/nml/agora_sigma20.nml rename to high/agora/nml/base.nml index c35062c..1ffdc0f 100644 --- a/med/agora/nml/agora_sigma20.nml +++ b/high/agora/nml/base.nml @@ -11,7 +11,7 @@ rt=.false. / &OUTPUT_PARAMS -delta_tout=0.25 +delta_tout=1 tend=50.0 / @@ -22,7 +22,7 @@ tendmov=50. movie=.true. nw_frame=1080 nh_frame=1080 -levelmax_frame=14 +levelmax_frame=16 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 @@ -36,7 +36,7 @@ deltaz_frame=30,0.0 30,0.0 &INIT_PARAMS filetype='dice' -initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out' +initfile(1)='/users/hlascomb/dice_ics/high/agora/out' / &UNITS_PARAMS @@ -46,9 +46,9 @@ units_length = 0.308567758128200E+22 ! kpc / &AMR_PARAMS -levelmin=7 -levelmax=14 -ngridmax=3000000 +levelmin=8 +levelmax=16 +ngridmax=1000000 npartmax=5000000 boxlen=320.0 / @@ -83,8 +83,8 @@ 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. +m_star=1.0 ! Star particle mass +sf_log_properties=.true. / &FEEDBACK_PARAMS @@ -107,14 +107,14 @@ riemann='hllc' / &REFINE_PARAMS -mass_sph=1e-6 +mass_sph=1e-7 m_refine=10*8. interpol_type=0 interpol_var=1 / &DICE_PARAMS -ic_file='agora_sigma20.g2' ! Name of the initial conditions file +ic_file='agora_sigma20.g2' 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' diff --git a/low/agora/params/agora_sigma10.params b/high/agora/params/agora_sigma00.params similarity index 96% rename from low/agora/params/agora_sigma10.params rename to high/agora/params/agora_sigma00.params index b65166c..aa476d1 100644 --- a/low/agora/params/agora_sigma10.params +++ b/high/agora/params/agora_sigma00.params @@ -49,11 +49,11 @@ hydro_eq_niter 3 # Fraction of the virial mass in the component 1 mass_frac1 0.9578424465115449 # Number of particles for the component 1 -npart1 100000 +npart1 1000000 # Target mass of individual particles [Msol] #part_mass1 2e6 # Number of particles for the potential computation -npart_pot1 200000 +npart_pot1 2000000 # Target mass of individual particles for potential computation [Msol] #part_mass_pot1 1e6 # Component 1 density model @@ -84,7 +84,7 @@ rcore1 0. # 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 +softening1 0.3 # Component 1 radial density cut [kpc] cut1 0.0 # Component 1 thickness parameters @@ -132,8 +132,8 @@ accept_max1 0.99 # Component 2: Thin stellar disk ################################ mass_frac2 0.03066166043860979 -npart2 100000 -npart_pot2 200000 +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) @@ -183,9 +183,9 @@ metal_gradient3 0 # Temperature of the gas particles [K] t_init3 1e4 # Turbulent velocity dispersion [km/s] -turb_sigma3 10.0 +turb_sigma3 0.0 # Turbulence injection scale [kpc] -turb_scale_inj3 0.34 +turb_scale_inj3 1.0 # Turbulence dissipation scale [kpc] turb_scale_diss3 0.01 # Seed for the turbulent gaussian field @@ -203,8 +203,8 @@ gamma_poly3 1.0 # Component 4: Stellar bulge ############################ mass_frac4 0.0038322616318995422 -npart4 12500 -npart_pot4 25000 +npart4 125000 +npart_pot4 250000 model4 4 scale_length4 0.343218 flatz4 0.80 diff --git a/med/agora/params/agora_sigma10.params b/high/agora/params/agora_sigma10.params similarity index 100% rename from med/agora/params/agora_sigma10.params rename to high/agora/params/agora_sigma10.params diff --git a/low/agora/params/agora_sigma30.params b/high/agora/params/agora_sigma20.params similarity index 96% rename from low/agora/params/agora_sigma30.params rename to high/agora/params/agora_sigma20.params index 345fa92..0e084b1 100644 --- a/low/agora/params/agora_sigma30.params +++ b/high/agora/params/agora_sigma20.params @@ -49,11 +49,11 @@ hydro_eq_niter 3 # Fraction of the virial mass in the component 1 mass_frac1 0.9578424465115449 # Number of particles for the component 1 -npart1 100000 +npart1 1000000 # Target mass of individual particles [Msol] #part_mass1 2e6 # Number of particles for the potential computation -npart_pot1 200000 +npart_pot1 2000000 # Target mass of individual particles for potential computation [Msol] #part_mass_pot1 1e6 # Component 1 density model @@ -84,7 +84,7 @@ rcore1 0. # 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 +softening1 0.3 # Component 1 radial density cut [kpc] cut1 0.0 # Component 1 thickness parameters @@ -132,8 +132,8 @@ accept_max1 0.99 # Component 2: Thin stellar disk ################################ mass_frac2 0.03066166043860979 -npart2 100000 -npart_pot2 200000 +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) @@ -183,11 +183,11 @@ metal_gradient3 0 # Temperature of the gas particles [K] t_init3 1e4 # Turbulent velocity dispersion [km/s] -turb_sigma3 30.0 +turb_sigma3 20.0 # Turbulence injection scale [kpc] turb_scale_inj3 0.34 # Turbulence dissipation scale [kpc] -turb_scale_diss3 0.01 +turb_scale_diss3 0.005 # Seed for the turbulent gaussian field turb_seed3 1234 # Compute hydrostatic equilibrium @@ -203,8 +203,8 @@ gamma_poly3 1.0 # Component 4: Stellar bulge ############################ mass_frac4 0.0038322616318995422 -npart4 12500 -npart_pot4 25000 +npart4 125000 +npart_pot4 250000 model4 4 scale_length4 0.343218 flatz4 0.80 diff --git a/med/agora/params/agora_sigma30.params b/high/agora/params/agora_sigma30.params similarity index 100% rename from med/agora/params/agora_sigma30.params rename to high/agora/params/agora_sigma30.params diff --git a/low/agora/nml/agora_sigma00.nml b/low/agora/nml/agora_sigma00.nml deleted file mode 100644 index c851ac8..0000000 --- a/low/agora/nml/agora_sigma00.nml +++ /dev/null @@ -1,145 +0,0 @@ -&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 deleted file mode 100644 index a9bdf93..0000000 --- a/low/agora/nml/agora_sigma10.nml +++ /dev/null @@ -1,145 +0,0 @@ -&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_sigma30.nml b/low/agora/nml/agora_sigma30.nml deleted file mode 100644 index c6d1f29..0000000 --- a/low/agora/nml/agora_sigma30.nml +++ /dev/null @@ -1,145 +0,0 @@ -&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/nml/agora_sigma20.nml b/low/agora/nml/base.nml similarity index 91% rename from low/agora/nml/agora_sigma20.nml rename to low/agora/nml/base.nml index e7df59c..b4f25ac 100644 --- a/low/agora/nml/agora_sigma20.nml +++ b/low/agora/nml/base.nml @@ -1,4 +1,4 @@ -&RUN_PARAMS +&RUN_PARAMS cosmo=.false. pic=.true. poisson=.true. @@ -11,7 +11,7 @@ rt=.false. / &OUTPUT_PARAMS -delta_tout=0.5 +delta_tout=1 tend=50.0 / @@ -19,7 +19,7 @@ tend=50.0 imovout=2000 imov=0 tendmov=50. -movie=.false. +movie=.true. nw_frame=1080 nh_frame=1080 levelmax_frame=12 @@ -36,7 +36,7 @@ deltaz_frame=30,0.0 30,0.0 &INIT_PARAMS filetype='dice' -initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out' +initfile(1)='/users/hlascomb/dice_ics/low/agora/out' / &UNITS_PARAMS @@ -48,12 +48,12 @@ units_length = 0.308567758128200E+22 ! kpc &AMR_PARAMS levelmin=7 levelmax=12 -ngridmax=3000000 +ngridmax=300000 npartmax=5000000 boxlen=320.0 / -&POISSON_PARAMS +&POISSON_PARAMS epsilon=1e-4 / @@ -78,13 +78,13 @@ self_shielding=.true. &SF_PARAMS sf_virial=.true. ! Enable turbulent SF -sf_model=5 ! See mom2/star_formation.f90 +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. +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. +m_star=1.0 ! Star particle mass +sf_log_properties=.true. ! do not set this to .true. / &FEEDBACK_PARAMS @@ -101,7 +101,7 @@ 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 +beta_fix=0.5 ! Controls threshold at which energy is ! truncated riemann='hllc' / @@ -114,12 +114,12 @@ interpol_var=1 / &DICE_PARAMS -ic_file='agora_sigma20.g2' ! Name of the initial conditions file +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_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 diff --git a/med/agora/config_agora_sigma10.dice b/med/agora/config_agora_sigma10.dice deleted file mode 100644 index 2b46073..0000000 --- a/med/agora/config_agora_sigma10.dice +++ /dev/null @@ -1,39 +0,0 @@ -# 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_sigma30.dice b/med/agora/config_agora_sigma30.dice deleted file mode 100644 index 32fdd4f..0000000 --- a/med/agora/config_agora_sigma30.dice +++ /dev/null @@ -1,39 +0,0 @@ -# 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 deleted file mode 100644 index 058dcb4..0000000 --- a/med/agora/nml/agora_sigma00.nml +++ /dev/null @@ -1,145 +0,0 @@ -&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 deleted file mode 100644 index 47b10d4..0000000 --- a/med/agora/nml/agora_sigma10.nml +++ /dev/null @@ -1,145 +0,0 @@ -&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_sigma30.nml b/med/agora/nml/base.nml similarity index 91% rename from med/agora/nml/agora_sigma30.nml rename to med/agora/nml/base.nml index 72834d6..c8938ff 100644 --- a/med/agora/nml/agora_sigma30.nml +++ b/med/agora/nml/base.nml @@ -1,4 +1,4 @@ -&RUN_PARAMS +&RUN_PARAMS cosmo=.false. pic=.true. poisson=.true. @@ -11,7 +11,7 @@ rt=.false. / &OUTPUT_PARAMS -delta_tout=0.25 +delta_tout=1 tend=50.0 / @@ -36,7 +36,7 @@ deltaz_frame=30,0.0 30,0.0 &INIT_PARAMS filetype='dice' -initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out' +initfile(1)='/users/hlascomb/dice_ics/med/agora/out' / &UNITS_PARAMS @@ -48,12 +48,12 @@ units_length = 0.308567758128200E+22 ! kpc &AMR_PARAMS levelmin=7 levelmax=14 -ngridmax=3000000 -npartmax=5000000 +ngridmax=500000 +npartmax=3000000 boxlen=320.0 / -&POISSON_PARAMS +&POISSON_PARAMS epsilon=1e-4 / @@ -78,13 +78,13 @@ self_shielding=.true. &SF_PARAMS sf_virial=.true. ! Enable turbulent SF -sf_model=5 ! See mom2/star_formation.f90 +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. +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. +m_star=1.0 ! Star particle mass +sf_log_properties=.true. ! do not set this to .true. / &FEEDBACK_PARAMS @@ -101,7 +101,7 @@ 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 +beta_fix=0.5 ! Controls threshold at which energy is ! truncated riemann='hllc' / @@ -114,12 +114,12 @@ interpol_var=1 / &DICE_PARAMS -ic_file='agora_sigma30.g2' ! Name of the initial conditions file +ic_file='agora_sigma20.g2' 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_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