北极振动とは.PPTVIP

①コメント ②定数 ③配列 ④入出力 ⑤文字の結合 ⑥DOループ ⑦IF文 ⑧副プログラム（サブルーチン） 相関係数とは？ プログラムでファイルに出力したデータをエクセルで描いてみよう。 同じようにB列、C列データの散布図を作成 散布図と相関の感覚的理解 直線からのズレが小さい 相関係数1とは？ 相関係数の数学的な意味 相関係数の計算では、この性質を使っている。 データの数が３個だったら？ 標本数（データ数）と相関係数の限界値 * 環境システム科学演習 ２回目 　ＦＯＲＴＲＡＮ　文法について ①コメント　②定数　③配列　④入出力文 　⑤文字の結合　⑥IF文　⑦IF文　 ⑧副プログラム（サブルーチン） その他???? 演習：相関 北極振動指数と日本各地の気温の関係を調べる 2008.12.2 北極海は国内旅行程度 の範囲にある 日本の気候は バランスで決まる 海 偏西風 陸 海 陸 偏西風 高温?低圧 低温?高圧 季節風 海 偏西風 高温?低圧 低温?高圧 季節風 陸 出典：探検しようみんなの地球 http://www.mofa.go.jp/mofaj/gaiko/oda/sanka/kyouiku/kaihatsu/chikyu/statistics/nature/graph_02.html 温暖 寒冷 寒冷 温暖 北極振動のパターン 白 黒 黒 灰 灰 北極振動とは？ 「北極海で気圧が下がれば、北大西洋、北太平洋で気圧が上がり」、「北極海で気圧が上がれば、北大西洋、北太平洋で気圧が下がる」というパターン。北極海で低気圧偏差のとき??＞北極振動指数は正北極海で高気圧偏差のとき??＞北極振動指数は負 Cから始まる行は、何もしない ?コメント行 ！以降は、何もしない ?コメント PARAMETER は定数を与える CHARACTER は文字であることを宣言 REAL　なら実数であることを宣言 INTEGERなら整数であることを宣言 暗黙の了解について： 何も宣言しないならば、 　　A-H,O-Zは実数 　　I-Nは整数 になる DIMENSION ?配列であることを宣言 装置（この場合はファイル）をOPENする ?ファイルを読み込んだり、 ファイルに書き込んだりする 暗黙の了解について： 5番はキーボード 6番はCRT（画面） READ(*,*) の最初の*はキーボード WRITE(*,*) の最初の*は画面 つまり、 READ(*,*)はREAD(5,*) と同じ WRITE(*,*) はWRITE(6,*)と同じ 5番、6番以外はファイルを指定する CHARACTERで宣言しないときは、 ’　’ で囲めば文字変数になる FNは3文字の文字であると宣言した。 ’A0’と’.TXT’は数値ではなく 文字「A0」と「TXT」である。 ’A0’ //FN// ’.TXT’ は、FNとして、 　　ATS を指定すると 「A0ATS.TXT」 という文字を意味する。 内側： Jが1～Nまでを繰返し行うという意味です。 上の表の矢印の順に作業を行います。 ここでは、11番のファイル「AO.TXT」から順番に1からM（＝12）まで読み込むことを意味します。 Jについて2つおき行いたい場合には、 　　J=1,M,2 とします。 外側： 外側も同じです。上の表の点線にのように作業します。 IFのブロック 判定を行う、IF文です。 ここでは、 もし、Jが11以上ならば、 IF(J.GE.11) THENと END IFに挟まれた部分の計算を行います。 実際には、Jが11よりも大きな値を取るのは、J=11かJ=12しかありません。 つまり、11月と12月だけ考えますよというものです。 ここでは、X(I)を計算しているものを説明します。 I=1（1950年）から始まります。 J＝11のとき初めて、IFブロックに入ります。 そして、 AOTMPにAO(1,11）の値を入れますよとしています。 そして、 配列X(1)にAOTMPとX(1)を足したものを与えます。 この演算を行う前には、X(1)は、まだ何も数値を与えていないので、X(1)は0です。 演算後、次に、J=12になります。これも、11以上ですから、演算を行います。 このとき、X(1)には既に11月の値が入っていることに注意してださい。 そして、AOTMPには、AO(1,12)の値を入れます。 X(1)=AOTMP+X(1) とは、AO(1,11)とAO(1,12)を足したものを意味しているのです。 Y(I)も基本は同じです。 11月と12月の和を計算したので、その0.5倍しています。これは11-12月の平均を求めているのです。 装置番号6