gmtset ANOT_FONT_SIZE 10
gmtset ANOT_OFFSET 0.2c
gmtset LABEL_FONT_SIZE 10
gmtset LABEL_OFFSET 0.1c
set fig=fig_gmt_cont.eps
psbasemap -R0.000/75.000/0.000/20.000 -JX15.000/4.000 -Ba5f1:"Distance (m)":/a5f1:"Altitude (m)":WS -P -K > %fig%
grdmask _mask.txt -G_mask.grd -I0.1 -R -NNaN/0/0
surface _temp1.txt -G_temp.grd -I0.1 -R -T0.0
grdmath _temp.grd _mask.grd OR = _masked_data.grd
grdimage _masked_data.grd -J -C_cp.cpt -K -O >> %fig%
grdcontour _masked_data.grd -C1 -A5 -J -K -O >> %fig%
surface _temp2.txt -G_temp.grd -I0.1 -R -T0.0
grdmath _temp.grd _mask.grd OR = _masked_data.grd
grdcontour _masked_data.grd -C1 -A5 -J -K -O >> %fig%
echo 0 0 | psxy -R -J -Sp -O >> %fig%
set fig=
del .gmt*