#!/bin/sh
# run from epsg5070/PERMANENT
g.region rast=ned_1_conus_clean -ap
r.mask conus_tl --o
r.mapcalc ex="ned_1_conus_clean_tl=ned_1_conus_clean"
r.watershed -sma ned_1_conus_clean_tl drain=drain_conus_tl memory=100000
r.mapcalc ex="fdr_conus_tl=2^(8-abs(drain_conus_tl))"
r.out.gdal fdr_conus_tl out=fdr_conus_tl.tif