Now PISMLakesHydrology can supply the 'bwp' diagnostic if requested.
This run script uses a function to avoid code duplication.  Seems to work.