初期水位分布 プロット メモ †コード例 †計算したディレクトリの直下で簡単にプロットするスクリプト. ## 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 ## move mv ${outps//.ps/.pdf} ${outps//.ps/.png} $figdir/ done rm $cpt 出力例 † |