=use same air_refractive_index formula as used by phosim_3.2.1