Takuya Miyashita

初期水位分布　プロット　メモ †

コード例 †

計算したディレクトリの直下で簡単にプロットするスクリプト．

## 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

出力例 †