Institute for Geography University of Augsburg Physical Geography and Quantitative Methods Energy Balance Model |
albedo = stdalbedo * ( 1.D0 - icealbedo * tanh( 1.548D0 * (temp-288.D0) ) ! albedo depending on temp
radiationoutput_shortwave = albedo * radiationinput_shortwave ! reflected short wave radiation
noise = random_normal() * noisefactor ! simulating small weather disturbances
radiationoutput_longwave = (emis * sigma * temp**4) * (transmiss + noise) ! long wave emission
radiationbalance = radiationinput_shortwave - radiationoutput_shortwave - radiationoutput_longwave ! radiation balance
delta_t = radiationbalance / heat_capacity ! temperature change
temp = temp + delta_t * timestep ! new temperature