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