#!/bin/bash

forests="Arapaho_and_Roosevelt_National_Forests Medicine_Bow_Routt_National_Forest Grand_Mesa_Uncompahgre_and_Gunnison_National_Forests Rio_Grande_National_Forest San_Juan_National_Forest"

# source="USFS_MVUM_Roads USFS_MVUM_Trails USFS_Trails USGS_Topo_Roads OSM_Roads"

conflator=~/projects/HOT/osm-merge.git/highways/osm_merge/conflator.py

if test x"$1" != x ; then
    forests="$1"
fi
   
for forest in ${forests}; do
    echo "Processing forest ${forest}, conflating ${primary} with ${secondary}..."
    for task in ${forest}_Tasks/*Tasks_[0-9].geojson; do
	base=$(echo ${task} | sed -e 's/Tasks_[0-9]*.geojson//')
	num=$(echo ${task} | grep -o "Tasks_[0-9]*")
	echo "    Processing task ${task}"
	${conflator} -v -s ${forest}_Tasks/${forest}_${num}_USGS_Topo_Roads.geojson -e ${forest}_Tasks/${forest}_${num}_USFS_MVUM_Roads.geojson -o ${forest}_Tasks/USGS_MVUM_${num}.geojson
	size=$(wc -l ${forest}_Tasks/USGS_MVUM_${num}-out.osm | cut -d ' ' -f 1)
	if test ${size} -eq 3; then
	    echo "No data in USGS_MVUM_${num}-out.osm)"
	    continue
	fi
	${conflator} -v -s ${forest}_Tasks/USGS_MVUM_${num}-out.geojson -e ${forest}_Tasks/${forest}_${num}_OSM_Roads.osm -o ${forest}_Tasks/OSM_${num}.geojson

    done
done