We present the results of a numerical simulation of the dependence of the J
osephson critical current on the parallel magnetic field in the presence of
a large number of Abrikosov 'monopole' vortices for the case of a circular
junction. The influence of the junction edge on the distribution of the ma
gnetic field of each vortex was taken into account by introducing image mon
opole vortices. Our simulations show that the dependence of the Josephson c
ritical current on the parallel magnetic field (the I-c vs H-par dependence
) critically depends on the existence of a vortex located near the center o
f the junction. The simulation data are compared with experiments made on c
ircular Nb/Al - AlOx/Nb junctions. A good agreement has been found between
calculated and I-c vs H-par dependences measured after the application of a
perpendicular magnetic field.