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:
print "ANEOS SiO2 not loaded"
pass
#-------------------------------------------------------------------------------
# A polytropic fluid EOS
#-------------------------------------------------------------------------------
K = 2e5 # polytropic constant
n = 1 # polytropic index
mu = 2.2e-3 # mean molecular weight
poly = sph.PolytropicEquationOfStateMKS3d(K,n,mu)
del K, n, mu
#-------------------------------------------------------------------------------
# Available materials table
#-------------------------------------------------------------------------------
Materials = {'granite':'Granite solid (Tillotson)',
'basalt':'Basalt solid (Tillotson)',
'nylon':'Nylon solid (Tillotson)',
'h2oice':'Water ice solid (Tillotson)',
'dirtyice':'30% silicate in water ice (Tillotson)',
'water':'Liquid water (Tillotson)',
'SiO2':'Multi phase SiO2 (M/ANEOS)',
'poly':'Polytrope (n=1 K=2e5)',
}
print
