function [gstop,isterminal,direction]=long_jump_landing(t,y,p)

gstop=y(4);     %Measure when y=0--jumper has landed
isterminal=1;   %stop the numerical solving
direction=-1;   %jumper is falling