#!/usr/bin/python3
# -*- coding: utf-8 -*-

import logging
import pandas as pd
import os.path as path

from . import coastdat
from . import powerplants as pg_pp
from . import tools
from feedinlib import powerplants as pp


class Feedin:
    ''

    def __init__(self):
        ''
        pass

    def aggregate_cap_val(self, conn, **kwargs):
        '''
        Returns the normalised feedin profile and installed capacity for
        a given region.

        Parameters
        ----------
        region : Region instance
        region.geom : shapely.geometry object
            Geo-spatial data with information for location/region-shape. The
            geometry can be a polygon/multi-polygon or a point.

        Returns
        -------
        feedin_df : pandas dataframe
            Dataframe containing the normalised feedin profile of the given
            region. Index of the dataframe is the hour of the year; columns
            are 'pv_pwr' and 'wind_pwr'.
        cap : pandas series
            Series containing the installed capacity (in W) of PV and wind
            turbines of the given region.
        '''
        region = kwargs['region']
        [pv_df, wind_df, cap] = self.get_timeseries(
            conn,
            geometry=region.geom,
            **kwargs)

        if kwargs.get('store', False):
            self.store_full_df(pv_df, wind_df, **kwargs)

        # Summerize the results to one column for pv and one for wind
        df = pd.concat([pv_df.sum(axis=1), wind_df.sum(axis=1)], axis=1)
        feedin_df = df.rename(columns={0: 'pv_pwr', 1: 'wind_pwr'})
        cap = cap.sum()

        return feedin_df, cap

    def get_timeseries(self, conn, **kwargs):
        ''
        weather = coastdat.get_weather(
            conn, kwargs['geometry'], kwargs['year'])

        pv_df = 0
        pv_cap = {}
        wind_df = 0
        wind_cap = {}

        if not isinstance(weather, list):
            weather = [weather]

        for w_cell in weather:
            ee_pps = pg_pp.get_energymap_pps(
                conn, geometry1=w_cell.geometry, geometry2=kwargs['geometry'])

            # Find type of wind turbine and its parameters according to the
            # windzone.
            wz = tools.get_windzone(conn, w_cell.geometry)

            kwargs['wind_conv_type'] = (kwargs['wka_model_dc'].get(
                wz, kwargs['wka_model']))
            kwargs['d_rotor'] = (kwargs['d_rotor_dc'].get(
                wz, kwargs['d_rotor']))
            kwargs['h_hub'] = (kwargs['h_hub_dc'].get(wz, kwargs['h_hub']))

            # Determine the feedin time series for the weather cell
            # Wind energy
            wind_peak_power = ee_pps[ee_pps.type == 'wind_power'].cap.sum()
            wind_power_plant = pp.WindPowerPlant(**kwargs)
            wind_series = wind_power_plant.feedin(
                weather=w_cell, installed_capacity=wind_peak_power)
            wind_series.name = w_cell.name
            wind_cap[w_cell.name] = wind_peak_power

            # PV
            pv_peak_power = ee_pps[ee_pps.type == 'solar_power'].cap.sum()
            pv_plant = pp.Photovoltaic(**kwargs)
            pv_series = pv_plant.feedin(
                weather=w_cell, peak_power=pv_peak_power)
            pv_series.name = w_cell.name
            pv_cap[w_cell.name] = pv_peak_power

            # Combine the results to a DataFrame
            try:
                pv_df = pd.concat([pv_df, pv_series], axis=1)
                wind_df = pd.concat([wind_df, wind_series], axis=1)
            except:
                pv_df = pv_series.to_frame()
                wind_df = wind_series.to_frame()

        # Write capacity into a dataframe
        capw = pd.Series(pd.DataFrame.from_dict(wind_cap, orient='index')[0])
        capw.name = 'wind_pwr'
        cappv = pd.Series(pd.DataFrame.from_dict(pv_cap, orient='index')[0])
        cappv.name = 'pv_pwr'
        cap = pd.concat([capw, cappv], axis=1)

        return pv_df, wind_df, cap

    def store_full_df(self, pv_df, wind_df, **kwargs):
        ''
        dpath = kwargs.get(
            'dpath', path.join(path.expanduser("~"), '.oemof'))
        filename = kwargs.get('filename', 'feedin_' + kwargs['region'].name)
        fullpath = path.join(dpath, filename)

        if kwargs['store'] == 'hf5':
            pv_df.to_hdf(fullpath + '.hf5', 'pv_pwr')
            wind_df.to_hdf(fullpath + '.hf5', 'wind_pwr')

        if kwargs['store'] == 'csv':
            pv_df.to_csv(fullpath + '_pv.csv')
            wind_df.to_csv(fullpath + '_wind.csv')

@@ -1449,24 +1449,48 @@
 ne for wind%0A
+        cap = cap.sum()%0A
         df =
@@ -1518,17 +1518,45 @@
 (axis=1)
-,
+ / cap%5B'pv_pwr'%5D,%0A           
  wind_df
@@ -1567,16 +1567,34 @@
 (axis=1)
+ / cap%5B'wind_pwr'%5D
 %5D, axis=
@@ -1667,32 +1667,8 @@
 r'%7D)
-%0A        cap = cap.sum()
 %0A%0A  

Scale feedin_df to installed capacity of 1