“…With this procedure it is necessary to apply the column approximation to the outer boundary of the mesh, i.e., in equation (3) ag = 1, b = 0, and cg = cp,.is calculated at the boundary using fhe column approximation. In order to avoid distortion of the image the step-size, Az, must be chosen carefully and the distance between columns (meshsize) must be sufficiently small (Anstis and Cockayne, 1979). The other method employs "free" boundary conditions defined by ag = 0, b = 1, and dc dz = 0, and which provides automatic control of the integration (e.g., the choice of Az) and provides error-estimates.…”