In this paper, Liebmann method (point iterative) and Local Green Function method (local point iterative) which are usually used in traditional atmospheric/oceanic calculating problems are modified to parallel algorithms adapting to the dual-CPU in IBM-4381 (P03) to solve the Poisson equation for zero-boundary value. The numerical experimental results show that the calculation speed of parallel algorithm is 1.8 or 1.9 times as much as serial algorithm for the same calculating problems. However, when expanding the areas, parallel iterative speed of Local Green Function can reach 2.8 time which exceeds the best speedup (1.7 times for P01) declared by IBM Company. The paper also points out that the data array can be shared by the IBM-4381 (P03) CPU when running parallel calculations