how to do AEP analysis for single turbine data ?

[bharathrajad]

Currently the plant data is working for farm with multiple turbines. can flag be added so that it can work in a focussed manner for AEP gap calculation alone without Wake losses and yaw losses. when i try with single turbine data, i get error during the calculate_asset_distance_matrix function invoking.

[RHammond2]

@bharathrajad could you share either a screenshot or the traceback itself in a python block like the following?

```python
<insert traceback here>
```

[bharathrajad]

@RHammond2 , attaching the traceback. i currently bypassed the error by adding a dummy turbine in asset table data.

WARNING:py.warnings:/anaconda3/envs/tsi_env/lib/python3.10/site-packages/pyproj/transformer.py:817: DeprecationWarning: Conversion of an array with ndim > 0 to a scalar is deprecated, and will error in future. Ensure you extract a single element from your array before performing this operation. (Deprecated NumPy 1.25.)
  return self._transformer._transform_point(

---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/pandas/core/indexes/range.py:413, in RangeIndex.get_loc(self, key)
    412 try:
--> 413     return self._range.index(new_key)
    414 except ValueError as err:

ValueError: 0 is not in range

The above exception was the direct cause of the following exception:

KeyError                                  Traceback (most recent call last)
Cell In[1], line 12
      9 get_ipython().run_line_magic('autoreload', '2')
     11 # Load plant object and process plant data
---> 12 project = project_GREENTECH.prepare('./data/greentech/', use_cleansed=False)
     14 # Add the analysis workflow validations needed below and re-validate
     15 project.analysis_type.extend(["TurbineLongTermGrossEnergy", "ElectricalLosses"])

File ~/home_bharath/windmill/22_LEAP/OpenOA-main/examples/project_GREENTECH.py:299, in prepare(path, return_value, use_cleansed)
    290     return (
    291         scada_df,
    292         meter_df,
   (...)
    296         dict(era5=reanalysis_era5_df),
    297     )
    298 elif return_value == "plantdata":
--> 299     engie_plantdata = PlantData(
    300         analysis_type=CONFIG["plant_data"]["analysis_type"],
    301         metadata=path.parent / files["metadata_file"],
    302         scada=scada_df,
    303         meter=meter_df,
    304         curtail=curtail_df,
    305         asset=asset_df,
    306         #reanalysis={"era5": reanalysis_era5_df, "merra2": reanalysis_merra2_df},
    307         reanalysis={"era5": reanalysis_era5_df},
    308     )
    309     return engie_plantdata
    310 else:

File <attrs generated methods openoa.plant.PlantData>:43, in __init__(self, log_level, metadata, analysis_type, scada, meter, tower, status, curtail, asset, reanalysis)
     41     __attr_validator_asset(self, __attr_asset, self.asset)
     42     __attr_validator_reanalysis(self, __attr_reanalysis, self.reanalysis)
---> 43 self.__attrs_post_init__()

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/openoa/plant.py:467, in PlantData.__attrs_post_init__(self)
    464     raise ValueError(error_message)
    466 # Post-validation data manipulations
--> 467 self.calculate_asset_geometries()
    468 if self.asset is not None:
    469     self.parse_asset_geometry()

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/openoa/plant.py:1384, in PlantData.calculate_asset_geometries(self)
   1382 if self.asset is not None:
   1383     self.parse_asset_geometry()
-> 1384     self.calculate_asset_distance_matrix()
   1385     self.calculate_asset_direction_matrix()

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/openoa/logging.py:33, in logged_method_call.<locals>._wrapper(self, *args, **kwargs)
     31 logger = logging.getLogger(the_method.__module__)
     32 logger.debug(f"{self.__class__.__name__}#{id(self)}.{the_method.__name__}: {msg}")
---> 33 return the_method(self, *args, **kwargs)

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/openoa/plant.py:1239, in PlantData.calculate_asset_distance_matrix(self)
   1231 """Calculates the distance between all assets on the site with `np.inf` for the distance
   1232 between an asset and itself.
   1233 
   1234 Returns:
   1235     pd.DataFrame: Dataframe containing distances between each pair of assets
   1236 """
   1237 ix = self.asset.index.values
   1238 distance = (
-> 1239     pd.DataFrame(
   1240         [i, j, self.asset.loc[i, "geometry"].distance(self.asset.loc[j, "geometry"])]
   1241         for i, j in itertools.combinations(ix, 2)
   1242     )
   1243     .pivot(index=0, columns=1, values=2)
   1244     .rename_axis(index={0: None}, columns={1: None})
   1245     .fillna(0)
   1246     .loc[ix[:-1], ix[1:]]
   1247 )
   1249 # Insert the first column and last row because the self-self combinations are not produced in the above
   1250 distance.insert(0, ix[0], 0.0)

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/pandas/core/frame.py:9346, in DataFrame.pivot(self, columns, index, values)
   9339 @Substitution("")
   9340 @Appender(_shared_docs["pivot"])
   9341 def pivot(
   9342     self, *, columns, index=lib.no_default, values=lib.no_default
   9343 ) -> DataFrame:
   9344     from pandas.core.reshape.pivot import pivot
-> 9346     return pivot(self, index=index, columns=columns, values=values)

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/pandas/core/reshape/pivot.py:553, in pivot(data, columns, index, values)
    549         index_list = [
    550             data._constructor_sliced(data.index, name=data.index.name)
    551         ]
    552 else:
--> 553     index_list = [data[idx] for idx in com.convert_to_list_like(index)]
    555 data_columns = [data[col] for col in columns_listlike]
    556 index_list.extend(data_columns)

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/pandas/core/reshape/pivot.py:553, in <listcomp>(.0)
    549         index_list = [
    550             data._constructor_sliced(data.index, name=data.index.name)
    551         ]
    552 else:
--> 553     index_list = [data[idx] for idx in com.convert_to_list_like(index)]
    555 data_columns = [data[col] for col in columns_listlike]
    556 index_list.extend(data_columns)

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/pandas/core/frame.py:4107, in DataFrame.__getitem__(self, key)
   4105 if self.columns.nlevels > 1:
   4106     return self._getitem_multilevel(key)
-> 4107 indexer = self.columns.get_loc(key)
   4108 if is_integer(indexer):
   4109     indexer = [indexer]

File /anaconda3/envs/tsi_env/lib/python3.10/site-packages/pandas/core/indexes/range.py:415, in RangeIndex.get_loc(self, key)
    413         return self._range.index(new_key)
    414     except ValueError as err:
--> 415         raise KeyError(key) from err
    416 if isinstance(key, Hashable):
    417     raise KeyError(key)

KeyError: 0