Fortran Coder

查看: 6170|回复: 1

[求助] 空间插值 程序代码有问题

[复制链接]

29

帖子

10

主题

0

精华

熟手

F 币
152 元
贡献
98 点
发表于 2021-9-13 20:11:06 | 显示全部楼层 |阅读模式
[Fortran] 纯文本查看 复制代码
C     空间插值
C     --------------------------------------------------
DOWN = 0.0
      UP = 0.0
      MU = 2
C
      DO I=1,10125
        Distance = (COORDS(1) - P_time(I,1))**2
     &           + (COORDS(2) - P_time(I,2))**2
     &           + (COORDS(3) - P_time(I,3))**2
        Distance = Distance**0.5
        IF(Distance == 0)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

1948

帖子

12

主题

5

精华

论坛跑堂

臭石头雪球

F 币
1298 元
贡献
547 点

美女勋章热心勋章星光勋章新人勋章贡献勋章管理勋章帅哥勋章爱心勋章规矩勋章元老勋章水王勋章

发表于 2021-9-13 21:40:46 | 显示全部楼层
请详细描述问题,而不是只给代码。

试试下面的代码?

[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
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-3-29 20:15

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表