bankfull: determine bankfull condition from yearly peak data

• changed bankfull calculation to use yearly peak events
• return-period is now really given in years
• vectorized calculations to prevent python-loops
• update call signature to synchronize with cli:

  bankfull_discharge(in_file, out_file, return_period, wetted_perimeter, var)