import math
from ClimateUtilities import *


def albedo(tau,g=0.,mu = math.cos(math.pi/4.)):
    beta = 1.-math.exp(-tau/mu)
    num = beta*(.5-gamma*mu)+ (1-g)*gamma*tau
    den = 1. + (1-g)*gamma*tau
    return num/den

def plotAlb(mu,g=0.):
	alb = [albedo(tau,g,mu) for tau in tauList]
	c = Curve()
	c.addCurve(tauList)
	c.addCurve(alb)
	plot(c)

tauList = [.01*i for i in range(1000)]

gamma = 1.