Sorry for the delay in getting back to you. Here is the response from our income team.
For complete households you can sum the individual incomes in the benefit unit to get the benefit unit total. The Stata code for this is something like this:
use pidp l_hidp l_buno_dv using l_indall, clear
merge 1:1 pidp using l_indresp, nogen keep(3) keepus(l_fimngrs_dv)
mvdecode _all, mv(-9/-1)
bys l_hidp l_buno_dv: egen l_fimngrs_dv_bu=sum(l_fimngrs_dv)
bys l_hidp l_buno_dv: keep if _n==1
drop l_fimngrs_dv pidp
isid l_hidp l_buno_dv
su
Similarly, you can do this for subcomponents of income that are provided like total benefit income or total earnings. Note, you will not be able to recreate benefit unit totals of individual benefits (eg. benefit unit child benefit income or universal credit income) as these amounts are not on the released datasets (at least not in a cleaned up form with imputed values. Only the raw survey reports are available on the benefit.dta file). Also note that you will also be ignoring council tax deductions which occur at the household level and not deducted from individual incomes (and are deducted in the variable fihhmnnet3_dv measured at the household level).
One thing to consider is that there can be joint reporting of benefits and unearned sources (eg. both partners report child benefit, we identify and give it only to one) which possibly could introduce an issue. When summing at the household level it doesnt matter – the household total will be correct. But at the benefit unit level it will only work if both of the joint recipients are in the same benefit unit – or if not it will be slightly random which benefit unit gets the amount. [In principle this should not be an issue as both recipients should always be in the same benefit unit. But there MAYBE small violations where there are reporting errors and respondents in different benefit units but in the same household report the same benefit (and so it will be a bit random which benefit unit gets assigned the amount). That should really be checked.]
For incomplete households (ie. households where not all members completed an interview) things are more tricky. See variable hhresp_dv. You will not observe the imputed income amounts for non-respondents as these are not available in the data. And so you will have to drop such households which is a kind of selection which they might want to create some weights to adjust depending on their research question.
Hope this helps. If you have further questions please let us know.