#author("2022-07-08T13:50:37+09:00","default:Miyashita","Miyashita") #author("2022-07-08T13:51:19+09:00","default:Miyashita","Miyashita") *初期水位分布 プロット メモ [#q7ba62e6] **コード例 [#ibe5b33d] 計算したディレクトリの直下で簡単にプロットするスクリプト. #codeprettify{{ ## directory figdir="_initl" if [ ! -e $figdir ]; then mkdir $figdir fi ## makecpt cpt="tmpinitl.cpt" gmt makecpt -Cpolar -T-1/1 -D > $cpt ## plot for each for f in SD0[0-9].initl_disp*.grd do ## filename outps=${f//\.grd/\.ps} ## plot proj="X"$(gmt grdinfo $f -Cn -o0,1,2,3 | awk '{print 10"/"10*($4-$3)/($2-$1)}') #echo $proj gmt grdimage $f -J$proj -Baf -BneSW+t$f -R$f -C$cpt -K > $outps gmt psscale -C$cpt -Bxa0.5f0.1 -By+lm -DJMR+w5.0/0.3+o1.0/0.0+e -J$proj -R$f -O >> $outps ## convert gmt psconvert -A -Tf $outps # PS -> PDF gmt psconvert -A -TG $outps # PS -> PNG gmt psconvert -A -TG $outps # PS -> PNG ## move mv ${outps//.ps/.pdf} ${outps//.ps/.png} $figdir/ done rm $cpt }} **出力例 [#e34d107e] &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD01.initl_disp001.png,438x263); &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD01.initl_disp003.png,438x263); &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD01.initl_disp006.png,438x263);