Topへ戻る　

10：データ観察（時空間解析1） 各海域における海面温度の上昇傾向 引き続き、NOAAが配布するOISST ver2を用いる。 各グリッドでトレンドを計算するだけ。 つまり前回示したスクリプトがほぼ答えで、改変するだけ！ バイナリ出力は以下の通り。 open(80,file='slope.bin',form='unformatted', & access="direct",recl=360*180*4) write(80,rec=1) slope close(80)

GMTで傾き。つまり上記の変数slopeをマッピングする。 以下は9月のデータのみで計算したもの。

GMTスクリプト gmt begin SSTslope gmt xyz2grd slope.bin -R0.5/359.5/-89.5/89.5 -I1 -ZTLf -Gsst -di-999 gmt makecpt -Chaxby -T-0.5/0.5/0.01 -N gmt grdimage sst -R0.5/359.5/-89.5/89.5 -JR16c gmt grdcontour sst -C0.05 -A0.1+f10p -Gd8c -S4 gmt coast -B+t"Longterm trend of SST" -B -Gwhite -Wthinnest gmt colorbar -DjTR+o0c+w3c/0.4c -Bx -By+loC/decade -I -F+gwhite+p1p gmt end show