%-------------------------------------------------------------------------
%----  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  /work/01799/phopkins/sim_ICs/ics_z10n192_B1hr 
OutputDir     /scratch/01799/phopkins/B1_lr_fixg_51 

%---- 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          output_times.txt  % list of times for snaps 
NumFilesPerSnapshot         1 
NumFilesWrittenInParallel   16  % must be < N_processors & power of 2

%---- Output frequency 
TimeOfFirstSnapshot     0.1 
TimeBetSnapshot         1.1 
TimeBetStatistics       0.05 

%---- CPU-time limits 
TimeLimitCPU            100000  % in seconds 
CpuTimeBetRestartFile   7200  	% in seconds 
ResubmitOn              0 
ResubmitCommand         my-scriptfile 

%----- Memory allocation
MaxMemSize          1800    % sets maximum MPI process memory use in MByte 
PartAllocFactor     5.0     % memory load allowed for better cpu balance 
BufferSize          100     % in MByte 

%---- Characteristics of run 
TimeBegin   0.01    % Beginning of the simulation 
TimeMax     1.0     % End of the simulation 

%---- Cosmological parameters 
ComovingIntegrationOn   1       % is it cosmological? (yes=1, no=0)
BoxSize                 10000.  % in code units
Omega0                  0.27    % =0 for non-cosmological
OmegaLambda             0.73    % =0 for non-cosmological
OmegaBaryon             0.044   % =0 for non-cosmological
HubbleParam             0.7     % little 'h'; =1 for non-cosmological runs

%---- Accuracy of time integration 
MaxSizeTimestep         0.005   % in code units, set for your problem
MinSizeTimestep         1.0e-12 % set this very low, or risk stability

%---- Tree algorithm, force accuracy, domain update frequency 
TreeDomainUpdateFrequency   0.005	% 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     300	    % 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.0e10  % minimum gas kernel length (some very large value to prevent errors)
MinGasHsmlFractional    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)
SofteningGas    0.001    % gas (type=0) (in units above, =1 pc softening)
SofteningHalo   0.020    % dark matter/collisionless particles (type=1)
SofteningDisk   0.150    % collisionless particles (type=2)
SofteningBulge  0.500    % collisionless particles (type=3)
SofteningStars  0.001    % stars spawned from gas (type=4)
SofteningBndry  0.001    % 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.0005    % switch to 0.5pc physical below z=1 
SofteningHaloMaxPhys    0.010 
SofteningDiskMaxPhys    0.075 
SofteningBulgeMaxPhys   0.250 
SofteningStarsMaxPhys   0.0005 
SofteningBndryMaxPhys   0.0005 
%----- 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-8    % initial B_x, in code units
BiniY   1.0e-8    % initial B_y, in code units
BiniZ   1.0e-8    % 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     100.    %  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.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     % 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  
MinFoFMassForNewSeed         10.	    % log(M_stars) for on-the-fly FoF seeding 
TimeBetOnTheFlyFoF           200000     % time (in sec) between FoF searches 
SeedBlackHoleMass            1.0e-7 	% initial mass (on-the-fly or single galaxy) 
%--- 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
