/* Make a separate mask for each group of unsampled basins
/* Renumber them, and group them, way more complicated than mk_unsampled_groups.aml
&sys touch subgroupxy
&do name &list dabob_bay east_canal east_canal2 great_bend_no great_bend_s northeast northwest s_shr_lynch_c west_canal
  &typ Centers for %name%
  &do bas &list [listfile %name%_* -workspace]
    &ty %name% %bas%
    &if [exist tmpgrid -grid] &then ; kill tmpgrid
    &if [exist tmpgrid2 -grid] &then ; kill tmpgrid2
    tmpgrid = zonalcentroid(%bas%/mask)
    tmpgrid2 = combine(tmpgrid,master/utmx,master/utmy)
    &sys echo %bas% [listunique tmpgrid2.vat -info TMPGRID] [listunique tmpgrid2.vat -info utmx] [listunique tmpgrid2.vat -info utmy] >> subgroupxy
  &end
&end