gmt set FONT_ANNOT_PRIMARY 12
gmt set MAP_ANNOT_OFFSET_PRIMARY 0.3c
gmt set FONT_LABEL 12p
gmt set MAP_LABEL_OFFSET 0.5c
gmt set MAP_TICK_LENGTH_PRIMARY -0.2c

cat << EOT > _awk.awk
BEGIN{
    printf "%.3f %.3f\n",5.5,-3.75
    printf "%.3f %.3f\n",0.0,-3.75
    printf "%.3f %.3f\n",0.0,-2.75
    a=4.5
    b=1.5
    x0=4.5
    y0=-2.75
    for(i=0;i<=450;i++){
        x=int(i)*0.01
        y=y0+b/a*sqrt(a*a-(x-x0)*(x-x0))
        printf "%.3f %.3f\n",x,y
    }
    printf "%.3f %.3f\n",5.0,-1.25
    printf "%.3f %.3f\n",5.0,1.25
    y0=2.75
    for(i=0;i<=450;i++){
        x=4.5-int(i)*0.01
        y=y0-b/a*sqrt(a*a-(x-x0)*(x-x0))
        printf "%.3f %.3f\n",x,y
    }
    printf "%.3f %.3f\n",0.0,3.75
    printf "%.3f %.3f\n",5.5,3.75
    printf "%.3f %.3f\n",5.5,3.25
    printf "%.3f %.3f\n",0.5,3.25
    a=4.25
    b=1.25
    x0=4.75
    y0=3.00
    for(i=0;i<=400;i++){
        x=0.5+int(i)*0.01
        y=y0-b/a*sqrt(a*a-(x-x0)*(x-x0))
        printf "%.3f %.3f\n",x,y
    }
    printf "%.3f %.3f\n",5.5,1.75
    printf "%.3f %.3f\n",5.5,-1.75
    y0=-3.00
    for(i=0;i<=400;i++){
        x=4.5-int(i)*0.01
        y=y0+b/a*sqrt(a*a-(x-x0)*(x-x0))
        printf "%.3f %.3f\n",x,y
    }
    printf "%.3f %.3f\n",0.5,-3.25
    printf "%.3f %.3f\n",5.5,-3.25
}
EOT
awk -f _awk.awk > _inp.txt
rangeSW=-3/9/-6/6
scale=12/12
afgS=a1f0.5
afgW=a1f0.5
labelS='Distance in x-direction (m)'
labelW='Distance in y-direction (m)'
fig=fig_plane.eps
gmt psbasemap -R$rangeSW -JX$scale -Bx$afgS+l"$labelS" -By$afgW+l"$labelW" -BWSen -P -K > $fig
gmt psxy _inp.txt -R -J -G#cccccc -W2 -K -O >> $fig
gmt psxy -R -J -W2 -K -O << EOT >> $fig
>
0.0  2.750
0.0 -2.750
>
0.5  2.063
0.5 -2.063
>
4.5  1.250
4.5 -1.250
>
 0.0  4.25
-0.5  4.25
-0.5 -4.25
 0.0 -4.25
EOT
gmt psxy -R -J -W2,#000000,3_3:0 -K -O << EOT >> $fig
>
0.0 -4.25
5.5 -4.25
5.5  4.25
0.0  4.25
>
5.5 1.75
7.0 1.75
>
5.5 1.25
7.0 1.25
>
5.5 -1.25
7.0 -1.25
>
5.5 -1.75
7.0 -1.75
EOT
gmt psxy -R -J -W0.5,#000000 -K -O << EOT >> $fig
>
1.5  1.632
1.5 -1.632
>
2.5  1.406
2.5 -1.406
>
3.5  1.288
3.5 -1.288
EOT
gmt psxy -R -J -G#ffffff -W1,#000000 -K -O << EOT >> $fig
2.0 -0.1
2.5 -0.1
2.5 -0.3
2.8  0.0
2.5  0.3
2.5  0.1
2.0  0.1
EOT
# Streight line
gmt psxy -R -JX -Sv0.2+s -W0.5 -N -K -O << EOT >> $fig
-0.5  4.25 -1.5  4.25
 0.0  3.75 -1.5  3.75
 0.0  2.75 -1.5  2.75
 0.0 -2.75 -1.5 -2.75
 0.0 -3.75 -1.5 -3.75
-0.5 -4.25 -1.5 -4.25
-0.5 -4.25 -0.5 -5.0
 0.0 -3.75  0.0 -5.0
 0.5 -2.063 0.5 -5.0
 4.5 -1.25  4.5 -5.0
 5.0 -1.25  5.0 -5.0
 5.5 -1.75  5.5 -5.0
5.5  4.25 7.5  4.25
5.5  3.75 7.5  3.75
5.5  3.25 7.5  3.25
7.0  1.75 7.5  1.75
7.0  1.25 7.5  1.25
7.0 -1.25 7.5 -1.25
7.0 -1.75 7.5 -1.75
5.5 -3.25 7.5 -3.25
5.5 -3.75 7.5 -3.75
5.5 -4.25 7.5 -4.25
EOT
# Both side arrow
gmt set MAP_VECTOR_SHAPE 1
gmt psxy -R -JX -Sv0.2+s+b+e+a45 -W0.5 -G0 -N -O -K << EOT >> $fig
-1.4  4.25 -1.4  3.75
-1.4  3.75 -1.4  2.75
-1.4  2.75 -1.4 -2.75
-1.4 -2.75 -1.4 -3.75
-1.4 -3.75 -1.4 -4.25
-0.5 -4.9 0.0 -4.9
 0.0 -4.9 0.5 -4.9
 0.5 -4.9 4.5 -4.9
 4.5 -4.9 5.0 -4.9
 5.0 -4.9 5.5 -4.9
7.4  4.25 7.4  3.75
7.4  3.75 7.4  3.25
7.4  3.25 7.4  1.75
7.4  1.75 7.4  1.25
7.4  1.25 7.4 -1.25
7.4 -1.25 7.4 -1.75
7.4 -1.75 7.4 -3.25
7.4 -3.25 7.4 -3.75
7.4 -3.75 7.4 -4.25
EOT
gmt psxy -R -JX -Sv0.2+s+e+a60 -W1 -G0 -N -O -K << EOT >> $fig
2.0 1.0 2.0 1.503
EOT
# Text
gmt pstext -R -J -F+f+a+j -N -K -O << EOT >> $fig
-0.25 0 12p,Helvetica 90 MC EL.804.0
 0.25 0 12p,Helvetica 90 MC EL.805.0
 4.75 0 12p,Helvetica 90 MC EL.804.0
 5.25 0 12p,Helvetica 90 MC EL.808.0
 3.00 0 12p,Helvetica 0 ML 1:4.0
 2.50  3.50 12p,Helvetica 0 MC EL.808.0
 2.50 -3.50 12p,Helvetica 0 MC EL.808.0
 2.50  4.00 12p,Helvetica 0 MC EL.804.0
 2.50 -4.00 12p,Helvetica 0 MC EL.804.0
-1.7  4.00 12p,Helvetica 90 MC 0.500
-1.1  3.25 12p,Helvetica 90 MC 1.000
-1.7  0.00 12p,Helvetica 90 MC 5.500
-1.1 -3.25 12p,Helvetica 90 MC 1.000
-1.7 -4.00 12p,Helvetica 90 MC 0.500
-0.25 -5.00 12p,Helvetica 0 TC 0.500
 0.25 -4.80 12p,Helvetica 0 BC 0.500
 2.50 -5.00 12p,Helvetica 0 TC 4.000
 4.75 -4.80 12p,Helvetica 0 BC 0.500
 5.25 -5.00 12p,Helvetica 0 TC 0.500
 3.50  2.50 12p,Helvetica 0 MC Backfill
 3.50 -2.50 12p,Helvetica 0 MC Backfill
7.7  4.00 12p,Helvetica 90 MC 0.500
7.1  3.50 12p,Helvetica 90 MC 0.500
7.7  2.50 12p,Helvetica 90 MC 1.500
7.1  1.50 12p,Helvetica 90 MC 0.500
7.7  0.00 12p,Helvetica 90 MC 2.500
7.1 -1.50 12p,Helvetica 90 MC 0.500
7.7 -2.50 12p,Helvetica 90 MC 1.500
7.1 -3.50 12p,Helvetica 90 MC 0.500
7.7 -4.00 12p,Helvetica 90 MC 0.500
2.5 0.75 10p,Helvetica 0 MC (X/4.5)@+2@++(Y/1.5)@+2@+=1
EOT
echo '0 0' | gmt psxy -R -J -Sp -O >> $fig