#Makefile for evanston
#https://web.mst.edu/sgao/gmtexs/grd/ex01/plot.gmt
v:
	gmt grdinfo tmp.grd
#	false
	gmt math -T1/3/1 T = -o1,1,1 |\
	gmt xyz2grd -R-126.000/-100.000/30.500/50.000 -I0.5 -Gtmp.grd 

export CampusDr_SheridanRd=-87.67440 42.05014
export CentralSt_AsburyAv=-87.6893661 42.0642361
export CentralSt_GreeleyAv=-87.73243 42.06466 #NW corner of area that needs separate chunk with oblique axes!
export ChicagoAv_DavisSt=-87.6796573 42.0461939
export ChicagoAv_SheridanRd=-87.6772906 42.0511062
export Howard_StAsburyAv=-87.6901786 42.0195323
export LakeSt_SheridanRd=-87.67314 42.04398
export LakeSt_ShermanAv=-87.68192 42.04409
export UniversityPl_ShermanAv=-87.68177 42.05126
rid=123419 #OSM relation id
# Tilted downtown area.
evanston.view:
%.view:downtown.gmt east.gmt southwest.gmt
	./$* > /tmp/$*.kml
	viking --map 13 --external /tmp/$*.kml
%.gmt:%.cvs
%.gmt:%.kml
	ogr2ogr -nlt POLYGON $@ $?
downtown.csv: #counter clockwise
	/bin/echo -e "\
	$(LakeSt_ShermanAv)\n\
	$(LakeSt_SheridanRd)\n\
	$(CampusDr_SheridanRd)\n\
	$(ChicagoAv_SheridanRd)\n\
	$(UniversityPl_ShermanAv)\n\
	$(LakeSt_ShermanAv)" > $@
$(rid).kml:; ../fetch_boundary $(rid)
#Whole city
boundary_file=$rid.gmt