 The discrete pressure equation is solved numerically with a fast elliptic solver in conjunction with a generalized conjugate gradient method.The fast elliptic solver is based on fast Fourier analysis in both horizontal directions and Gaussian elimination in the vertical direction (Moussiopoulos and Flassak, 1989).The second deviation from the explicit treatment is related to the turbulent diffusion in vertical direction.