{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# MHKiT WEC-Sim Example\n",
    "\n",
    "This notebook demonstrates using the MHKiT wave module and [WEC-Sim](http://wec-sim.github.io/WEC-Sim/index.html) together to perform a resource characterization study in MHKiT, simulate representative cases with WEC-Sim, and visualize the results in MHKiT to estimate MAEP (Mean Annual Energy Production).\n",
    "\n",
    "  1. Characterize the available resource at a location\n",
    "     - See the PacWave example notebook\n",
    "  2. Write a WEC-Sim batch file for the given clusters\n",
    "  3. Simulate the device *in WEC-Sim*.\n",
    "     - Ensure that the spectra used in WEC-Sim is identical to the one used in MHKiT.\n",
    "  4. Load WEC-Sim batch results\n",
    "  5. Assess results and visualize quantities of interest\n",
    "\n",
    "This example uses WEC-Sim to simulate the [Oscillating Surge Wave Energy Converter (OSWEC)](https://wec-sim.github.io/WEC-Sim/main/user/tutorials.html#oscillating-surge-wec-oswec), a flap-type device.\n",
    "\n",
    "Start by importing MHKiT and the necessary python packages (e.g.`scipy.io`, `matplotlib.pyplot`, `pandas`, `numpy`)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from mhkit import wave\n",
    "import scipy.io as sio\n",
    "import matplotlib.pyplot as plt\n",
    "import pandas as pd\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 1. Characterize the available resource at a location\n",
    "\n",
    "This example will use an abbreviated version of `PacWave_resource_characterization_example.ipynb`. \n",
    "For full details on downloading, calculating, and visualizing the k-means clusters representation of the site's wave resouce, see that example.\n",
    "\n",
    "We select the N=32 cluster as it's total energy flux is closet to the total energy flux of the site considering all wave conditions. We will load the PacWave example output, which can be easily saved after running the example with the command `results[32].to_csv(\"pacwave_cluster_32.csv\")`. We will start this example by reading in that csv output and formatting it for WEC-Sim."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Te</th>\n",
       "      <th>Hm0</th>\n",
       "      <th>weights</th>\n",
       "      <th>Tp</th>\n",
       "      <th>J</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>7.974491</td>\n",
       "      <td>1.253970</td>\n",
       "      <td>0.058861</td>\n",
       "      <td>9.294279</td>\n",
       "      <td>6031.115427</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>10.794533</td>\n",
       "      <td>2.641403</td>\n",
       "      <td>0.035216</td>\n",
       "      <td>12.581041</td>\n",
       "      <td>37004.325098</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.901979</td>\n",
       "      <td>1.953122</td>\n",
       "      <td>0.052001</td>\n",
       "      <td>8.044264</td>\n",
       "      <td>12516.214519</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>12.667628</td>\n",
       "      <td>7.310116</td>\n",
       "      <td>0.005070</td>\n",
       "      <td>14.764135</td>\n",
       "      <td>367451.945581</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>12.893701</td>\n",
       "      <td>2.262294</td>\n",
       "      <td>0.016046</td>\n",
       "      <td>15.027624</td>\n",
       "      <td>36455.136139</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>10.557621</td>\n",
       "      <td>4.754297</td>\n",
       "      <td>0.017311</td>\n",
       "      <td>12.304920</td>\n",
       "      <td>116784.361789</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>8.766664</td>\n",
       "      <td>2.739380</td>\n",
       "      <td>0.043646</td>\n",
       "      <td>10.217557</td>\n",
       "      <td>31825.667989</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>6.537403</td>\n",
       "      <td>1.305578</td>\n",
       "      <td>0.050746</td>\n",
       "      <td>7.619350</td>\n",
       "      <td>5272.441394</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>9.666291</td>\n",
       "      <td>1.340694</td>\n",
       "      <td>0.037070</td>\n",
       "      <td>11.266073</td>\n",
       "      <td>8443.649821</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>12.787307</td>\n",
       "      <td>3.920397</td>\n",
       "      <td>0.016464</td>\n",
       "      <td>14.903621</td>\n",
       "      <td>107680.986511</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>11.605879</td>\n",
       "      <td>1.821016</td>\n",
       "      <td>0.022323</td>\n",
       "      <td>13.526666</td>\n",
       "      <td>19446.169168</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>7.584082</td>\n",
       "      <td>1.878735</td>\n",
       "      <td>0.054624</td>\n",
       "      <td>8.839256</td>\n",
       "      <td>12827.143861</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>10.175411</td>\n",
       "      <td>6.133932</td>\n",
       "      <td>0.008836</td>\n",
       "      <td>11.859453</td>\n",
       "      <td>188189.689187</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>9.319357</td>\n",
       "      <td>4.587432</td>\n",
       "      <td>0.018817</td>\n",
       "      <td>10.861722</td>\n",
       "      <td>95155.068368</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>7.228996</td>\n",
       "      <td>1.256691</td>\n",
       "      <td>0.057692</td>\n",
       "      <td>8.425403</td>\n",
       "      <td>5449.034309</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>5.646967</td>\n",
       "      <td>1.339107</td>\n",
       "      <td>0.032118</td>\n",
       "      <td>6.581547</td>\n",
       "      <td>4750.825175</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>7.615980</td>\n",
       "      <td>2.634620</td>\n",
       "      <td>0.036497</td>\n",
       "      <td>8.876433</td>\n",
       "      <td>25339.645267</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>9.460406</td>\n",
       "      <td>3.381147</td>\n",
       "      <td>0.033824</td>\n",
       "      <td>11.026114</td>\n",
       "      <td>52508.670941</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>16.000441</td>\n",
       "      <td>3.044223</td>\n",
       "      <td>0.004295</td>\n",
       "      <td>18.648532</td>\n",
       "      <td>87446.895732</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>10.550343</td>\n",
       "      <td>1.563634</td>\n",
       "      <td>0.030265</td>\n",
       "      <td>12.296437</td>\n",
       "      <td>12622.321616</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>11.817436</td>\n",
       "      <td>2.982923</td>\n",
       "      <td>0.021717</td>\n",
       "      <td>13.773236</td>\n",
       "      <td>53748.089829</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>12.101122</td>\n",
       "      <td>5.305727</td>\n",
       "      <td>0.014540</td>\n",
       "      <td>14.103872</td>\n",
       "      <td>177305.220432</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>10.400035</td>\n",
       "      <td>3.588297</td>\n",
       "      <td>0.032312</td>\n",
       "      <td>12.121253</td>\n",
       "      <td>65405.272028</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>9.132540</td>\n",
       "      <td>2.011568</td>\n",
       "      <td>0.048130</td>\n",
       "      <td>10.643986</td>\n",
       "      <td>17913.129344</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>9.880296</td>\n",
       "      <td>2.462908</td>\n",
       "      <td>0.050634</td>\n",
       "      <td>11.515496</td>\n",
       "      <td>29157.316025</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>8.321803</td>\n",
       "      <td>2.003080</td>\n",
       "      <td>0.054171</td>\n",
       "      <td>9.699071</td>\n",
       "      <td>16104.920042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>6.131352</td>\n",
       "      <td>1.794449</td>\n",
       "      <td>0.035429</td>\n",
       "      <td>7.146098</td>\n",
       "      <td>9295.960216</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>11.430727</td>\n",
       "      <td>3.979891</td>\n",
       "      <td>0.025331</td>\n",
       "      <td>13.322525</td>\n",
       "      <td>90734.821382</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>14.263809</td>\n",
       "      <td>2.781733</td>\n",
       "      <td>0.009359</td>\n",
       "      <td>16.624486</td>\n",
       "      <td>66183.179571</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>13.744161</td>\n",
       "      <td>5.465225</td>\n",
       "      <td>0.007518</td>\n",
       "      <td>16.018835</td>\n",
       "      <td>240475.837506</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>8.735251</td>\n",
       "      <td>1.270630</td>\n",
       "      <td>0.047066</td>\n",
       "      <td>10.180945</td>\n",
       "      <td>6821.344537</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>8.301748</td>\n",
       "      <td>3.676767</td>\n",
       "      <td>0.022070</td>\n",
       "      <td>9.675697</td>\n",
       "      <td>54122.886646</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           Te       Hm0   weights         Tp              J\n",
       "0    7.974491  1.253970  0.058861   9.294279    6031.115427\n",
       "1   10.794533  2.641403  0.035216  12.581041   37004.325098\n",
       "2    6.901979  1.953122  0.052001   8.044264   12516.214519\n",
       "3   12.667628  7.310116  0.005070  14.764135  367451.945581\n",
       "4   12.893701  2.262294  0.016046  15.027624   36455.136139\n",
       "5   10.557621  4.754297  0.017311  12.304920  116784.361789\n",
       "6    8.766664  2.739380  0.043646  10.217557   31825.667989\n",
       "7    6.537403  1.305578  0.050746   7.619350    5272.441394\n",
       "8    9.666291  1.340694  0.037070  11.266073    8443.649821\n",
       "9   12.787307  3.920397  0.016464  14.903621  107680.986511\n",
       "10  11.605879  1.821016  0.022323  13.526666   19446.169168\n",
       "11   7.584082  1.878735  0.054624   8.839256   12827.143861\n",
       "12  10.175411  6.133932  0.008836  11.859453  188189.689187\n",
       "13   9.319357  4.587432  0.018817  10.861722   95155.068368\n",
       "14   7.228996  1.256691  0.057692   8.425403    5449.034309\n",
       "15   5.646967  1.339107  0.032118   6.581547    4750.825175\n",
       "16   7.615980  2.634620  0.036497   8.876433   25339.645267\n",
       "17   9.460406  3.381147  0.033824  11.026114   52508.670941\n",
       "18  16.000441  3.044223  0.004295  18.648532   87446.895732\n",
       "19  10.550343  1.563634  0.030265  12.296437   12622.321616\n",
       "20  11.817436  2.982923  0.021717  13.773236   53748.089829\n",
       "21  12.101122  5.305727  0.014540  14.103872  177305.220432\n",
       "22  10.400035  3.588297  0.032312  12.121253   65405.272028\n",
       "23   9.132540  2.011568  0.048130  10.643986   17913.129344\n",
       "24   9.880296  2.462908  0.050634  11.515496   29157.316025\n",
       "25   8.321803  2.003080  0.054171   9.699071   16104.920042\n",
       "26   6.131352  1.794449  0.035429   7.146098    9295.960216\n",
       "27  11.430727  3.979891  0.025331  13.322525   90734.821382\n",
       "28  14.263809  2.781733  0.009359  16.624486   66183.179571\n",
       "29  13.744161  5.465225  0.007518  16.018835  240475.837506\n",
       "30   8.735251  1.270630  0.047066  10.180945    6821.344537\n",
       "31   8.301748  3.676767  0.022070   9.675697   54122.886646"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "results = pd.read_csv(\"data/wave/pacwave_cluster_32.csv\", index_col=0)\n",
    "results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 2. Write a WEC-Sim batch file for the given clusters\n",
    "\n",
    "[WEC-Sim MCR](https://wec-sim.github.io/WEC-Sim/main/user/advanced_features.html#multiple-condition-runs-mcr) (multiple condition run) files should contain a structure `mcr` that contains two variables: `header` and `cases`. Each column of `header` and `cases` denotes a variable and it's value respectively. Each row is another simulation. WEC-Sim defines waves using the significant wave height and peak period. We will isolate these values from the results of the cluster analysis and create a dictionary that is written to the `.mat` file."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "ws_mcr_cases = results[[\"Hm0\",\"Tp\"]]\n",
    "ws_mcr_header = np.array([\"waves.height\",\"waves.period\"], dtype='object')\n",
    "ws_mcr_out = {'mcr': {'header': ws_mcr_header, 'cases': ws_mcr_cases}}\n",
    "sio.savemat('mcr_mhkit.mat', ws_mcr_out)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 3. Simulate the device *in WEC-Sim*\n",
    "\n",
    "Now that the MCR file is created, we need to go simulate WEC performance in these wave conditions using WEC-Sim. To recreate the data used in the next step, use the created MCR file with WEC-Sim's [OSWEC example](https://github.com/WEC-Sim/WEC-Sim/tree/main/examples/OSWEC). For an accurate comparison to the power calculated in the resource characterization, we should ensure that the WEC-Sim cases use irregular JONSWAP wave spectra as in the PacWave example.\n",
    "\n",
    "For convenience in this demonstration, we enforce OSWEC model stability in the extreme wave conditions by arbitrarily applying a large PTO stiffness and damping:\n",
    "```\n",
    "pto(1).stiffness = 1e5;\n",
    "pto(1).damping = 5e7;\n",
    "```\n",
    "\n",
    "To reduce the amount of extranenous data saved for this example, we limit the WEC-Sim output to the PTO's power output in the `userDefinedFunctions.m` script:\n",
    "```\n",
    "if exist('imcr','var')\n",
    "    if imcr == 1\n",
    "        nmcr = size(mcr.cases,1);\n",
    "        power = nan(1, nmcr);\n",
    "    end\n",
    "\n",
    "    iRampEnd = simu.rampTime./simu.dtOut + 1;\n",
    "    power(imcr) = -mean(output.ptos(1).powerInternalMechanics(iRampEnd:end,5));\n",
    "\n",
    "    if imcr == nmcr\n",
    "        % Save output\n",
    "        save('mcr_mhkit_power.mat', 'power');\n",
    "    end\n",
    "end\n",
    "```"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 4. Load WEC-Sim batch results\n",
    "\n",
    "Note that in this example we do not save the entire WEC-Sim `output` structure for each case. See the `wecsim_example.ipynb` for information on loading that WEC-Sim data. Here the output is one array of average power output that we will load and compare to the resource characterization.\n",
    "\n",
    "Note that the power output \\[W\\] is significantly larger than the energy flux \\[W/m\\] due to the width of the OSWEC. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Te</th>\n",
       "      <th>Hm0</th>\n",
       "      <th>weights</th>\n",
       "      <th>Tp</th>\n",
       "      <th>J</th>\n",
       "      <th>P</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>7.974491</td>\n",
       "      <td>1.253970</td>\n",
       "      <td>0.058861</td>\n",
       "      <td>9.294279</td>\n",
       "      <td>6031.115427</td>\n",
       "      <td>6.861312e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>10.794533</td>\n",
       "      <td>2.641403</td>\n",
       "      <td>0.035216</td>\n",
       "      <td>12.581041</td>\n",
       "      <td>37004.325098</td>\n",
       "      <td>3.873519e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.901979</td>\n",
       "      <td>1.953122</td>\n",
       "      <td>0.052001</td>\n",
       "      <td>8.044264</td>\n",
       "      <td>12516.214519</td>\n",
       "      <td>1.923400e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>12.667628</td>\n",
       "      <td>7.310116</td>\n",
       "      <td>0.005070</td>\n",
       "      <td>14.764135</td>\n",
       "      <td>367451.945581</td>\n",
       "      <td>1.951187e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>12.893701</td>\n",
       "      <td>2.262294</td>\n",
       "      <td>0.016046</td>\n",
       "      <td>15.027624</td>\n",
       "      <td>36455.136139</td>\n",
       "      <td>3.115873e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>10.557621</td>\n",
       "      <td>4.754297</td>\n",
       "      <td>0.017311</td>\n",
       "      <td>12.304920</td>\n",
       "      <td>116784.361789</td>\n",
       "      <td>8.410281e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>8.766664</td>\n",
       "      <td>2.739380</td>\n",
       "      <td>0.043646</td>\n",
       "      <td>10.217557</td>\n",
       "      <td>31825.667989</td>\n",
       "      <td>3.398579e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>6.537403</td>\n",
       "      <td>1.305578</td>\n",
       "      <td>0.050746</td>\n",
       "      <td>7.619350</td>\n",
       "      <td>5272.441394</td>\n",
       "      <td>8.332614e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>9.666291</td>\n",
       "      <td>1.340694</td>\n",
       "      <td>0.037070</td>\n",
       "      <td>11.266073</td>\n",
       "      <td>8443.649821</td>\n",
       "      <td>6.951609e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>12.787307</td>\n",
       "      <td>3.920397</td>\n",
       "      <td>0.016464</td>\n",
       "      <td>14.903621</td>\n",
       "      <td>107680.986511</td>\n",
       "      <td>5.068824e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>11.605879</td>\n",
       "      <td>1.821016</td>\n",
       "      <td>0.022323</td>\n",
       "      <td>13.526666</td>\n",
       "      <td>19446.169168</td>\n",
       "      <td>1.601980e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>7.584082</td>\n",
       "      <td>1.878735</td>\n",
       "      <td>0.054624</td>\n",
       "      <td>8.839256</td>\n",
       "      <td>12827.143861</td>\n",
       "      <td>2.112478e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>10.175411</td>\n",
       "      <td>6.133932</td>\n",
       "      <td>0.008836</td>\n",
       "      <td>11.859453</td>\n",
       "      <td>188189.689187</td>\n",
       "      <td>2.402013e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>9.319357</td>\n",
       "      <td>4.587432</td>\n",
       "      <td>0.018817</td>\n",
       "      <td>10.861722</td>\n",
       "      <td>95155.068368</td>\n",
       "      <td>1.067714e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>7.228996</td>\n",
       "      <td>1.256691</td>\n",
       "      <td>0.057692</td>\n",
       "      <td>8.425403</td>\n",
       "      <td>5449.034309</td>\n",
       "      <td>6.414256e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>5.646967</td>\n",
       "      <td>1.339107</td>\n",
       "      <td>0.032118</td>\n",
       "      <td>6.581547</td>\n",
       "      <td>4750.825175</td>\n",
       "      <td>7.328911e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>7.615980</td>\n",
       "      <td>2.634620</td>\n",
       "      <td>0.036497</td>\n",
       "      <td>8.876433</td>\n",
       "      <td>25339.645267</td>\n",
       "      <td>3.671972e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>9.460406</td>\n",
       "      <td>3.381147</td>\n",
       "      <td>0.033824</td>\n",
       "      <td>11.026114</td>\n",
       "      <td>52508.670941</td>\n",
       "      <td>4.020208e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>16.000441</td>\n",
       "      <td>3.044223</td>\n",
       "      <td>0.004295</td>\n",
       "      <td>18.648532</td>\n",
       "      <td>87446.895732</td>\n",
       "      <td>2.414505e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>10.550343</td>\n",
       "      <td>1.563634</td>\n",
       "      <td>0.030265</td>\n",
       "      <td>12.296437</td>\n",
       "      <td>12622.321616</td>\n",
       "      <td>1.335059e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>11.817436</td>\n",
       "      <td>2.982923</td>\n",
       "      <td>0.021717</td>\n",
       "      <td>13.773236</td>\n",
       "      <td>53748.089829</td>\n",
       "      <td>3.552203e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>12.101122</td>\n",
       "      <td>5.305727</td>\n",
       "      <td>0.014540</td>\n",
       "      <td>14.103872</td>\n",
       "      <td>177305.220432</td>\n",
       "      <td>1.343228e+06</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>10.400035</td>\n",
       "      <td>3.588297</td>\n",
       "      <td>0.032312</td>\n",
       "      <td>12.121253</td>\n",
       "      <td>65405.272028</td>\n",
       "      <td>7.368317e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>9.132540</td>\n",
       "      <td>2.011568</td>\n",
       "      <td>0.048130</td>\n",
       "      <td>10.643986</td>\n",
       "      <td>17913.129344</td>\n",
       "      <td>1.998769e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>9.880296</td>\n",
       "      <td>2.462908</td>\n",
       "      <td>0.050634</td>\n",
       "      <td>11.515496</td>\n",
       "      <td>29157.316025</td>\n",
       "      <td>2.715734e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>8.321803</td>\n",
       "      <td>2.003080</td>\n",
       "      <td>0.054171</td>\n",
       "      <td>9.699071</td>\n",
       "      <td>16104.920042</td>\n",
       "      <td>2.319710e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>6.131352</td>\n",
       "      <td>1.794449</td>\n",
       "      <td>0.035429</td>\n",
       "      <td>7.146098</td>\n",
       "      <td>9295.960216</td>\n",
       "      <td>1.478045e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>11.430727</td>\n",
       "      <td>3.979891</td>\n",
       "      <td>0.025331</td>\n",
       "      <td>13.322525</td>\n",
       "      <td>90734.821382</td>\n",
       "      <td>9.391133e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>14.263809</td>\n",
       "      <td>2.781733</td>\n",
       "      <td>0.009359</td>\n",
       "      <td>16.624486</td>\n",
       "      <td>66183.179571</td>\n",
       "      <td>2.000335e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>13.744161</td>\n",
       "      <td>5.465225</td>\n",
       "      <td>0.007518</td>\n",
       "      <td>16.018835</td>\n",
       "      <td>240475.837506</td>\n",
       "      <td>8.686179e+05</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>8.735251</td>\n",
       "      <td>1.270630</td>\n",
       "      <td>0.047066</td>\n",
       "      <td>10.180945</td>\n",
       "      <td>6821.344537</td>\n",
       "      <td>8.550423e+04</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>8.301748</td>\n",
       "      <td>3.676767</td>\n",
       "      <td>0.022070</td>\n",
       "      <td>9.675697</td>\n",
       "      <td>54122.886646</td>\n",
       "      <td>7.578856e+05</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           Te       Hm0   weights         Tp              J             P\n",
       "0    7.974491  1.253970  0.058861   9.294279    6031.115427  6.861312e+04\n",
       "1   10.794533  2.641403  0.035216  12.581041   37004.325098  3.873519e+05\n",
       "2    6.901979  1.953122  0.052001   8.044264   12516.214519  1.923400e+05\n",
       "3   12.667628  7.310116  0.005070  14.764135  367451.945581  1.951187e+06\n",
       "4   12.893701  2.262294  0.016046  15.027624   36455.136139  3.115873e+05\n",
       "5   10.557621  4.754297  0.017311  12.304920  116784.361789  8.410281e+05\n",
       "6    8.766664  2.739380  0.043646  10.217557   31825.667989  3.398579e+05\n",
       "7    6.537403  1.305578  0.050746   7.619350    5272.441394  8.332614e+04\n",
       "8    9.666291  1.340694  0.037070  11.266073    8443.649821  6.951609e+04\n",
       "9   12.787307  3.920397  0.016464  14.903621  107680.986511  5.068824e+05\n",
       "10  11.605879  1.821016  0.022323  13.526666   19446.169168  1.601980e+05\n",
       "11   7.584082  1.878735  0.054624   8.839256   12827.143861  2.112478e+05\n",
       "12  10.175411  6.133932  0.008836  11.859453  188189.689187  2.402013e+06\n",
       "13   9.319357  4.587432  0.018817  10.861722   95155.068368  1.067714e+06\n",
       "14   7.228996  1.256691  0.057692   8.425403    5449.034309  6.414256e+04\n",
       "15   5.646967  1.339107  0.032118   6.581547    4750.825175  7.328911e+04\n",
       "16   7.615980  2.634620  0.036497   8.876433   25339.645267  3.671972e+05\n",
       "17   9.460406  3.381147  0.033824  11.026114   52508.670941  4.020208e+05\n",
       "18  16.000441  3.044223  0.004295  18.648532   87446.895732  2.414505e+05\n",
       "19  10.550343  1.563634  0.030265  12.296437   12622.321616  1.335059e+05\n",
       "20  11.817436  2.982923  0.021717  13.773236   53748.089829  3.552203e+05\n",
       "21  12.101122  5.305727  0.014540  14.103872  177305.220432  1.343228e+06\n",
       "22  10.400035  3.588297  0.032312  12.121253   65405.272028  7.368317e+05\n",
       "23   9.132540  2.011568  0.048130  10.643986   17913.129344  1.998769e+05\n",
       "24   9.880296  2.462908  0.050634  11.515496   29157.316025  2.715734e+05\n",
       "25   8.321803  2.003080  0.054171   9.699071   16104.920042  2.319710e+05\n",
       "26   6.131352  1.794449  0.035429   7.146098    9295.960216  1.478045e+05\n",
       "27  11.430727  3.979891  0.025331  13.322525   90734.821382  9.391133e+05\n",
       "28  14.263809  2.781733  0.009359  16.624486   66183.179571  2.000335e+05\n",
       "29  13.744161  5.465225  0.007518  16.018835  240475.837506  8.686179e+05\n",
       "30   8.735251  1.270630  0.047066  10.180945    6821.344537  8.550423e+04\n",
       "31   8.301748  3.676767  0.022070   9.675697   54122.886646  7.578856e+05"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Relative location and filename of simulated WEC-Sim data (run with mooring)\n",
    "filename = \"./data/wave/mcr_mhkit_power.mat\"\n",
    "\n",
    "# Load data using the `wecsim.read_output` function which returns a dictionary of Datasets.\n",
    "wecsim_data = sio.loadmat(filename)\n",
    "results['P'] = wecsim_data['power'][0]\n",
    "results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## 5. Assess results and visualize quantities of interest\n",
    "Now that we have loaded the OSWEC's modeled power, we can assess it's performance relative to the incoming wave and calculate the mean annual energy production (MAEP) using MHKiT."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Te</th>\n",
       "      <th>Hm0</th>\n",
       "      <th>weights</th>\n",
       "      <th>Tp</th>\n",
       "      <th>J</th>\n",
       "      <th>P</th>\n",
       "      <th>CW</th>\n",
       "      <th>CWR</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>7.974491</td>\n",
       "      <td>1.253970</td>\n",
       "      <td>0.058861</td>\n",
       "      <td>9.294279</td>\n",
       "      <td>6031.115427</td>\n",
       "      <td>6.861312e+04</td>\n",
       "      <td>11.376523</td>\n",
       "      <td>0.632029</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>10.794533</td>\n",
       "      <td>2.641403</td>\n",
       "      <td>0.035216</td>\n",
       "      <td>12.581041</td>\n",
       "      <td>37004.325098</td>\n",
       "      <td>3.873519e+05</td>\n",
       "      <td>10.467747</td>\n",
       "      <td>0.581542</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>6.901979</td>\n",
       "      <td>1.953122</td>\n",
       "      <td>0.052001</td>\n",
       "      <td>8.044264</td>\n",
       "      <td>12516.214519</td>\n",
       "      <td>1.923400e+05</td>\n",
       "      <td>15.367264</td>\n",
       "      <td>0.853737</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>12.667628</td>\n",
       "      <td>7.310116</td>\n",
       "      <td>0.005070</td>\n",
       "      <td>14.764135</td>\n",
       "      <td>367451.945581</td>\n",
       "      <td>1.951187e+06</td>\n",
       "      <td>5.310048</td>\n",
       "      <td>0.295003</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>12.893701</td>\n",
       "      <td>2.262294</td>\n",
       "      <td>0.016046</td>\n",
       "      <td>15.027624</td>\n",
       "      <td>36455.136139</td>\n",
       "      <td>3.115873e+05</td>\n",
       "      <td>8.547144</td>\n",
       "      <td>0.474841</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>10.557621</td>\n",
       "      <td>4.754297</td>\n",
       "      <td>0.017311</td>\n",
       "      <td>12.304920</td>\n",
       "      <td>116784.361789</td>\n",
       "      <td>8.410281e+05</td>\n",
       "      <td>7.201547</td>\n",
       "      <td>0.400086</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>8.766664</td>\n",
       "      <td>2.739380</td>\n",
       "      <td>0.043646</td>\n",
       "      <td>10.217557</td>\n",
       "      <td>31825.667989</td>\n",
       "      <td>3.398579e+05</td>\n",
       "      <td>10.678736</td>\n",
       "      <td>0.593263</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>6.537403</td>\n",
       "      <td>1.305578</td>\n",
       "      <td>0.050746</td>\n",
       "      <td>7.619350</td>\n",
       "      <td>5272.441394</td>\n",
       "      <td>8.332614e+04</td>\n",
       "      <td>15.804090</td>\n",
       "      <td>0.878005</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>9.666291</td>\n",
       "      <td>1.340694</td>\n",
       "      <td>0.037070</td>\n",
       "      <td>11.266073</td>\n",
       "      <td>8443.649821</td>\n",
       "      <td>6.951609e+04</td>\n",
       "      <td>8.232943</td>\n",
       "      <td>0.457386</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>12.787307</td>\n",
       "      <td>3.920397</td>\n",
       "      <td>0.016464</td>\n",
       "      <td>14.903621</td>\n",
       "      <td>107680.986511</td>\n",
       "      <td>5.068824e+05</td>\n",
       "      <td>4.707260</td>\n",
       "      <td>0.261514</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>11.605879</td>\n",
       "      <td>1.821016</td>\n",
       "      <td>0.022323</td>\n",
       "      <td>13.526666</td>\n",
       "      <td>19446.169168</td>\n",
       "      <td>1.601980e+05</td>\n",
       "      <td>8.238022</td>\n",
       "      <td>0.457668</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>7.584082</td>\n",
       "      <td>1.878735</td>\n",
       "      <td>0.054624</td>\n",
       "      <td>8.839256</td>\n",
       "      <td>12827.143861</td>\n",
       "      <td>2.112478e+05</td>\n",
       "      <td>16.468814</td>\n",
       "      <td>0.914934</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>10.175411</td>\n",
       "      <td>6.133932</td>\n",
       "      <td>0.008836</td>\n",
       "      <td>11.859453</td>\n",
       "      <td>188189.689187</td>\n",
       "      <td>2.402013e+06</td>\n",
       "      <td>12.763787</td>\n",
       "      <td>0.709099</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>9.319357</td>\n",
       "      <td>4.587432</td>\n",
       "      <td>0.018817</td>\n",
       "      <td>10.861722</td>\n",
       "      <td>95155.068368</td>\n",
       "      <td>1.067714e+06</td>\n",
       "      <td>11.220783</td>\n",
       "      <td>0.623377</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>7.228996</td>\n",
       "      <td>1.256691</td>\n",
       "      <td>0.057692</td>\n",
       "      <td>8.425403</td>\n",
       "      <td>5449.034309</td>\n",
       "      <td>6.414256e+04</td>\n",
       "      <td>11.771363</td>\n",
       "      <td>0.653965</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>5.646967</td>\n",
       "      <td>1.339107</td>\n",
       "      <td>0.032118</td>\n",
       "      <td>6.581547</td>\n",
       "      <td>4750.825175</td>\n",
       "      <td>7.328911e+04</td>\n",
       "      <td>15.426605</td>\n",
       "      <td>0.857034</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>7.615980</td>\n",
       "      <td>2.634620</td>\n",
       "      <td>0.036497</td>\n",
       "      <td>8.876433</td>\n",
       "      <td>25339.645267</td>\n",
       "      <td>3.671972e+05</td>\n",
       "      <td>14.491015</td>\n",
       "      <td>0.805056</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>9.460406</td>\n",
       "      <td>3.381147</td>\n",
       "      <td>0.033824</td>\n",
       "      <td>11.026114</td>\n",
       "      <td>52508.670941</td>\n",
       "      <td>4.020208e+05</td>\n",
       "      <td>7.656274</td>\n",
       "      <td>0.425349</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>16.000441</td>\n",
       "      <td>3.044223</td>\n",
       "      <td>0.004295</td>\n",
       "      <td>18.648532</td>\n",
       "      <td>87446.895732</td>\n",
       "      <td>2.414505e+05</td>\n",
       "      <td>2.761110</td>\n",
       "      <td>0.153395</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>10.550343</td>\n",
       "      <td>1.563634</td>\n",
       "      <td>0.030265</td>\n",
       "      <td>12.296437</td>\n",
       "      <td>12622.321616</td>\n",
       "      <td>1.335059e+05</td>\n",
       "      <td>10.576970</td>\n",
       "      <td>0.587609</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>11.817436</td>\n",
       "      <td>2.982923</td>\n",
       "      <td>0.021717</td>\n",
       "      <td>13.773236</td>\n",
       "      <td>53748.089829</td>\n",
       "      <td>3.552203e+05</td>\n",
       "      <td>6.608984</td>\n",
       "      <td>0.367166</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>12.101122</td>\n",
       "      <td>5.305727</td>\n",
       "      <td>0.014540</td>\n",
       "      <td>14.103872</td>\n",
       "      <td>177305.220432</td>\n",
       "      <td>1.343228e+06</td>\n",
       "      <td>7.575793</td>\n",
       "      <td>0.420877</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>10.400035</td>\n",
       "      <td>3.588297</td>\n",
       "      <td>0.032312</td>\n",
       "      <td>12.121253</td>\n",
       "      <td>65405.272028</td>\n",
       "      <td>7.368317e+05</td>\n",
       "      <td>11.265631</td>\n",
       "      <td>0.625868</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>9.132540</td>\n",
       "      <td>2.011568</td>\n",
       "      <td>0.048130</td>\n",
       "      <td>10.643986</td>\n",
       "      <td>17913.129344</td>\n",
       "      <td>1.998769e+05</td>\n",
       "      <td>11.158124</td>\n",
       "      <td>0.619896</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>9.880296</td>\n",
       "      <td>2.462908</td>\n",
       "      <td>0.050634</td>\n",
       "      <td>11.515496</td>\n",
       "      <td>29157.316025</td>\n",
       "      <td>2.715734e+05</td>\n",
       "      <td>9.314074</td>\n",
       "      <td>0.517449</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>8.321803</td>\n",
       "      <td>2.003080</td>\n",
       "      <td>0.054171</td>\n",
       "      <td>9.699071</td>\n",
       "      <td>16104.920042</td>\n",
       "      <td>2.319710e+05</td>\n",
       "      <td>14.403735</td>\n",
       "      <td>0.800207</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>6.131352</td>\n",
       "      <td>1.794449</td>\n",
       "      <td>0.035429</td>\n",
       "      <td>7.146098</td>\n",
       "      <td>9295.960216</td>\n",
       "      <td>1.478045e+05</td>\n",
       "      <td>15.899861</td>\n",
       "      <td>0.883326</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>11.430727</td>\n",
       "      <td>3.979891</td>\n",
       "      <td>0.025331</td>\n",
       "      <td>13.322525</td>\n",
       "      <td>90734.821382</td>\n",
       "      <td>9.391133e+05</td>\n",
       "      <td>10.350087</td>\n",
       "      <td>0.575005</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>14.263809</td>\n",
       "      <td>2.781733</td>\n",
       "      <td>0.009359</td>\n",
       "      <td>16.624486</td>\n",
       "      <td>66183.179571</td>\n",
       "      <td>2.000335e+05</td>\n",
       "      <td>3.022422</td>\n",
       "      <td>0.167912</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>29</th>\n",
       "      <td>13.744161</td>\n",
       "      <td>5.465225</td>\n",
       "      <td>0.007518</td>\n",
       "      <td>16.018835</td>\n",
       "      <td>240475.837506</td>\n",
       "      <td>8.686179e+05</td>\n",
       "      <td>3.612080</td>\n",
       "      <td>0.200671</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>30</th>\n",
       "      <td>8.735251</td>\n",
       "      <td>1.270630</td>\n",
       "      <td>0.047066</td>\n",
       "      <td>10.180945</td>\n",
       "      <td>6821.344537</td>\n",
       "      <td>8.550423e+04</td>\n",
       "      <td>12.534806</td>\n",
       "      <td>0.696378</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>31</th>\n",
       "      <td>8.301748</td>\n",
       "      <td>3.676767</td>\n",
       "      <td>0.022070</td>\n",
       "      <td>9.675697</td>\n",
       "      <td>54122.886646</td>\n",
       "      <td>7.578856e+05</td>\n",
       "      <td>14.003052</td>\n",
       "      <td>0.777947</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "           Te       Hm0   weights         Tp              J             P  \\\n",
       "0    7.974491  1.253970  0.058861   9.294279    6031.115427  6.861312e+04   \n",
       "1   10.794533  2.641403  0.035216  12.581041   37004.325098  3.873519e+05   \n",
       "2    6.901979  1.953122  0.052001   8.044264   12516.214519  1.923400e+05   \n",
       "3   12.667628  7.310116  0.005070  14.764135  367451.945581  1.951187e+06   \n",
       "4   12.893701  2.262294  0.016046  15.027624   36455.136139  3.115873e+05   \n",
       "5   10.557621  4.754297  0.017311  12.304920  116784.361789  8.410281e+05   \n",
       "6    8.766664  2.739380  0.043646  10.217557   31825.667989  3.398579e+05   \n",
       "7    6.537403  1.305578  0.050746   7.619350    5272.441394  8.332614e+04   \n",
       "8    9.666291  1.340694  0.037070  11.266073    8443.649821  6.951609e+04   \n",
       "9   12.787307  3.920397  0.016464  14.903621  107680.986511  5.068824e+05   \n",
       "10  11.605879  1.821016  0.022323  13.526666   19446.169168  1.601980e+05   \n",
       "11   7.584082  1.878735  0.054624   8.839256   12827.143861  2.112478e+05   \n",
       "12  10.175411  6.133932  0.008836  11.859453  188189.689187  2.402013e+06   \n",
       "13   9.319357  4.587432  0.018817  10.861722   95155.068368  1.067714e+06   \n",
       "14   7.228996  1.256691  0.057692   8.425403    5449.034309  6.414256e+04   \n",
       "15   5.646967  1.339107  0.032118   6.581547    4750.825175  7.328911e+04   \n",
       "16   7.615980  2.634620  0.036497   8.876433   25339.645267  3.671972e+05   \n",
       "17   9.460406  3.381147  0.033824  11.026114   52508.670941  4.020208e+05   \n",
       "18  16.000441  3.044223  0.004295  18.648532   87446.895732  2.414505e+05   \n",
       "19  10.550343  1.563634  0.030265  12.296437   12622.321616  1.335059e+05   \n",
       "20  11.817436  2.982923  0.021717  13.773236   53748.089829  3.552203e+05   \n",
       "21  12.101122  5.305727  0.014540  14.103872  177305.220432  1.343228e+06   \n",
       "22  10.400035  3.588297  0.032312  12.121253   65405.272028  7.368317e+05   \n",
       "23   9.132540  2.011568  0.048130  10.643986   17913.129344  1.998769e+05   \n",
       "24   9.880296  2.462908  0.050634  11.515496   29157.316025  2.715734e+05   \n",
       "25   8.321803  2.003080  0.054171   9.699071   16104.920042  2.319710e+05   \n",
       "26   6.131352  1.794449  0.035429   7.146098    9295.960216  1.478045e+05   \n",
       "27  11.430727  3.979891  0.025331  13.322525   90734.821382  9.391133e+05   \n",
       "28  14.263809  2.781733  0.009359  16.624486   66183.179571  2.000335e+05   \n",
       "29  13.744161  5.465225  0.007518  16.018835  240475.837506  8.686179e+05   \n",
       "30   8.735251  1.270630  0.047066  10.180945    6821.344537  8.550423e+04   \n",
       "31   8.301748  3.676767  0.022070   9.675697   54122.886646  7.578856e+05   \n",
       "\n",
       "           CW       CWR  \n",
       "0   11.376523  0.632029  \n",
       "1   10.467747  0.581542  \n",
       "2   15.367264  0.853737  \n",
       "3    5.310048  0.295003  \n",
       "4    8.547144  0.474841  \n",
       "5    7.201547  0.400086  \n",
       "6   10.678736  0.593263  \n",
       "7   15.804090  0.878005  \n",
       "8    8.232943  0.457386  \n",
       "9    4.707260  0.261514  \n",
       "10   8.238022  0.457668  \n",
       "11  16.468814  0.914934  \n",
       "12  12.763787  0.709099  \n",
       "13  11.220783  0.623377  \n",
       "14  11.771363  0.653965  \n",
       "15  15.426605  0.857034  \n",
       "16  14.491015  0.805056  \n",
       "17   7.656274  0.425349  \n",
       "18   2.761110  0.153395  \n",
       "19  10.576970  0.587609  \n",
       "20   6.608984  0.367166  \n",
       "21   7.575793  0.420877  \n",
       "22  11.265631  0.625868  \n",
       "23  11.158124  0.619896  \n",
       "24   9.314074  0.517449  \n",
       "25  14.403735  0.800207  \n",
       "26  15.899861  0.883326  \n",
       "27  10.350087  0.575005  \n",
       "28   3.022422  0.167912  \n",
       "29   3.612080  0.200671  \n",
       "30  12.534806  0.696378  \n",
       "31  14.003052  0.777947  "
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "results['CW'] = wave.performance.capture_width(results['P'], results['J'])\n",
    "oswec_width = 18\n",
    "results['CWR'] = results['CW'] / oswec_width\n",
    "results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "2865149.0"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "MAEP = wave.performance.mean_annual_energy_production_matrix(results['CW'], results['J'], results['weights']) / 1000 # kWh\n",
    "MAEP = np.round(MAEP, 0).item()\n",
    "MAEP"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}