Nbreen run is pretty much what I want, except for reasonably-well-understood issue #125.  (Another known issue is that the 'nbreen' mask is not being used the way it was in the hydrolakes Matlab scripts, wherein (nbreen==0) areas get W set to 0.  But this may be a feature too far.)