cat << EOT > _inp.txt
12 $ 12 $ tr $ 0.5 $ 0.5
14 $ 0.7 $ 4
Gauss   $ 2         $ C0.3 $ 0/0/0 $ ---
Simpson $ 2,3_3:0   $ S0.3 $ 255/255/255 $ 1
EOT
awk -f /Users/kk/DATA_GMT/awk_klegend.awk _inp.txt > _legend.txt

rangeSW=1/1000000/0.000000000000001/0.1
rangeNE=1/1000000/0.000000000000001/0.1
scale=12l/12l
afgS=g1a1p
afgW=g1a1p
afgN=
afgE=
labelS="Number of elements"
labelW="Calculation error"
labelN=""
labelE=""
gmt set FONT_ANNOT_PRIMARY 14
gmt set FONT_LABEL 14
gmt set MAP_TICK_LENGTH_PRIMARY 0.0c
gmt set MAP_ANNOT_OFFSET_PRIMARY 0.3c
fig=fig_di.eps
#===============================
#Plotting
#===============================
gmt psbasemap -R$rangeSW -JX$scale -Bx$afgS+l"$labelS" -By$afgW+l"$labelW" -BSWne -P -K > $fig
inp_1=_out_g.txt
inp_2=_out_s.txt
awk '{if($4<0)$4=-$4;if($4<1e-15)$4=1e-15;print $1*$2,$4}' $inp_1 | gmt psxy -R -J -B -W2       -O -K >> $fig
awk '{if($4<0)$4=-$4;if($4<1e-15)$4=1e-15;print $1*$2,$4}' $inp_2 | gmt psxy -R -J -B -W2,3_3:0 -O -K >> $fig
awk '{if($4<0)$4=-$4;if($4<1e-15)$4=1e-15;print $1*$2,$4}' $inp_1 | gmt psxy -R -J -B -SC0.3 -G0/0/0 -N -O -K >> $fig
awk '{if($4<0)$4=-$4;if($4<1e-15)$4=1e-15;print $1*$2,$4}' $inp_2 | gmt psxy -R -J -B -SS0.3 -G255/255/255 -W1 -N -O -K >> $fig
#===============================
#Drawing of legend
#===============================
source _legend.txt


#********************************************
#Delete of work files
#********************************************
rm _*.*