We study analytic solutions to the Bardeen-Cooper-Schrieffer (BCS) gap equation for isotropic superconductors with finite-range interaction potentials over the full range of temperatures from absolute zero to the superconducting critical temperature 0≤T≤Tc. Using these solutions Δ(ϵ,T), we provide a proof of the universality of the temperature dependence of the BCS gap ratio at the Fermi level Δ(ϵ=0,T)/Tc. Moreover, by examining the behavior of this ratio as a function of energy ϵ, we find that nonuniversal features emerge away from the Fermi level, and these features take the form of a temperature-independent multiplicative factor F(ϵ), which is equal to Δ(ϵ,T)/Δ(ϵ=0,T) up to exponentially small corrections, i.e., the error terms vanish like e-1/λ in the weak-coupling limit λ→0. We discuss the model-dependent features of both F(ϵ) and Tc, and we illustrate their behavior focusing on several concrete examples of physically relevant finite-range potentials. Comparing these cases for fixed coupling constants, we highlight the importance of the functional form of the interaction potential in determining the size of the critical temperature and provide guidelines for choosing potentials which lead to higher values of Tc. We also propose experimental signatures which could be used to probe the energy dependence of the gap and potentially shed light on the underlying mechanisms giving rise to superconductivity.