空间点到曲线距离计算fortran程序

!DOI=1,IM !DOJ=1,JM !r_min=10000.0 !DOk=2,IS !TEMP=(xs(k)-xs(k-1))/(ys(k)-ys(k-1)) !TEMPP=(ys(k)-ys(k-1))/(xs(k)-xs(k-1)) !x0=1.0/(TEMP+TEMPP)*(TEMP*x_heng(I,J)+TEMPP*xs(k)+y_heng(I,J)-ys(k)) !y0=