Faster calculations in general by simplifying code.
Checks are not implemented by default but can be triggered.
Bug fixes for high level function.
For flux_bulk, now kw must be explicitly given due to the considerations required in calculating kw.
If input arrays are xarray.DataArrays, then attributes are written to output (e.g. units).