I am using the combinef() function of the R package hts to make my incoherent hierarchical / grouped time series forecasts coherent.
combinef()
hts
The