import math
import sys
sys.path.append('/home/www/claudev.cz/orb3')

from solar_simulation_astronomical import BifacialSolarSimulation
from datetime import datetime

# Create simulation instance
sim = BifacialSolarSimulation()

# Test June 15, 10:00
date = datetime(2024, 6, 15)
solar_elev, solar_az = sim.get_solar_position(date, 10, 0)
print(f'June 15, 10:00: {solar_elev:.1f}°/{solar_az:.0f}°')
print()

configs = [
    ('Vertical North-South', 90, 0, 'ns_vertical'),
    ('Vertical East-West', 90, 90, 'ew_vertical'), 
    ('Tilt 20°', 20, 0, 'tilted'),
    ('Tilt 30°', 30, 0, 'tilted'),
    ('Tilt 45°', 45, 0, 'tilted'),
    ('Tilt 60°', 60, 0, 'tilted'),
    ('Tilt 70°', 70, 0, 'tilted')
]

for name, tilt, azimuth, panel_type in configs:
    front_irr, back_irr = sim.calculate_panel_irradiance(date, 10, 0, tilt, azimuth, panel_type)
    total_irradiance = front_irr + back_irr
    power = (total_irradiance / 1000) * 700
    bifacial_contribution = (back_irr / total_irradiance * 100) if total_irradiance > 0 else 0
    
    print(f'{name}: {power:.1f} W ({bifacial_contribution:.1f}%)')