加速度波形を積分して速度・変位波形をつくる
地震観測によって得られた加速度記録(0.01秒間隔,データ数7100個)を積分して速度および変位波形を描きたいのですが,これをFORTRANによってプログラミングして描きたいと思ってますが,主プログラムをどう書けばよいのか分かりません.サブルーチンは分かりますが...
私が思うに,FORMATの設定がおかしいのではないのかと考えています.
下に読み込む加速度記録および主プログラムを載せていますので,これを参考に分かる方は教えて頂ければ幸いです.
加速度記録(既に生成されたデータ)は,テキストファイル形式で保存しています.
ちなみにFORTRANのバージョンは95です.
加速度記録↓
0.105
0.069
0.041
0.051
-0.028
-0.188
-0.256
-0.243
-0.259
-0.233
-0.129
-0.054
0.005
0.098
0.141
0.103
・・・
主プログラム
DIMENSION DDY(800),DY(800),Y(800)
C
OPEN(5,FILE='ACC01.txt',STATUS='OLD')
WRITE(5,501)
READ(5,501) DT,NN,(DDY(M),M=1,NN)
CALL IACC(DT,NN,DDY,DY,Y,800,DYMAX,YMAX)
STOP
C
501 FORMAT(T51,F10.0,I10/(8F10.0))
CLOSE(5)
END