GMT 海洋データの時系列鉛直断面図②（演習９-２,３,４）
Rangkaan topik
-
-
問題が生じたら、まず、こちらを確認してください。
-
質問と解決策を投稿してください。皆さんのギモンは、みんなのギモンかもしれません。
（北大構成員のみ。北大Googleアカウントでログイン）
-
-
-
演習９－２
- BAT-9-1.bat をコピーして、作業フォルダ（c:\GMT_exe\Exercise-9\）にBAT-9-2.bat として保存してください。
- バッチファイルの、１～２行目（フォルダパス、ファイル名の定義）と最後の行（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の２行は、この順に記す必要があります。
できましたか？
まず、
psbasemap の命令文の中身を理解しましょう。
できなかったら、以下をダウンロードしてください。
-
-
-
-
注意：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
くらいにカラーパレットの範囲を調整するとよいでしょう。
-
-
演習９－４コンター（等値線）に、その値を注釈（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 の説明です。
各自、見やすい図にアレンジしましょう。
-
-