Asymptotic approach to the theory of large (

) two-dimensional Josephson junction is developed. As a result, the usual dc sine-Gordon equation for the rapidly changing Josephson phase difference φ is reduced to the much more simple hydrodynamic-type equations

, curl

, for the slowly changing wave vector

(A is a simple function of

). The reduced equations are applied for analysis of a square-shaped Josephson junction and are solved using the rigorous boundary conditions. The obtained dependence of the junction critical current on the applied magnetic field is discussed and compared with the recent experimental data.