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

/* We originally wrote this to create new masks with a buffer of 1.5 cells.
/* However, row&col numbers were inconsistent with earlier results,
/* so we are now copying earlier masks.
/*&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

/* NOT THIS VERSION &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 24 21 43 /* 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%
  setmask off
  newid = scalar(%num%)
  &if ^ [exist %basin% -directory] &then
    arc createworkspace %basin%
  &work %basin%
/*  mask = con(..\..\mission\master\newids == %num%,%num%)
  &if [exist mask -grid] &then ; kill mask
  copy ..\..\harvey-mission\test_basins\%basin%\mask mask
  copy ..\..\harvey-mission\test_basins\%basin%\groundwater groundwater
  /* Let it calculate soil/soildepth
  copy ..\..\harvey-mission\test_basins\%basin%\lulc lulc
  copy ..\..\harvey-mission\test_basins\%basin%\soiltexture soiltexture
  setwindow mask
  &if [exist dem -grid] &then
    kill dem
  dem = ..\..\mission\master\demff * 0.1 /* standard version needs meters
  &if %num% eq 5 &then /* union_store, but how do I tweak the threshold to get a network?
    flowd = ..\..\mission\master\tweakdir 
  &else &if  %num% eq 6 &then  /* enati
    flowd = ..\..\mission\master\tweakdir6
  &else
    flowd = ..\..\mission\master\flowdir 
  &work ..
  &label CONTINUE
  &typ On to the next
&end