Takuya Miyashita
This site
Web
Start:
*到達時間 プロット メモ [#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...
fmin=${f//\.grd/\-min.grd}
gmt grdmath $f 60 DIV = $fmin=cf
gmt grdimage $fmin -J$proj -Baf -BneSW -R$f -C$cpt -K...
gmt grdcontour $fmin -J$proj -R$f -C30 -Wthinnest -K ...
gmt psscale -C$cpt -Np -Bxa60f15g30 -By+lmin -DJMR+w5...
## 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
}}
&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
End:
*到達時間 プロット メモ [#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...
fmin=${f//\.grd/\-min.grd}
gmt grdmath $f 60 DIV = $fmin=cf
gmt grdimage $fmin -J$proj -Baf -BneSW -R$f -C$cpt -K...
gmt grdcontour $fmin -J$proj -R$f -C30 -Wthinnest -K ...
gmt psscale -C$cpt -Np -Bxa60f15g30 -By+lmin -DJMR+w5...
## 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
}}
&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/i...
Page:
Edit with a page name which already exists