# 1 Layer problem for solar absorption and infrared emissivity driven
# effects on atmospheric and surface temperature. 

import math

A=0.3         # Albedo.
eps=0.7       # infrared emissivity.
sigma=5.67e-8 # Stefan Boltzmann constant.
Io=342.0      # Solar irradiance.
MAE=10.0      # Mass absorption efficiency of black carbon aerosol in m^2/gram.
L=5000.0      # Soot layer thickness in meters.

print "SolarTransmissivity SolarAbsorptivity Tsurface_K Tatmosphere_K BC_Concentration_ug/m3"
for i in range (100):
    t=1.0 - i/100.0 # Solar transmissivity.
    Ts=(Io*(1.0+t)*(1.0-A*t)/(sigma*(2.0-eps)))**0.25 # Surface Temperature.
    Ta=(Io*((1.0-t)*(1.0+A*t)+eps*t*(1.0-A))/(sigma*eps*(2.0-eps)))**0.25 # Atmosphere temperature.
    BC=-math.log(t)*1.0e6/(MAE*L) # BC concentration in micrograms/m^3
    print t,1.0-t,Ts,Ta,BC