Homec4science
Diffusion solar_potential (master) Browse

solarPOT/EW_scenarios_MMH/baseline_scenario/compute_PV_output_all_V2.py
29f993e4425amaster

Last Change
Raw File

compute_PV_output_all_V2.py
View Options

import pandas as pd
import xarray as xr
import numpy as np
import util
from pvlib.tools import cosd
import os
import time
import sys
## for more details, check reference in Total_radiation_spielwiese
TILTED_IRRAD_FP = sys.argv[1]
START = int(sys.argv[2])
N_LINES = int(sys.argv[3])
OUTNAME = sys.argv[4]
## ===== CONSTANTS ===== ##
CHUNKSIZE = 10000
N_DF_UIDs = 9639231
workdir = '/scratch/walch/workdir_solar_MMH_base'
PARTLY_SHADING_FP = '/work/hyenergy/output/solar_potential/geographic_potential/available_area/shading_stats/CH_partly_shading_replaced_debiased.nc'
SVF_FP = '/work/hyenergy/output/solar_potential/geographic_potential/tilted_irradiance/SVF_CH_replaced.csv'
AVAIL_AREA_FP = '/work/hyenergy/output/solar_potential/geographic_potential/available_area/PRED_available_area_unc_filled.csv'
MUTUAL_SHADE_FP = '/work/hyenergy/output/solar_potential/technical_potential/flat_roofs/mutual_shading.nc'
ROOFS_FP = '/work/hyenergy/output/solar_potential/geographic_potential/CH_ROOFS_all_replaced.csv'
EFFICIENCY_FP = '/work/hyenergy/output/solar_potential/technical_potential/pv_efficiencies_proj_roofs.nc'
OUTFILE_FP = OUTNAME
OUT_BASE = os.path.splitext(OUTFILE_FP)[0]
END = min(START + N_LINES, N_DF_UIDs)
FLAT_TILT = 30
FLAT_DIST = 1
print('\nShading file path: %s' %PARTLY_SHADING_FP)
print('Tilted irradiance file path: %s' %TILTED_IRRAD_FP)
print('Skyview factor file path: %s' %SVF_FP)
print('Available area file path: %s' %AVAIL_AREA_FP)
print('Mutual shading file path: %s' %MUTUAL_SHADE_FP)
print('Roof surfaces file path: %s' %ROOFS_FP)
print('Efficiency file path: %s' %EFFICIENCY_FP)
print('Output file path: %s' %OUTFILE_FP)
## =========================================================================================== ##
## ===== LOAD DATA ===== ##
tt = time.time()
t0 = time.time()
irrad = xr.open_dataset( TILTED_IRRAD_FP )
print(irrad)
shading = xr.open_dataset( PARTLY_SHADING_FP ).sel(DF_UID = irrad.DF_UID.values).compute()
svf = util.read_csv_rows( SVF_FP, START, N_LINES )
av_area = util.read_csv_rows( AVAIL_AREA_FP, START, N_LINES, usecols = ['DF_UID', 'tilted_area', 'available_area'])
mutual_shade = xr.open_dataset( MUTUAL_SHADE_FP )
roofs = util.read_csv_rows( ROOFS_FP, START, N_LINES, usecols = ['DF_UID', 'NEIGUNG', 'panel_tilt'])
eff = xr.open_dataset( EFFICIENCY_FP ).sel(DF_UID = irrad.DF_UID.values).compute()
print('\n')
print(shading)
print('\n')
print(irrad)
print('\n')
print(svf.head(2))
print(svf.tail(2))
print('\n')
print(mutual_shade)
print('\n')
print(av_area.head(2))
print(av_area.tail(2))
print('\n')
print(roofs.head(2))
print(roofs.tail(2))
print('\n')
print(eff)
print('Loaded all data in %.3f seconds' %(time.time()-tt))
flat_roofs = roofs[roofs.NEIGUNG == 0]
av_area = av_area.merge( roofs, on = 'DF_UID', how = 'left' )
av_area['panelled_area'] = av_area.available_area / cosd( av_area.panel_tilt - av_area.NEIGUNG )
print('Merged roof info with available area')
## ===== COMPUTE MEAN IRRADIANCE ON ROOFTOPS ===== ##
print('\nComputing mean irradiance on rooftops ...')
tt = time.time()
tilted_irrad = xr.merge( [ shading, irrad ] ).fillna(0)
mutual_shade_flat = mutual_shade.sel(tilt = FLAT_TILT, distance = FLAT_DIST, timestamp = irrad.timestamp).drop(['tilt', 'distance'])
tilted_irrad['svf'] = svf.set_index('DF_UID')['SVF_unbiased'].to_xarray()
tilted_irrad['is_flat'] = tilted_irrad.DF_UID.isin( flat_roofs.DF_UID.values ).astype(int)
print(tilted_irrad)
# compute tilted irradiance for tilted and flat surfaces
# by multiplying mutual shading with is_flat (0, 1)
tilted_irrad['tilted_irradiance'] = tilted_irrad.partly_shading * tilted_irrad.poa_direct * (1 - mutual_shade_flat.shading * tilted_irrad.is_flat) + tilted_irrad.svf * tilted_irrad.poa_sky_diffuse * (1 - mutual_shade_flat.SVF_shade * tilted_irrad.is_flat) + tilted_irrad.poa_ground_diffuse
print(tilted_irrad)
print('Finished in %.3f seconds' %(time.time()-tt))
## ===== PREPARE PV OUTPUT ===== ##
tt = time.time()
print('\nComputing PV potential ...')
# First get all required variables for computation of PV output
pv_output = tilted_irrad.tilted_irradiance.to_dataset()
pv_output['panelled_area'] = av_area.set_index('DF_UID').to_xarray().panelled_area
pv_output['tilted_area'] = av_area.set_index('DF_UID').to_xarray().tilted_area
pv_output['available_area'] = av_area.set_index('DF_UID').to_xarray().available_area
pv_output['efficiency'] = eff['module_eff']
pv_output['performance_factor'] = eff['PF_mean']
# Compute total PV output
pv_output['pv_potential'] = pv_output.tilted_irradiance * pv_output.panelled_area * pv_output.efficiency * pv_output.performance_factor
print('Finished in %.3f seconds' %(time.time()-tt))
tt = time.time()
print('\nComputing yearly sums ...')
# Compute yearly sums
pv_output['yearly_kWh_m2'] = (util.get_yearly_sum(pv_output, 'tilted_irradiance')) / 1000
pv_output['yearly_PV_kWh'] = (util.get_yearly_sum(pv_output, 'pv_potential' )) / 1000
# Compute also yearly PV output for entire rooftops (without considering available area ratio)
# pv_output['yearly_PV_all_area'] = pv_output.yearly_kWh_m2 * pv_output.tilted_area * pv_output.efficiency * pv_output.performance_factor
print(pv_output)
print('Finished in %.3f seconds' %(time.time()-tt))
tt = time.time()
print('\nSaving output ...')
pv_save = pv_output[['tilted_irradiance', 'pv_potential', 'yearly_kWh_m2', 'yearly_PV_kWh']]
pv_save.to_netcdf( OUTFILE_FP )
print('Saved output in %s (%.3fs)' %(OUTFILE_FP, time.time()-tt))
print('\nDONE - total execution time: %.3f seconds' %(time.time()-t0))
c4science · Help