def init_BEM_obj(self):
if not os.path.exists(self.readDOE_file_path):<EOL><INDENT>raise Exception("<STR_LIT>".format(readDOE_file_path))<EOL><DEDENT>readDOE_file = open(self.readDOE_file_path, '<STR_LIT:rb>')  <EOL>refDOE = pickle.load(readDOE_file)<EOL>refBEM = pickle.load(readDOE_file)<EOL>refSchedule = pickle.load(readDOE_file)<EOL>readDOE_file.close()<EOL>k = <NUM_LIT:0><EOL>self.r_glaze_total = <NUM_LIT:0.>             <EOL>self.SHGC_total = <NUM_LIT:0.>                <EOL>self.alb_wall_total = <NUM_LIT:0.>            <EOL>h_floor = self.flr_h or <NUM_LIT>  <EOL>total_urban_bld_area = math.pow(self.charLength, <NUM_LIT:2>)*self.bldDensity *self.bldHeight/h_floor  <EOL>area_matrix = utilities.zeros(<NUM_LIT:16>, <NUM_LIT:3>)<EOL>self.BEM = []           <EOL>self.Sch = []           <EOL>for i in range(<NUM_LIT:16>):    <EOL><INDENT>for j in range(<NUM_LIT:3>):  <EOL><INDENT>if self.bld[i][j] > <NUM_LIT:0.>:<EOL><INDENT>self.BEM.append(refBEM[i][j][self.zone])<EOL>self.BEM[k].frac = self.bld[i][j]<EOL>self.BEM[k].fl_area = self.bld[i][j] * total_urban_bld_area<EOL>if self.glzR:<EOL><INDENT>self.BEM[k].building.glazingRatio = self.glzR<EOL><DEDENT>if self.albRoof:<EOL><INDENT>self.BEM[k].roof.albedo = self.albRoof<EOL><DEDENT>if self.vegRoof:<EOL><INDENT>self.BEM[k].roof.vegCoverage = self.vegRoof<EOL><DEDENT>if self.SHGC:<EOL><INDENT>self.BEM[k].building.shgc = self.SHGC<EOL><DEDENT>if self.albWall:<EOL><INDENT>self.BEM[k].wall.albedo = self.albWall<EOL><DEDENT>if self.flr_h:<EOL><INDENT>self.BEM[k].building.floorHeight = self.flr_h<EOL><DEDENT>self.r_glaze_total += self.BEM[k].frac * self.BEM[k].building.glazingRatio<EOL>self.SHGC_total += self.BEM[k].frac * self.BEM[k].building.shgc<EOL>self.alb_wall_total += self.BEM[k].frac * self.BEM[k].wall.albedo<EOL>self.Sch.append(refSchedule[i][j][self.zone])<EOL>k += <NUM_LIT:1><EOL><DEDENT><DEDENT><DEDENT>
Define BEM for each DOE type (read the fraction)
self.BEM                # list of BEMDef objects
self.r_glaze            # Glazing ratio for total building stock
self.SHGC               # SHGC addition for total building stock
self.alb_wall           # albedo wall addition for total building stock