Earth_Energy(t) = Earth_Energy(t - dt) + (Solar_to_Earth - Infrared_to_Space) * dt
INIT Earth_Energy = 0.0 {We do not know to put here yet. Let the model compute it for us.}

Solar_to_Earth = Solar_Constant {J/m2 yr} * (1-Earth_Albedo) * Earth_Crossection {m2}
Infrared_to_Space = Earth_Surface_Area {m2} * Stefan_Boltzmann {J/m2 yr k4} * Temperature^4 {K4}
Density_Water = 1000. {kg/m3}
Earth_Albedo = 0.30 {30% as a fraction}
Earth_Crossection = PI*Earth_Diameter^2/4
Earth_Diameter = 12742e3 {m}
Earth_Surface_Area = PI*Earth_Diameter^2
Heat_Capacity = Water_Depth * Earth_Surface_Area * Density_Water * Specific_Heat_Water
Solar_Constant = 1376 {J/m2 s} * 3.15576E7 {s/yr}
Specific_Heat_Water = 4218. {J/kg K}
Stefan_Boltzmann = 5.67E-8 {J/m2 s K4} * 3.15576E7 {s/yr}
Temperature = Earth_Energy / Heat_Capacity {1st Law of Thermodynamics}
Water_Depth = 1.0 {m, temporary assumption}