def calc_electricity(scenario, gp, avg_photoperiod, light, days_in_year, years):
    """Calculate Electricity costs Function

    Args:
       scenario (object): The farm scenario
       gp (object): Growth plan for Farm
       avg_photoperiod (float): Average photoperiod of lights considering crop requirements
       light (object): The light selected by the user
       days_in_year (float): The number of days in a year
       years (int): The no. of years the simulation will analyse

    Returns:
        cogs_electricity (list): Cost of Goods Sold expenditure on Electricity as a time series for each year
        electricity_consumption (list): The amount of eletricity consumed each year

    To Do:
        Add HVAC model and remove HVAC multiplier of 1.25 for enhanced accuracy... leave multiplier as a default?
    """

    electricity_consumption =[0]
    HVAC_multiplier = 1.25

    for y in range(years+1):
        if y == 1:
            electricity_consumption.append(avg_photoperiod * light.max_power * scenario.no_lights_pilot * days_in_year/ 1000)
        elif y > 1: # Multiplied by 1.25 to accomodate HVAC upgrade
            electricity_consumption.append(avg_photoperiod * light.max_power * gp.no_lights_full * days_in_year * HVAC_multiplier/ 1000)

    cogs_electricity = [i * scenario.electricity_price for i in electricity_consumption]
    return cogs_electricity, electricity_consumption