[1.1] 授業の目標
この授業では，プログラミング言語の習得を目標とします．習得すると何ができるようになるのでしょうか？ 大きくわけると二つのことができるようになります．
- 数値シミュレーション
- データ解析
私が最近行った数値シミュレーションの例と研究の背景を下に示します．
研究の背景：地球に落下してくる隕石のうち8割はコンドライト隕石と呼ばれるタイプです．この隕石の中にはコンドリュールとよばれる大きさ0.1-2mm程度の球状の粒子が多数含まれています．形状が球であることから，何らかの加熱イベントによって鉱物が溶融し，表面張力によって球状になったと考えられています．体積の8割以上がコンドリュールである隕石も多く，また形成年代は太陽系形成以後100-300万年程度と測定されています．これらのことから，太陽系形成時において発生した加熱イベントによってコンドリュールは形成されたはずです．しかしその加熱イベントが何であったかは未だに謎です．
太陽系の形成時は，微惑星とよばれる小天体の相互衝突が頻繁におこっていました．衝突すると温度が上昇するので，シリケイトが溶融する可能性があります．そこで微惑星の相互衝突を数値シミュレーションし，どの程度溶融するのかを調べました．それが上の動画です．
衝突のシミュレーション結果を元に，溶融したシリケイトがどの大きさまで分裂するかを流体力学シミュレーションにより求めました．同時に温度低下の速度も算出しました．実験室ではそれぞれ0.1-2mm および 10 - 1000 度/時 と求められています．シミュレーションの結果は上図に示されており，幅広くパラメータを変化させても計測値とちょうどおさまる結果となりました．
[1.2] 今日の目標
- この授業の意義を知ること
- エディタの使い方を覚える
- プログラムの実行のしかたを覚える
- メールの出し方を覚える
[2.1] テキストエディタ
プログラムは「エディタ」とよばれるアプリケーションで作成します。エディタとは、文字列のみ作成できるもので、ワード等のワープロとは根本的に異なります。メモ帳（windows)，テキストエディット（Mac）が標準で付属しているエディタです．さまざまなフリーのエディタがありますので好みに応じて使ってください．標準のものでも十分です．以下の解説は，サテライトラボのiMacにインストールされている emacs エディタの解説です．
画面下のアイコンが並んでいる中に，バックがうすむらさき色の上にペンの絵があるアイコンがあります．一度クリックすると起動します．
[2.1.1] Emacsの基本操作
以下で C-x とは CTRLキーを押しながら「x」キーを押すことを意味します. 次の3つは最低限覚えましょう！
- ファイルを開く: C-x C-f ファイル名入力
- ファイルをセーブする: C-x C-s
- emacsの終了: C-x C-c
- コマンドを中断する（訳がわからなくなった時の回復手
段）: C-g
[2.1.2] 練習
ここで練習をしてみましょう. test.txtというファイルを開いてください. 今の場合, test.txtはまだ存在しません. その場合は新たにファイルが作成されます. 開くことができたら, ローマ字，かな漢字を入力してみましょう. 一通り書いたらセーブしてみましょう.
作成したプログラムは，「コンパイル」という作業をすることで実際に動くようになります．この作業は「ターミナル」という画面で行います．この画面での操作は，マウスでクリックといった直感的なものではなく，コマンドを入力することで実行します．windowsの方は gfortran を起動してください．macの方は「ターミナル」を起動してください．
windowsにおけるコマンド
| cd フォルダ名 |
フォルダ間の移動 |
| cd Mail |
Mail フォルダに移動 |
| cd ../ |
一つ上の階層のフォルダに移動 |
| mkdir フォルダ名 |
フォルダの作成 |
| mkdir hoge |
hogeフォルダを作成 |
| del ファイル名 |
ファイルの削除 |
| del hogehoge |
ファイルhogehogeを削除 |
MacOSにおけるコマンド
| cd フォルダ名 |
フォルダ間の移動 |
| cd Mail |
Mail フォルダに移動 |
| cd ../ |
一つ上の階層のフォルダに移動 |
| mkdir フォルダ名 |
フォルダの作成 |
| mkdir hoge |
hogeフォルダを作成 |
| rm ファイル名 |
ファイルの削除 |
| rm hogehoge |
ファイルhogehogeを削除 |
図1: フォルダの階層構造（macの場合）
[4.1]ディレクトリを作る
毎週必ず行なう手順を示します.
- 日付を名前にしたディレクトリを作る：mkdir 201008
- そのディレクトリに移動：cd 201008
- 毎回ディレクトリを作成し, その中にプログラムを保存すること.
こうしておくと後で見返すことが容易にできます.
[4.2]プログラムの作成, コンパイル
次のプログラムを書いてみましょう．example1.f90という名前でさきほど作成したフォルダの中に保存してください. もし間違えたら再度 201008 フォルダに保存してください．保存したら，dir (ls：今後，特記しなかぎりwindowsのコマンドをかっこ無しで明記し，かっこの中にmac でのコマンドを書きます) コマンドでファイルが存在していることを確認してください．
program example1 ! プログラムの開始
implicit none ! 暗黙の型宣言を使わない
double precision a,b,c ! 倍精度型実数a,b,cを宣言
a = 1.D0 ! aに倍精度の1を代入
b = 2.D0 ! bに倍精度の2を代入
c = a + b ! cにaとbの和を代入
write(*,*) a,b,c ! a,b,cに格納されている数値をモニタに出力
end program example1 ! プログラム修了
まず覚えてほしいフォートランのきまり：
- 上の行から下の行へと順番に処理される．
- 1-132文字目までで記述する．
- 大文字と小文字の区別は無い．WRITEとwrite，a と Aは同じ意味となる．
- 次の行へ継続するためには，その行の末尾に&をつける．
- !の後はコメント．メモに使う. できるだけコメントは書き込むこと. 後で役にたつ.
- プログラムはprogram プログラム名 で開始する．
- implicit none は，暗黙の型宣言を抑制することを示す．
- double precision は, 続く変数（数値が格納される文字）が倍精度であることを示す.
- 1.D0は倍精度実数の1.0を表す. 今後定数を書くときには倍精度実数にすること．
- +, -, = 等の演算子の前後にはスペースを入れると可読性が増す．
- write(*,*) a,b,c はa,b,cの数値を標準出力(モニタ)に出す．
- プログラムの最後は end program プログラム名 とする．
さて作成したプログラムはこのままでは計算機は理解できません. 実行形式（バイナリ）に変換する必要があります. この操作をコンパイルと呼びます. コンパイルはエディタ上ではできません. ターミナル上で入力します.
gfortran example1.f90
で行われます．すると a.exe (a.out) と言うファイルが作成されます．以下のコマンドで実行してみましょう.
a.exe (./a.out)
- a + bをa * b, a / bに変更して実行してみましょう.
[4.3]定数, 変数の型
FORTRANでは定数（1.D0などの具体的な数）および変数（a, bなど定数を格納できる文字）には型が定められています. 常に型を意識して変数および定数を使用してください. 当面登場する型は次の三つです. そのうち使用するのは始めの二つとなります.
- 倍精度型実数型：実数を用いる時は常にこの型を使います.
変数の宣言の仕方：double precision a
定数の書き方：1.D0（1かける10のゼロ乗）, 0.124D-3（0.124かける10のマイナス3乗）などDを用いる.
- 整数型：後述する配列, DOループで主に用います.
変数の宣言の仕方：integer a
定数の書き方：1, -10など． 小数点を書かない.
- 単精度型実数型：
変数の宣言の仕方：real a
定数の書き方：1.E0, 0.124E-3などEを用いる.
注意：宣言文はプログラムの先頭に置いてください！
次のプログラムを実行してください.
program example2
implicit none
double precision a,b
real c
integer n
a = 2.1D0
b = 2.1
c = 2.1
n = 2.1
write(*,'(D24.16,1x,D24.16,1x,D24.16)') a,b,c
write(*,'(I3)') n
end program example2
モニタ上の表示はどうなったでしょうか？ 何桁目まで正しく表示されていますか？
nの表示はどうなっていますか？
再度注意点をまとめると：
- 用いる変数は必ず宣言する
- 定数にはD0（ディーゼロ）を必ずつける
以上二点をいつも守るように心がけて下さい.
[4.4]書式の指定
write(*,*) の右側*の部分には本来「書式指定子」と呼ばれるものが入る. 書式を指定すると,計算結果の出力を自分の都合に合わせてカスタマイズできる. 使用頻度が高い形式の例を二つ挙げておきます.
- Dn.m：全体でn文字をとり, そのうちm文字を小数部にあてて実数を出力する. n は m+8以上. 8文字が別に使われるため．
- nX: n個の空白を空ける
- In（大文字アイ）: 全体でn文字をとり, 右に詰めて整数を出力する
以上の記号を, '( )'で囲んで記述してください. 複数ある場合は ,（コンマ）で続けて下さい.
[4.5]read文とwrite文
[4.5.1]モニタに出力, キーボードで入力
前節では画面に結果を出力させました．キーボードから入力もしたいと思います．入力は read文を使って行われます．先ほどのプログラムを
program example2
implicit none
double precision a,b,c,d
a = 1.D0
b = 2.D0
c = a + b
read(*,*) d
write(*,*) d
write(*,*) a,b,c
end program example2
と変更して実行してみよう．
- read文の入力として, キーボードから e をいれてみましょう．
- dの宣言をdouble precisionからintegerに書き替えて実行してみましょう.
[4.5.2]ファイルに入出力
計算させた結果は別のファイルに書き出すことがよくあります．そのためにはプログラムの宣言文（double precision 等）の次に
open(1,FILE='test.dat')
を追加してください. こうするとその後write(1,*)とすることによりtest.datに結果が出力されます．
- open(1,FILE='test.dat')を追加し, write(*,*) をwrite(1,*)に書き換えてプログラムを実行してください. ファイルtest.datが生成されていることを確認しましょう. write(1,*)の1はopen文で開いたファイル番号で，一方今まで使っていた*は標準出力（モニタ）に対応します．
[5.1] 今日の宿題
メールを送る．メールに含まれる内容は：
自己紹介文
数値計算を将来何に使うつもりか？
宿題の提出先: 城野(sirono@eps.nagoya-u.ac.jp) 濵本(hamamoto.miki@k.mbox.nagoya-u.ac.jp)
宿題のしめきり： 10月13日（火曜日）
[5.2] 提出先
レポートは毎回とも城野，濵本にメール
で送ること．二人以上に同時にメールを送るときは
sirono@eps.nagoya-u.ac.jp,hamamoto.miki@k.mbox.nagoya-u.ac.jpとコンマをつけて並べれば良い．
日程表へ戻る <<