Takuya Miyashita

Start:
*GMT -J オプション [#q288969a] #contents **はじめに [#jac8759d] 共通オプション -J は projection を指し，描画をするコマン... 図の種類，形状，スケールがこのオプションによって決定され... J の次に必ずアルファベット1〜2文字が入り，描く図法の種類... **projection指定子の大文字・小文字の違い [#d57d020d] J に続く英字で描く図法が決まるが，その英字が大文字か小文... メルカトル図法を例にすると，おおよそ下記のような書き方が... - -JM<width> - -Jm1:<scale> - -Jm</degree> ここで，<width> は出力する図の幅，<scale> は比，</degree>... 大文字 M は <width> インチなどの出力サイズを指定する方法，~ 小文字 m は実スケールとの比を指定する(出力される紙面にと... この原則は -J オプションのその他(メルカトル図法を表す M, ... ~ それぞれの使用例と出力結果は下記の通り．~ このページ3つの画像の表示サイズを統一しているが，実際の出... #codeprettify(lang-bash){{ #!/bin/bash }} #codeprettify(lang-bash){{ # filename output="JMm" proj1=M10 region1=125/150/25/50 axes1="-Ba5f5g5 -B+t-JM<wid>,-R125/50/25/50" proj2=m1:2.75e7 region2=125/150/25/50 axes2="-Ba5f5g5 -B+t-Jm1:<scale>,-R125/150/25/50" proj3=m0.4 region3=125/150/25/50 axes3="-Ba5f5g5 -B+t-Jm<wid_perdeg>,-R125/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt coast -J$proj3 -R$region3 $axes3 -W -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ** projection, 各図法の記述例 [#d78dd727] 以降では，projection ごとに pscoast または psbasemap を用... 描ける図法の種類が多いので，一部に限定する．~ また，煩雑さを避けるため大文字，-J<uppercase> で統一し，... 図法によってオプション -R や -B も変更する必要があるが，... ~ **-JM, Mercator, メルカトル図法 [#ufe55c9a] 見慣れたメルカトル図法．記法は下記の通り．~ - -JM<width> - -JM<lon0>/<width> - -JM<lon0>/<lat0>/<width> #codeprettify(lang-bash){{ # filename output="JM" proj=M10 region=125/150/25/50 axes="-Ba5f5g5 -B+t-JM<wid>,-R125/50/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JB, Albers conic equal area projection [#e1384514] Conic projection の一つ． - -JB<lon0>/<lat0>/<lat1>/<lat2>/<width> ここで，<lat1>,<lat2> は Alberts projection の two standa... 参考: [[Wikipedia>https://en.wikipedia.org/wiki/Albers_pr... #codeprettify(lang-bash){{ # output file output="JB" # options proj=B135/35/30/45/12 region=120/150/25/50 axes="-Ba5f5g5 -B+t-JB135/35/30/45/<wid>,-R120/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P -V > "${output}... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JU, UTM, ユニバーサル横メルカトル図法 [#j529cb9a] 数値計算でよく UTM 座標を用いるため．計算結果に対して -JU... - -JU<zone>/<width> ここで，<zone> は UTM のゾーン(1〜60の整数)． #codeprettify(lang-bash){{ # output file output="JU" # options proj=U53/10 region=125/150/25/50 axes="-Ba5f5g5 -B+t-JU<utmzone>/<wid>,-R125/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JC, Cylindrical projection, 円筒図法 [#hc899061] UTM に似ている．UTM でいう zone の基準点がないので，<lon0... - -JC<lon0>/<lat0>/<width> #codeprettify(lang-bash){{ # output file output="JC" # options proj=C135/35/10 region=125/150/25/50 axes="-Ba5f5g5 -B+t-JC135/35/<wid>,-R125/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JA, Lambert Azimuthal projection, ランベルト正積方位図... ランベルト図法．図の特性上,描ける最大範囲は半球．~ オプション -R の標記が上記のタイプの指定方法と異なるので... - -JA<lon0>/<lat0>/<width> -Rg - -JA<lon0>/<lat0>/<width> -R<lonSW>/<latSW>/<lonNE>/<lat... Rgは<lon0>,<lat0> を中心とする半球を描く．~ "-R<...>r" は南西，北東の端を指定する．最後の【r】をつけ... g も r も付けず，R<lonmin>/<lonmax>/<latmin>/<latmax> の... #codeprettify(lang-bash){{ # output file output="JARgARr" # options proj1=A136/35/10 region1=g axes1="-Ba15f15g15 -B+t-JA135/35/<wid>,-Rg" # options proj2=A120/30/10 region2=90/5/180/65r axes2="-Ba15f15g15 -B+t-JA120/30/<wid>,-R90/5/180/65r" # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JG, Azimuthal orthographic projection, 正射図法 [#db01... 地球表面を正射影して描かれる地図．オプションの表記方法は ... - -JG<lon0>/<lat0>/<width> -Rg さらに詳細な指定もできるが，ここでは省略． #codeprettify(lang-bash){{ # output file output="JG" # options proj=G135/35/10 region=g axes="-Bf30g30 -B+t-JG135/35/<wid>,-Rg" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JX, Non-geographical projection (linear) [#s1cb788c] 地図でない X-Y 座標には -JX を使う． - -JX<width> - -JX<width>/<height> JX<width> の場合は，<height> は <width> と同じ値．~ 範囲指定 -R も基本的には -JM などの場合と同じ． #codeprettify(lang-bash){{ # output file output="JXwXwh" # options proj1=X10 region1=0/100/0/50 axes1="-Bxa20g20 -By10g10 -BSWne+t-JX<wid>,-R0/100/0/50" # options proj2=X10/5 region2=0/100/0/50 axes2="-Bxa20g20 -By10g10 -BSWne+t-JX<wid>/<hei>,-R0/100/... # gmt6 gmt begin $output png gmt basemap -J$proj1 -R$region1 $axes1 gmt basemap -J$proj2 -R$region2 $axes2 -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JXl, Non-geographical projection (log-log, semilog), ... 対数軸にしたい場合は，横軸なら <width> の後に，縦軸なら <... - -JX<width>l/<height> : 片対数 X - -JX<width>/<height>l : 片対数 Y - -JX<width>l/<height>l : 両対数 ここでは B オプションに記述する数値の意味が他と異なる点に... a, f, g の値が大きいほど間隔が短くなる．p をつけることで ... #codeprettify(lang-bash){{ # output file output="JXl" # options proj=X10l/5 region=1e0/1e5/0/50 axes="-Bxa1pg3p -By10g10 -BSWne+t-JX<wid>l/<hei>,-R1e0/1e... # gmt6 gmt begin $output png gmt basemap -J$proj -R$region $axes gmt end # gmt5 gmt psbasemap -J$proj -R$region $axes -P > "${output}_gmt... gmt psconvert -Tg -A "${output}_gmt5.ps" }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JX,-JZ, Non-geographical projection (3-dimensional) [#... 3D図は，-JX と -JZ で2回 J オプションを使う． - -JX<width> -JZ<height> - -JX<width>/<length> -JZ<height> ここで，<length> は水平平面で <width> と垂直な方向の長さ．~ 回転すると入力した素直な(?)長さにならないので，最終的に出... #codeprettify(lang-bash){{ # output file output="JXJZ.ps" # output file output="JXJZ" # options proj=X10 projZ=Z10 pers=-120/25 region=0/10/0/10/0/10 axes="-Bxa2g2 -Bya2g2 -Bza2f2 -BSWZ+b+t-JX<wid>,-JZ<hei>,... # gmt6 gmt begin $output png gmt basemap -J$proj -R$region $axes -J$projZ -p$pers gmt end # gmt5 gmt psbasemap -J$proj -R$region $axes -P -J$projZ -p$pers... gmt psconvert -Tg -A "${output}_gmt5.ps" }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JX, -Rd or -Rg, 地球全体図 [#m9a91360] メルカトル図法 -JM では，R オプションで極の -90,90 を含め... JX の projection で，地球全体を描く方法がある．-Rg または... - -JX<width> -Rd - -JX<width> -Rd<lonmin>/<lonmax>/<latmin>/<latmax> - -JX<width> -Rg - -JX<width> -Rg<lonmin>/<lonmax>/<latmin>/<latmax> 経度方向の中心は，Rd で本初子午線，Rg で日付変更線になる．~ それ以外の経度を中心にしたい場合は -R で中心に来るように ... #codeprettify(lang-bash){{ # output file # output file output="JXRd_JXRg" # options proj1=X10/5 region1=d axes1="-Bxa60f60g60 -Bya30f30g30 -B+t-JX<wid>/<hei>,-Rd" # options proj2=X10/5 region2=g-45/315/-90/90 axes2="-Bxa60f60g60 -Bya30f30g30 -B+t-JX<wid>/<hei>,-Rg-4... # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt end }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JKf, -JKs, Eckert projection, エケルト図法 [#s2c7a506] - -JKf<lon0>/<width> :エケルト第4図法 - -JKs<lon0>/<width> :エケルト第6図法 lon0 は経度方向の中心である．緯度の annotations を入れる... #codeprettify(lang-bash){{ # output file output="JKf_JKs" # options proj1=Kf180/10 region1=d0/360/-89/89 axes1="-Bxa60f60g60 -Bya30f30g30 -B+t-JKf180/<wid>,-Rd0/3... # options proj2=Ks180/10 region2=d0/360/-89/89 axes2="-Bxa60f60g60 -Bya30f30g30 -B+t-JKs180/<wid>,-Rd0/3... # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~

End:
*GMT -J オプション [#q288969a] #contents **はじめに [#jac8759d] 共通オプション -J は projection を指し，描画をするコマン... 図の種類，形状，スケールがこのオプションによって決定され... J の次に必ずアルファベット1〜2文字が入り，描く図法の種類... **projection指定子の大文字・小文字の違い [#d57d020d] J に続く英字で描く図法が決まるが，その英字が大文字か小文... メルカトル図法を例にすると，おおよそ下記のような書き方が... - -JM<width> - -Jm1:<scale> - -Jm</degree> ここで，<width> は出力する図の幅，<scale> は比，</degree>... 大文字 M は <width> インチなどの出力サイズを指定する方法，~ 小文字 m は実スケールとの比を指定する(出力される紙面にと... この原則は -J オプションのその他(メルカトル図法を表す M, ... ~ それぞれの使用例と出力結果は下記の通り．~ このページ3つの画像の表示サイズを統一しているが，実際の出... #codeprettify(lang-bash){{ #!/bin/bash }} #codeprettify(lang-bash){{ # filename output="JMm" proj1=M10 region1=125/150/25/50 axes1="-Ba5f5g5 -B+t-JM<wid>,-R125/50/25/50" proj2=m1:2.75e7 region2=125/150/25/50 axes2="-Ba5f5g5 -B+t-Jm1:<scale>,-R125/150/25/50" proj3=m0.4 region3=125/150/25/50 axes3="-Ba5f5g5 -B+t-Jm<wid_perdeg>,-R125/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt coast -J$proj3 -R$region3 $axes3 -W -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ** projection, 各図法の記述例 [#d78dd727] 以降では，projection ごとに pscoast または psbasemap を用... 描ける図法の種類が多いので，一部に限定する．~ また，煩雑さを避けるため大文字，-J<uppercase> で統一し，... 図法によってオプション -R や -B も変更する必要があるが，... ~ **-JM, Mercator, メルカトル図法 [#ufe55c9a] 見慣れたメルカトル図法．記法は下記の通り．~ - -JM<width> - -JM<lon0>/<width> - -JM<lon0>/<lat0>/<width> #codeprettify(lang-bash){{ # filename output="JM" proj=M10 region=125/150/25/50 axes="-Ba5f5g5 -B+t-JM<wid>,-R125/50/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JB, Albers conic equal area projection [#e1384514] Conic projection の一つ． - -JB<lon0>/<lat0>/<lat1>/<lat2>/<width> ここで，<lat1>,<lat2> は Alberts projection の two standa... 参考: [[Wikipedia>https://en.wikipedia.org/wiki/Albers_pr... #codeprettify(lang-bash){{ # output file output="JB" # options proj=B135/35/30/45/12 region=120/150/25/50 axes="-Ba5f5g5 -B+t-JB135/35/30/45/<wid>,-R120/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P -V > "${output}... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JU, UTM, ユニバーサル横メルカトル図法 [#j529cb9a] 数値計算でよく UTM 座標を用いるため．計算結果に対して -JU... - -JU<zone>/<width> ここで，<zone> は UTM のゾーン(1〜60の整数)． #codeprettify(lang-bash){{ # output file output="JU" # options proj=U53/10 region=125/150/25/50 axes="-Ba5f5g5 -B+t-JU<utmzone>/<wid>,-R125/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JC, Cylindrical projection, 円筒図法 [#hc899061] UTM に似ている．UTM でいう zone の基準点がないので，<lon0... - -JC<lon0>/<lat0>/<width> #codeprettify(lang-bash){{ # output file output="JC" # options proj=C135/35/10 region=125/150/25/50 axes="-Ba5f5g5 -B+t-JC135/35/<wid>,-R125/150/25/50" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JA, Lambert Azimuthal projection, ランベルト正積方位図... ランベルト図法．図の特性上,描ける最大範囲は半球．~ オプション -R の標記が上記のタイプの指定方法と異なるので... - -JA<lon0>/<lat0>/<width> -Rg - -JA<lon0>/<lat0>/<width> -R<lonSW>/<latSW>/<lonNE>/<lat... Rgは<lon0>,<lat0> を中心とする半球を描く．~ "-R<...>r" は南西，北東の端を指定する．最後の【r】をつけ... g も r も付けず，R<lonmin>/<lonmax>/<latmin>/<latmax> の... #codeprettify(lang-bash){{ # output file output="JARgARr" # options proj1=A136/35/10 region1=g axes1="-Ba15f15g15 -B+t-JA135/35/<wid>,-Rg" # options proj2=A120/30/10 region2=90/5/180/65r axes2="-Ba15f15g15 -B+t-JA120/30/<wid>,-R90/5/180/65r" # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JG, Azimuthal orthographic projection, 正射図法 [#db01... 地球表面を正射影して描かれる地図．オプションの表記方法は ... - -JG<lon0>/<lat0>/<width> -Rg さらに詳細な指定もできるが，ここでは省略． #codeprettify(lang-bash){{ # output file output="JG" # options proj=G135/35/10 region=g axes="-Bf30g30 -B+t-JG135/35/<wid>,-Rg" # gmt6 gmt begin $output png gmt coast -J$proj -R$region $axes -W gmt end # gmt5 gmt pscoast -J$proj -R$region $axes -W -P > "${output}_gm... gmt psconvert -Tg -A "${output}_gmt5.ps" }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JX, Non-geographical projection (linear) [#s1cb788c] 地図でない X-Y 座標には -JX を使う． - -JX<width> - -JX<width>/<height> JX<width> の場合は，<height> は <width> と同じ値．~ 範囲指定 -R も基本的には -JM などの場合と同じ． #codeprettify(lang-bash){{ # output file output="JXwXwh" # options proj1=X10 region1=0/100/0/50 axes1="-Bxa20g20 -By10g10 -BSWne+t-JX<wid>,-R0/100/0/50" # options proj2=X10/5 region2=0/100/0/50 axes2="-Bxa20g20 -By10g10 -BSWne+t-JX<wid>/<hei>,-R0/100/... # gmt6 gmt begin $output png gmt basemap -J$proj1 -R$region1 $axes1 gmt basemap -J$proj2 -R$region2 $axes2 -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JXl, Non-geographical projection (log-log, semilog), ... 対数軸にしたい場合は，横軸なら <width> の後に，縦軸なら <... - -JX<width>l/<height> : 片対数 X - -JX<width>/<height>l : 片対数 Y - -JX<width>l/<height>l : 両対数 ここでは B オプションに記述する数値の意味が他と異なる点に... a, f, g の値が大きいほど間隔が短くなる．p をつけることで ... #codeprettify(lang-bash){{ # output file output="JXl" # options proj=X10l/5 region=1e0/1e5/0/50 axes="-Bxa1pg3p -By10g10 -BSWne+t-JX<wid>l/<hei>,-R1e0/1e... # gmt6 gmt begin $output png gmt basemap -J$proj -R$region $axes gmt end # gmt5 gmt psbasemap -J$proj -R$region $axes -P > "${output}_gmt... gmt psconvert -Tg -A "${output}_gmt5.ps" }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JX,-JZ, Non-geographical projection (3-dimensional) [#... 3D図は，-JX と -JZ で2回 J オプションを使う． - -JX<width> -JZ<height> - -JX<width>/<length> -JZ<height> ここで，<length> は水平平面で <width> と垂直な方向の長さ．~ 回転すると入力した素直な(?)長さにならないので，最終的に出... #codeprettify(lang-bash){{ # output file output="JXJZ.ps" # output file output="JXJZ" # options proj=X10 projZ=Z10 pers=-120/25 region=0/10/0/10/0/10 axes="-Bxa2g2 -Bya2g2 -Bza2f2 -BSWZ+b+t-JX<wid>,-JZ<hei>,... # gmt6 gmt begin $output png gmt basemap -J$proj -R$region $axes -J$projZ -p$pers gmt end # gmt5 gmt psbasemap -J$proj -R$region $axes -P -J$projZ -p$pers... gmt psconvert -Tg -A "${output}_gmt5.ps" }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ ~ **-JX, -Rd or -Rg, 地球全体図 [#m9a91360] メルカトル図法 -JM では，R オプションで極の -90,90 を含め... JX の projection で，地球全体を描く方法がある．-Rg または... - -JX<width> -Rd - -JX<width> -Rd<lonmin>/<lonmax>/<latmin>/<latmax> - -JX<width> -Rg - -JX<width> -Rg<lonmin>/<lonmax>/<latmin>/<latmax> 経度方向の中心は，Rd で本初子午線，Rg で日付変更線になる．~ それ以外の経度を中心にしたい場合は -R で中心に来るように ... #codeprettify(lang-bash){{ # output file # output file output="JXRd_JXRg" # options proj1=X10/5 region1=d axes1="-Bxa60f60g60 -Bya30f30g30 -B+t-JX<wid>/<hei>,-Rd" # options proj2=X10/5 region2=g-45/315/-90/90 axes2="-Bxa60f60g60 -Bya30f30g30 -B+t-JX<wid>/<hei>,-Rg-4... # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt end }} #ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~ **-JKf, -JKs, Eckert projection, エケルト図法 [#s2c7a506] - -JKf<lon0>/<width> :エケルト第4図法 - -JKs<lon0>/<width> :エケルト第6図法 lon0 は経度方向の中心である．緯度の annotations を入れる... #codeprettify(lang-bash){{ # output file output="JKf_JKs" # options proj1=Kf180/10 region1=d0/360/-89/89 axes1="-Bxa60f60g60 -Bya30f30g30 -B+t-JKf180/<wid>,-Rd0/3... # options proj2=Ks180/10 region2=d0/360/-89/89 axes2="-Bxa60f60g60 -Bya30f30g30 -B+t-JKs180/<wid>,-Rd0/3... # gmt6 gmt begin $output png gmt coast -J$proj1 -R$region1 $axes1 -W gmt coast -J$proj2 -R$region2 $axes2 -W -X11.5 gmt end }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i... ~

Page: Edit with a page name which already exists