hlasco
146 lines
4.7 KiB
Plaintext
146 lines
4.7 KiB
Plaintext
&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
|
|
/
|