Topic outline

  • 　

  • 時系列のベースマップ　psbasemap　（演習９－２）

    • 演習９－２


      • BAT-9-1.bat  をコピーして、作業フォルダ（c:\GMT_exe\Exercise-9\）にBAT-9-2.bat  として保存してください。
      • バッチファイルの、１～２行目（フォルダパス、ファイル名の定義）と最後の行（psconvert）だけを残して、ほかの行を削除してください。書き出しファイル名を、Fig-9-2  としてください。
      • 以下、デフォルト設定を変更を加えてください。
        （コピー＆ペーストすれば楽です）
      gmt set FORMAT_DATE_MAP=-o

      gmt 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  の命令文の中身を理解しましょう。



      できなかったら、以下をダウンロードしてください。


    • Mac用命令文（GMT_script_Ex9-2.sh）ダウンロード File 227 bytes application/x-sh Uploaded 26/11/20, 16:01


    • つぎは、xyDATA（時系列データ）を読み込み、このベースマップに色を塗ってゆきます。

  • 時系列のコンター図作成　（演習９－３）


    • 演習９－３

      • 以下のデータファイルをダウンロードして、作業フォルダに保存してください。

    • xyDATA-9-3.txt をダウンロード Folder

      注意：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

      くらいにカラーパレットの範囲を調整するとよいでしょう。






    • できなかったら、以下をダウンロードしてください。



    • 演習９－３　（つづき）


      先ほど作った図（ Fig-9-3 ）を見やすいように、

          • 左と下だけに、軸の数値注釈（annotation）をつけ、軸ラベルをつける
      （コンター図でプロットするデータ列の種類は、データフォルダの先頭行で確認してください）

          • 適した軸範囲に調整する
          • 他のデータ列のパラメタのカラーコンター図も作成する（書き出しファイル名を自由に付けてください）



  • コンターに数値注釈をつける　（演習９－４）


    • 演習９－４


      コンター（等値線）に、その値を注釈（annotation）をつけます。

      • BAT-9-3.bat  を同じ作業フォルダに、BAT-9-4.bat として別名保存してください。書き出しファイル名を、Fig-9-4  としてください。


      • 2番目の pscontour  にあるオプション  -A-  （数値注釈ない）を、-A+f15,22+a0+g255 -Gd2+r2c  に置き換えてください。
        （　-A+f15,22+a0+g255  と、-Gd2+r2c　の間にスペースが入っています）

      実行しましょう。




      pscontour オプションの -A  と  -G  の説明です。



      各自、見やすい図にアレンジしましょう。