GMT on Mac

JAGURS のファイルが GMT グリッドファイルということなので，いよいよちゃんと GMT を使うべく手元にインストール．MacPorts にあるようなのでそのまま使います．新しいの入れておけば良いだろうと思って， gmt5 をインストールしたけど， grdinfo コマンドが無い．．．

        % port contents gmt5 | grep grdinfo
          /opt/local/share/doc/gmt/html/_sources/grdinfo.rst.txt
          /opt/local/share/doc/gmt/html/grdinfo.html
          /opt/local/share/doc/gmt/man/man1/grdinfo.1.gz

じゃぁ，と思って gmt4 をインストールしたら，変わったディレクトリにコマンドが入ってました．

        % port contents gmt4 | grep grdinfo
          /opt/local/lib/gmt4/bin/grdinfo
          /opt/local/share/doc/gmt4/html/man/grdinfo.html
          /opt/local/share/doc/gmt4/man/man1/grdinfo.1

なので，最後の仕上げに
        export PATH=$PATH:/opt/local/lib/gmt4/bin
を追加して完成．

第一歩として試したかった，JAGURS 付属の海底地形データの情報を取得．zmax を使って，計算ステップ間隔を決定します．
% grdinfo bathy.SD01.grd
bathy.SD01.grd: Title:
bathy.SD01.grd: Command: grdsample -V -R140/147:30/34/43 -I18c /home/babat/data/bathymetry_thk/bathy.thk.27s.grd -Gbathy.whole.grd=cf
bathy.SD01.grd: Remark:
bathy.SD01.grd: Gridline node registration used
bathy.SD01.grd: Grid file format: cf (# 10) GMT netCDF format (float) (deprecated)
bathy.SD01.grd: x_min: 140 x_max: 147.5 x_inc: 0.005 name: x nx: 1501
bathy.SD01.grd: y_min: 34 y_max: 43 y_inc: 0.005 name: latitude [degrees_north] ny: 1801
bathy.SD01.grd: z_min: -2096.24316406 z_max: 9788.53125 name: z
bathy.SD01.grd: scale_factor: 1 add_offset: 0