Advertisement
Guest User

Gizmo Parameter File

a guest
Mar 25th, 2018
319
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
INI file 10.37 KB | None | 0 0
  1. %-------------------------------------------------------------------------
  2. %----  This file contains the input parameters needed at run-time for
  3. %       simulations. It is based on and closely resembles the GADGET-3
  4. %       parameterfile (format of which and parsing routines written by
  5. %       Volker Springel [volker.springel@h-its.org]). It has been updated
  6. %       with new naming conventions and additional variables as needed by
  7. %       Phil Hopkins [phopkins@caltech.edu] for GIZMO.
  8. %-------------------------------------------------------------------------
  9.  
  10. %----  Relevant files
  11. InitCondFile  gizmo_ch1_mill.dat
  12. %%%InitCondFile snapshot_149
  13. OutputDir     .
  14.  
  15. %---- File formats
  16. ICFormat    1  % 1=binary, 3=hdf5, 4=cluster
  17. %%%ICFormat    3  % 1=binary, 3=hdf5, 4=cluster
  18. SnapFormat  3  % 1=binary, 3=hdf5
  19.  
  20. %---- Output parameters
  21. RestartFile                 restart
  22. SnapshotFileBase            snapshot
  23. OutputListOn                1  % =1 to use list in "OutputListFilename"
  24. OutputListFilename          output_times.txt  % list of times for snaps
  25. NumFilesPerSnapshot         1
  26. NumFilesWrittenInParallel   1  % must be < N_processors & power of 2
  27.  
  28. %---- Output frequency
  29. TimeOfFirstSnapshot     0.02
  30. TimeBetSnapshot         1.5
  31.  
  32. TimeBetStatistics       0.01
  33.  
  34. %---- CPU-time limits
  35. TimeLimitCPU            4500000  % in seconds
  36. CpuTimeBetRestartFile   1250    % in seconds
  37. %%%CpuTimeBetRestartFile   43200    % in seconds
  38. ResubmitOn              0
  39. ResubmitCommand         my-scriptfile
  40.  
  41. %----- Memory alloction
  42. MaxMemSize          3000    % sets maximum MPI process memory use in MByte
  43. %%%MaxMemSize          3500    % sets maximum MPI process memory use in MByte
  44. PartAllocFactor     10.0     % memory load allowed for better cpu balance
  45. BufferSize          10     % in MByte
  46.  
  47. %---- Characteristics of run
  48. TimeBegin   0.0049751243    % Beginning of the simulation
  49. %%%TimeBegin   0.04690762567    % Beginning of the simulation
  50. TimeMax     1.0    % End of the simulation
  51.  
  52. %---- Cosmological parameters
  53. ComovingIntegrationOn   1       % is it cosmological? (yes=1, no=0)
  54. BoxSize                 500000.  % in code units
  55. Omega0                  0.3089    % =0 for non-cosmological
  56. OmegaLambda             0.6911    % =0 for non-cosmological
  57. OmegaBaryon             0.04860   % =0 for non-cosmological
  58. HubbleParam             0.6774     % little 'h'; =1 for non-cosmological runs
  59.  
  60. %---- Accuracy of time integration
  61. ErrTolIntAccuracy       0.01   % <0.02
  62. CourantFac              0.1     % <0.20
  63. MaxRMSDisplacementFac   0.125   % <0.25
  64. MaxSizeTimestep         0.5
  65. MinSizeTimestep         1.0e-12
  66.  
  67. %---- Tree algorithm, force accuracy, domain update frequency
  68. ErrTolTheta                 0.7     % 0.7=standard
  69. ErrTolForceAcc              0.0025  % 0.0025=standard
  70. TreeDomainUpdateFrequency   0.005   % 0.0005-0.05, dept on core+particle number  
  71.  
  72. %---- System of units
  73. UnitLength_in_cm            3.085678e21     % 1.0 kpc/h
  74. UnitMass_in_g               1.989e43        % 1.0e10 solar masses/h
  75. UnitVelocity_in_cm_per_s    1e5             % 1 km/sec
  76. GravityConstantInternal     0               % calculated by code if =0
  77.  
  78. %---- Initial temperature & temperature floor
  79. InitGasTemp     300     % set by IC file if =0
  80. MinGasTemp      10.     % don't set <10 in explicit feedback runs, otherwise 0
  81.  
  82. %---- Density/volume estimation (kernel)
  83. DesNumNgb               32      % 32 for standard kernel, 60-114 for quintic
  84. MaxNumNgbDeviation      0.25    % <<DesNumNgb (values<1 are fine)
  85. MaxHsml                 1.0e10  % some very large value (here to prevent errors)
  86. MinGasHsmlFractional    1.0     % minimum gas smoothing in terms of softening  (=1)
  87.  
  88. %---- Gravitational softening lengths
  89. %----- Softening lengths per particle type. If ADAPTIVE_GRAVSOFT is set, these  
  90. %-------- are the minimum softening allowed for each type -------
  91. %-------- (units are co-moving for cosmological integrations)
  92. %%%SofteningGas    1.0e-10    % gas (type=0) (in units above, =1 pc softening)
  93. %%%SofteningHalo   1.0e-10    % dark matter/collisionless particles (type=1)
  94. SofteningGas    2500   % gas (type=0) (in units above, =1 pc softening)
  95. SofteningHalo   5000    % dark matter/collisionless particles (type=1)
  96. SofteningDisk   2500    % collisionless particles (type=2)
  97. SofteningBulge  2500    % collisionless particles (type=3)
  98. SofteningStars  2500    % stars spawned from gas (type=4)
  99. SofteningBndry  2500    % black holes (if active), or collisionless (type=5)
  100. %---- if these are set in cosmo runs, SofteningX switches from comoving to physical
  101. %------- units when the comoving value exceeds the choice here
  102. SofteningGasMaxPhys     2500   % switch to 0.5pc physical below z=1
  103. SofteningHaloMaxPhys    5000
  104. SofteningDiskMaxPhys    2500
  105. SofteningBulgeMaxPhys   2500
  106. SofteningStarsMaxPhys   2500
  107. SofteningBndryMaxPhys   2500
  108. %----- parameters for adaptive gravitational softening
  109. AGS_DesNumNgb           32  % neighbor number for calculating adaptive gravsoft
  110. AGS_MaxNumNgbDeviation  8   % tolerance in Nngb (make larger than gas)
  111.  
  112.  
  113. %---------- SPH-Specific Parameters ---------------------------------
  114. %---- Artificial viscosity  
  115. ArtBulkViscConst    1.0     % multiplies 'standard' AV (use 1.0)
  116. %---- P&M artificial conductivity (if present); normalized to Alpha_Visc:
  117. ArtCondConstant     0.50    % multiplies 'standard' (use 0.25-0.5)
  118. %---- Cullen & Dehnen viscosity suppression
  119. ViscosityAMin       0.05    % minimum viscosity away from shocks (>0.025)
  120. ViscosityAMax       2.00    % maximum viscosity in shocks (>1)
  121.  
  122.  
  123. %------------------------------------------------------------
  124. %---------- Physics Modules ---------------------------------
  125. %------------------------------------------------------------
  126.  
  127. %---- Turbulent Diffusion Master Switch (TURB_DIFFUSION on)
  128. TurbDiffusionCoefficient    1.0     % Normalizes diffusion rates, should be ~1
  129.  
  130. %---- Magneto-Hydrodynamics Parameters (MAGNETIC on)
  131. %--- Dedner 2002 div-cleaning parameters
  132. DivBcleaningParabolicSigma      1.0     % (3D~0.7-1.0,2D~0.2-0.3)
  133. DivBcleaningHyperbolicSigma     1.0     % (~1)
  134. DivBcleaningQ                   0.5     % (~0.5-1)
  135. %--- Initial B-Field Strengths (BINISET on) (in Gauss)
  136. BiniX   1.0e-8    % B_x, in Gauss
  137. BiniY   1.0e-8
  138. BiniZ   1.0e-8
  139. %-------------------- SPH-Specific Parameters----------------------
  140. ArtificialResistivityMax    1.  % maximum alpha_B (~1-2) for art. res. (like art. visc)
  141.  
  142.  
  143. %---- Star Formation parameters (GALSF on)
  144. CritPhysDensity     100.    %  critical physical density for star formation (cm^(-3))
  145. SfEffPerFreeFall    1.0     %  SFR/(Mgas/tfreefall) for gas which meets SF criteria
  146.  
  147.  
  148. %---- sub-grid (Springel+Hernquist/GADGET/AREPO) "effective equation of state"
  149. %------- star formation+feedback model (GALSF_EFFECTIVE_EQS on)
  150. MaxSfrTimescale     4.0         % code units (SF timescale at 2-phase threshold)
  151. TempSupernova       3.0e8       % in Kelvin (temp of hot gas in 2-phase model)
  152. TempClouds          1000.0      % in Kelvin (temp of cold gas in 2-phase model)
  153. FactorSN            0.1         % SNe coupling frac (frac of egy retained in hot)
  154. FactorEVP           3000.0      % controls Kennicutt normalization
  155. FactorForSofterEQS  1.0         % interpolate between 'stiff' and isothermal EOS
  156. %------- the sub-grid "decoupled winds" model (GALSF_SUBGRID_WINDS on)
  157. WindEfficiency          2.0     % mass-loading (Mdot_wind = SFR * WindEfficiency)
  158. WindEnergyFraction      0.06    % fraction of SNe energy in winds (sets velocity)
  159. WindFreeTravelMaxTime   0.1     % 'free-stream time' in units of t_Hubble(z)
  160. WindFreeTravelDensFac   0.1     % 'free-stream' until density < this * CritPhysDensity
  161.  
  162.  
  163. %-------------- FIRE (PFH) explicit star formation & feedback model
  164. %--- initial metallicity of gas & stars in simulation
  165. InitMetallicity             0.0001  % (in solar)
  166. InitStellarAge              0.001   % mean age (in Gyr)
  167. %--- local radiation-pressure driven winds (GALSF_FB_RPWIND_LOCAL)
  168. WindMomentumLoading         1.0     % boost momentum-loading (relative to ssp)
  169. WindInitialVelocityBoost    1.0     % boost velocity coupled (fixed momentum)
  170. %--- SneII Heating Model (GALSF_FB_SNE_HEATING)
  171. SNeIIEnergyFrac             1.0     % boost fraction of mechanical energy to couple
  172. SNeIIBW_Radius_Factor       1.0     % boost cooling radius for resolution-check
  173. %--- HII region photo-heating model (GALSF_FB_HII_HEATING)  
  174. HIIRegion_Temp              10000.  % temperature (in K) of heated gas
  175. HIIRegion_fLum_Coupled      1.0     % boost fraction of ionizing photons
  176. %--- long-range radiation pressure acceleration (GALSF_FB_RT_PHOTONMOMENTUM)
  177. PhotonMomentum_Coupled_Fraction  1.0  % boost fraction of L to allow incident
  178. PhotonMomentum_fUV          0.01    % incident SED f(L) in UV (minimum scattering)
  179. PhotonMomentum_fOPT         0.01    % incident SED f(L) in optical/near-IR
  180. %--- gas return/recycling (GALSF_FB_GASRETURN)
  181. GasReturnFraction           1.0     % boost fraction of gas returns (relative to ssp)
  182. GasReturnEnergy             1.0     % boost temperature of gas returns (relative to ssp)
  183.  
  184.  
  185. %-------------- Black Hole accretion & formation (BLACK_HOLES on)
  186. %--- formation/seeding  
  187. MinFoFMassForNewSeed         8.     % log(M_stars) for on-the-fly FoF seeding
  188. TimeBetOnTheFlyFoF           800     % time (in sec) between FoF searches
  189. SeedBlackHoleMass            1.0e-7     % initial mass (on-the-fly or single galaxy)
  190. %--- accretion  
  191. BlackHoleAccretionFactor     1.0        % multiplier for mdot
  192. BlackHoleEddingtonFactor     1.0        % fraction of eddington to cap (can be >1)
  193. BlackHoleNgbFactor           3.0        % multiplier for kernel neighbors for BH
  194. BlackHoleMaxAccretionRadius  0.02       % max radius for BH neighbor search/accretion
  195. BlackHoleRadiativeEfficiency 0.1        % radiative efficiency
  196. SeedAlphaDiskMass            100.
  197. SeedBlackHoleMassSigma       2.0
  198. SeedBlackHoleMinRedshift     0.1
  199. massDMpart                   1.0
  200. %--- feedback
  201. BlackHoleFeedbackFactor      1.0        % generic feedback strength multiplier
  202. BH_FluxMomentumFactor        1.0        % multiply radiation pressure (BH_PHOTONMOMENTUM)
  203. BAL_f_accretion              0.5        % fraction of gas swallowed by BH (BH_BAL_WINDS)
  204. BAL_v_outflow                30000      % velocity (km/s) of BAL outflow (BH_BAL_WINDS)
  205.  
  206. % Sink
  207. SinkHsml 0.02
  208. SinkDensThresh 10000
  209.  
  210. %-------------- Grackle UVB file (GRACKLE on)
  211. GrackleDataFile          CloudyData_UVB=HM2012.h5
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement