Adding the blackbody radiation term was simple. We added a fixed differential for the blackbody radiation equation at the top-most level of cells in the i+1 direction.
We assumed the sink temperature (T0) was the average temperature of space, which means there are no clouds above our model.