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: