2006年度計算力学第10回プログラミング実習（2） 2006.06.12 ■ レポート提出について提出するもの • プログラムリスト • 結果のグラフ（Excelを使うとラク） • 課題に対する解答 Wordなどで一つのファイルにまとめ，以下のWebサイトから提出する（６月３０日締め切り） https://www.e-skillup.jp/shibaura-it/ ■ プログラムの流れ変数の設定，型宣言必要な変数は書き出しておくクーラン数の設定初期条件の設定クーラン数は変更できるように差分法による計算配列を使って繰り返し計算計算結果の出力まずは最終結果を！ ■ 繰り返し処理（DO 文）のプログラム例【問題】Ｓ＝１＋２＋３＋・・・＋５０を計算せよ program goukei integer sum, i sum=0 do 10 i=1,50,1 sum=sum+i 10 continue write(*,*) ’S=‘,sum stop end do N i=i1,i2,i3 文１文２・・ N continue i を i1 から i2 まで i3 ずつ増やしながら do と continue にはさまれた文を繰り返し実行する ■ 配列を使ったプログラムの例【問題】３個のデータ a1，a2，a3 を読み込んで，その平均値を求めよ． a1，a2，a3 はいずれも実数とする． ■ 配列を使ったプログラムの例配列を使うと・・・ ▲ 平均値を求めるプログラム ■ 配列について（補足１）添え字が１つの変数 a1, a2, a3, ・・ , ai , ・・に対して a(i), i=1,2,3, ・・・・１次元配列添え字が２つの変数 f11, f21, f31, ・・ , fi j, ・・に対して f(i,j), i=1,2,3, ・・・・ j=1,2,3, ・・・・２次元配列 ■ 配列について（補足２）配列の番号を０から使いたければ， a(0:100), f(0:100,0:200) のように，上限と下限を指定する． ■ 風上差分法のプログラム n：時間ステップ， j：空間刻み， v：クーラン数として， do 10 n=1, nmax do 20 j=1, jmax f(j)=f(j)-v*(f(j)-f(j-1)) 20 continue 10 continue このプログラムで正しく計算できるか？ ■ 境界条件について ex. 風上差分 j=1 j=1 j=1 j=0 計算領域の左端（ｘ軸上の最初の格子点）における f の値は計算しない． → 境界条件として与えてしまう． ※初期分布として設定した f1 の値を使えばよい． ■ 差分スキームの補足 ▼ Friedrichs-Lax 法 ▼ Fromm 法