/* merge DEM-based basins which retain station ids
/* This will be edited as more basins are added.
setwindow basins
/* newids = merge(..\alderbrook\mask,..\big_bend\mask,..\unnamed\mask,..\union_store\mask,..\enati\mask,..\twanoh\mask,..\twanoh_falls\mask,..\mulberg\mask,..\happy_hollow\mask,..\trails_end\mask,..\hill\mask,..\holyoke\mask,..\finch\mask,..\devereaux\mask,..\stimson\mask,..\miller\mask,..\little_mission\mask,..\eagle\mask,..\lilliwaup\mask,..\skokomish\mask,..\jorsted\mask,..\mission\mask,..\dewatto\mask,..\union\mask,..\big_anderson\mask,..\tahuya\mask,..\stavis\mask,..\wacketickeh\mask,..\seabeck\mask,..\big_beef\mask,..\hamma_hamma\mask,..\fulton\mask,..\duckabush\mask,..\dosewallips\mask,..\big_quilcene\mask,..\thorndyke\mask,..\little_quilcene\mask,..\tarboo\mask)


allgroups = merge(..\alderbrook\mask,..\big_anderson\mask,..\big_beef\mask,..\big_bend\mask,..\big_quilcene\mask,..\dabob_bay\mask,..\devereaux\mask,..\dewatto\mask,..\dosewallips\mask,..\duckabush\mask,..\eagle\mask,..\east_canal\mask,..\east_canal_2\mask,..\enati\mask,..\finch\mask,..\fulton\mask,..\great_bend_no\mask,..\great_bend_s\mask,..\hamma_hamma\mask,..\happy_hollow\mask,..\hill\mask,..\holyoke\mask,..\jorsted\mask,..\lilliwaup\mask,..\little_mission\mask,..\little_quilcene\mask,..\miller\mask,..\mission\mask,..\mulberg\mask,..\northeast\mask,..\northwest\mask,..\s_shr_lynch_c\mask,..\seabeck\mask,..\skokomish\mask,..\stavis\mask,..\stimson\mask,..\tahuya\mask,..\tarboo\mask,..\thorndyke\mask,..\trails_end\mask,..\twanoh\mask,..\twanoh_falls\mask,..\union\mask,..\union_store\mask,..\unnamed\mask,..\wacketickeh\mask,..\west_canal\mask)