gmt begin Map png
# size of map
map_radii=5i
# horizon of the map
map_horizon=160
# longitude and latitude of the projection center
centerN=39	 			
centerE=141
# plot map with azimuthal equidistant projection
gmt pscoast -Rg -JE$centerE/$centerN/$map_horizon/$map_radii -A100000 -Gdarkgrey -Swhite -Baf
# add earthquakes data on map
awk '{print ($2,$1)}' earthquakes | gmt psxy -Rg -JE -Sc0.16c -Gbrown
# plot circles at 70,140 degrees
gmt psxy -Rg -JE -SE- -W0.6p,blue << EOF
$centerE $centerN 15540
$centerE $centerN 31080
EOF
gmt end show