A technique is presented for modeling two‐dimensional horizontal flow in a system of two infinite aquifers. The aquifers are separated by a thin impermeable layer and are hydraulically connected through openings in the separating bed. The upper aquifer is unconfined, and the flow in the lower aquifer is confined in some regions and unconfined in others. The model includes cross flow, inflow from rainfall, and flow toward creeks, canals, lake boundaries, and wells. Line sinks are used for incorporation of the creeks, canals, and lake boundaries. The method is validated against various exact solutions, and some examples of applications of the method are discussed.