When you have an equation with more unknowns that constraints and you want to add some assumptions to solve it you may add a relative propensity calculation. Best served by an example.

Given:

 households[uss,usdt,th]
 population[uss,th]

Want to get:

 peoplePerHousehold[uss,usdt,th]

constraints:

 population[uss,th] = sum (households[uss,usdt,th] * peoplePerHousehold[uss,usdt,th]; dim=usdt)

makeup a relative propensity:

 relPopPersPerHH[usdt] = where each value is related to one of the entries 

Finally the math to calculate peoplePerHousehold[uss,usdt,th] is:

  !Equation A:
  local populationEst[uss,th] = sum (households[uss,usdt,th] * relPopPersPerHH[usdt]; dim=usdt)
  !Equation B:
  peoplePerHousehold[uss,usdt,th] = (population[uss,th] / populationEst[uss,th]) * relPopPersPerHH[usdt]