Question

我正在帮助一位有一些问题的老科学家，他有一个他发现的程序修改了大约20年。以前，并运行良好作为32位os \ 2可执行文件，但我需要它在DOS下运行！我不是程序员，而是一个好的硬件和软件男，所以我关于这个问题非常愚蠢，但这里去了我已经下载了6个不同的编译器watcom77，silverfrost ftn95，gfortran，g77和f80的2个版本。 WATCOM 说它是旧的程序，找到旧的编译器，silverfrost打开它，调试等。但是正在将所有子程序从“真实”改为“复杂”和副本 g77似乎安装完美（库链接等），但甚至不编译 test.f程序。我的问题是1;重新编译“按原样”或“升级”代码？ PROGRAM xconvlv INTEGER N，N2，M

PARAMETER  (N=2048,N2=2048,M=128)
INTEGER  i,isign
REAL data(n),respns(m),resp(n),ans(n2),t3(n),DUMMY
OPEN(UNIT=1, FILE='C:\QKBAS20\FDATA1.DAT')
DO 1 i=1,N
  READ(1,*) T3(i), data(i), DUMMY
continue
  CLOSE(UNIT-1)
do  12 i=1,N
  respns(i)=data(i)
  resp(i)=respns(i)
continue
isign=-1
call  convlv(data,N,resp,M,isign,ans)
  OPEN(UNIT=1,FILE='C:\QKBAS20\FDATA9.DAT')
  DO 14 i=1,N
  WRITE(1,*)  T3(i), ans(i)
   continue
END
SUBROUTINE CONVLV(data,n,respns,m,isign,ans)
INTEGER  isign,m,n,NMAX
REAL data(n),respns(n)
COMPLEX ans(n)
PARAMETER (NMAX=4096)
* uses realft, twofft
INTEGER i,no2
COMPLEX fft (NMAX)
do 11 i=1, (m-1)/2
  respns(n+1-i)=respns(m+1-i)
continue
do 12 i=(m+3)/2,n-(m-1)/2
  respns(i)=0.0
continue
call twofft (data,respns,fft,ans,n)
no2=n/2
do 13 i=1,no2+1
  if (isign.eq.1) then
    ans(i)=fft(i)*ans(i)/no2
  else if (isign.eq.-1) then
    if (abs(ans(i)) .eq.0.0) pause
    ans(i)=fft(i)/ans(i)/no2
  else
    pause 'no meaning for isign in convlv'
  endif
continue
ans(1)=cmplx(real (ans(1)),real (ans(no2+1)))
call realft(ans,n,-1)
return
END

SUBROUTINE realft(data,n,isign)
INTEGER isign,n
REAL data(n)
* uses four1
INTEGER i,i1,i2,i3,i4,n2p3
REAL c1,c2,hli,hir,h2i,h2r,wis,wrs
DOUBLE PRECISION theta,wi,wpi,wpr,wr,wtemp
theta=3.141592653589793d0/dble(n/2)
cl=0.5
if  (isign.eq.1) then
  c2=-0.5
  call four1(data,n/2,+1)
else
  c2=0.5
  theta=-theta
endif
 (etc.,etc., etc.)


SUBROUTINE twofft(data,data2,fft1,fft2,n)
INTEGER n
REAL data1(n,data2(n)
COMPLEX fft1(n), fft2(n)
* uses four1
INTEGER j,n2
COMPLEX h1,h2,c1,c2
c1=cmplx(0.5,0.0)
c2=cmplx(0.0,-0.5)
do 11 j=1,n
  fft1(j)=cmplx(data1(j),data2(j)
continue
call four1 (fft1,n,1)
fft2(1)=cmplx(aimag(fft1(1)),0.0)
fft1(1)=cmplx(real(fft1(1)),0.0)
n2=n+2
do 12 j=2,n/2+1
  h1=c1*(fft1(j)+conjg(fft1(n2-j)))
  h2=c2*(fft1(j)-conjg(fft1(n2-j)))
  fft1(j)=h1
  fft1(n2-j)=conjg(h1)
  fft2(j)=h2
  fft2(n2-j)=conjg(h2)
continue
return
END

SUBROUTINE four1(data,nn,isign)
INTEGER isign,nn
REAL data(2*nn)
INTEGER i,istep,j,m,mmax,n
REAL tempi,tempr
DOUBLE PRECISION theta, wi,wpi,wpr,wr,wtemp
n=2*nn
j=1
do 11 i=1,n,2
  if(j.gt.i)then
    tempr=data(j)
    tempi=data(j+1)
(etc.,etc.,etc.,)
continue
mmax=istep
goto 2
endif
return
END

有4个子程序，大约有3页代码，应该很多如果他们能够帮助我，那么更容易通过电子邮件发送给某人。我的电子邮件是 mastermikie4356@yahoo.com，或者如果有人能告诉我哪里可以“工作” 编译器可以重新编译吗？谢谢，谢谢，谢谢你对此的任何帮助！我得到的错误是; 1.在从另一个过程调用CONVLV时，第一个参数是一个类型 REAL（kind = 1），它现在是COMPLEX（kind = 1） 2.在从另一个程序调用REALFT时，... COMPLEX（kind = 1）现在是真实的（种类= 1） 3.在来自... COMPLEX（kind-1）的TWOFFT调用中它现在是REAL（kind = 1） 4.在之前的FOUR1调用中，第一个参数是REAL类型（kind = 1）它现在是一个复合体（种类= 1）。

Answer 1

在快速阅读该程序时，我没有看到任何重大问题。文件名被“硬编码”到打开的语句中 - 您可能需要更改这些文件名以匹配您将程序移植到的任何计算机上的位置。

对于FORTRAN 77，行号以外的源代码应从第7列开始。

如果您需要更具体的建议，请向我们显示错误消息。

我建议尝试使用一两个编译并找出代码所需的更改，而不是安装许多编译器，其中的想法是以某种方式使程序正常工作。

Answer 2

错误消息表明调用例程和子例程之间存在参数不匹配（发送类型REAL，期望类型为COMPLEX）。这可能是一个错误，或者更可能的是，某些非标准技巧适用于较旧的编译器，但不适用于较新的编译器。

你可以发布每个子程序的标题吗？（或整个程序的链接？）

重新编译为os \ 2编译的旧fortran 2/4 \ 66程序需要它在dos中运行

2 个答案: