#!/bin/sh
# run from epsg5070/PERMANENT
for co in $(cat $(dirname $0)/conus_countyns.txt); do
	[ -f fdr_county_$co.tif ] && continue
	g.region rast=ned_1_conus_clean -ap
	r.mask conus_countyns mask=$(echo $co | sed 's/^0\+//') --o
	g.region rast=ned_1_conus_clean zoom=MASK -ap
	r.mapcalc ex="ned_1_county_$co=ned_1_conus_clean"
	r.watershed -sa ned_1_county_$co drain=drain_county_$co
	r.mapcalc ex="fdr_county_$co=2^(8-abs(drain_county_$co))"
	r.out.gdal fdr_county_$co out=fdr_county_$co.tif
done