プログラミング初心者です

Question

前提・実現したいこと

fortranで移流方程式の解を求めようとしています。

発生している問題・エラーメッセージ

解が表示されません

該当のソースコード

      dimension un(100),uf(100)
      
      open(unit=6,file='result.dat')

      dx=1000.*100.
      dt=60.
      imax=30
      nmax=200
     
      do 20 i=1,imax
       un(i)=0.
       uf(i)=0.
20    continue

      n=1
1     continue

      un(1)=un(2)
      un(imax)=0.

      do 30 i=2,imax-1
       uf(i)=un(i)-dt*c*(un(i)-un(i-1))/dx
30    continue

      do 40 i=1,imax
       un(i)=uf(i)
40    continue

      n=n+1
      
c      write(6,500)n,(un(i),i=1,imax)

      if(n.le.nmax)go to 1

500   format(i3,30f7.2)
      close(6)

      stop
      end

試したこと

ここに問題に対して試したことを記載してください。

補足情報（FW/ツールのバージョンなど）

ここにより詳細な情報を記載してください。

Answer 1

変数Cの初期値を与えましょう。
Cの初期値をゼロにするコンパイラだと

uf(i)=un(i)-dtc(un(i)-un(i-1))/dx

がいつもゼロになってしまいます。

Answer 2

c      write(6,500)n,(un(i),i=1,imax)

先頭のCは注釈行を意味しますから、Cを外しましょう。そうすれぱwrite文が働きます。
Cと同様に先頭の*も注釈を意味します。
なお無造作に削除して、wが6文字目に来ないように気をつけましょう。

Comments

[PT239]

write文が働きました。ありがとうございます。
しかし、結果の数字が全て0と表示されてしまいます。

[seastar3]

nも0ですか。
正確に出力結果の文字列を示してもらえますか。
デバッグ用の各変数の途中経過を調べるwrite文をくどいぐらい割り込ましてみましょう。