Add permanent description of the changes made last commit
The shortwave radiation was being ignored when the function was given precomputed surface heat flux. This fixes that error and has been tested on both precomputed surface flux and component flux (i.e. longwave + shortwave + latent + sensible)