CosmicRayDiffusionCoeff 69.0 % multiplies anisotropic diffusion/streaming coefficients CosmicRay_SNeFraction 0.1 % fraction of SNe ejecta kinetic energy into cosmic rays (~10%) %------------------------------------------------------------------------- %---- This file contains the input parameters needed at run-time for % simulations. It is based on and closely resembles the GADGET-3 % parameterfile (format of which and parsing routines written by % Volker Springel [volker.springel@h-its.org]). It has been updated % with new naming conventions and additional variables as needed by % Phil Hopkins [phopkins@caltech.edu] for GIZMO. %------------------------------------------------------------------------- %---- Relevant files InitCondFile ICs/m09_cube_ref12.ics OutputDir output %---- File formats ICFormat 1 % 1=binary, 3=hdf5, 4=cluster SnapFormat 3 % 1=binary, 3=hdf5 %---- Output parameters RestartFile restart SnapshotFileBase snapshot OutputListOn 1 % =1 to use list in "OutputListFilename" OutputListFilename snapshot_scale-factors.txt % list of scale factors for snapshots NumFilesPerSnapshot 8 NumFilesWrittenInParallel 32 % must be < N_processors & power of 2 %---- Output frequency TimeOfFirstSnapshot 0.1 TimeBetSnapshot 1.1 TimeBetStatistics 1.e10 %0.05 %---- CPU-time limits TimeLimitCPU 999999999 % in seconds CpuTimeBetRestartFile 3450 % in seconds ResubmitOn 0 ResubmitCommand sbatch submit_gizmo_slurm.py %----- Memory allocation MaxMemSize 7800 % sets maximum MPI process memory use in MByte PartAllocFactor 2.8 %2.2 %5.0 % memory load allowed for better cpu balance BufferSize 30 %50 %200 % in MByte %---- Characteristics of run TimeBegin 0.0079365079 % Beginning of the simulation TimeMax 1.0 % End of the simulation %---- Cosmological parameters ComovingIntegrationOn 1 % is it cosmological? (yes=1, no=0) BoxSize 5000.0 % in code units Omega0 0.2660 % =0 for non-cosmological OmegaLambda 0.7340 % =0 for non-cosmological OmegaBaryon 0.0440 % =0 for non-cosmological HubbleParam 0.7100 % little 'h'; =1 for non-cosmological runs %---- Accuracy of time integration MaxSizeTimestep 0.005 % in code units, set for your problem MinSizeTimestep 1.0e-9 % set this very low, or risk stability %---- Tree algorithm, force accuracy, domain update frequency TreeDomainUpdateFrequency 5.e-7 %5.e-6 %0.05 % 0.0005-0.05, dept on core+particle number %---- System of units UnitLength_in_cm 3.085678e21 % 1.0 kpc/h UnitMass_in_g 1.989e43 % 1.0e10 solar masses/h UnitVelocity_in_cm_per_s 1.0e5 % 1 km/sec UnitMagneticField_in_gauss 1.0 % 1 gauss GravityConstantInternal 0 % calculated by code if =0 %---- Initial temperature & temperature floor InitGasTemp 0 % set by IC file if =0 MinGasTemp 10. % don't set <10 in explicit feedback runs, otherwise 0 %---- Density/volume estimation (kernel) DesNumNgb 32 % 32 for standard kernel, 60-114 for quintic MaxHsml 1.0e6 % minimum gas kernel length (some very large value to prevent errors) MinGasHsmlFractional 1 % 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) SofteningGas 0.0005 % gas (type=0) (in units above, =1 pc softening) SofteningHalo 0.100 % dark matter/collisionless particles (type=1) SofteningDisk 1.000 % collisionless particles (type=2) SofteningBulge 2.000 % collisionless particles (type=3) SofteningStars 0.001 % stars spawned from gas (type=4) SofteningBndry 2.000 % black holes (if active), or collisionless (type=5) %---- if these are set in cosmo runs, SofteningX switches from comoving to physical %------- units when the comoving value exceeds the choice here %------- (these are ignored, and *only* the above are used, for non-cosmo runs) SofteningGasMaxPhys 0.00005 % switch to 0.5pc physical below z=10 SofteningHaloMaxPhys 0.010 SofteningDiskMaxPhys 0.100 SofteningBulgeMaxPhys 0.200 SofteningStarsMaxPhys 0.0001 SofteningBndryMaxPhys 0.200 %----- parameters for adaptive gravitational softening AGS_DesNumNgb 32 % neighbor number for calculating adaptive gravsoft %------------------------------------------------------------------------- %------------------------------------------------------------------------- %---------- Physics Modules ---------------------------------------------- %------------------------------------------------------------------------- %------------------------------------------------------------------------- %------------------------------------------------------------ %------------------ Additional Fluid Physics ---------------- %------------------------------------------------------------ %---- 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 %------------------------------------------------------------------------- %------------------ Star, Black Hole, and Galaxy Formation --------------- %------------------------------------------------------------------------- %---- Star Formation parameters (GALSF on) CritPhysDensity 1.0e5 % critical physical density for star formation (cm^(-3)) SfEffPerFreeFall 1.0 % SFR/(Mgas/tfreefall) for gas which meets SF criteria %---- sub-grid (Springel+Hernquist/GADGET/AREPO) "effective equation of state" %------- star formation+feedback model (GALSF_EFFECTIVE_EQS on) MaxSfrTimescale 4.0 % code units (SF timescale at 2-phase threshold) TempSupernova 3.0e8 % in Kelvin (temp of hot gas in 2-phase model) TempClouds 1000.0 % in Kelvin (temp of cold gas in 2-phase model) FactorSN 0.1 % SNe coupling frac (frac of egy retained in hot) FactorEVP 3000.0 % controls Kennicutt normalization FactorForSofterEQS 1.0 % interpolate between 'stiff' and isothermal EOS %------- the sub-grid "decoupled winds" model (GALSF_SUBGRID_WINDS on) WindEfficiency 2.0 % mass-loading (Mdot_wind = SFR * WindEfficiency) WindEnergyFraction 0.06 % fraction of SNe energy in winds (sets velocity) WindFreeTravelMaxTime 0.1 % 'free-stream time' in units of t_Hubble(z) WindFreeTravelDensFac 0.1 % 'free-stream' until density < this * CritPhysDensity %------- alternative winds (GALSF_SUBGRID_DMDISPERSION or GALSF_SUBGRID_VARIABLEVELOCITY) %------- (scaling with local dark matter properties, as Dave/Oppenheimer/Mannucci/Illustris) VariableWindVelFactor 1.0 % wind velocity relative to estimated halo v_escape VariableWindSpecMomentum 5000. % wind momentum per unit stellar mass (code velocity units) %-------------- FIRE (PFH) explicit star formation & feedback model %--- initial metallicity of gas & stars in simulation InitMetallicity 0.0001 % initial gas+stellar metallicity (in solar) InitStellarAge 0.001 % initial mean age (in Gyr; for stars in sim ICs) %--- local radiation-pressure driven winds (GALSF_FB_RPWIND_LOCAL) WindMomentumLoading 1.0 % fraction of photon momentum to couple %--- SneII Heating Model (GALSF_FB_SNE_HEATING) SNeIIEnergyFrac 1.0 % fraction of mechanical energy to couple %--- HII region photo-heating model (GALSF_FB_HII_HEATING) HIIRegion_fLum_Coupled 1.0 % fraction of ionizing photons allowed to see gas %--- long-range radiation pressure acceleration (GALSF_FB_RT_PHOTONMOMENTUM) 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 (GALSF_FB_GASRETURN) 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 1.0e-7 % initial mass (on-the-fly or single galaxy) SeedAlphaDiskMass 0.0 % initial mass in the alpha disk (for idealized tests: if BH_ALPHADISK_ACCRETION on) %----- (specific options for on-the-fly friends-of-friends based BH seeding: FOF on) MinFoFMassForNewSeed 10. % log(M_stars) for on-the-fly FoF seeding TimeBetOnTheFlyFoF 200000 % time (in sec) between FoF searches SeedBlackHoleMinRedshift 2.0 % this minimum redshift where new BH particles are seeded via FoF SeedBlackHoleMassSigma 0.5 % log normal standard deviation in the init BH seed masses from FoF %--- accretion BlackHoleAccretionFactor 1.0 % multiplier for mdot BlackHoleEddingtonFactor 1.0 % fraction of eddington to cap (can be >1) BlackHoleNgbFactor 3.0 % multiplier for kernel neighbors for BH BlackHoleMaxAccretionRadius 0.02 % max radius for BH neighbor search/accretion BlackHoleRadiativeEfficiency 0.1 % radiative efficiency %--- feedback BlackHoleFeedbackFactor 1.0 % generic feedback strength multiplier BH_FluxMomentumFactor 1.0 % multiply radiation pressure (BH_PHOTONMOMENTUM) BAL_f_accretion 0.5 % fraction of gas swallowed by BH (BH_BAL_WINDS) BAL_v_outflow 30000 % velocity (km/s) of BAL outflow (BH_BAL_WINDS) %-------------- Grackle UVB file (GRACKLE on) GrackleDataFile CloudyData_UVB=HM2012.h5 %------------------------------------------------------------------------- %------------------ Driven Turbulence (Large-Eddy boxes) ----------------- %------------------------------------------------------------------------- %-------------- Turbulent stirring parameters (TURB_DRIVING on) ST_decay 1. % decay time for driving-mode phase correlations ST_energy 0.0002 % energy of driving-scale modes: sets norm of turb ST_DtFreq 0.005 % time interval for driving updates (set by hand) ST_Kmin 6.27 % minimum driving-k: should be >=2.*M_PI/All.BoxSize ST_Kmax 12.57 % maximum driving-k: set to couple times Kmin or more if cascade desired ST_SolWeight 1.0 % fractional wt of solenoidal modes (wt*curl + (1-wt)*div) ST_AmplFac 1.0 % multiplies turb amplitudes ST_SpectForm 2 % driving pwr-spec: 0=Ek~const; 1=sharp-peak at kc; 2=Ek~k^(-5/3); 3=Ek~k^-2 ST_Seed 42 % random number seed for modes (so you can reproduce it) IsoSoundSpeed 1.0 % initializes gas sound speed in box to this value TimeBetTurbSpectrum 0.5 % time (code units) between evaluations of turb pwrspec