/* Create workspaces for the Hood Canal basins, so that we can repeat
/* all processing witht the new software to guarantee consistency.

/*&do bas &list alderbrook big_anderson big_beef big_bend big_quilcene devereaux dewatto dosewallips duckabush eagle enati finch fulton hamma_hamma happy_hollow hill holyoke jorsted lilliwaup little_mission little_quilcene miller mission mulberg seabeck skokomish stavis stimson tahuya tarboo thorndyke trails_end twanoh twanoh_falls union union_store unnamed wacketickeh

&s buffer 225
&sys touch basinslabs.gen
&set name3 alderbrook
&set name4 big_bend
&set name5 unnamed
&set name6 union_store
&set name8 enati
&set name9 twanoh
&set name10 twanoh_falls
&set name11 mulberg
&set name12 happy_hollow
&set name13 trails_end
&set name14 hill
&set name15 holyoke
&set name16 finch
&set name17 devereaux
&set name19 stimson
&set name20 miller
&set name21 little_mission
&set name22 eagle
&set name24 lilliwaup
&set name25 skokomish
&set name26 jorsted
&set name27 mission
&set name28 dewatto
&set name29 union
&set name30 big_anderson
&set name31 tahuya
&set name32 stavis
&set name33 wacketickeh
&set name34 seabeck
&set name36 big_beef
&set name37 hamma_hamma
&set name38 fulton
&set name39 duckabush
&set name40 dosewallips
&set name41 big_quilcene
&set name42 thorndyke
&set name43 little_quilcene
&set name44 tarboo

&do num &list 3 4 5 6 8 9 10 11 12 13 14 15 16 17 19 20 21 22 24 25 26 27 28 29 30 31 32 33 34 36 37 38 39 40 41 42 43 44
  &s basin [value name%num%]
    &ty %num% %basin%
  setwindow ../mission/master/newids
  setmask off
  newid = scalar(%num%)
  xcenterutm = scalar(0)
  ycenterutm = scalar(0)
  xmin = scalar(999999)
  xmax = scalar(0)
  ymin = scalar (9999999)
  ymax = scalar(0)
  docell
    if (..\mission\master\newids == newid)
     begin
      if (xmin > $$xmap)
        xmin = $$xmap
      if (xmax < $$xmap)
        xmax = $$xmap
      if (ymax == 0)
        ymax = $$ymap
      ymin = $$ymap
    end
    if (..\mission\master\centroids == newid)
     begin
      xcenterutm = $$xmap
      ycenterutm = $$ymap
    end
  end
  xm = scalar(xmin - %buffer%)
  xx = scalar(xmax + %buffer%)
  ym = scalar(ymin - %buffer%)
  yx = scalar(ymax + %buffer%)
/*print %num%,xcenterutm,ycenterutm,xm,xx,ym,yx
  &sys echo %num%,xcenterutm,ycenterutm >> ../basinslabs.gen
  setwindow [show xm] [show ym] [show xx] [show yx]
  &if ^ [exist %basin% -directory] &then
    arc createworkspace %basin%
  &work %basin%
  mask = con(..\..\mission\master\newids == %num%,%num%)
  if [exist dem -grid] &then
    kill dem
  dem = ..\..\mission\master\demff * 0.1 /* standard version needs meters
  &if %num% eq 5 &then
    flowd = ..\..\mission\master\tweakdir 
  &else &if  %num% eq 6 &then
    flowd = ..\..\mission\master\tweakdir6
  &else
    flowd = ..\..\mission\master\flowdir 
  &work ..
  &label CONTINUE
  &typ On to the next
&end