お二人とも、アドバイスありがとうございます。

>>970
1から書いています。
昨日、パラメータと基本物理量以外は全てローカルに直すことに成功しました。
あとは、基本物理量をローカルにするだけです。

構造体は使ったことがないのですが、
おっしゃっていることは、このような感じのことでしょうか？
3つの流体変数を時間発展させていくサンプルプログラムです。(若干、端折っています
構造体の要素の代入文のせいで、サブルーチンが
長くなってしまうのは、仕方がないですよね？
まだ、削れるところありますか？