# Show signs of life.
print "Loading spheral equations of state..."
# EOS constructors take a units object. I usually work in MKS.
units = sph.PhysicalConstants(1.0, # Unit length in meters
1.0, # Unit mass in kg
1.0) # Unit time in seconds
#-------------------------------------------------------------------------------
# Tillotson EOS for common materials
#-------------------------------------------------------------------------------
mats = ['Granite', 'Basalt', 'Nylon', 'Pure Ice', '30% Silicate Ice', 'Water']
etamin, etamax = 0.94, 10.0
pext, pmin, pmax = 0.0, -1e200, 1e200 # these are actually the defaults
EOSes = [sph.TillotsonEquationOfState(mat, 1e-20, 1e20, units,
etamin_solid = etamin,
externalPressure = pext, minimumPressure = pmin, maximumPressure = pmax)
for mat in mats]
granite = EOSes[0]
basalt = EOSes[1]
nylon = EOSes[2]
h2oice = EOSes[3]
dirtyice = EOSes[4]
water = EOSes[5]
del EOSes, mats, etamin, etamax
#-------------------------------------------------------------------------------
# M/ANEOS improved SiO2
#-------------------------------------------------------------------------------
try:
izetl = sph.vector_of_int(1, -1)
