#Script for generating the figure illustrating
#how the balance between OLR and absorbed solar radiation
#determines the surface temperature
#
#Data on section of text which this script is associated with
Chapter = '3.**'
Figure = '**'
#

from ClimateUtilities import *
import phys

#Define albedo and OLR functions here

def OLR(T,pRad):
    return phys.sigma * (T**4.)*(pRad/1000.)**(4.*2./7.) 



#Main body of script starts here
S = 1370./4. #Incident solar per unit surface area
albedo = .3


Tlist = [200.+ 2.*i for i in range(60)]
Slist = [(1.-albedo)*S for T in Tlist] # Generates a list of constant Sabs

#Plot the results
c = Curve()
c.addCurve(Tlist)
c.addCurve(Slist,'Sabs','Absorbed Solar')
c.addCurve([OLR(T,1000.) for T in Tlist],'OLR1','Outgoing Longwave,p_rad = 1000mb')
c.addCurve([OLR(T,600.) for T in Tlist],'OLR2','Outgoing Longwave,p_rad = 600mb')
w = plot(c)