ENH prepare_bem_model:  add sol option (output filename)