这段程序是磨损过程中，计算磨损的深度，由于磨损，部件的网格节点的坐标会发生变化，但是这段程序我不是很明白是怎么回事，有没有老师帮忙讲解一下

IF(NODE.GE.18383) THEN
  LOCNUM1 = NODE + 344
  CALL GETVRN(LOCNUM1,'COORD',ARRAY,JRCD,JGVBLOCK,LTRN)
  current_x1=ARRAY(1)
  current_y1=ARRAY(2)
  current_z1=ARRAY(3)
  dist_n1=SQRT((current_x1-current_x)**2+(current_y1-current_y)**2+
     &     (current_z1-current_z)**2)
! WRITE(22,103),LOCNUM1
  WVGLOBAL(1) = w_dist*(current_x1-current_x)/dist_n1
  WVGLOBAL(2) = w_dist*(current_y1-current_y)/dist_n1
  WVGLOBAL(3) = w_dist*(current_z1-current_z)/dist_n1
  do k1=1,NDIM
  WVLOCAL(k1) = 0
  do k2=1,NDIM
  WVLOCAL(k1)=WVLOCAL(k1)+WVGLOBAL(k2)*ALOCAL(k2,k1)
  end do
  end do
  do k1=1,NDIM
  ULOCAL(k1) = ULOCAL(k1) + WVLOCAL(k1)
  end do