decompose main bayes_metab function

there's a subsection for munging discharge data if needed, a section for deciding whether the split-dates option is acceptable, and different handling for split_dates == T vs F. this makes the function quite long.

consider removing the option to split dates, and with it the bayes_1ply function