gmtset ANOT_FONT_SIZE 10 gmtset ANOT_OFFSET 0.2c gmtset LABEL_FONT_SIZE 10 gmtset LABEL_OFFSET 0.1c set range=-1.000/2.000/-1.000/5.000 set scale=6.000/12.000 set ba=a1.000f0.500 set grdspc=0.050 set sradi=0.0354 set conint=100.000 set annint=100.000 set strx="x-distance (mm)" set stry="y-distance (mm)" set scld=6.500/10.200/3.000/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%:%stry%: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% rem *** Drawing of contour of ps1 *** set fig=fig_gmt_ps1.eps grdmask _mask_data.txt -G_mask.grd -I%grdspc% -R%range% -S%sradi% -NNan/0/0 surface _ps1.txt -G_ps1.grd -I%grdspc% -R -T0.0 grdmath _ps1.grd _mask.grd OR = _ps1m.grd psbasemap -R%range% -JX%scale% -B%ba%:%strx%:/%ba%:%stry%:WS -P -K > %fig% grdimage _ps1m.grd -J -C0_cpt.cpt -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 contour of ps1 *** set fig=fig_gmt_ps2.eps surface _ps2.txt -G_ps2.grd -I%grdspc% -R -T0.0 grdmath _ps2.grd _mask.grd OR = _ps2m.grd psbasemap -R%range% -JX%scale% -B%ba%:%strx%:/%ba%:%stry%:WS -P -K > %fig% grdimage _ps2m.grd -J -C0_cpt.cpt -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*