EBM STELLA Equations

Atmosphere(t) = Atmosphere(t - dt) + (SW_to_Atm + Land_to_Atm_LW + Ocean_to_ATM_LW - Atm_to_Ocean_LW - Atm_to_Land_LW - Atm_LW_loss) * dt
INIT Atmosphere = 16.5

INFLOWS:
SW_to_Atm = (Solar_Input*(.18/.44)*Cloud_Cover)*(1+Percent_increase_in_solar_input)
Land_to_Atm_LW = 119*Land_fraction*Land_Temp/288
Ocean_to_ATM_LW = 119*Ocean_fraction*Ocean_Temp/288
OUTFLOWS:
Atm_to_Ocean_LW = 77*Ocean_fraction*Atm_Temp/255
Atm_to_Land_LW = 77*Land_fraction*Atm_Temp/255
Atm_LW_loss = 60*Atm_Temp/255
Land(t) = Land(t - dt) + (Atm_to_Land_LW + SW_to_Land - Land_to_Atm_LW - Land_LW_loss) * dt
INIT Land = Initial_Land

INFLOWS:
Atm_to_Land_LW = 77*Land_fraction*Atm_Temp/255
SW_to_Land = (Solar_Input*(1+Percent_increase_in_solar_input))*Land_fraction*((Cloud_Cover*(3.24/13.2))+(1-Cloud_Cover)*(1-Land_albedo))
OUTFLOWS:
Land_to_Atm_LW = 119*Land_fraction*Land_Temp/288
Land_LW_loss = 8*Land_fraction*Land_Temp/288
Ocean(t) = Ocean(t - dt) + (Atm_to_Ocean_LW + SW_to_Ocean - Ocean_LW_loss - Ocean_to_ATM_LW) * dt
INIT Ocean = 270*(Ocean_fraction/0.7)*(Mixed_Layer/35)

INFLOWS:
Atm_to_Ocean_LW = 77*Ocean_fraction*Atm_Temp/255
SW_to_Ocean = (Solar_Input*(1+Percent_increase_in_solar_input))*Ocean_fraction*((Cloud_Cover*(1.288/30.8))+(1-Cloud_Cover)*(1-ocean_albedo))
OUTFLOWS:
Ocean_LW_loss = 8*Ocean_fraction*Ocean_Temp/288
Ocean_to_ATM_LW = 119*Ocean_fraction*Ocean_Temp/288
Atm_Temp = 255*Atmosphere/16.5
Cloud_Cover = .44
GLOBAL_TEMP = (Land_Temp*Land_fraction)+(Ocean_Temp*Ocean_fraction)
Global_temp_in_Celsius = GLOBAL_TEMP-273.15
Initial_Land = 0.6*(Land_fraction/0.3)*Soil_thickness
Initial_Ocean = 270*(Ocean_fraction/0.7)*(Mixed_Layer/35)
Land_albedo = .3
Land_fraction = 1-Ocean_fraction
Land_Temp = 288*Land/Initial_Land
Mixed_Layer = 35
ocean_albedo = .14
Ocean_fraction = .7
Ocean_Temp = 288*Ocean/Initial_Ocean
Percent_increase_in_solar_input = 0
Soil_thickness = 1
Solar_Input = 100
ĘT_Atmosphere = (Atm_Temp-255)+18
ĘT_Global = (GLOBAL_TEMP-288)-15
ĘT_Land = (Land_Temp-288)-15
ĘT_Ocean = (Ocean_Temp-288)-15