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)
sph.initializeANEOS('/proj/nmovshov_hindmost/collisions/ANEOS/ANEOS.INPUT',
'ANEOS.barf', izetl)
etamin, etamax = 0.94, 10
rho0 = 2650
pext, pmin, pmax = 0.0, -1e200, 1e200 # these are actually the defaults
SiO2 = sph.ANEOS(0, # Material number
1000, # num rho vals
1000, # num T vals
etamin*rho0, # minimum density (kg/m^3)
etamax*rho0, # maximum density (kg/m^3)
1.0, # minimum temperature (K)
1.0e4, # maximum temperature (K)
units,
pext, pmin, pmax)
os.system('rm -f ANEOS.barf')
del izetl, etamin, etamax, rho0
except:
