High-resolution atmospheric simulations can become numerically unstable over steep terrain. Traditional terrain smoothing approaches modify the domain globally, reducing terrain detail. We developed a local smoothing method that improves simulation stability while only modifying steep-slope points, helping to retain the benefits of high-resolution modeling. Tested in a mesoscale and a microscale atmospheric model, it is computationally efficient, easy to implement, and adaptable to other models.