#!/bin/sh
# run from epsg5070/PERMANENT
for st in $(cat $(dirname $0)/conus_statefp.txt); do
	g.region rast=ned_1_conus_clean -ap
	r.mask conus_statefp mask=$(echo $st | sed 's/^0\+//') --o
	g.region rast=ned_1_conus_clean zoom=MASK -ap
	st=$(v.db.select -c tl_2022_us_state col=STUSPS where="STATEFP='$st'" | tr A-Z a-z)
	r.mapcalc ex="ned_1_state_$st=ned_1_conus_clean"
	r.watershed -sa ned_1_state_$st drain=drain_state_$st
	r.mapcalc ex="fdr_state_$st=2^(8-abs(drain_state_$st))"
	r.out.gdal fdr_state_$st out=fdr_state_$st.tif
done