Refactor bulk formulas for use of xarray Datasets and DataArrays

I needed to update what the bulk formulas take as arguments, which became possible due to the CF name convention.