This model uses new coupling correction alternation direction implicit iterative solution, parallel calculation program design, alternative direction grid delineation and parallel data transference.