%---- Relevant files (filenames and directories) % if restarting from snapshot file: % (1) rename InitCondFile to snapshot file name base, excluding block number and '.hdf5' % (2) change ICFormat to hdf5 = 3 % (3) change submission script to pass restart argument = 2 InitCondFile initial_condition/ic_agora_m11b_ref13.ics % IC file OutputDir output % name (full path or sub-folder name) for outputs %---- File formats (input and output) ICFormat 1 % format of IC file: 1=unformatted (gadget) binary, 3=hdf5, 4=cluster SnapFormat 3 % format of snapshots: 1=unformatted (gadget) binary, 3=hdf5 [Default=3] %---- Output parameters RestartFile restart % base name for restartfiles [Default='restart'] SnapshotFileBase snapshot % base name for snapshot files [Default='snapshot'] OutputListOn 1 % =1 to use list in "OutputListFilename", =0 to specify OutputListFilename snapshot_scale-factors.txt % list of times (in code units) for snaps if OutputListOn=1 NumFilesPerSnapshot 1 % number of sub-files to break each snapshot into [Default=1] NumFilesWrittenInParallel 16 % must be < N_processors & power of 2 [Default=1] %---- Output frequency (use if OutputListOn = 0) TimeOfFirstSnapshot 0.1 % time (code units) of first snapshot TimeBetSnapshot 1.1 % time between snapshot writes (if OutputListOn=0), code units TimeBetStatistics 0.05 % code time between additional statistics (e.g. energy) [Default=infinity] %---- CPU run-time and checkpointing time-limits % 5 day -> 500000, 2 day -> 202000, 1 day -> 100000, 8 hr -> 33000, 2hr -> 8000, 1 hr -> 3900 TimeLimitCPU 202000 % wall-clock time limit in seconds, before self-stop at 85% of this CpuTimeBetRestartFile 3400 % wall-clock time between restartfile writes, in seconds [Default=1 hour] %---- Desired simulation beginning and end times (in code units) for run TimeBegin 0.01 % Code-unit time at the beginning of the simulation TimeMax 1.0 % Code-unit time at the end of the simulation (run will stop) %---- Maximum and minimum timesteps allowed MaxSizeTimestep 0.005 % maximum timestep in code units, set for your problem. in units of delta(ln(a)) MinSizeTimestep 1e-99 % minimum timestep (code units). set this very low, or get the wrong answer [Default=0] %---- Rebuild domains when >this fraction of particles active, set lower to rebuild more frequently TreeDomainUpdateFrequency 0.0005 % 0.0005-0.05, dept on core+particle number %----- Memory allocation MaxMemSize 3500 % sets maximum memory cap per MPI task in MByte (to prevent segfaults/memory overreach) PartAllocFactor 2.5 % memory load allowed (ratio of element number to mean) for better cpu balance BufferSize 100 % size of communications buffer per MPI sub-domain, in MByte %---- Cosmological parameters ComovingIntegrationOn 1 % is the simulation a cosmological integration? (yes=1, no=0) BoxSize 40000. % box side-length in code units [default kpc/h] Omega0 0.272 % Omega_Matter, set = 0 for non-cosmological OmegaLambda 0.728 % Omega_Lambda, set = 0 for non-cosmological OmegaBaryon 0.0455 % Omega_Baryon, set = 0 for non-cosmological HubbleParam 0.702 % Hubble Parameter little 'h', set = 1 for non-cosmological runs %---- System of units. This defines what one 'code unit' is in CGS, so it can be converted for subroutines like cooling UnitLength_in_cm 3.085678e21 % sets to 1.0 kpc / h UnitMass_in_g 1.989e43 % sets to 1e10 M_sun / h UnitVelocity_in_cm_per_s 1.0e5 % sets to 1 km / sec UnitMagneticField_in_gauss 1.0 % sets to 1 gauss GravityConstantInternal 0 % Set the value of 'G' in code units, if desired. calculated by code if =0 [Default=0] %---- (Optional) Initial hydro temperature & temperature floor (in Kelvin) InitGasTemp 0 % Initial temperature: set by IC file if =0 [Default=0] MinGasTemp 10. % Temperature floor. [Default=0] %---- Hydro reconstruction (kernel) parameters DesNumNgb 32 % Cell-reconstruction kernel number: 32-60 standard in 3D, 60-114 for quintic MaxHsml 1e10 % minimum gas kernel length (some very large value to prevent errors) [Default=infinity] MinGasHsmlFractional 1.0 % minimum kernel length relative to gas force softening (<= 1) [Default=0] %----- adaptive grav softening or other kernel-interactions for non-gas elements AGS_DesNumNgb 32 % Neighbor number for calculating interaction kernel for non-gas neighbor maps. %---- Gravitational softening lengths, in code units %----- Softening lengths per element type. If ADAPTIVE_GRAVSOFT is set, these %-------- are the minimum softening allowed for each type ------- %-------- (units are co-moving for cosmological integrations) Softening_Type2 2 % low-resolution dark matter, 'SofteningDisk' Softening_Type1 0.18 % high-resolution dark matter, 'SofteningHalo' Softening_Type4 0.018 % stars spawned from gas, 'SofteningStars' Softening_Type0 0.007 % gas, 'SofteningGas' Softening_Type5 1 % black holes/sinks, 'SofteningBndry' Softening_Type3 1 % grains/PIC particles, 'SofteningBulge' %--- If these are set in cosmological runs, SofteningX switches from comoving to physical %------- units when the comoving value exceeds the choice here %------- (these are ignored/do not need to be set, and *only* the above are used, for non-cosmo runs) Softening_Type2_MaxPhysLimit 0.2 % 'SofteningDiskMaxPhys' Softening_Type1_MaxPhysLimit 0.018 % 'SofteningHaloMaxPhys' Softening_Type4_MaxPhysLimit 0.0018 % 'SofteningStarsMaxPhys' Softening_Type0_MaxPhysLimit 0.0007 % 'SofteningGasMaxPhys' Softening_Type5_MaxPhysLimit 1 % 'SofteningBndryMaxPhys' Softening_Type3_MaxPhysLimit 1 % 'SofteningBulgeMaxPhys' %------------------------------------------------------------ %------------------ Additional Fluid Physics ---------------- %------------------------------------------------------------ %---- Magneto-Hydrodynamics Parameters (MAGNETIC on) %----- Initial B-Field Strengths (if MHD_B_SET_IN_PARAMS on, otherwise read from IC file) BiniX 1e-10 % initial B_x, in code units BiniY 1e-10 % initial B_y, in code units BiniZ 1e-10 % initial B_z, in code units %---- Thermal Conduction (CONDUCTION on) %----- set coefficient kappa [code units] or, if CONDUCTION_SPITZER on, multiplies value ConductionCoeff 1.0 % set/multiply conduction coefficient %---- Navier-Stokes Viscosity (VISCOSITY on) %--- set coefficients eta,zeta [code units] or, if VISCOSITY_BRAGINSKII on, multiplies value ShearViscosityCoeff 1.0 % set/multiply shear viscosity coefficient BulkViscosityCoeff 1.0 % set/multiply bulk viscosity coefficient %---- Turbulent Diffusion Master Switch (TURB_DIFFUSION on) TurbDiffusionCoefficient 1.0 % Normalizes diffusion rates relative to Smagorinsky-Lilly theory [best calibration] (~0.5-2) (~0.5-2) %---- Cosmic Ray + Gas Fluids (COSMIC_RAYS on) CosmicRayDiffusionCoeff 1.0 % multiplies anisotropic diffusion/streaming coefficients %------------------------------------------------------------------------- %------------------ Star, Black Hole, and Galaxy Formation --------------- %------------------------------------------------------------------------- %---- Star Formation parameters (GALSF on) CritPhysDensity 1000. % critical physical density for star formation (cm^(-3)) SfEffPerFreeFall 1.0 % SFR/(Mgas/tfreefall) for gas which meets SF criteria %-------------- FIRE (PFH) explicit star formation & feedback model (FIRE on) %--- initial metallicity of gas & stars in simulation InitMetallicity 1e-4 % initial gas+stellar metallicity (in solar) InitStellarAge 1e-3 % initial mean age (in Gyr; for stars in sim ICs) %--- local radiation-pressure driven winds (GALSF_FB_FIRE_RT_LOCALRP) WindMomentumLoading 1.0 % fraction of photon momentum to couple %--- SneII Heating Model (GALSF_FB_MECHANICAL) SNeIIEnergyFrac 1.0 % fraction of mechanical energy to couple %--- HII region photo-heating model (GALSF_FB_FIRE_RT_HIIHEATING) HIIRegion_fLum_Coupled 1.0 % fraction of ionizing photons allowed to see gas %--- long-range radiation pressure acceleration (GALSF_FB_FIRE_RT_LONGRANGE) PhotonMomentum_Coupled_Fraction 1.0 % fraction of L to allow incident PhotonMomentum_fUV 0.01 % incident SED f(L) in UV (minimum scattering) PhotonMomentum_fOPT 0.01 % incident SED f(L) in optical/near-IR %--- gas return/recycling GasReturnFraction 1.0 % fraction of gas mass returned (relative to ssp) GasReturnEnergy 1.0 % fraction of returned gas energy+momentum (relative to ssp) %--- cosmic rays (COSMIC_RAYS) CosmicRay_SNeFraction 0.1 % fraction of SNe ejecta kinetic energy into cosmic rays (~10%) %-------------- Black Hole accretion & formation (BLACK_HOLES on) %--- formation/seeding SeedBlackHoleMass 0.702e-6 % initial mass (on-the-fly or single galaxy) SeedAlphaDiskMass 0.0 % initial mass in the alpha disk (BH_ALPHADISK_ACCRETION) SeedBlackHoleMinRedshift 2.0 % minimum redshift where new BH particles are seeded (lower-z ceases seeding) SeedBlackHoleMassSigma 0 % lognormal standard deviation (in dex) in initial BH seed masses SeedBlackHolePerUnitMass 1.e-4 % stellar mass (code unit) per BH seed (BH_SEED_FROM_LOCALGAS) %----- (specific options for on-the-fly friends-of-friends based BH seeding: FOF on) MinFoFMassForNewSeed 0.702e-3 % minimum mass of FOF group (stars or DM) to get seed, in code units TimeBetOnTheFlyFoF 1.01 % time (in code units, e.g. scale-factor) between on-the-fly FOF searches %--- accretion BlackHoleAccretionFactor 0.5 % multiplier for mdot (relative to model) BlackHoleEddingtonFactor 10 % fraction of eddington to cap (can be >1) BlackHoleNgbFactor 8 % multiplier for kernel neighbors for BH BlackHoleMaxAccretionRadius 0.1 % max radius for BH neighbor search/accretion (code units) BlackHoleRadiativeEfficiency 0.1 % radiative efficiency (for accretion and feedback) %--- feedback BlackHoleFeedbackFactor 1.0 % generic feedback strength multiplier BH_FluxMomentumFactor 1.0 % multiply radiation pressure (BH_PHOTONMOMENTUM) BAL_f_accretion 1.0 % fraction of gas swallowed by BH (BH_WIND options) BAL_v_outflow 30000 % velocity (km/s) of BAL outflow (BH_WIND options) BAL_internal_temperature 1e4 % internal temperature (K) of BAL outflow (BH_WIND_SPAWN) BAL_wind_particle_mass 1e-10 % mass of 'virtual wind particles' in code units (BH_WIND_SPAWN) BH_CosmicRay_Injection_Efficiency 0.1 % fraction of accretion energy into cosmic rays (BH_COSMIC_RAYS) %------------------------------------------------------------------------- %------------------ GRACKLE cooling module ----------------- %------------------------------------------------------------------------- %-------------- Grackle UVB file (GRACKLE on) GrackleDataFile CloudyData_UVB=HM2012.h5 % filename for UV background info for GRACKLE %------------------------------------------------------------------------------------------------- %------------------ Non-Standard Dark Matter, Dark Energy, Gravity, or Expansion ----------------- %------------------------------------------------------------------------------------------------- %-------------- Parameters for non-standard or time-dependent Gravity/Dark Energy/Expansion (GR_TABULATED_COSMOLOGY on) DarkEnergyConstantW -1 % time-independent DE parameter w, used only if no table TabulatedCosmologyFile CosmoTbl % table with cosmological parameters %-------------- Parameters for self-interacting dark matter (DM_SIDM on) DM_InteractionCrossSection 1 % cross-section per unit mass in cm^2/g (normalized at 1 km/s, if vel-dependent) DM_InteractionVelocityScale 0 % velocity [code units]: if>0, cross-section scales ~1/(1+(v/v0)^4) above this "v0" DM_DissipationFactor 0 % 0=elastic, 1=pure dissipative (fractional dissipation of kinetic energy in event) DM_KickPerCollision 0 % velocity 'kick' [code units] per collision (this^2=specific energy/mass released) %-------------- Parameters for scalar-field dark matter (DM_FUZZY on) FuzzyDM_Mass_in_eV 1e-22 % boson mass in eV for SFDM/Fuzzy DM module