Pc = 2.477e17
Tc = 1.571e7
m = 1e-30 * Ms
print ""
print "Inner Boundary Conditions"
print "Testing Delta Movement of Temperature & Pressure from Boundary"
print " Initial Guess: (polytrope)"
print " P = %10.6e Dyne/cm^2" % Pc
print " T = %10.6e K" % Tc
print " Initial Step:"
print " m = %10.6e g" % m
(rs, ls, Ps, Ts) = inner_boundary(Pc=Pc,
Tc=Tc,
M=Ms,
mu=mmw(X=X, Y=Y),
m=m,
epsilon=epsilon,
optable=Opacity)
(rc, lc, Pc, Tc) = (0, 0, Pc, Tc)
print "At m=m:"
print " r = %10.6e cm" % rs
print " l = %10.6e erg/s" % ls
print " P = %10.6e Dyne/cm^2" % Ps
print " T = %10.6e K" % Ts
print " -- Initial ---------------- | -- 1st Step -------------- | -- Delta --------| "
print u" R = %10.6e cm, | R = %10.6e cm | ΔR = %11.4e |" % (
rs, rc, abs(rs - rc) / rs)
print u" L = %10.6e erg/s, | L = %10.6e erg/s | ΔL = %11.4e |" % (
ls, lc, abs(ls - lc) / ls)
print u" P = %10.6e Dyne/cm^2, | P = %10.6e Dyne/cm^2 | ΔP = %11.4e |" % (
# Tc = 1.4e7
Pc = 2.477e17
Tc = 1.571e7
m = 1e-10 * Ms
print "Inner Boundary Conditions"
print " Initial Guess: (polytrope)"
print " P = %10.6e Dyne/cm^2" % Pc
print " T = %10.6e K" % Tc
print " Initial Step:"
print " m = %10.6e g" % m
(r, l, P, T) = inner_boundary(Pc=Pc,
Tc=Tc,
M=Ms,
mu=mmw(X=X, Y=Y),
m=m,
optable=Opacity,
X=star.X,
XCNO=star.XCNO,
cfg=star.config["Data.Energy"])
print "At m=m:"
print " r = %10.6e cm" % r
print " l = %10.6e erg/s" % l
print " P = %10.6e Dyne/cm^2" % P
print " T = %10.6e K" % T
# print "Derivatives"
# (dr,dl,dP,dT) = derivatives(np.array([m]),np.atleast_2d([r,l,P,T]),mu=mmw(X=X,Y=Y),optable=Opacity,epsilon=epsilon)
# print " dr = %10.6e" % dr
# print " dl = %10.6e" % dl
# print " dP = %10.6e" % dP
# print " dT = %10.6e" % dT
print u" P = %10.6e Dyne/cm^2, | P = %10.6e Dyne/cm^2 | ΔP = %11.4e |" % (P,Pi,(P-Pi)/P) print u" T = %10.6e K, | T = %10.6e K | ΔT = %11.4e |" % (T,Ti,(T-Ti)/T) Pc = 2.477e17 Tc = 1.571e7 m = 1e-30*Ms print "" print "Inner Boundary Conditions" print "Testing Delta Movement of Temperature & Pressure from Boundary" print " Initial Guess: (polytrope)" print " P = %10.6e Dyne/cm^2" % Pc print " T = %10.6e K" % Tc print " Initial Step:" print " m = %10.6e g" % m (rs,ls,Ps,Ts) = inner_boundary(Pc=Pc,Tc=Tc,M=Ms,mu=mmw(X=X,Y=Y),m=m,epsilon=epsilon,optable=Opacity) (rc,lc,Pc,Tc) = (0,0,Pc,Tc) print "At m=m:" print " r = %10.6e cm" % rs print " l = %10.6e erg/s" % ls print " P = %10.6e Dyne/cm^2" % Ps print " T = %10.6e K" % Ts print " -- Initial ---------------- | -- 1st Step -------------- | -- Delta --------| " print u" R = %10.6e cm, | R = %10.6e cm | ΔR = %11.4e |" % (rs,rc,abs(rs-rc)/rs) print u" L = %10.6e erg/s, | L = %10.6e erg/s | ΔL = %11.4e |" % (ls,lc,abs(ls-lc)/ls) print u" P = %10.6e Dyne/cm^2, | P = %10.6e Dyne/cm^2 | ΔP = %11.4e |" % (Ps,Pc,abs(Ps-Pc)/Ps) print u" T = %10.6e K, | T = %10.6e K | ΔT = %11.4e |" % (Ts,Tc,abs(Ts-Tc)/Ts)
