#!/bin/csh

/bin/rm tmp.ps

## The center of the map is lon = 5.52838, and lat = 22.79149 (Algeria)


#########Basemap

pscoast -R0/360/-90/90 -JE5.52838/22.79149/180/4 \
-P -A10000 -W0.5 -K -G250/250/0 -X2 -Y2 > tmp.ps

######### circles
psxy -R -J -O -K  -Sc -L -P -W2 -V <<END>> tmp.ps
5.52838 22.79149 1.0
5.52838 22.79149 2.0
5.52838 22.79149 3.0
5.52838 22.79149 4.0
END

## Center (blue triangle)

psxy -R -J -St0.1 -G0/0/255 -K -O <<END>>tmp.ps
5.52838 22.79149
END

######### Add event locations #########################################
awk '{print $3, $4}' pie.dat | \
psxy -R -J -O -Sc0.04 -L -W6/255/0/0 -K >> tmp.ps

pstext -R-10/10/-10/10 -JX4 -O -N <<END>>tmp.ps
0 -2.5 13 0 5 6 45
0 -5   13 0 5 6 90
0 -7.5 13 0 5 6 135
0 -10  13 0 5 6 180
END


ps2jpeg150 tmp.ps