空间插值 程序代码有问题

fcode

请详细描述问题，而不是只给代码。

试试下面的代码？

[Fortran] 纯文本查看 复制代码

C     空间插值
C     --------------------------------------------------
      DOWN = 0.0
      UP = 0.0
      MU = 2
C
      DO I=1,size(P_time,1)
        Distance = sqrt(dot_product(coords(1:3),P_time(I,1:3)))
        IF(Distance < 1.e-6)THEN
          DOWN = 1
          UP = P_time(I,4)
          EXIT
        ELSE
          DOWN = DOWN + 1 / Distance**MU
          UP = UP + P_time(I,4) / Distance**MU
        END IF
      END DO
C
      F = UP / down