問題が生じたら、まず、こちらを確認してください。
質問と解決策を投稿してください。皆さんのギモンは、みんなのギモンかもしれません。
(北大構成員のみ。北大Googleアカウントでログイン)
先のコース(GMT演習9-1)のダウンロードファイルです。
Topic outline
時系列のベースマップ psbasemap (演習9-2)
演習9-2
- BAT-9-1.bat をコピーして、作業フォルダ(c:\GMT_exe\Exercise-9\)にBAT-9-2.bat として保存してください。
- バッチファイルの、1~2行目(フォルダパス、ファイル名の定義)と最後の行(psconvert)だけを残して、ほかの行を削除してください。書き出しファイル名を、Fig-9-2 としてください。
- 以下、デフォルト設定を変更を加えてください。
(コピー&ペーストすれば楽です)
gmt set FORMAT_DATE_MAP=-ogmt set FORMAT_TIME_PRIMARY_MAP=abbreviated
- 以下の命令文(psbasemap )を定義文のあとに加えて、実行してください。
gmt psbasemap -R2019-2-14T00:00:00/2019-8-30T00:00:00/0/95 -JX12/-8 -Bxa1O -Bya20f10 -P > %PathName%%FileName%.eps
※ -Bxa1O の ' O ' は、大文字の 'オー' です。
※ gmt set でデフォルトの変更が保存されているので、以後、その変更が反映されます。
※ 上の、gmt setの2行は、この順に記す必要があります。
できましたか?
まず、
psbasemap の命令文の中身を理解しましょう。
できなかったら、以下をダウンロードしてください。
つぎは、xyDATA(時系列データ)を読み込み、このベースマップに色を塗ってゆきます。
時系列のコンター図作成 (演習9-3)
演習9-3- 以下のデータファイルをダウンロードして、作業フォルダに保存してください。
注意:xyDATA9_3.txt はMac用です。
- BAT-9-1.bat をコピーして、作業フォルダに BAT-9-3.bat として保存してください。
- 書き出しファイル名を Fig-9-3、読み込みファイル名を xyDATA-9-3.txt としてください。
- BAT-9-2.bat の、psbasemap の命令文をコピーして、BAT-9-3.bat の、psbasemap の命令文と置き換えてください。
- psbasemap の命令文のオプションに、追記予告オプション -K を加えてください。
変更箇所を赤字、置き換え箇所を紫枠で記しました。
実行しましょう。
このような図が描けましたか?(makecpt のオプション -T の設定値が違うと違う表示になります)
このように色が単調で、コンター線が無い図になってしまいましたか?
カラーパレット作成の命令文: gmt makecpt -Crainbow -T0/70/5 -Do > %PathName%contour.cpt
数値範囲が0~70、幅が5 に指定されていますね。
この作図で使用したデータファイル(xyDATA-9-3.txt)を確認してください。
gmt pscontour %xyDATA% -i0,1,2 ~~ にて、データを読み込み、
X軸を0列目(時間)、Y軸を1列目(水深)、データを2列目(NO2+NO3)に指定しました。
データ列(i=2)の値を見てください。10未満が多く、最大でも17です。
gmt makecpt -Crainbow -T0/20/2 -Do > %PathName%contour.cpt
くらいにカラーパレットの範囲を調整するとよいでしょう。
できなかったら、以下をダウンロードしてください。
演習9-3 (つづき)
先ほど作った図( Fig-9-3 )を見やすいように、
- 左と下だけに、軸の数値注釈(annotation)をつけ、軸ラベルをつける
(コンター図でプロットするデータ列の種類は、データフォルダの先頭行で確認してください)- 適した軸範囲に調整する
- 他のデータ列のパラメタのカラーコンター図も作成する(書き出しファイル名を自由に付けてください)
コンターに数値注釈をつける (演習9-4)
演習9-4コンター(等値線)に、その値を注釈(annotation)をつけます。
- BAT-9-3.bat を同じ作業フォルダに、BAT-9-4.bat として別名保存してください。書き出しファイル名を、Fig-9-4 としてください。
- 2番目の pscontour にあるオプション -A- (数値注釈ない)を、-A+f15,22+a0+g255 -Gd2+r2c に置き換えてください。
実行しましょう。pscontour オプションの -A と -G の説明です。
各自、見やすい図にアレンジしましょう。