- 610 名前：デフォルトの名無しさん [2007/06/20(水) 11:28:28 ]
- 質問させていただきます。
以下のプログラムで、y1とy3は答えとともにオープンしてくれるのですが、y2のオープンしたファイルが何も表示されません。
サブルーチン文はy1y2y3ともに引数は表示してくれてます。何か問題がありましたら教えて下さい。
お願いします。
program IFFT
implicit real*8(a-h,o-z)
dimension y1(1000000),y2(1000000),y3(1000000)
complex*8 w1(1000000),w2(1000000),w3(1000000)
dimension H(1000000),f(1000000)
c
open(114,file='d.dat')
open(115,file='e.dat')
open(116,file='f.dat')
c
m=15.d0
dt=0.01
n=2.0**m
tp=n*dt
df=1.0/tp
c
call SPECFILT(w1,w2,w3)
do 21 j=n/2+2,n
w1(j)=conjg(w1(n-j+2))
w2(j)=conjg(w2(n-j+2))
w3(j)=conjg(w3(n-j+2))
21 continue
c
|
 |
