check_comm_all.py
View Options

			import xarray as xr
			import sys
			import os
			import pandas as pd

			for year in [2010, 2016, 2017, 2018, 2019]:
			print(year)
			fp = 'files_%d' %year

			# CHECK FOR CONSISTENCE WITH ROOFS:
			roof_pv = pd.read_csv(os.path.join(fp, '%d_annual_roofs.csv' %year), index_col = 0)
			comm_pv = pd.read_csv(os.path.join(fp, '%d_commune_hourly_pv.csv' %year), index_col = 0)

			roof_pv['PV_roof'] = roof_pv.suitable * roof_pv.EPV_kWh
			print(roof_pv['PV_roof'].sum() / 1e9 )

			PV_roof_aggr = roof_pv.groupby('BFS_NUMMER').PV_roof.sum()
			PV_roof_aggr.index = PV_roof_aggr.index.astype('int')

			PV_comm_aggr = comm_pv.sum().rename('PV_comm')
			PV_comm_aggr.index = PV_comm_aggr.index.astype('int')

			max_diff = (abs(PV_comm_aggr-PV_roof_aggr ) / PV_roof_aggr).max()

			if max_diff > 1e-5:
			print('Significant deviation (%f)' %(max_diff))
			else:
			print('Execution OK')

solarPOT/Potential_hourly/check_comm_all.py29f993e4425aaw_solar