這段程序是磨損過程中，計算磨損的深度，由于磨損，部件的網(wǎng)格節(jié)點的坐標會發(fā)生變化，但是這段程序我不是很明白是怎么回事，有沒有老師幫忙講解一下

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