Atmos and surface models update radiation in different frequencies

[szy21]

Currently, radiation is called every dt_rad (typically one hour) in the atmosphere, but is called every timestep in surface models (land or sea ice). This is inconsistent and energy is not conserved. To fix this, we can still calculate radiation at every timestep for land and pass the averaged quantities (albedo, T_sfc, emissivity) to the atmosphere. For sea ice we can do the same thing or move radiation to a callback.