The first results of a planned continuing investigation of astronomical refraction are reported in this paper. Pure and anomalous refraction are investigated using meteorological data and environmental models. Numerical methods based upon Chebyshev polynomials are developed. Applications are made to Winslow, Arizona (which is somewhat similar to El Leoncito, Argentina), and to results of PZT observations made in Washington, D.C.