Such a feature is common in many atmospheric models. Physically, it also makes sense to use different advection schemes for the vertical and horizontal directions because of stratification. This capability might help address the spurious cooling at the surface caused by WENO seen in #21 and #215
