%------------------------------------------------------------------------- %------------------------------------------------------------------------- %---------- Generic Parameters Required in Any Simulation ---------------- %------------------------------------------------------------------------- %------------------------------------------------------------------------- %---- Relevant files InitCondFile output/snapdir_179/snapshot_179 OutputDir output % output directory %---- File formats ICFormat 3 % 1 = binary, 3 = hdf5: set appropriately SnapFormat 3 % 3 = hdf5: ALWAYS USE THIS for FIRE simulations %---- Output parameters (snapshots, etc) RestartFile restart % prefix for restartfiles, always use 'restart' for FIRE SnapshotFileBase snapshot % prefix for snapshots, always use 'snapshot' for FIRE OutputListOn 1 % = 1 to use list in "OutputListFilename": for FIRE zoom-in simulations [see below] OutputListFilename snapshot_scale-factors.txt % list of scale factors for snapshots NumFilesPerSnapshot 8 % generally set this =1 for snapshots <~ 1 GB, but ~4-8 for larger snapshots NumFilesWrittenInParallel 512 % must be <= N_processors & power of 2; set to >1 for efficiency %---- Output frequency (use if OutputListOn = 0) TimeOfFirstSnapshot 0.0 % Time for snapshot 0 (if OutputListOn=1, this parameter is irrelevant) TimeBetSnapshot 1.0 % Time (code units) between snapshots TimeBetStatistics 0.05 % Time (code units) between internal statistics of certain types: set to large fraction of total time %---- CPU-time limits TimeLimitCPU 432000 % in seconds. stops at 85%. CpuTimeBetRestartFile 7200 % in seconds ResubmitOn 0 ResubmitCommand qsub resub_script.pbs %----- Memory allocation MaxMemSize 6250 % in MByte. maximum memory per MPI task (increase if multi-threading) PartAllocFactor 2.5 % memory load (x average particles per MPI task) allowed for better cpu balance BufferSize 25 % in MByte %---- Characteristics of run TimeBegin 0.00793650793 % beginning time (scale factor) of simulation: this is 1/(1+z_initial) for cosmological runs, =0 for non-cosmological TimeMax 1.0 % ending time (scale factor) of simulation (=1/(1+z_final) for cosmological runs) %---- Cosmological parameters ComovingIntegrationOn 1 % is it cosmological? (yes = 1, no = 0) BoxSize 50000. % in code units [default in kpc/h] Omega0 0.266 % = 0 for non-cosmological OmegaLambda 0.734 % = 0 for non-cosmological OmegaBaryon 0.0455 % = 0 for non-cosmological HubbleParam 0.71 % little 'h'; = 1 for non-cosmological runs %---- Accuracy of time integration MaxSizeTimestep 0.01 % for cosmological, in units of delta(ln(a)) MinSizeTimestep 1.0e-11 %---- Tree algorithm, force accuracy, domain update frequency TreeDomainUpdateFrequency 0.0025 % 0.0005 - 0.05, depending on core+particle number %---- System of units -- ALWAYS USE THESE UNITS for FIRE simulations, our analysis software and ICs assume them 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 % calculated by code if = 0 %---- Initial temperature & temperature floor InitGasTemp 0 % set by IC file if = 0. if IC file value != 0, use that instead MinGasTemp 10. % =10 K in FIRE simulations, or else you will get errors from the cooling tables %---- Density/volume estimation (kernel) DesNumNgb 32 % 32 for standard kernel, 60 - 114 for quintic MaxHsml 1.0e10 % maximum gas kernel length (some large value to prevent errors) MinGasHsmlFractional 1.0 % minimum kernel length relative to gas force softening (<= 1) %---- Gravitational softening lengths %----- Softening lengths per particle type. If ADAPTIVE_GRAVSOFT is set, these %-------- are the minimum softening allowed for each type ------- %-------- (units are co-moving for cosmological integrations) SofteningHalo 0.22 % dark matter/collisionless particles (type = 1) SofteningDisk 1.0 % collisionless particles (type = 2) SofteningGas 0.0027 % gas (type = 0) SofteningStars 0.019 % stars spawned from gas (type = 4) SofteningBulge 2.0 % collisionless particles (type = 3) SofteningBndry 2.0 % black holes (if active), or collisionless (type = 5) %---- if these are set in cosmological runs, SofteningX switches from comoving to physical %------- units when the comoving value exceeds the choice here SofteningHaloMaxPhys 0.022 SofteningDiskMaxPhys 0.1 SofteningGasMaxPhys 0.00027 SofteningStarsMaxPhys 0.0019 % ~7 * gas softening SofteningBulgeMaxPhys 0.2 SofteningBndryMaxPhys 0.2 %----- parameters for adaptive gravitational softening AGS_DesNumNgb 32 % neighbor number for calculating adaptive gravsoft %------------------------------------------------------------------------- %------------------------------------------------------------------------- %---------- FIRE-Specific Physics Parameters ----------------------------- %------------------------------------------------------------------------- %------------------------------------------------------------------------- %---- Star Formation parameters (GALSF on) CritPhysDensity 1000. % critical physical density for star formation (cm^(-3)): our default is 1000 in FIRE cosmological sims SfEffPerFreeFall 1.0 % SFR/(Mgas/tfreefall) for gas which meets SF criteria: always set=1 in FIRE simulations %-------------- FIRE (PFH) explicit star formation & feedback model %----- DO NOT change these values from those below in FIRE runs. %----- unless you want to turn on/off feedback or artificially make it stronger/weaker %----- for some reason. Then use these (rather than the compile-time options) to, say, %----- turn off SNe by setting "SNeIIEnergyFrac" = 0.0 %----- The InitMetallicity = 0.0001 sets cosmological simulations to initial %----- metallicities of 10^-4 solar -- this is a floor to prevent bad behavior in %----- some FIRE lookup tables and fitting functions. For non-cosmological simulations %----- it must be set to the correct metallicity of the initial conditions. %----- Likewise InitStellarAge (irrelevant in cosmological simulations) must be set for %----- non-cosmological simulations to the correct mean age of the stars which are already present in the simulation ICs %--- initial metallicity of gas & stars in simulation InitMetallicity 0.0001 % in solar InitStellarAge 0.001 % mean age (in Gyr) %--- local radiation-pressure driven winds (GALSF_FB_RPWIND_LOCAL) WindMomentumLoading 1.0 % boost momentum-loading (relative to ssp) %--- SneII Heating Model (GALSF_FB_SNE_HEATING) SNeIIEnergyFrac 1.0 % boost fraction of mechanical energy to couple %--- HII region photo-heating model (GALSF_FB_HII_HEATING) HIIRegion_fLum_Coupled 1.0 % boost fraction of ionizing photons %--- long-range radiation pressure acceleration (GALSF_FB_RT_PHOTONMOMENTUM) PhotonMomentum_Coupled_Fraction 1.0 % boost fraction of L to allow incident PhotonMomentum_fUV 0.0 % incident SED f(L) in UV (minimum scattering) PhotonMomentum_fOPT 0.0 % incident SED f(L) in optical/near-IR %--- gas return/recycling (GALSF_FB_GASRETURN) GasReturnFraction 1.0 % boost fraction of gas returns (relative to ssp) GasReturnEnergy 1.0 % boost temperature of gas returns (relative to ssp) %------------------------------------------------------------------------- %------------------------------------------------------------------------- %---------- NON-FIRE Physics Modules ------------------------------------- %------------------------------------------------------------------------- %------------------------------------------------------------------------- %---- Magneto-Hydrodynamics Parameters (MAGNETIC on) %--- Initial B-Field Strengths (if B_SET_IN_PARAMS on, otherwise read from IC file) BiniX 1.0e-10 % initial B_x, in code units BiniY 1.0e-10 % initial B_y, in code units BiniZ 1.0e-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 (~0.5-2) %--- Cosmic Ray + Gas Fluids (COSMIC_RAYS) CosmicRayDiffusionCoeff 1.0 % multiplies anisotropic diffusion/streaming coefficients CosmicRay_SNeFraction 0.1 % fraction of SNe ejecta kinetic energy into cosmic rays (~10%) %-------------- Grackle UVB file (GRACKLE on) GrackleDataFile CloudyData_UVB=HM2012.h5 % data file for GRACKLE cooling tables