subroutineでの計算

Question

fcの値がコンパイルしても０という値しか出てきません。
どう改善すればいいでしょうか？

C2345678
C
E = 2.05E+5
F = 235.0
c
CLAMBDA = SQRT(3.143.14E/0.6/F)
write(,) 'CLAMBDA=',CLAMBDA
c
call wx( DLAMBDA, RLAMBDA, r, fc)
c
stop
end
c************************************************************
subroutine wx( DLAMBDA, RLAMBDA, r, fc)
DLAMBDA = 1.0
RLAMBDA = 0.0
c
do i = 1, 250
RLAMBDA = RLAMBDA + DLAMBDA
c
r = RLAMBDA/CLAMBDA
if ( r.le.1.0 ) then
fc = (1.0-rr2.0/5.0)/((3.0/2.0)+rr2.0/3.0)F
else
fc = 18.0/65.0/(rr)F
end if
c
write(,*) RLAMBDA, fc
c
end do
return
end

Answer 1

メインプログラムで変数 r に値が与えられておらず、FORTRAN は初期値を与えないのでゴミが入っている。たまたま今は大きな値になっていて、サブルーチンの if 文では else の方へ行き、分母で大きな r で割るので 0 に近い値が出ているのだと思います。

FORTRAN77 に近い形式で書いていますが、もう少し Fortran90/95 に寄せて、implicit none を付けて変数は全部宣言するようにすれば、この様なミスは減ると思います。余力があれば更に一歩進んで、サブルーチンを module の中にしまうと引数のチェックなどが厳しくなっていいことあるかと思います。。

Comments

[kawasann]

回答ありがとうございます。
再度質問なのですが、メインにはどのように書いた方がいいのでしょうか？

[curehoney]

すみません。私の回答が間違っていました。
結果的に r がおかしいのですが、r に変な値を渡しているからではなく、そもそもサブルーチン wx に CLAMBDA を渡していないためでした。

いまは CLAMBDA が未定義のまま使われているので結果的に r が変な値となって昨日の回答へ続きます。

メインルーチンでやるべきことは、サブルーチンに CLAMBDA を渡すことです。それはサブルーチンの引数に CLAMBDA を加えることとセットでやります。

[kawasann]

サブルーチンの引数にCLAMDAを加えたのですがうまく機能しません。どうすればいいですか？

[curehoney]

質問が漠然としているので、漠然としか答えられませんが、慎重に一つ一つ確かめるしかないでしょう。