Takuya Miyashita
This site
Web
Start:
*最大水位・最大流速のプロット メモ [#tc2c1cfa]
**例 [#f4fb8307]
***最大水位 [#s2e2b9da]
浸水深をプロットするにはもともとの地盤高と,(設定次第で...
#codeprettify{{
#!/bin/bash
## directory
figdir="_zmax"
if [ ! -e $figdir ]; then
mkdir $figdir
fi
## makecpt
cpt="tmpzmax.cpt"
gmt makecpt -Cseis -I -T0/5 -D > $cpt
## plot for each
for f in SD*.zmax.grd
do
## filename
outps=${f//\.grd/\.ps}
## plot
proj="X"$(gmt grdinfo $f -Cn -o0,1,2,3 | awk '{print...
gmt grdimage $f -J$proj -Baf -BneSW -R$f -C$cpt -K > ...
gmt psscale -C$cpt -Bxa1f0.5 -By+lm -DJMR+w5.0/0.3+o1...
## 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...
~
~
***最大流速 [#fec3e36f]
#codeprettify{{
#!/bin/bash
## directory
figdir="_vmax"
if [ ! -e $figdir ]; then
mkdir $figdir
fi
## makecpt
cpt="tmpvmax.cpt"
echo "0 gray80 0.0001 blue L" > $cpt
gmt makecpt -Cjet -T0.0001/5 -D >> $cpt
## plot for each
for f in SD*.vmax.grd
do
## filename
outps=${f//\.grd/\.ps}
## plot
proj="X"$(gmt grdinfo $f -Cn -o0,1,2,3 | awk '{print...
gmt grdimage $f -J$proj -Baf -BneSW -R$f -C$cpt -K > ...
gmt psscale -C$cpt -Bxa1f0.5 -By+l"m/s" -DJMR+w5.0/0....
## 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:
*最大水位・最大流速のプロット メモ [#tc2c1cfa]
**例 [#f4fb8307]
***最大水位 [#s2e2b9da]
浸水深をプロットするにはもともとの地盤高と,(設定次第で...
#codeprettify{{
#!/bin/bash
## directory
figdir="_zmax"
if [ ! -e $figdir ]; then
mkdir $figdir
fi
## makecpt
cpt="tmpzmax.cpt"
gmt makecpt -Cseis -I -T0/5 -D > $cpt
## plot for each
for f in SD*.zmax.grd
do
## filename
outps=${f//\.grd/\.ps}
## plot
proj="X"$(gmt grdinfo $f -Cn -o0,1,2,3 | awk '{print...
gmt grdimage $f -J$proj -Baf -BneSW -R$f -C$cpt -K > ...
gmt psscale -C$cpt -Bxa1f0.5 -By+lm -DJMR+w5.0/0.3+o1...
## 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...
~
~
***最大流速 [#fec3e36f]
#codeprettify{{
#!/bin/bash
## directory
figdir="_vmax"
if [ ! -e $figdir ]; then
mkdir $figdir
fi
## makecpt
cpt="tmpvmax.cpt"
echo "0 gray80 0.0001 blue L" > $cpt
gmt makecpt -Cjet -T0.0001/5 -D >> $cpt
## plot for each
for f in SD*.vmax.grd
do
## filename
outps=${f//\.grd/\.ps}
## plot
proj="X"$(gmt grdinfo $f -Cn -o0,1,2,3 | awk '{print...
gmt grdimage $f -J$proj -Baf -BneSW -R$f -C$cpt -K > ...
gmt psscale -C$cpt -Bxa1f0.5 -By+l"m/s" -DJMR+w5.0/0....
## 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