#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