Currently matrix_multiply_grouped_table called under compute_biomass has ops cross between using pandas dataframes and xarray dataset/dataarray. If all inputs uses xr.dataarray as containers, the dimensions would align automatically and this special function is not needed.