From c5429a8f798c004f89b50e6fcde669519338c2ee Mon Sep 17 00:00:00 2001 From: ssolson Date: Wed, 3 Jul 2024 09:55:25 -0600 Subject: [PATCH] calc maep_matrix manually --- examples/wave_example.ipynb | 528 +++++++++++++++++++++++------------- 1 file changed, 332 insertions(+), 196 deletions(-) diff --git a/examples/wave_example.ipynb b/examples/wave_example.ipynb index 02680f530..8c5dd63a6 100644 --- a/examples/wave_example.ipynb +++ b/examples/wave_example.ipynb @@ -35,6 +35,14 @@ "execution_count": 2, "metadata": {}, "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "c:\\users\\sterl\\codes\\mhkit-python\\mhkit\\wave\\io\\ndbc.py:99: UserWarning: Could not infer format, so each element will be parsed individually, falling back to `dateutil`. To ensure parsing is consistent and as-expected, please specify a format.\n", + " data = pd.read_csv(\n" + ] + }, { "data": { "text/html": [ @@ -527,6 +535,10 @@ " \n", " Te\n", " \n", + " \n", + " index\n", + " \n", + " \n", " \n", " \n", " \n", @@ -555,6 +567,7 @@ ], "text/plain": [ " Te\n", + "index \n", "2018-01-01 00:40:00 7.458731\n", "2018-01-01 01:40:00 7.682413\n", "2018-01-01 02:40:00 7.498263\n", @@ -601,6 +614,10 @@ " \n", " Hm0\n", " \n", + " \n", + " index\n", + " \n", + " \n", " \n", " \n", " \n", @@ -629,6 +646,7 @@ ], "text/plain": [ " Hm0\n", + "index \n", "2018-01-01 00:40:00 0.939574\n", "2018-01-01 01:40:00 1.001399\n", "2018-01-01 02:40:00 0.924770\n", @@ -675,6 +693,10 @@ " \n", " J\n", " \n", + " \n", + " index\n", + " \n", + " \n", " \n", " \n", " \n", @@ -703,6 +725,7 @@ ], "text/plain": [ " J\n", + "index \n", "2018-01-01 00:40:00 3354.825613\n", "2018-01-01 01:40:00 3916.541523\n", "2018-01-01 02:40:00 3278.298930\n", @@ -742,6 +765,7 @@ { "data": { "text/plain": [ + "index\n", "2018-01-01 00:40:00 7.458731\n", "2018-01-01 01:40:00 7.682413\n", "2018-01-01 02:40:00 7.498263\n", @@ -770,6 +794,7 @@ "name": "stdout", "output_type": "stream", "text": [ + "index\n", "2018-01-01 00:40:00 0.939574\n", "2018-01-01 01:40:00 1.001399\n", "2018-01-01 02:40:00 0.924770\n", @@ -782,6 +807,7 @@ "2018-01-31 22:40:00 2.941428\n", "2018-01-31 23:40:00 2.895928\n", "Name: Hm0, Length: 743, dtype: float64\n", + "index\n", "2018-01-01 00:40:00 3354.825613\n", "2018-01-01 01:40:00 3916.541523\n", "2018-01-01 02:40:00 3278.298930\n", @@ -861,7 +887,7 @@ "\n", " \n", " \n", - " \n", + " \n", " \n", " \n", " \n", @@ -880,6 +906,26 @@ " \n", " \n", " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", " \n", " \n", " \n", @@ -1327,77 +1373,80 @@ "" ], "text/plain": [ - " 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 \\\n", - "0.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "0.5 NaN NaN NaN NaN NaN NaN NaN 0.120286 0.053376 \n", - "1.0 NaN NaN NaN NaN NaN NaN 0.110686 0.068070 0.049452 \n", - "1.5 NaN NaN NaN NaN NaN NaN NaN 0.019749 0.018673 \n", - "2.0 NaN NaN NaN NaN NaN NaN NaN 0.013882 0.012547 \n", - "2.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.007244 \n", - "3.0 NaN NaN NaN NaN NaN NaN NaN 0.004500 0.005660 \n", - "3.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.003924 \n", - "4.0 NaN NaN NaN NaN NaN NaN NaN NaN 0.003185 \n", - "4.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.002343 \n", - "5.0 NaN NaN NaN NaN NaN NaN NaN NaN 0.001913 \n", - "5.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.002101 \n", - "6.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "6.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "7.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "7.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "8.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "8.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "9.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "9.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "10.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "10.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "y_centers 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 \\\n", + "x_centers \n", + "0.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "0.5 NaN NaN NaN NaN NaN NaN NaN 0.120286 0.053376 \n", + "1.0 NaN NaN NaN NaN NaN NaN 0.110686 0.068070 0.049452 \n", + "1.5 NaN NaN NaN NaN NaN NaN NaN 0.019749 0.018673 \n", + "2.0 NaN NaN NaN NaN NaN NaN NaN 0.013882 0.012547 \n", + "2.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.007244 \n", + "3.0 NaN NaN NaN NaN NaN NaN NaN 0.004500 0.005660 \n", + "3.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.003924 \n", + "4.0 NaN NaN NaN NaN NaN NaN NaN NaN 0.003185 \n", + "4.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.002343 \n", + "5.0 NaN NaN NaN NaN NaN NaN NaN NaN 0.001913 \n", + "5.5 NaN NaN NaN NaN NaN NaN NaN NaN 0.002101 \n", + "6.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "6.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "7.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "7.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "8.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "8.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "9.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "9.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "10.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "10.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", "\n", - " 9.0 10.0 11.0 12.0 13.0 14.0 15.0 \\\n", - "0.0 NaN NaN NaN NaN NaN NaN NaN \n", - "0.5 NaN NaN NaN NaN NaN NaN NaN \n", - "1.0 0.065912 NaN 0.056593 0.029950 0.017234 NaN NaN \n", - "1.5 NaN NaN 0.012473 0.011205 0.012307 0.010432 NaN \n", - "2.0 0.009672 0.008770 0.008585 0.007525 0.005272 0.007809 NaN \n", - "2.5 0.006488 0.005788 0.005652 0.005180 0.004260 0.003623 0.004509 \n", - "3.0 0.004691 0.004109 0.003952 0.003104 0.003408 0.002291 0.001792 \n", - "3.5 0.003674 0.003020 0.002746 0.002247 0.002000 0.002257 0.002033 \n", - "4.0 0.002513 0.002386 0.002147 0.002246 0.001605 0.001730 NaN \n", - "4.5 0.002087 0.001919 0.001590 0.001438 NaN NaN NaN \n", - "5.0 0.001720 0.001716 0.001411 0.001219 0.001345 NaN NaN \n", - "5.5 0.001516 0.001331 0.000902 0.001033 NaN NaN NaN \n", - "6.0 0.001097 0.000895 NaN 0.000858 0.000987 NaN NaN \n", - "6.5 0.000837 0.001024 0.000419 NaN 0.000688 NaN NaN \n", - "7.0 NaN NaN NaN 0.000461 0.000633 NaN NaN \n", - "7.5 NaN 0.000553 NaN NaN 0.000312 0.000437 NaN \n", - "8.0 NaN NaN NaN NaN NaN 0.000443 0.000351 \n", - "8.5 NaN NaN NaN NaN NaN NaN 0.000418 \n", - "9.0 NaN NaN NaN NaN NaN NaN NaN \n", - "9.5 NaN NaN NaN NaN NaN NaN 0.000153 \n", - "10.0 NaN NaN NaN NaN NaN NaN NaN \n", - "10.5 NaN NaN NaN NaN NaN NaN 0.000204 \n", + "y_centers 9.0 10.0 11.0 12.0 13.0 14.0 \\\n", + "x_centers \n", + "0.0 NaN NaN NaN NaN NaN NaN \n", + "0.5 NaN NaN NaN NaN NaN NaN \n", + "1.0 0.065912 NaN 0.056593 0.029950 0.017234 NaN \n", + "1.5 NaN NaN 0.012473 0.011205 0.012307 0.010432 \n", + "2.0 0.009672 0.008770 0.008585 0.007525 0.005272 0.007809 \n", + "2.5 0.006488 0.005788 0.005652 0.005180 0.004260 0.003623 \n", + "3.0 0.004691 0.004109 0.003952 0.003104 0.003408 0.002291 \n", + "3.5 0.003674 0.003020 0.002746 0.002247 0.002000 0.002257 \n", + "4.0 0.002513 0.002386 0.002147 0.002246 0.001605 0.001730 \n", + "4.5 0.002087 0.001919 0.001590 0.001438 NaN NaN \n", + "5.0 0.001720 0.001716 0.001411 0.001219 0.001345 NaN \n", + "5.5 0.001516 0.001331 0.000902 0.001033 NaN NaN \n", + "6.0 0.001097 0.000895 NaN 0.000858 0.000987 NaN \n", + "6.5 0.000837 0.001024 0.000419 NaN 0.000688 NaN \n", + "7.0 NaN NaN NaN 0.000461 0.000633 NaN \n", + "7.5 NaN 0.000553 NaN NaN 0.000312 0.000437 \n", + "8.0 NaN NaN NaN NaN NaN 0.000443 \n", + "8.5 NaN NaN NaN NaN NaN NaN \n", + "9.0 NaN NaN NaN NaN NaN NaN \n", + "9.5 NaN NaN NaN NaN NaN NaN \n", + "10.0 NaN NaN NaN NaN NaN NaN \n", + "10.5 NaN NaN NaN NaN NaN NaN \n", "\n", - " 16.0 \n", - "0.0 NaN \n", - "0.5 NaN \n", - "1.0 NaN \n", - "1.5 NaN \n", - "2.0 NaN \n", - "2.5 NaN \n", - "3.0 NaN \n", - "3.5 NaN \n", - "4.0 NaN \n", - "4.5 NaN \n", - "5.0 NaN \n", - "5.5 NaN \n", - "6.0 NaN \n", - "6.5 NaN \n", - "7.0 NaN \n", - "7.5 NaN \n", - "8.0 NaN \n", - "8.5 0.000405 \n", - "9.0 NaN \n", - "9.5 NaN \n", - "10.0 0.000281 \n", - "10.5 0.000225 " + "y_centers 15.0 16.0 \n", + "x_centers \n", + "0.0 NaN NaN \n", + "0.5 NaN NaN \n", + "1.0 NaN NaN \n", + "1.5 NaN NaN \n", + "2.0 NaN NaN \n", + "2.5 0.004509 NaN \n", + "3.0 0.001792 NaN \n", + "3.5 0.002033 NaN \n", + "4.0 NaN NaN \n", + "4.5 NaN NaN \n", + "5.0 NaN NaN \n", + "5.5 NaN NaN \n", + "6.0 NaN NaN \n", + "6.5 NaN NaN \n", + "7.0 NaN NaN \n", + "7.5 NaN NaN \n", + "8.0 0.000351 NaN \n", + "8.5 0.000418 0.000405 \n", + "9.0 NaN NaN \n", + "9.5 0.000153 NaN \n", + "10.0 NaN 0.000281 \n", + "10.5 0.000204 0.000225 " ] }, "execution_count": 10, @@ -1458,7 +1507,7 @@ "
y_centers0.01.02.015.016.0
x_centers
\n", " \n", " \n", - " \n", + " \n", " \n", " \n", " \n", @@ -1477,6 +1526,26 @@ " \n", " \n", " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", " \n", " \n", " \n", @@ -1924,77 +1993,80 @@ "" ], "text/plain": [ - " 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 \\\n", - "0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "0.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.002692 0.001346 \n", - "1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.001346 0.006729 0.004038 \n", - "1.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.005384 0.002692 \n", - "2.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.002692 0.005384 \n", - "2.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.013459 \n", - "3.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.001346 0.021534 \n", - "3.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.006729 \n", - "4.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.009421 \n", - "4.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.016151 \n", - "5.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.002692 \n", - "5.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.001346 \n", - "6.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "6.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "7.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "7.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "8.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "8.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "9.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "9.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", - "10.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "y_centers 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 \\\n", + "x_centers \n", + "0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "0.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.002692 0.001346 \n", + "1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.001346 0.006729 0.004038 \n", + "1.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.005384 0.002692 \n", + "2.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.002692 0.005384 \n", + "2.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.013459 \n", + "3.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.001346 0.021534 \n", + "3.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.006729 \n", + "4.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.009421 \n", + "4.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.016151 \n", + "5.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.002692 \n", + "5.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.001346 \n", + "6.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "6.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "7.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "7.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "8.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "8.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "9.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "9.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "10.0 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", + "10.5 0.0 0.0 0.0 0.0 0.0 0.0 0.000000 0.000000 0.000000 \n", "\n", - " 9.0 10.0 11.0 12.0 13.0 14.0 15.0 \\\n", - "0.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", - "0.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", - "1.0 0.001346 0.000000 0.002692 0.002692 0.001346 0.000000 0.000000 \n", - "1.5 0.000000 0.000000 0.009421 0.004038 0.006729 0.005384 0.000000 \n", - "2.0 0.018843 0.018843 0.029610 0.021534 0.001346 0.002692 0.000000 \n", - "2.5 0.052490 0.055182 0.018843 0.025572 0.022880 0.005384 0.001346 \n", - "3.0 0.044415 0.047106 0.020188 0.012113 0.010767 0.010767 0.001346 \n", - "3.5 0.040377 0.029610 0.047106 0.004038 0.008075 0.004038 0.001346 \n", - "4.0 0.017497 0.029610 0.040377 0.002692 0.004038 0.005384 0.000000 \n", - "4.5 0.013459 0.017497 0.022880 0.012113 0.000000 0.000000 0.000000 \n", - "5.0 0.008075 0.008075 0.010767 0.022880 0.001346 0.000000 0.000000 \n", - "5.5 0.012113 0.006729 0.004038 0.014805 0.000000 0.000000 0.000000 \n", - "6.0 0.002692 0.002692 0.000000 0.005384 0.001346 0.000000 0.000000 \n", - "6.5 0.002692 0.002692 0.001346 0.000000 0.002692 0.000000 0.000000 \n", - "7.0 0.000000 0.000000 0.000000 0.001346 0.004038 0.000000 0.000000 \n", - "7.5 0.000000 0.001346 0.000000 0.000000 0.001346 0.008075 0.000000 \n", - "8.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.002692 0.002692 \n", - "8.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.001346 \n", - "9.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", - "9.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.001346 \n", - "10.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", - "10.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.001346 \n", + "y_centers 9.0 10.0 11.0 12.0 13.0 14.0 \\\n", + "x_centers \n", + "0.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "0.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "1.0 0.001346 0.000000 0.002692 0.002692 0.001346 0.000000 \n", + "1.5 0.000000 0.000000 0.009421 0.004038 0.006729 0.005384 \n", + "2.0 0.018843 0.018843 0.029610 0.021534 0.001346 0.002692 \n", + "2.5 0.052490 0.055182 0.018843 0.025572 0.022880 0.005384 \n", + "3.0 0.044415 0.047106 0.020188 0.012113 0.010767 0.010767 \n", + "3.5 0.040377 0.029610 0.047106 0.004038 0.008075 0.004038 \n", + "4.0 0.017497 0.029610 0.040377 0.002692 0.004038 0.005384 \n", + "4.5 0.013459 0.017497 0.022880 0.012113 0.000000 0.000000 \n", + "5.0 0.008075 0.008075 0.010767 0.022880 0.001346 0.000000 \n", + "5.5 0.012113 0.006729 0.004038 0.014805 0.000000 0.000000 \n", + "6.0 0.002692 0.002692 0.000000 0.005384 0.001346 0.000000 \n", + "6.5 0.002692 0.002692 0.001346 0.000000 0.002692 0.000000 \n", + "7.0 0.000000 0.000000 0.000000 0.001346 0.004038 0.000000 \n", + "7.5 0.000000 0.001346 0.000000 0.000000 0.001346 0.008075 \n", + "8.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.002692 \n", + "8.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "9.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "9.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "10.0 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", + "10.5 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 \n", "\n", - " 16.0 \n", - "0.0 0.000000 \n", - "0.5 0.000000 \n", - "1.0 0.000000 \n", - "1.5 0.000000 \n", - "2.0 0.000000 \n", - "2.5 0.000000 \n", - "3.0 0.000000 \n", - "3.5 0.000000 \n", - "4.0 0.000000 \n", - "4.5 0.000000 \n", - "5.0 0.000000 \n", - "5.5 0.000000 \n", - "6.0 0.000000 \n", - "6.5 0.000000 \n", - "7.0 0.000000 \n", - "7.5 0.000000 \n", - "8.0 0.000000 \n", - "8.5 0.001346 \n", - "9.0 0.000000 \n", - "9.5 0.000000 \n", - "10.0 0.001346 \n", - "10.5 0.001346 " + "y_centers 15.0 16.0 \n", + "x_centers \n", + "0.0 0.000000 0.000000 \n", + "0.5 0.000000 0.000000 \n", + "1.0 0.000000 0.000000 \n", + "1.5 0.000000 0.000000 \n", + "2.0 0.000000 0.000000 \n", + "2.5 0.001346 0.000000 \n", + "3.0 0.001346 0.000000 \n", + "3.5 0.001346 0.000000 \n", + "4.0 0.000000 0.000000 \n", + "4.5 0.000000 0.000000 \n", + "5.0 0.000000 0.000000 \n", + "5.5 0.000000 0.000000 \n", + "6.0 0.000000 0.000000 \n", + "6.5 0.000000 0.000000 \n", + "7.0 0.000000 0.000000 \n", + "7.5 0.000000 0.000000 \n", + "8.0 0.002692 0.000000 \n", + "8.5 0.001346 0.001346 \n", + "9.0 0.000000 0.000000 \n", + "9.5 0.001346 0.000000 \n", + "10.0 0.000000 0.001346 \n", + "10.5 0.001346 0.001346 " ] }, "execution_count": 11, @@ -2085,7 +2157,7 @@ "
y_centers0.01.02.015.016.0
x_centers
\n", " \n", " \n", - " \n", + " \n", " \n", " \n", " \n", @@ -2104,6 +2176,26 @@ " \n", " \n", " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", + " \n", " \n", " \n", " \n", @@ -2551,53 +2643,80 @@ "" ], "text/plain": [ - " 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 9.0 \\\n", - "0.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "0.5 NaN NaN NaN NaN NaN NaN NaN 224.996 117.594 NaN \n", - "1.0 NaN NaN NaN NaN NaN NaN 212.762 202.713 188.707 187.103 \n", - "1.5 NaN NaN NaN NaN NaN NaN NaN 176.402 199.802 NaN \n", - "2.0 NaN NaN NaN NaN NaN NaN NaN 203.667 216.857 192.965 \n", - "2.5 NaN NaN NaN NaN NaN NaN NaN NaN 193.397 203.529 \n", - "3.0 NaN NaN NaN NaN NaN NaN NaN 170.739 216.459 197.484 \n", - "3.5 NaN NaN NaN NaN NaN NaN NaN NaN 194.894 214.108 \n", - "4.0 NaN NaN NaN NaN NaN NaN NaN NaN 217.289 189.403 \n", - "4.5 NaN NaN NaN NaN NaN NaN NaN NaN 197.994 194.238 \n", - "5.0 NaN NaN NaN NaN NaN NaN NaN NaN 198.149 196.527 \n", - "5.5 NaN NaN NaN NaN NaN NaN NaN NaN 249.158 212.561 \n", - "6.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN 182.314 \n", - "6.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN 164.712 \n", - "7.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "7.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "8.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "8.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "9.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "9.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "10.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", - "10.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "y_centers 0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 \\\n", + "x_centers \n", + "0.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "0.5 NaN NaN NaN NaN NaN NaN NaN 224.996 117.594 \n", + "1.0 NaN NaN NaN NaN NaN NaN 212.762 202.713 188.707 \n", + "1.5 NaN NaN NaN NaN NaN NaN NaN 176.402 199.802 \n", + "2.0 NaN NaN NaN NaN NaN NaN NaN 203.667 216.857 \n", + "2.5 NaN NaN NaN NaN NaN NaN NaN NaN 193.397 \n", + "3.0 NaN NaN NaN NaN NaN NaN NaN 170.739 216.459 \n", + "3.5 NaN NaN NaN NaN NaN NaN NaN NaN 194.894 \n", + "4.0 NaN NaN NaN NaN NaN NaN NaN NaN 217.289 \n", + "4.5 NaN NaN NaN NaN NaN NaN NaN NaN 197.994 \n", + "5.0 NaN NaN NaN NaN NaN NaN NaN NaN 198.149 \n", + "5.5 NaN NaN NaN NaN NaN NaN NaN NaN 249.158 \n", + "6.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "6.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "7.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "7.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "8.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "8.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "9.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "9.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "10.0 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", + "10.5 NaN NaN NaN NaN NaN NaN NaN NaN NaN \n", "\n", - " 10.0 11.0 12.0 13.0 14.0 15.0 16.0 \n", - "0.0 NaN NaN NaN NaN NaN NaN NaN \n", - "0.5 NaN NaN NaN NaN NaN NaN NaN \n", - "1.0 NaN 213.926 174.154 164.886 NaN NaN NaN \n", - "1.5 NaN 201.883 191.598 221.705 190.124 NaN NaN \n", - "2.0 201.633 216.268 209.634 162.569 232.530 NaN NaN \n", - "2.5 196.907 212.883 211.277 202.760 199.263 272.421 NaN \n", - "3.0 200.895 212.107 193.837 222.185 169.497 122.296 NaN \n", - "3.5 202.725 206.901 184.099 186.077 221.659 186.201 NaN \n", - "4.0 201.362 207.532 207.971 172.771 213.854 NaN NaN \n", - "4.5 205.559 203.195 197.980 NaN NaN NaN NaN \n", - "5.0 222.219 215.221 204.002 254.004 NaN NaN NaN \n", - "5.5 212.734 168.655 208.220 NaN NaN NaN NaN \n", - "6.0 159.418 NaN 208.418 241.347 NaN NaN NaN \n", - "6.5 233.890 110.517 NaN 207.919 NaN NaN NaN \n", - "7.0 NaN NaN 155.691 229.022 NaN NaN NaN \n", - "7.5 166.855 NaN NaN 128.897 198.053 NaN NaN \n", - "8.0 NaN NaN NaN NaN 230.281 184.510 NaN \n", - "8.5 NaN NaN NaN NaN NaN 248.338 264.534 \n", - "9.0 NaN NaN NaN NaN NaN NaN NaN \n", - "9.5 NaN NaN NaN NaN NaN 116.230 NaN \n", - "10.0 NaN NaN NaN NaN NaN NaN 244.634 \n", - "10.5 NaN NaN NaN NaN NaN 190.849 212.411 " + "y_centers 9.0 10.0 11.0 12.0 13.0 14.0 15.0 \\\n", + "x_centers \n", + "0.0 NaN NaN NaN NaN NaN NaN NaN \n", + "0.5 NaN NaN NaN NaN NaN NaN NaN \n", + "1.0 187.103 NaN 213.926 174.154 164.886 NaN NaN \n", + "1.5 NaN NaN 201.883 191.598 221.705 190.124 NaN \n", + "2.0 192.965 201.633 216.268 209.634 162.569 232.530 NaN \n", + "2.5 203.529 196.907 212.883 211.277 202.760 199.263 272.421 \n", + "3.0 197.484 200.895 212.107 193.837 222.185 169.497 122.296 \n", + "3.5 214.108 202.725 206.901 184.099 186.077 221.659 186.201 \n", + "4.0 189.403 201.362 207.532 207.971 172.771 213.854 NaN \n", + "4.5 194.238 205.559 203.195 197.980 NaN NaN NaN \n", + "5.0 196.527 222.219 215.221 204.002 254.004 NaN NaN \n", + "5.5 212.561 212.734 168.655 208.220 NaN NaN NaN \n", + "6.0 182.314 159.418 NaN 208.418 241.347 NaN NaN \n", + "6.5 164.712 233.890 110.517 NaN 207.919 NaN NaN \n", + "7.0 NaN NaN NaN 155.691 229.022 NaN NaN \n", + "7.5 NaN 166.855 NaN NaN 128.897 198.053 NaN \n", + "8.0 NaN NaN NaN NaN NaN 230.281 184.510 \n", + "8.5 NaN NaN NaN NaN NaN NaN 248.338 \n", + "9.0 NaN NaN NaN NaN NaN NaN NaN \n", + "9.5 NaN NaN NaN NaN NaN NaN 116.230 \n", + "10.0 NaN NaN NaN NaN NaN NaN NaN \n", + "10.5 NaN NaN NaN NaN NaN NaN 190.849 \n", + "\n", + "y_centers 16.0 \n", + "x_centers \n", + "0.0 NaN \n", + "0.5 NaN \n", + "1.0 NaN \n", + "1.5 NaN \n", + "2.0 NaN \n", + "2.5 NaN \n", + "3.0 NaN \n", + "3.5 NaN \n", + "4.0 NaN \n", + "4.5 NaN \n", + "5.0 NaN \n", + "5.5 NaN \n", + "6.0 NaN \n", + "6.5 NaN \n", + "7.0 NaN \n", + "7.5 NaN \n", + "8.0 NaN \n", + "8.5 264.534 \n", + "9.0 NaN \n", + "9.5 NaN \n", + "10.0 244.634 \n", + "10.5 212.411 " ] }, "execution_count": 14, @@ -2636,8 +2755,20 @@ "name": "stdout", "output_type": "stream", "text": [ - "MAEP from timeseries = 1767087.527586333\n", - "MAEP from matrices = 1781210.8652839188\n" + "MAEP from timeseries = 1767087.5275863332\n" + ] + }, + { + "ename": "ValueError", + "evalue": "The truth value of a Series is ambiguous. Use a.empty, a.bool(), a.item(), a.any() or a.all().", + "output_type": "error", + "traceback": [ + "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m", + "\u001b[1;31mValueError\u001b[0m Traceback (most recent call last)", + "\u001b[1;32m~\\AppData\\Local\\Temp\\ipykernel_22536\\832768017.py\u001b[0m in \u001b[0;36m?\u001b[1;34m()\u001b[0m\n\u001b[0;32m 2\u001b[0m \u001b[0mmaep_timeseries\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mwave\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mperformance\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmean_annual_energy_production_timeseries\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mL\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mJ\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 3\u001b[0m \u001b[0mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"MAEP from timeseries = \"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mmaep_timeseries\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 4\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 5\u001b[0m \u001b[1;31m# Calcaulte maep from matrix\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 6\u001b[1;33m maep_matrix = wave.performance.mean_annual_energy_production_matrix(\n\u001b[0m\u001b[0;32m 7\u001b[0m \u001b[0mLM_mean\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mJM\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mLM_freq\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 8\u001b[0m )\n\u001b[0;32m 9\u001b[0m \u001b[0mprint\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"MAEP from matrices = \"\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mmaep_matrix\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32mc:\\users\\sterl\\codes\\mhkit-python\\mhkit\\wave\\performance.py\u001b[0m in \u001b[0;36m?\u001b[1;34m(LM, JM, frequency)\u001b[0m\n\u001b[0;32m 327\u001b[0m \u001b[0mfrequency\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mconvert_to_dataarray\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfrequency\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 328\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 329\u001b[0m \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mLM\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mshape\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0mJM\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mshape\u001b[0m \u001b[1;33m==\u001b[0m \u001b[0mfrequency\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mshape\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 330\u001b[0m \u001b[1;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"LM, JM, and frequency must be of the same size\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 331\u001b[1;33m \u001b[1;32mif\u001b[0m \u001b[1;32mnot\u001b[0m \u001b[0mnp\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mabs\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mfrequency\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0msum\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m-\u001b[0m \u001b[1;36m1\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m<\u001b[0m \u001b[1;36m1e-6\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m 332\u001b[0m \u001b[1;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34m\"Frequency components must sum to one.\"\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 333\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 334\u001b[0m \u001b[0mT\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;36m8766\u001b[0m \u001b[1;31m# Average length of a year (h)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n", + "\u001b[1;32mc:\\Users\\sterl\\anaconda3\\Lib\\site-packages\\pandas\\core\\generic.py\u001b[0m in \u001b[0;36m?\u001b[1;34m(self)\u001b[0m\n\u001b[0;32m 1517\u001b[0m \u001b[1;33m@\u001b[0m\u001b[0mfinal\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1518\u001b[0m \u001b[1;32mdef\u001b[0m \u001b[0m__nonzero__\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;33m->\u001b[0m \u001b[0mNoReturn\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 1519\u001b[1;33m raise ValueError(\n\u001b[0m\u001b[0;32m 1520\u001b[0m \u001b[1;34mf\"The truth value of a {type(self).__name__} is ambiguous. \"\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1521\u001b[0m \u001b[1;34m\"Use a.empty, a.bool(), a.item(), a.any() or a.all().\"\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m 1522\u001b[0m )\n", + "\u001b[1;31mValueError\u001b[0m: The truth value of a Series is ambiguous. Use a.empty, a.bool(), a.item(), a.any() or a.all()." ] } ], @@ -2646,10 +2777,15 @@ "maep_timeseries = wave.performance.mean_annual_energy_production_timeseries(L, J)\n", "print(\"MAEP from timeseries = \", maep_timeseries)\n", "\n", - "# Calcaulte maep from matrix\n", - "maep_matrix = wave.performance.mean_annual_energy_production_matrix(\n", - " LM_mean, JM, LM_freq\n", - ")\n", + "# Calcaulte maep from matrix \n", + "# See Issue #339\n", + "# maep_matrix = wave.performance.mean_annual_energy_production_matrix(\n", + "# LM_mean, JM, LM_freq\n", + "# )\n", + "\n", + "T = 8766 # Average length of a year (h)\n", + "maep_matrix = T * np.nansum(LM_mean * JM * LM_freq)\n", + "\n", "print(\"MAEP from matrices = \", maep_matrix)" ] }, @@ -2663,7 +2799,7 @@ }, { "cell_type": "code", - "execution_count": 16, + "execution_count": null, "metadata": {}, "outputs": [ { @@ -2694,7 +2830,7 @@ }, { "cell_type": "code", - "execution_count": 17, + "execution_count": null, "metadata": { "scrolled": true }, @@ -2755,7 +2891,7 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.8.8" + "version": "3.11.7" } }, "nbformat": 4,
y_centers0.01.02.015.016.0
x_centers