R='70/200/0/73'
J='S140/90/15i'

#RM=117/128/5/20
#JM=M122/30/10i

RP='60/80/22/50'
JP='S70/90/5i'

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 -R$R -Bsg10/g10 -J$J -Gbeige -Wsandybrown -Dl -A1000 -K > map.ps
	awk '{print $$1, $$2, "4 0 4 BL", $$3}' stns.tsv | pstext -R -J -Gbrown -Dj2p/-1p -K -O >> map.ps
	psxy -R -J -Sc2p -Wbrown -O >> map.ps < stns.tsv 

mappk.png: stns.tsv
	gmtset FRAME_WIDTH = 0.25p GRID_PEN_PRIMARY = 0.5p,sandybrown GRID_PEN_SECONDARY = 0.25p,sandybrown
	pscoast -R$(RP) -Bsg10/g10 -J$(JP) -Gbeige -Wsandybrown -Dl -A1000 -K > mappk.ps
	awk '{print $$1, $$2, "4 0 4 BL", $$3}' stns.tsv | pstext -R -J -Gbrown -Dj2p/-1p -K -O >> mappk.ps
	psxy -R -J -Sc2p -Wbrown -O >> mappk.ps < stns.tsv 
	ps2raster -P -Tg mappk.ps
	rm -f mappk.ps

stns.tsv: upr.tsv
	ruby -pe 'f=$$_.split(/\s/); next if f.first.empty?; $$_=f.values_at(1, 0, 2).join("\t") + "\n"' upr.tsv > stns.tsv