        dhi = self.tmy["DHI"]
        albedo = 0.16  # Urban environement is 0.14 - 0.22

        # POA Beam component
        df_poa["E_b_poa"] = dni * np.cos(np.radians(aoi))

        # POA Ground component
        df_poa["E_g_poa"] = (
            ghi * albedo * ((1 - np.cos(np.radians(self.surface_tilt))) / 2)
        )