#!/bin/sh
r.mask -r
(
echo "stusps,cells"
for st in $(cat conus_statefp.txt); do
	ST=$(v.db.select -c tl_2022_us_state col=STUSPS where="STATEFP='$st'")
	st=$(echo $ST | tr A-Z a-z)
	g.region rast=fdr_state_$st -a
	r.stats -cn fdr_state_$st | awk -v st=$ST '{cells+=$2}END{printf "%s,%d\n", st, cells}'
done
) > fdr_states.cells
