#author("2021-05-13T02:38:07+09:00","default:Miyashita","Miyashita") #author("2022-07-08T14:02:07+09:00","default:Miyashita","Miyashita") *到達時間 プロット メモ [#ddeafb5e] **例 [#j102dbdd] #codeprettify{{ #!/bin/bash ## directory figdir="_arrival" if [ ! -e $figdir ]; then mkdir $figdir fi ## makecpt cpt="tmparrival.cpt" gmt makecpt -Chaxby -I -T0/180/10 -D > $cpt ## plot for each for f in SD*.arrival_time.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)}') fmin=${f//\.grd/\-min.grd} gmt grdmath $f 60 DIV = $fmin=cf gmt grdimage $fmin -J$proj -Baf -BneSW -R$f -C$cpt -K > $outps gmt grdcontour $fmin -J$proj -R$f -C30 -Wthinnest -K -O >> $outps && rm -f $fmin gmt psscale -C$cpt -Np -Bxa60f15g30 -By+lmin -DJMR+w5.0/0.3+o1.0/0.0+ef0.2 -J -R -O >> $outps ## convert ps2pdf $outps ${outps//.ps/.pdf} && rm $outps # PDF pdfcrop --margins '5 5 5 5' ${outps//.ps/.pdf} | grep -iv pdf # PDF crop convert -density 300 ${outps//.ps/-crop.pdf} ${outps//.ps/.png} && rm ${outps//.ps/.pdf} # PDF to PNG ## convert gmt psconvert -A -Tf $outps # PS -> PDF gmt psconvert -A -TG $outps # PS -> PNG ## move mv ${outps//.ps/-crop.pdf} ${outps//.ps/.png} $figdir/ mv ${outps//.ps/.pdf} ${outps//.ps/.png} $figdir/ done rm $cpt }} &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD02.arrival_time.png,579x403); &ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD03.arrival_time.png,579x425);