A stable and accurate convective modeling procedure based on quadratic upstream interpolation