rem ****************************************************************************
rem カラーパレットファイル定義
rem ****************************************************************************
echo -20 255 255 190 0 255 255 125 > cp.cpt
echo 0 0 255 255 20 0 125 255 >> cp.cpt
echo N 255 255 255 N 255 255 255 >> cp.cpt

rem ****************************************************************************
rem GMTでの文字の大きさ等定義
rem ****************************************************************************
gmtset ANOT_FONT_SIZE 10
gmtset ANOT_OFFSET 0.2c
gmtset LABEL_FONT_SIZE 10
gmtset LABEL_OFFSET 0.1c

rem ****************************************************************************
rem 出力ファイル：fig_cont1.eps 作図
rem ****************************************************************************
rem ---------------------------------------------------
rem 作図用データ準備
rem ---------------------------------------------------
gawk "BEGIN{FS=\",\"}NR==1715,NR==2314{print $2,$3,$6}" out_core1_2.csv > temp1.txt
gawk "BEGIN{FS=\",\"}NR==1715,NR==2314{print $2,$3,$4}" out_core1_2.csv > temp2.txt
echo  0  0 > mask.txt
echo 10 20 >> mask.txt
echo 12 20 >> mask.txt
echo 22  0 >> mask.txt
set fig=fig_cont1.eps
rem ---------------------------------------------------
rem グラフ軸設定とマスクファイル作成
rem ---------------------------------------------------
psbasemap -R0/26/0/20 -JX5.2/4 -Ba5f1:"Distance (m)":/a5f1:"Altitude (m)":WS -P -K > %fig%
grdmask mask.txt -Gmask.grd -I0.1 -R -NNaN/0/0
rem ---------------------------------------------------
rem 圧力水頭コンター作成（コンター線と塗りつぶし）
rem ---------------------------------------------------
del temp.txt
ren temp1.txt temp.txt
surface temp.txt -Gtemp.grd -I0.1 -R -T0.0
grdmath temp.grd mask.grd OR = masked_data.grd
grdimage masked_data.grd -J -Ccp.cpt -K -O >> %fig%
grdcontour masked_data.grd -C1 -A5 -J -K -O >> %fig%
rem ---------------------------------------------------
rem テキスト挿入（上下流水位）
rem ---------------------------------------------------
echo 1 18 > temp.txt
echo 9 18 >> temp.txt
psxy temp.txt -R -J -W2 -O -K >> %fig%
echo 1 18.5 10 0 0 BL H@#1@#=18m | pstext -R -J -N -K -O >> %fig%
rem ===================================================
echo 20.5 3 > temp.txt
echo 26 3 >> temp.txt
psxy temp.txt -R -J -W2 -O -K >> %fig%
echo 26 3.5 10 0 0 BR H@#2@#=3m | pstext -R -J -N -K -O >> %fig%
rem ---------------------------------------------------
rem 全水頭コンター作成（コンター線）
rem ---------------------------------------------------
del temp.txt
ren temp2.txt temp.txt
surface temp.txt -Gtemp.grd -I0.1 -R -T0.0
grdmath temp.grd mask.grd OR = masked_data.grd
grdcontour masked_data.grd -C1 -A5 -J -O >> %fig%

rem ****************************************************************************
rem 出力ファイル：fig_cont2.eps 作図
rem ****************************************************************************
rem ---------------------------------------------------
rem 作図用データ準備
rem ---------------------------------------------------
gawk "BEGIN{FS=\",\"}NR==2005,NR==2704{print $2,$3,$6}" out_core1_3.csv > temp1.txt
gawk "BEGIN{FS=\",\"}NR==2005,NR==2704{print $2,$3,$4}" out_core1_3.csv > temp2.txt
echo  0  0 > mask.txt
echo 20 20 >> mask.txt
echo 22 20 >> mask.txt
echo 42  0 >> mask.txt
set fig=fig_cont2.eps
rem ---------------------------------------------------
rem グラフ軸設定とマスクファイル作成
rem ---------------------------------------------------
psbasemap -R0/46/0/20 -JX9.2/4 -Ba5f1:"Distance (m)":/a5f1:"Altitude (m)":WS -P -K > %fig%
grdmask mask.txt -Gmask.grd -I0.1 -R -NNaN/0/0
rem ---------------------------------------------------
rem 圧力水頭コンター作成（コンター線と塗りつぶし）
rem ---------------------------------------------------
del temp.txt
ren temp1.txt temp.txt
surface temp.txt -Gtemp.grd -I0.1 -R -T0.0
grdmath temp.grd mask.grd OR = masked_data.grd
grdimage masked_data.grd -J -Ccp.cpt -K -O >> %fig%
grdcontour masked_data.grd -C1 -A5 -J -K -O >> %fig%
rem ---------------------------------------------------
rem テキスト挿入（上下流水位）
rem ---------------------------------------------------
echo  1 18 > temp.txt
echo 18 18 >> temp.txt
psxy temp.txt -R -J -W2 -O -K >> %fig%
echo 1 18.5 10 0 0 BL H@#1@#=18m | pstext -R -J -N -K -O >> %fig%
rem ===================================================
echo 39 3 > temp.txt
echo 46 3 >> temp.txt
psxy temp.txt -R -J -W2 -O -K >> %fig%
echo 46 3.5 10 0 0 BR H@#2@#=3m | pstext -R -J -N -K -O >> %fig%
rem ---------------------------------------------------
rem 全水頭コンター作成（コンター線）
rem ---------------------------------------------------
del temp.txt
ren temp2.txt temp.txt
surface temp.txt -Gtemp.grd -I0.1 -R -T0.0
grdmath temp.grd mask.grd OR = masked_data.grd
grdcontour masked_data.grd -C1 -A5 -J -O >> %fig%

rem ****************************************************************************
rem 出力ファイル：fig_cont3.eps 作図
rem ****************************************************************************
rem ---------------------------------------------------
rem 作図用データ準備
rem ---------------------------------------------------
gawk "BEGIN{FS=\",\"}NR==2554,NR==3870{print $2,$3,$6}" out_fem_zone4m.csv > temp1.txt
gawk "BEGIN{FS=\",\"}NR==2554,NR==3870{print $2,$3,$4}" out_fem_zone4m.csv > temp2.txt
echo  0  0 > mask.txt
echo 40 20 >> mask.txt
echo 45 20 >> mask.txt
echo 75  0 >> mask.txt
set fig=fig_cont3.eps
rem ---------------------------------------------------
rem グラフ軸設定とマスクファイル作成
rem ---------------------------------------------------
psbasemap -R0/79/0/20 -JX15.8/4 -Ba5f1:"Distance (m)":/a5f1:"Altitude (m)":WS -P -K > %fig%
grdmask mask.txt -Gmask.grd -I0.1 -R -NNaN/0/0
rem ---------------------------------------------------
rem 圧力水頭コンター作成（コンター線と塗りつぶし）
rem ---------------------------------------------------
del temp.txt
ren temp1.txt temp.txt
surface temp.txt -Gtemp.grd -I0.1 -R -T0.0
grdmath temp.grd mask.grd OR = masked_data.grd
grdimage masked_data.grd -J -Ccp.cpt -K -O >> %fig%
grdcontour masked_data.grd -C1 -A5 -J -K -O >> %fig%
rem ---------------------------------------------------
rem テキスト挿入（上下流水位）
rem ---------------------------------------------------
echo  1 18 > temp.txt
echo 36 18 >> temp.txt
psxy temp.txt -R -J -W2 -O -K >> %fig%
echo 1 18.5 10 0 0 BL H@#1@#=18m | pstext -R -J -N -K -O >> %fig%
rem ===================================================
echo 69 4 > temp.txt
echo 79 4 >> temp.txt
psxy temp.txt -R -J -W2 -O -K >> %fig%
echo 79 4.5 10 0 0 BR H@#2@#=4m | pstext -R -J -N -K -O >> %fig%
rem ---------------------------------------------------
rem 全水頭コンター作成（コンター線）
rem ---------------------------------------------------
del temp.txt
ren temp2.txt temp.txt
surface temp.txt -Gtemp.grd -I0.1 -R -T0.0
grdmath temp.grd mask.grd OR = masked_data.grd
grdcontour masked_data.grd -C1 -A5 -J -O >> %fig%

rem ****************************************************************************
rem 作業ファイル削除
rem ****************************************************************************
set fig=
del cp.cpt
del mask.*
del temp.*
del masked_data.*
del .gmt*