gmtset ANOT_FONT_SIZE 10
gmtset ANOT_OFFSET 0.2c
gmtset LABEL_FONT_SIZE 10
gmtset LABEL_OFFSET 0.1c
set range=-1.000/10.000/-1.000/4.000
set scale=12.000/5.455
set ba=a2.000f1.000
set grdspc=0.020
set sradi=0.0141
set conint=1.000
set annint=5.000
set strx="x-distance (m)"
set strz="y-distance (m)"
set scld=12.500/3.791/2.727/0.3
set psize=0.180
set csize=0.720

rem *** Drawing of mesh ***
set fig=fig_gmt_mesh.eps
psbasemap -R%range% -JX%scale% -B%ba%:%strx%:/%ba%:%strz%:WS -P -K > %fig%
set num=01
set col=255/255/255
psxy _mesh%num%.txt -R -J -W0.5/black -G%col% -K -O >> %fig%
psxy _ne_circ.txt -R -J -SC%csize% -N -K -O >> %fig%
pstext _text_ne.txt -R -J -N -K -O >> %fig%
pstext _text_nd.txt -R -J -N -K -O >> %fig%
echo 0 0 | psxy -R -J -Sp -O >> %fig%

grdmask _mask_data.txt -G_mask.grd -I%grdspc% -R%range% -S%sradi% -NNan/0/0
surface _pvec.txt -G_pvec.grd -I%grdspc% -R%range% -T0.0
surface _hvec.txt -G_hvec.grd -I%grdspc% -R%range% -T0.0
grdmath _pvec.grd _mask.grd OR = _pvecm.grd
grdmath _hvec.grd _mask.grd OR = _hvecm.grd
rem *** Drawing of contour ***
set fig=fig_gmt_conp.eps
psbasemap -R%range% -JX%scale% -B%ba%:%strx%:/%ba%:%strz%:WS -P -K > %fig%
grdimage _pvecm.grd -J -C0_cpt.cpt -K -O >> %fig%
grdcontour _pvecm.grd -C%conint% -A%annint% -J -K -O >> %fig%
grdcontour _hvecm.grd -C%conint% -A%annint% -J -K -O >> %fig%
psscale -Ba%annint% -C0_cpt.cpt -D%scld% -E0.3 -K -O >> %fig%
echo 0 0 | psxy -R -J -Sp -O >> %fig%
set fig=fig_gmt_conh.eps
psbasemap -R%range% -JX%scale% -B%ba%:%strx%:/%ba%:%strz%:WS -P -K > %fig%
grdimage _hvecm.grd -J -C0_cpt.cpt -K -O >> %fig%
grdcontour _hvecm.grd -C%conint% -A%annint% -J -K -O >> %fig%
grdcontour _pvecm.grd -C%conint% -A%annint% -J -K -O >> %fig%
psscale -Ba%annint% -C0_cpt.cpt -D%scld% -E0.3 -K -O >> %fig%
echo 0 0 | psxy -R -J -Sp -O >> %fig%

rem *** Drawing of vector ***
set fig=fig_gmt_vect.eps
psbasemap -R%range% -JX%scale% -B%ba%:%strx%:/%ba%:%strz%:WS -P -K > %fig%
grdimage _pvecm.grd -J -C0_cpt.cpt -K -O >> %fig%
psxy _vect.txt -R -J -Svs0.005/0.1/0.03 -G0 -N -K -O >> %fig%
psscale -Ba%annint% -C0_cpt.cpt -D%scld% -E0.3 -K -O >> %fig%
echo 0 0 | psxy -R -J -Sp -O >> %fig%

set fig=
set range=
set scale=
set ba=
set grdspc=
set sradi=
set conint=
set annint=
set strx=
set strz=
set scld=
set psize=
set csize=
del .gmt*