Diff of JAGURS/最大水位・最大流速のプロット

#author("2021-05-13T02:48:42+09:00","default:Miyashita","Miyashita")
#author("2022-07-08T13:57:16+09:00","default:Miyashita","Miyashita")
*最大水位・最大流速のプロット　メモ [#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 10"/"10*($4-$3)/($2-$1)}')

    gmt grdimage $f -J$proj -Baf -BneSW -R$f -C$cpt -K > $outps
    gmt psscale -C$cpt -Bxa1f0.5 -By+lm -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_SD01.zmax.png,631x347)
#ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD05.zmax.png,468x369)
~
~
***最大流速 [#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 10"/"10*($4-$3)/($2-$1)}')

    gmt grdimage $f -J$proj -Baf -BneSW -R$f -C$cpt -K > $outps
    gmt psscale -C$cpt -Bxa1f0.5 -By+l"m/s" -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_SD03.vmax.png,421x302);
&ref(https://main-t-miyashita.ssl-lolipop.jp/hydrocoast/image/JAGURS/nankai03_SD04.vmax.png,421x348);
Takuya Miyashita
