/* Harvey Greenberg  Fri Apr 18 11:05:18 PDT 2008
/* Creates a watershed mask with minimum extent for each of these watersheds:
&s name3  alderbrook  /* Not named 1:24K, build by hand?
&s name4  big_bend    /* tweaked to distinguish it from alderbrook
&s name5  unnamed     /* flowtweak5* needed to build basin
&s name6  union_store /* Are you joking?  One cell in x:/riverhistory/amir/hood_canal/hood_shed.  Never mind, tweakdir6 made it work.
&s name8  enati
&s name9  twanoh
&s name10  twanoh_falls
&s name11  mulberg
&s name12  happy_hollow
&s name13  trails_end
&s name14  hill
&s name15  holyoke
&s name16  finch
&s name17  devereaux
&s name19  stimson
&s name20  miller
&s name21  little_mission
&s name22  eagle
&s name24  lilliwaup
&s name25  skokomish
&s name26  jorsted
&s name27  mission
&s name28  dewatto
&s name29  union
&s name30  big_anderson
&s name31  tahuya
&s name32  stavis
&s name33  wacketickeh
&s name34  seabeck
&s name36  big_beef
&s name37  hamma_hamma
&s name38  fulton
&s name39  duckabush
&s name40  dosewallips
&s name41  big_quilcene
&s name42  thorndyke
&s name43  little_quilcene
&s name44  tarboo
/* Following list edited after examining grid and topos
&s newid44     44
&s newid43     90
&s newid42     95
&s newid41    175
&s newid40    201
&s newid39    394
&s newid38    538
&s newid34    552
&s newid36    521
&s newid33    586
&s newid32    585
&s newid37    537
&s newid30    725
&s newid26    803
&s newid22    832
&s newid29    732
&s newid27    771
&s newid28    774
&s newid24    824
&s newid31    726
&s newid25    808
&s newid21    866
&s newid20    888
&s newid19    894
&s newid16    902
&s newid17    912
&s newid13    986
&s newid14    982
&s newid31    970
&s newid15    954
&s newid11   1003
&s newid12    999
&s newid10   1022
&s newid9   1033
&s newid8   1102
&s newid4   1133

/*&do n &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
/* skip those that are done: 11 22 26 27 37 40 42 , 4 8 9 3 5 6
&do n &list 10 12 13 14 15 16 17 19 20 21 24 25 28 29 30 31 32 33 34 36 38 39 41 43 44
  setwindow basins
  &if ^ [var newid%n%] &then ; &do
    &typ Cannot identify %n%
    xmin = scalar(999999)
    xmax = scalar(0)
    ymin = scalar (9999999)
    ymax = scalar(0)
    n = scalar(%n%)
    docell
      if (..\ws27\allold == n)
       begin
        if (xmin > $$xmap)
          xmin = $$xmap
        if (xmax < $$xmap)
          xmax = $$xmap
        if (ymax == 0)
          ymax = $$ymap
        ymin = $$ymap
      end
    end
    &type Look for %n% in [show xmin] [show ymin] [show xmax] [show ymax] --------------
   &end
  &else ; &do
    &type %n% [value newid%n%] [value name%n%]
    newid = scalar([value newid%n%])
    xcenterutm = scalar(0)
    ycenterutm = scalar(0)
    xmin = scalar(999999)
    xmax = scalar(0)
    ymin = scalar (9999999)
    ymax = scalar(0)
    docell
      if (basins == newid)
       begin
        if (xmin > $$xmap)
          xmin = $$xmap
        if (xmax < $$xmap)
          xmax = $$xmap
        if (ymax == 0)
          ymax = $$ymap
        ymin = $$ymap
      end
      if (centroids == newid)
       begin
        xcenterutm = $$xmap
        ycenterutm = $$ymap
      end
    end
    xm = scalar(xmin - 75)
    xx = scalar(xmax + 75)
    ym = scalar(ymin - 75)
    yx = scalar(ymax + 75)
  /*print %n%,xcenterutm,ycenterutm,xm,xx,ym,yx
    setwindow [show xm] [show ym] [show xx] [show yx]
    &work ..\[value name%n%]
    mask = con(..\master\basins == [value newid%n%],%n%)
    &work ../master
    &label CONTINUE
    &typ On to the next
  &end
&end
/*  &if [exist tmpmask -grid] &then ; kill tmpmask
/*    tmpmask = con(basins eq [value newid%n%],$$xmap)