aux_df[col]=aux_df['buildings'].values# one dwelling per building assumed when no dwellings column available (the assumption is consistent with the Com and Ind models, for Res the dwellings column exists)
elifcol=='area_per_dwelling_sqm':
aux_df[col]=aux_df['area_per_building_sqm'].values# one dwelling per building assumed