Crack tip opening displacement of a crack lying parallel to a free surface is calculated by counting the number of dislocations emitted into the plastic zone from a crack tip. A discrete dislocation model was used to simulate the crack, while varying the strength of dislocations so as to satisfy the boundary condition. The result coincides numerically with the predictions made in a previous paper, in which the stress intensity factor appearing in a theory of bulk materials was replaced with the one which includes the surface correction.