# Parameter file for Initial Condition Generation for a generic pion run #################### ndim 1 # number of dimensions in grid coordinates spherical # coordinate system eqn euler # equations to solve solver 9 # Riemann solver GAMMA 1.666666666666666666666 # Adiabatic index of gas CFL 0.4 # Courant number for simulation ArtificialViscosity 1 # Include artificial viscosity following Falle+(1998) EtaViscosity 0.15 # viscosity parameter OrderOfAccSpace 2 # spatial order of accuracy OrderOfAccTime 2 # temporal order of accuracy # Which ICs to set up ics Uniform # which problem to set up smooth -1 # smooth initial conditions? (<0 -> no) noise -1 # add noise to initial conditions? (<0 -> no) ########################## ## Output File Properties OutputFile Pion_output_z004_r000_M1.000 # Name of snapshot !change OutputPath ./ # path to save data OutputFileType silo # format of snapshot to save OutputFrequency 4096 # frequency of outputs (if OutputCriterion==0) OutputCriterion 1 # Save every N steps (0) or D time intervals (1) OPfreqTime 1e11 # frequency of outputs (s) (if OutputCriterion==1) StartTime 1.000e-20 # initial time of simulation (s) FinishTime 6e+15 # Finishing time of simulation (if not at end of stellar evolution model) # Which physics to use chem_code MPv3 # which microphysics solver to use ntracer 9 # number of passive tracers InitIons LEAVE # How to initialise the tracers Tracer000 X_H # H mass fraction Tracer001 X_He # He mass fraction Tracer002 X_C # C mass fraction Tracer003 X_N # N mass fraction Tracer004 X_O # O mass fraction Tracer005 X_Z # other metal mass fraction (not in dust) Tracer006 X_D # Dust mass fraction Tracer007 H1+ # H+ ionization fraction Tracer008 WIND # passive tracer for wind EP_dynamics 1 # do hydrodynamics EP_raytracing 1 # do raytracing EP_phot_ionisation 1 # do photoionization EP_cooling 1 # do radiative cooling EP_chemistry 1 # integrate microphysics EP_coll_ionisation 1 # include collisional ionization EP_rad_recombination 1 # include radiative recombination EP_update_erg 1 # update internal energy in microphysics EP_MP_timestep_limit 1 # limit hydro timestep by microphysics timescales EP_Min_Temperature 1.0e1 # imposed temperature floor EP_Max_Temperature 1.0e10 # imposed temperature ceiling EP_Hydrogen_MassFrac 0.7114286 # H Mass fraction of ISM !change EP_Helium_MassFrac 0.2845714 # He Mass fraction of ISM !change EP_Metal_MassFrac 0.0040000 # Metallicity of ISM !change ################### ## Grid Parameters NGridX 256 # Number of gridpoints in X-direction on each level NGridY 256 # Number of gridpoints in Y-direction (if >1D) NGridZ 256 # Number of gridpoints in Z-direction (if >2D) Xmin 0.0e18 # negative X grid boundary (level 0) Ymin 0.0e18 # negative Y grid boundary (level 0) Zmin 0.0e18 # negative Z grid boundary (level 0) Xmax 2.048e20 # positive X grid boundary (level 0) Ymax 0.0e18 # positive Y grid boundary (level 0) Zmax 0.0e18 # positive Z grid boundary (level 0) ############################################ # Boundary Conditions BC_XN reflecting # boundary condition at negative X boundary BC_XP outflow # boundary condition at positive X boundary BC_YN outflow # boundary condition at negative Y boundary BC_YP outflow # boundary condition at positive Y boundary BC_ZN outflow # boundary condition at negative Z boundary BC_ZP outflow # boundary condition at positive Z boundary BC_Ninternal 1 # number of internal boundary conditions BC_INTERNAL_000 stellar-wind # internal boundary is stellar wind ######################## # Nested grid grid_nlevels 6 # Number of grid levels for static mesh-refinement !change grid_aspect_ratio_XX 1 # Aspect ratio of grid in each dimension grid_aspect_ratio_YY 1 # Aspect ratio of grid in each dimension grid_aspect_ratio_ZZ 1 # Aspect ratio of grid in each dimension NG_centre_XX 0.0e18 # focus of nested grid refinement in X NG_centre_YY 0.0e18 # focus of nested grid refinement in Y NG_centre_ZZ 0.0e18 # focus of nested grid refinement in Z NG_refine_XX 1 # refine in x-direction? (1==yes) NG_refine_YY 1 # refine in y-direction? (1==yes) NG_refine_ZZ 1 # refine in z-direction? (1==yes) ######################## ############### ## Code Units # Reference systems of units are MKS (i.e. SI), CGS, others? units cgs # if left blank, default to CGS units # One code unit is this number of reference units. rhoval 1 # unused for CGS lenval 1 # unused for CGS velval 1 # unused for CGS magval 1 # unused for CGS ## Ambient Medium Parameters # ############################## # This core part is not used here, can set up an overdense cloud with it. UNIFORM_radial_slope 0.0 # power-law slope of density gradient (0-->uniform) UNIFORM_radial_velocity 0.0 # expansion velocity of ambient medium (cm/s) UNIFORM_core_radius 0.0 # if radial_slope != 0, this is radius of constant-density core (cm) UNIFORM_core_centre_XX 1.1e200 # Location of core, in cm. UNIFORM_core_centre_YY 1.1e200 # Location of core, in cm. UNIFORM_core_centre_ZZ 1.1e200 # Location of core, in cm. # Ambient medium density, pressure, velocity, B-field, etc. # We set up a simulation where the star is at rest in the ISM. UNIFORM_ambRO 2.185e-22 # initial ISM density !change UNIFORM_ambPG 7.086e-13 # initial ISM pressure !change UNIFORM_ambVX 0.0e5 # initial ISM velocity (cm/s) (X) UNIFORM_ambVY 0.0 # initial ISM velocity (cm/s) (Y) UNIFORM_ambVZ 0.0 # initial ISM velocity (cm/s) (Z) UNIFORM_ambBX 0.0 # initial ISM magnetic field (Gauss) (X) UNIFORM_ambBY 0.0 # initial ISM magnetic field (Gauss) (Y) UNIFORM_ambBZ 0.0 # initial ISM magnetic field (Gauss) (Z) UNIFORM_ambTR0 0.7114286 # initial ISM H mass fraction !change; UNIFORM_ambTR1 0.2845714 # initial ISM He mass fraction !change UNIFORM_ambTR2 0.0006883 # initial ISM C mass fraction !change UNIFORM_ambTR3 0.0002016 # initial ISM N mass fraction !change UNIFORM_ambTR4 0.0018721 # initial ISM O mass fraction !change UNIFORM_ambTR5 0.0012380 # initial ISM metal mass fraction remaining UNIFORM_ambTR6 0.0000000 # initial ISM Dust mass fraction (unused here) UNIFORM_ambTR7 1.0e-4 # initial ISM H+ fraction UNIFORM_ambTR8 0.0 # initial ISM Wind tracer mass fraction ######################## # radiation source list ######################## RT_Nsources 1 # number of sources for raytracing # This source is time varying, so the values are read from the RT_EVO_FILE_0 and WIND_0_evofile. RT_position_0_0 0.0e18 # source position X RT_position_0_1 0.0e18 # source position Y RT_position_0_2 0.0e18 # source position Z RT_strength_0 1.0 # luminosity in erg/s for multifrequency source (if constant) RT_Rstar____0 # Radius of star (R_sun) (if constant) RT_Tstar____0 # Stellar effective temperature (K) (if constant) RT_src_type_0 1 # 1=single src, 2=diffuse-radiation. RT_at_infty_0 0 # 0=no, 1=yes. RT_update___0 2 # 1=C2Ray update, 2=Mackey(2012) update. RT_effect___0 3 # 1=UV heating, 2=mono-photoionisation, 3=multifrequency-photoion RT_Tau_src__0 10 # 1=total column, 2=nH*(1-y_i), 3=nH*(y_i), 10=request from microphysics RT_Tau_var__0 7 # opacity variably y_i (this is its index in tracers) RT_Nbins____0 1 # Number of frequency bins (only MPv10) RT_EVO_FILE_0 StellarModel_z004_r000_M1.000.pion # file with time-dependent stellar properties (with path, if needed) !change WIND_NSRC 1 # Number of stellar wind sources WIND_0_pos0 0.0e18 # source position (X) WIND_0_pos1 0.0e18 # source position (Y) WIND_0_pos2 0.0e18 # source position (Z) WIND_0_radius 0.5e17 # radius of wind-injection region !change WIND_0_type 1 # type of wind source (1==evolving) WIND_0_mdot # not needed for evolving wind. WIND_0_vinf # not needed for evolving wind. WIND_0_vrot 0.0 # not needed for evolving wind. WIND_0_temp # not needed for evolving wind. WIND_0_Rstr # not needed for evolving wind. WIND_0_Bsrf 0.0 # not needed for evolving wind. WIND_0_TR0 # H # not needed for evolving wind. WIND_0_TR1 # He # not needed for evolving wind. WIND_0_TR2 # C # not needed for evolving wind. WIND_0_TR3 # N # not needed for evolving wind. WIND_0_TR4 # O # not needed for evolving wind. WIND_0_TR5 # Z (all other metals) # not needed for evolving wind. WIND_0_TR6 # Dust # not needed for evolving wind. WIND_0_TR7 0.0 # H+ fraction # not needed for evolving wind. WIND_0_TR8 1.0 # Wind fraction # not needed for evolving wind. WIND_0_evofile StellarModel_z004_r000_M1.000.pion # file with time-dependent stellar properties (with path, if needed) !change WIND_0_t_offset 0.0 # offset of simulation vs. stellar evolution time (s) WIND_0_t_scalefac 1.0 # Accelerates evolution by this factor (for MS phase). WIND_0_updatefreq 1.0 # no longer used (seconds) WIND_0_enhance_mdot 0 # artificially increase Mdot if !=0 WIND_0_xi -2.0 # only used for non-spherical winds # Reference Vector used by some Riemann solvers to set floor values (not used here) refvec0 1.0e-22 # ~10 particles /cc refvec1 1.0e-12 # T ~ 10^3K refvec2 1.0e6 # v ~ 10 km/s = 1.e7 cm/s refvec3 1.0e6 # v refvec4 1.0e6 # v refvec5 1.0e-6 # B (gauss) refvec6 1.0e-6 refvec7 1.0e-6 refvec8 1.0 refvec9 1.0 refvec10 1.0 refvec11 1.0 refvec12 1.0 refvec13 1.0 refvec14 1.0 refvec15 1.0 refvec16 1.0 refvec17 1.0 refvec18 1.0 refvec19 1.0 refvec20 1.0