R=117/128/5/20
J=M122/30/10i

map.png: map.ps
	ps2raster -Tg map.ps

map.pdf: map.ps
	ps2raster -Tf map.ps

map.ps: stns.tsv
	gmtset FRAME_WIDTH = 0.25p GRID_PEN_PRIMARY = 0.5p,sandybrown GRID_PEN_SECONDARY = 0.25p,sandybrown
	pscoast -P -R$R -Bsg1/g1 -J$J -Gbeige -Wsandybrown -Dl -A1000 -K > map.ps
	awk '{sub(/^98/, "", $$3); print $$1, $$2, "8 0 4 BL", $$3}' stns.tsv | pstext -R -J -Gsandybrown -Dj5p/0p -K -O >> map.ps
	psxy -R -J -Sc9p -Wsandybrown -O >> map.ps < stns.tsv 

stns.tsv:
	wget http://toyoda-eizi.net/tmp/wmo9.sqlite
	sqlite3 -separator "	" -noheader wmo9.sqlite 'select lon, lat, idxnum from vola where lon between 105 and 130 and lat between 5 and 20' > stns.tsv
	rm -f wmo9.sqlite