IMSL编译错误，This name has already been used as an external function name.

楚香饭

use     NUMERICAL
use     CONSTANT
USE     VARIABLE
USE     IMSL
这些模块里，已经有定义过 Sxyz 了，可能冲突

楚香饭

试试
use IMSL
改为
use IMSL , sxyz_1 => sxyz

可能 IMSL 里面有一个叫 sxyz 的函数，用来算三维向量的内积。所以给他改个名字。

楚香饭

参数类型不一致。ivprk 用于处理单精度，而 T 是双精度。
尝试使用 DIVPRK 函数

楚香饭

VARIABLE 模块中，
real*8  sk3(3,0:2,0:22)
而在  shouli.f90 第 50 行：
do i=1,20
do j=0,4
    sk3(1,i,j)=sXYZ(1,i,j+1)-sXYZ(1,i,j)  !// 此处 sk3 越界，i 从1 到 20，但第二维是 0：2
    sk3(2,i,j)=sXYZ(2,i,j+1)-sXYZ(2,i,j)

楚香饭

楼主的完整代码在 7 楼。

但是我不知道他之前的越界怎么改的，所以也不知道问题在哪儿。

sqrt domian 常见于负数开方。

楚香饭

依旧是则徐 发表于 2015-4-7 14:39
越界问题 是把variable模块中的sk3（3,0:2,0:22）改成了sk3（3,0：22,0:5） 改过之后，在运行就出现了sqr ...

sqrt domain 不一定就是负数，有时候本身就是 NaN 也会出错。
至少在我这里，发现很多容易出负数的地方。

例如 130.f90 文件
YITABOU3(i)=Amp*DCOS(WK*X3(i)-W1*T)&
                  -Amp**2*WK/(2*DSINH(2*WK*H)) &
这里的 WK 为 0，导致分母为0

wk 应该是在更早的代码计算：
wk=newton(a,func,dfunc)
这里计算结果为0

可能是由此导致的 NaN

正如 pasuka 所言，没有注释的代码比较难看懂。我想你学习一下如何 Debug 单步运行，或许会让你自己更容易检查错误。

本站的文章：http://debug.w.fcode.cn 会指导你如何学习 Debug