■Program: BEM8QUDR.FOR■
このプログラムは、BEM8QUDQ.FOR をベースにして作成しました。追加した変数を下表にまとめてあります。

上の変数のMXWELLは、main program の PARAMETER 文で定義してあります。その他の変数は、Subroutine INPUT で読み込まれ、Subroutine WELL で計算に用いられています。Subroutine WELL は、境界積分式の右辺 を計算します。そして、Subroutine WELL は、Source point (ξ) 毎に、Subroutine MTXFORM から call されます。

Subroutine WELL では、まず、下図に示す様に、point source 毎に ξ とx_i の間の距離 r_i を計算します。そして、 Q_iG( r_i) を計算しています。

全ての point source についてQ_iG( r_i) を計算し、その合計を 変数 WELLS に格納しています。