fortran将txt文件中的一列数据赋值给一维数组

yangpeinian · 发表于 2020-5-4 16:00:27

vvt 发表于 2020-5-4 15:09
我试了一下，读出来的数据画图如你给出的图片2一样。

那么这张图，是如何得到的呢？

首先表示感谢您的有效回复。这个是我把文件中的数据赋值给一个一维数组，然后让这个数组是声场计算中的声源。计算后提取声源处的波形，就得到是第一幅图的波形。图1和图2的起跳点大概一致，有没有可能是数据计算的有效位数问题？但是我定义的是real(8),精度应该够吧

yangpeinian · 发表于 2020-5-4 16:02:08

vvt 发表于 2020-5-4 15:09
我试了一下，读出来的数据画图如你给出的图片2一样。

那么这张图，是如何得到的呢？

所以可以肯定的是，数据赋值应该没有问题。但是在计算过程中可能有什么有效位数的问题。

vvt · 发表于 2020-5-4 17:09:29

这需要你提供计算过程的代码，相关逻辑对应的公式，才能判断。

yangpeinian · 发表于 2020-5-4 18:40:26

本帖最后由 yangpeinian 于 2020-5-15 10:45 编辑

vvt 发表于 2020-5-4 17:09
这需要你提供计算过程的代码，相关逻辑对应的公式，才能判断。

麻烦请指教！不胜感激！

yangpeinian · 发表于 2020-5-4 18:48:23

vvt 发表于 2020-5-4 17:09
这需要你提供计算过程的代码，相关逻辑对应的公式，才能判断。

程序是在原来有限差分的基础上改的。之前的程序声源直接是一个声源函数，现在的是把声源波形的txt文件数据赋值，设为新的声源。程序可能不好看懂，模型就是包括轴向和径向两个方向的声场计算，需要得到不同时刻不同位置的跟声场相关的物理量。总之就是时间循环里嵌套了一个轴向网格循环，里面再嵌套了一个径向网格循环。