#Script for plotting the no-atmosphere equilibrium temperature
#of an isothermal planet, as a function of the distance
#from the Sun
#
#Data on section of text which this script is associated with
Chapter = '3.**'
Figure = '**'
#
#
#ToDo: Compute temperatures corresponding to actual planetary
# albedo, and plot as points on graph.
import phys,math
from ClimateUtilities import *
from planets import *
Tstar = 5800. #Photospheric temperature of the star
albedo = 0. #Planetary albedo (held constant for all planets)
#ToDo: Move solar radius into planet database
SolarRadius = 695.e6 #Solar radius in m
#List of orbital radii relative to stellar radius
#NormOrbit = [r*EarthOrbit/SolarRadius for r in Orbit]
rList = [5.*i for i in range(1,2000)]
TList = []
for r in rList:
TList.append(Tstar*(1.-albedo)**.25*math.sqrt(1./(2.*r)))
c = Curve()
c.addCurve(rList)
c.addCurve(TList)
c.XlogAxis = c.YlogAxis = True #Plot results on log-log axis
plot(c)
#Print out table of planetary "photosphere" temperatures based
#on actual observed albedo
for planet in [Mercury,Venus,Earth,Mars,Jupiter,Saturn,Uranus,Neptune,Pluto]:
r = planet.rsm/SolarRadius
T = Tstar*(1.-planet.albedo)**.25*math.sqrt(1./(2.*r))
print planet.name,r,T