def test_pvengine_float_inputs_iso(params):
"""Test that PV engine works for float inputs"""
irradiance_model = IsotropicOrdered()
pvarray = OrderedPVArray.init_from_dict(params)
eng = PVEngine(pvarray, irradiance_model=irradiance_model)
# Irradiance inputs
timestamps = dt.datetime(2019, 6, 11, 11)
DNI = 1000.
DHI = 100.
# Fit engine
eng.fit(timestamps, DNI, DHI, params['solar_zenith'],
params['solar_azimuth'], params['surface_tilt'],
params['surface_azimuth'], params['rho_ground'])
# Checks
np.testing.assert_almost_equal(eng.irradiance.direct['front_illum_pvrow'],
DNI)
# Run timestep
pvarray = eng.run_full_mode_timestep(0)
# Checks
assert isinstance(pvarray, OrderedPVArray)
np.testing.assert_almost_equal(
pvarray.pvrows[0].front.get_param_weighted('qinc'), 1099.22245374)
np.testing.assert_almost_equal(
pvarray.pvrows[1].front.get_param_weighted('qinc'), 1099.6948573)
np.testing.assert_almost_equal(
pvarray.pvrows[2].front.get_param_weighted('qinc'), 1102.76149246)
def test_pvengine_float_inputs_perez(params):
"""Test that PV engine works for float inputs"""
irradiance_model = HybridPerezOrdered()
pvarray = OrderedPVArray.init_from_dict(params)
eng = PVEngine(pvarray, irradiance_model=irradiance_model)
# Irradiance inputs
timestamps = dt.datetime(2019, 6, 11, 11)
DNI = 1000.
DHI = 100.
# Fit engine
eng.fit(timestamps, DNI, DHI, params['solar_zenith'],
params['solar_azimuth'], params['surface_tilt'],
params['surface_azimuth'], params['rho_ground'])
# Checks
np.testing.assert_almost_equal(eng.irradiance.direct['front_illum_pvrow'],
DNI)
# Run timestep
pvarray = eng.run_full_mode_timestep(0)
# Checks
assert isinstance(pvarray, OrderedPVArray)
np.testing.assert_almost_equal(
pvarray.pvrows[0].front.get_param_weighted('qinc'), 1110.1164773159298)
np.testing.assert_almost_equal(
pvarray.pvrows[1].front.get_param_weighted('qinc'), 1110.595903991)
np.testing.assert_almost_equal(
pvarray.pvrows[2].front.get_param_weighted('qinc'), 1112.37717553)
np.testing.assert_almost_equal(
pvarray.pvrows[1].back.get_param_weighted('qinc'), 116.49050349491208)
def test_pvengine_float_inputs_perez_transparency_spacing_full(params):
"""Test that module transparency and spacing are having the
expected effect to calculated PV back side irradiance"""
# Irradiance inputs
timestamps = dt.datetime(2019, 6, 11, 11)
DNI = 1000.
DHI = 100.
# --- with 0 transparency and spacing
# Create models
irr_params = {'module_transparency': 0., 'module_spacing_ratio': 0.}
irradiance_model = HybridPerezOrdered(**irr_params)
pvarray = OrderedPVArray.init_from_dict(params)
eng = PVEngine(pvarray, irradiance_model=irradiance_model)
# Fit engine
eng.fit(timestamps, DNI, DHI, params['solar_zenith'],
params['solar_azimuth'], params['surface_tilt'],
params['surface_azimuth'], params['rho_ground'])
# Run timestep
pvarray = eng.run_full_mode_timestep(0)
no_spacing_transparency_back_qinc = (
pvarray.pvrows[1].back.get_param_weighted('qinc'))
# --- with non-0 transparency and spacing
# Create models
irr_params = {'module_transparency': 0.1, 'module_spacing_ratio': 0.1}
irradiance_model = HybridPerezOrdered(**irr_params)
pvarray = OrderedPVArray.init_from_dict(params)
eng = PVEngine(pvarray, irradiance_model=irradiance_model)
# Fit engine
eng.fit(timestamps, DNI, DHI, params['solar_zenith'],
params['solar_azimuth'], params['surface_tilt'],
params['surface_azimuth'], params['rho_ground'])
# Run timestep
pvarray = eng.run_full_mode_timestep(0)
# Checks
expected_back_qinc = 132.13881181118185 # higher than when params are 0
w_spacing_transparency_back_qinc = (
pvarray.pvrows[1].back.get_param_weighted('qinc'))
np.testing.assert_almost_equal(w_spacing_transparency_back_qinc,
expected_back_qinc)
assert no_spacing_transparency_back_qinc < w_spacing_transparency_back_qinc
def test_vf_matrix_subset_calculation(params):
"""Check that the vf matrix subset is calculated correctly"""
# Run in fast mode
irradiance_model = HybridPerezOrdered()
pvarray = OrderedPVArray.init_from_dict(params)
fast_mode_pvrow_index = 1
eng = PVEngine(pvarray,
irradiance_model=irradiance_model,
fast_mode_pvrow_index=fast_mode_pvrow_index)
timestamps = dt.datetime(2019, 6, 11, 11)
DNI = 1000.
DHI = 100.
eng.fit(timestamps, DNI, DHI, params['solar_zenith'],
params['solar_azimuth'], params['surface_tilt'],
params['surface_azimuth'], params['rho_ground'])
pvarray_fast = _fast_mode_with_loop(eng.pvarray, eng.irradiance,
eng.vf_calculator,
fast_mode_pvrow_index, 0)
vf_mat_fast = pvarray_fast.vf_matrix
# Run in full mode
irradiance_model = HybridPerezOrdered()
pvarray = OrderedPVArray.init_from_dict(params)
eng = PVEngine(pvarray, irradiance_model=irradiance_model)
timestamps = dt.datetime(2019, 6, 11, 11)
DNI = 1000.
DHI = 100.
eng.fit(timestamps, DNI, DHI, params['solar_zenith'],
params['solar_azimuth'], params['surface_tilt'],
params['surface_azimuth'], params['rho_ground'])
pvarray_full = eng.run_full_mode_timestep(0)
vf_mat_full = pvarray_full.vf_matrix
index_of_back_side_surface = 13
np.testing.assert_array_almost_equal(
vf_mat_fast, vf_mat_full[[index_of_back_side_surface], :])