From 65f5f2b2c0d76a4d85228c4d6ef8c108d167dee5 Mon Sep 17 00:00:00 2001 From: Chris Leaman Date: Fri, 14 Dec 2018 12:02:35 +1100 Subject: [PATCH] Update notebook --- notebooks/01_exploration.ipynb | 1925 ++++++++++++++++++-------------- 1 file changed, 1091 insertions(+), 834 deletions(-) diff --git a/notebooks/01_exploration.ipynb b/notebooks/01_exploration.ipynb index deb0348..630f192 100644 --- a/notebooks/01_exploration.ipynb +++ b/notebooks/01_exploration.ipynb @@ -8,14 +8,24 @@ "This notebook provides an example how the data has been loaded and accessed for further analysis." ] }, + { + "cell_type": "markdown", + "metadata": { + "heading_collapsed": true + }, + "source": [ + "## Setup notebook" + ] + }, { "cell_type": "code", - "execution_count": 1, + "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2018-12-03T23:02:22.927101Z", - "start_time": "2018-12-03T23:02:22.612233Z" - } + "end_time": "2018-12-11T22:16:57.761856Z", + "start_time": "2018-12-11T22:16:55.841860Z" + }, + "hidden": true }, "outputs": [], "source": [ @@ -29,12 +39,13 @@ }, { "cell_type": "code", - "execution_count": 2, + "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2018-12-03T23:02:24.527369Z", - "start_time": "2018-12-03T23:02:22.929088Z" + "end_time": "2018-12-11T22:17:15.644076Z", + "start_time": "2018-12-11T22:16:57.763862Z" }, + "hidden": true, "scrolled": true }, "outputs": [], @@ -67,16 +78,17 @@ "cell_type": "markdown", "metadata": {}, "source": [ + "## Import data\n", "Import our data into pandas Dataframes for the analysis. Data files are `.csv` files which are stored in the `./data/interim/` folder." ] }, { "cell_type": "code", - "execution_count": 3, + "execution_count": 56, "metadata": { "ExecuteTime": { - "end_time": "2018-12-03T23:02:39.868010Z", - "start_time": "2018-12-03T23:02:24.529339Z" + "end_time": "2018-12-13T02:08:52.148080Z", + "start_time": "2018-12-13T02:08:29.576649Z" }, "pixiedust": { "displayParams": {} @@ -149,16 +161,16 @@ ] }, { - "cell_type": "code", - "execution_count": null, + "cell_type": "markdown", "metadata": { "ExecuteTime": { "end_time": "2018-11-27T23:02:57.631306Z", "start_time": "2018-11-27T23:02:57.615263Z" } }, - "outputs": [], - "source": [] + "source": [ + "## Profile/timeseries dashboard" + ] }, { "cell_type": "markdown", @@ -171,23 +183,22 @@ }, { "cell_type": "code", - "execution_count": 15, + "execution_count": 79, "metadata": { "ExecuteTime": { - "end_time": "2018-12-05T03:57:14.533063Z", - "start_time": "2018-12-05T03:57:13.745017Z" + "end_time": "2018-12-13T22:22:40.283477Z", + "start_time": "2018-12-13T22:22:39.430903Z" }, "code_folding": [ 408 ], - "hide_input": false, - "scrolled": false + "hide_input": false }, "outputs": [ { "data": { "application/vnd.jupyter.widget-view+json": { - "model_id": "eae0e9440a5f45599b2c9b43352d3d13", + "model_id": "47ba792667fd42f9b19c2fa574d8667c", "version_major": 2, "version_minor": 0 }, @@ -216,10 +227,10 @@ "for forecast in impacts['forecasted']:\n", " selectboxes.append(\n", " widgets.SelectMultiple(\n", - " options=impacts['forecasted'][\n", - " forecast].storm_regime.dropna().unique().tolist(),\n", - " value=impacts['forecasted'][forecast].storm_regime.dropna()\n", - " .unique().tolist(),\n", + " options=impacts['forecasted'][forecast].storm_regime.dropna().\n", + " unique().tolist(),\n", + " value=impacts['forecasted'][forecast].storm_regime.dropna().\n", + " unique().tolist(),\n", " disabled=False))\n", " titles.append('Forecasted: {}'.format(forecast))\n", "\n", @@ -237,8 +248,10 @@ "site_id_select = widgets.Dropdown(\n", " description='site_id: ',\n", " value='NARRA0001',\n", - " options=df_profiles.index.get_level_values('site_id').unique()\n", - " .sort_values().tolist())\n", + " # options=df_profiles.index.get_level_values('site_id').unique().\n", + " # sort_values().tolist(),\n", + " options=df_no_crests.index.tolist(),\n", + ")\n", "\n", "site_id_impacts = widgets.HTML(value=\"\", )\n", "\n", @@ -308,7 +321,6 @@ " symbol='square'),\n", ")\n", "\n", - "\n", "forecast_traces = []\n", "for forecast, color in zip(impacts['forecasted'], colors):\n", " forecast_traces.append(\n", @@ -343,7 +355,8 @@ " range=[-1, 20]))\n", "\n", "g_profiles = go.FigureWidget(\n", - " data=[trace1, trace2, trace3, trace4, trace5, trace6] + forecast_traces, layout=layout)\n", + " data=[trace1, trace2, trace3, trace4, trace5, trace6] + forecast_traces,\n", + " layout=layout)\n", "\n", "# Add panel for google maps\n", "mapbox_access_token = 'pk.eyJ1IjoiY2hyaXNsZWFtYW4iLCJhIjoiY2pvNTY1MzZpMDc2OTN2bmw5MGsycHp5bCJ9.U2dwFg2c7RFjUNSayERUiw'\n", @@ -537,14 +550,13 @@ " prestorm_dune_toe_x = prestorm_site_features.dune_toe_x\n", " prestorm_dune_toe_z = prestorm_site_features.dune_toe_z\n", "\n", - " \n", " poststorm_site_features = df_profile_features.query(\n", " 'site_id == \"{}\" and profile_type==\"poststorm\"'.format(site_id))\n", " poststorm_dune_crest_x = poststorm_site_features.dune_crest_x\n", " poststorm_dune_crest_z = poststorm_site_features.dune_crest_z\n", " poststorm_dune_toe_x = poststorm_site_features.dune_toe_x\n", " poststorm_dune_toe_z = poststorm_site_features.dune_toe_z\n", - " \n", + "\n", " # Update beach profile section plots\n", " with g_profiles.batch_update():\n", " g_profiles.data[0].x = prestorm_x\n", @@ -559,7 +571,7 @@ " g_profiles.data[4].y = poststorm_dune_crest_z\n", " g_profiles.data[5].x = poststorm_dune_toe_x\n", " g_profiles.data[5].y = poststorm_dune_toe_z\n", - " \n", + "\n", " for n, forecast in enumerate(impacts['forecasted']):\n", " R_high = max(impacts['forecasted'][forecast].query(\n", " \"site_id=='{}'\".format(site_id)).R_high)\n", @@ -612,10 +624,10 @@ " g_timeseries.data[3].x = [min(times), max(times)]\n", " g_timeseries.data[4].x = df_tide_site.index.get_level_values(\n", " 'datetime')\n", - " g_timeseries.data[2].y = prestorm_dune_crest_z.tolist()[\n", - " 0], prestorm_dune_crest_z.tolist()[0],\n", - " g_timeseries.data[3].y = prestorm_dune_toe_z.tolist()[0], prestorm_dune_toe_z.tolist()[\n", - " 0],\n", + " g_timeseries.data[2].y = prestorm_dune_crest_z.tolist(\n", + " )[0], prestorm_dune_crest_z.tolist()[0],\n", + " g_timeseries.data[3].y = prestorm_dune_toe_z.tolist(\n", + " )[0], prestorm_dune_toe_z.tolist()[0],\n", " g_timeseries.data[4].y = df_tide_site.tide.tolist()\n", "\n", " # Update rhigh values\n", @@ -633,14 +645,12 @@ " # Update site id impacts\n", " observed_regime = impacts['observed'].query(\n", " \"site_id=='{}'\".format(site_id)).storm_regime.values[0]\n", - " site_id_impacts.value = \"Observed: {}
\".format(\n", - " observed_regime)\n", + " site_id_impacts.value = \"Observed: {}
\".format(observed_regime)\n", "\n", " for forecast in impacts['forecasted']:\n", " regime = impacts['forecasted'][forecast].query(\n", " \"site_id=='{}'\".format(site_id)).storm_regime.values[0]\n", - " site_id_impacts.value += '{}: {}
'.format(\n", - " forecast, regime)\n", + " site_id_impacts.value += '{}: {}
'.format(forecast, regime)\n", "\n", "\n", "site_id_select.observe(update_profile, names=\"value\")\n", @@ -692,41 +702,29 @@ "end_time": "2018-11-22T22:52:36.039701Z", "start_time": "2018-11-22T22:52:36.035189Z" }, + "heading_collapsed": true, "hide_input": true, "scrolled": true }, "source": [ + "## Confusion matrix\n", "This visualization looks at how well the storm impact predictions performed. " ] }, { "cell_type": "code", - "execution_count": 5, + "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2018-12-03T23:02:42.021445Z", - "start_time": "2018-12-03T23:02:41.468637Z" + "end_time": "2018-12-11T22:41:46.514495Z", + "start_time": "2018-12-11T22:41:27.050029Z" }, "code_folding": [], + "hidden": true, "hide_input": false, "scrolled": false }, - "outputs": [ - { - "data": { - "application/vnd.jupyter.widget-view+json": { - "model_id": "a2b52c5ad861454db1756a427f13b55d", - "version_major": 2, - "version_minor": 0 - }, - "text/plain": [ - "VBox(children=(VBox(children=(HTML(value='Filter by beach:'), SelectMultiple(index=(0, 1, 2, 3, 4, 5, 6…" - ] - }, - "metadata": {}, - "output_type": "display_data" - } - ], + "outputs": [], "source": [ "# Create colorscale\n", "rdylgr_cmap = matplotlib.cm.get_cmap('RdYlGn')\n", @@ -829,12 +827,13 @@ }, { "cell_type": "code", - "execution_count": 8, + "execution_count": null, "metadata": { "ExecuteTime": { - "end_time": "2018-12-03T23:03:48.457193Z", - "start_time": "2018-12-03T23:03:48.274709Z" - } + "end_time": "2018-12-11T22:17:29.215527Z", + "start_time": "2018-12-11T22:16:56.023Z" + }, + "hidden": true }, "outputs": [], "source": [ @@ -855,168 +854,69 @@ } }, "source": [ - "### Look at time dependance" + "## Identify sites with no results" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### Check forecast TWL\n", + "Most probable explanation for TWL's which are NaN'ed is that the prestorm dune toe are not defined." ] }, { "cell_type": "code", - "execution_count": 11, + "execution_count": 57, "metadata": { "ExecuteTime": { - "end_time": "2018-12-03T23:49:16.581105Z", - "start_time": "2018-12-03T23:49:16.274275Z" + "end_time": "2018-12-13T02:09:10.914191Z", + "start_time": "2018-12-13T02:09:09.913622Z" } }, "outputs": [ { - "data": { - "application/vnd.jupyter.widget-view+json": { - "model_id": "8a366d02e4564347a5950b0f24c86363", - "version_major": 2, - "version_minor": 0 - }, - "text/plain": [ - "VBox(children=(HBox(children=(VBox(children=(HTML(value='Filter by site_id:'), HBox(children=(Dropdown(…" - ] - }, - "metadata": {}, - "output_type": "display_data" + "name": "stdout", + "output_type": "stream", + "text": [ + "The following sites have no slope defined in the twl csv file:\n", + "['ENTRA0078', 'ENTRA0079', 'MANNING0109']\n", + "\n", + "The following sites have no R_high defined in the twl csv file:\n", + "['ENTRA0078', 'ENTRA0079', 'MANNING0109']\n" + ] } ], "source": [ - "# Create widgets for selecting site_id\n", - "site_id_title = widgets.HTML(value=\"Filter by site_id:\", )\n", - "\n", - "site_id_select = widgets.Dropdown(\n", - " description='site_id: ',\n", - " value='NARRA0001',\n", - " options=df_profiles.index.get_level_values('site_id').unique()\n", - " .sort_values().tolist())\n", - "\n", - "site_id_impacts = widgets.HTML(value=\"\", )\n", - "\n", - "site_id_container = widgets.HBox(children=[\n", - " widgets.VBox(\n", - " children=[site_id_title,\n", - " widgets.HBox(children=[site_id_select])]), site_id_impacts\n", - "])\n", - "\n", - "\n", - "# Plot profiles\n", - "\n", - "# Add panel for pre/post storm profiles\n", - "trace1 = go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Pre Storm Profile',\n", - " line=dict(color=('rgb(51,160,44)'), width=2))\n", - "trace2 = go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Post Storm Profile',\n", - " line=dict(color=('rgb(255,127,0)'), width=2))\n", - "trace3 = go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Pre-storm dune crest',\n", - " mode='markers',\n", - " marker=dict(\n", - " color='rgba(255,255,255,0)',\n", - " size=10,\n", - " line=dict(color='rgba(106,61,154, 1)', width=2)),\n", - ")\n", - "trace4 = go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Pre-storm dune toe',\n", - " mode='markers',\n", - " marker=dict(\n", - " color='rgba(255,255,255,0)',\n", - " size=10,\n", - " line=dict(color='rgba(202,178,214,1)', width=2)),\n", - ")\n", - "\n", - "trace5 = go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Post-storm dune crest',\n", - " mode='markers',\n", - " marker=dict(\n", - " color='rgba(255,255,255,0)',\n", - " size=10,\n", - " line=dict(color='rgba(106,61,154, 1)', width=2),\n", - " symbol='square'),\n", - ")\n", - "trace6 = go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Post-storm dune toe',\n", - " mode='markers',\n", - " marker=dict(\n", - " color='rgba(255,255,255,0)',\n", - " size=10,\n", - " line=dict(color='rgba(202,178,214,1)', width=2),\n", - " symbol='square'),\n", - ")\n", - "\n", - "\n", - "forecast_traces = []\n", - "for forecast, color in zip(impacts['forecasted'], colors):\n", - " forecast_traces.append(\n", - " go.Scatter(\n", - " x=[0],\n", - " y=[0],\n", - " name='Peak R_high: {}'.format(forecast),\n", - " mode='lines',\n", - " line=dict(\n", - " color=color,\n", - " width=4,\n", - " )))\n", - "\n", - "layout = go.Layout(\n", - " title='Bed Profiles',\n", - " height=300,\n", - " legend=dict(font={'size': 10}),\n", - " margin=dict(t=50, b=50, l=50, r=20),\n", - " xaxis=dict(\n", - " title='x (m)',\n", - " autorange=True,\n", - " showgrid=True,\n", - " zeroline=True,\n", - " showline=True,\n", - " range=[0, 200]),\n", - " yaxis=dict(\n", - " title='z (m)',\n", - " autorange=False,\n", - " showgrid=True,\n", - " zeroline=True,\n", - " showline=True,\n", - " range=[-1, 20]))\n", + "df_twls = twls['forecasted']['mean_slope_sto06']\n", "\n", - "g_profiles = go.FigureWidget(\n", - " data=[trace1, trace2, trace3, trace4, trace5, trace6] + forecast_traces, layout=layout)\n", + "slope_mask = df_twls.groupby('site_id').agg({'beta': lambda x: x.isnull().sum() == len(x)}).beta\n", + "print('The following sites have no slope defined in the twl csv file:')\n", + "print(slope_mask.index[slope_mask].tolist())\n", + "print()\n", "\n", - "widgets.VBox([\n", - " site_id_container,\n", - " g_profiles\n", - "])" + "R_high_mask = df_twls.groupby('site_id').agg({'R_high': lambda x: x.isnull().sum() == len(x)}).R_high\n", + "print('The following sites have no R_high defined in the twl csv file:')\n", + "print(slope_mask.index[slope_mask].tolist())\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ - "## Scatter plot" + "### Check observed impacts\n", + "Find sites which have no observed impacts. If we do not identify an observed storm regime, the site cannot be included when we're trying to compare predicted and observed impacts." ] }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 58, "metadata": { "ExecuteTime": { - "end_time": "2018-12-07T05:08:12.117885Z", - "start_time": "2018-12-07T05:08:12.078780Z" - } + "end_time": "2018-12-13T02:09:14.785919Z", + "start_time": "2018-12-13T02:09:14.520986Z" + }, + "scrolled": false }, "outputs": [ { @@ -1065,107 +965,83 @@ " \n", " \n", " \n", - " AVOCAn0001\n", - " 113.9091\n", - " 77.9830\n", - " 35.6103\n", - " 31.2620\n", - " 165.4760\n", - " 166.7296\n", - " 0.0000\n", - " 0.0000\n", - " swash\n", - " \n", - " \n", - " AVOCAn0002\n", - " 106.8959\n", - " 67.0913\n", - " 39.6373\n", - " 37.0803\n", - " 256.4137\n", - " 258.9174\n", - " 0.0000\n", - " 0.0000\n", - " swash\n", - " \n", - " \n", - " AVOCAn0003\n", - " 99.0484\n", - " 53.6563\n", - " 45.1621\n", - " 45.5960\n", - " 372.7031\n", - " 373.9198\n", - " -0.3147\n", - " -0.0844\n", - " swash\n", + " AVOCAn0009\n", + " 4.5783\n", + " 0.1110\n", + " 4.4673\n", + " 97.5750\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " AVOCAn0004\n", - " 74.7543\n", - " 29.3280\n", - " 45.4262\n", - " 60.7674\n", - " 275.1689\n", - " 276.0476\n", - " 0.4104\n", - " 0.1492\n", - " swash\n", + " AVOCAs0001\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " AVOCAn0005\n", - " 70.5968\n", - " 24.1071\n", - " 46.4897\n", - " 65.8524\n", - " 268.5194\n", - " 263.4262\n", - " 7.4196\n", - " 2.7631\n", - " collision\n", + " AVOCAs0002\n", + " 97.9463\n", + " 26.6638\n", + " 71.2825\n", + " 72.7771\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " AVOCAn0006\n", - " 68.7583\n", - " 23.9665\n", - " 44.7918\n", - " 65.1438\n", - " 202.6770\n", - " 198.2397\n", - " 4.7944\n", - " 2.3655\n", - " collision\n", + " AVOCAs0003\n", + " 70.7306\n", + " 40.2020\n", + " 30.7232\n", + " 43.4369\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " AVOCAn0007\n", - " 75.7895\n", - " 27.2715\n", - " 48.5180\n", - " 64.0168\n", - " 149.8479\n", - " 143.1312\n", - " 7.1323\n", - " 4.7597\n", - " collision\n", + " AVOCAs0004\n", + " 98.2859\n", + " 45.4986\n", + " 52.6330\n", + " 53.5509\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " AVOCAn0008\n", - " 93.3107\n", - " 42.4968\n", - " 50.8139\n", - " 54.4567\n", - " 187.9201\n", - " 187.3459\n", - " 2.8297\n", - " 1.5058\n", - " collision\n", + " AVOCAs0005\n", + " 95.5841\n", + " 54.9753\n", + " 40.5733\n", + " 42.4478\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " AVOCAn0009\n", - " 3.6955\n", - " 0.1038\n", - " 3.5917\n", - " 97.1908\n", + " AVOCAs0006\n", + " 113.0441\n", + " 67.8912\n", + " 45.2582\n", + " 40.0359\n", " NaN\n", " NaN\n", " NaN\n", @@ -1173,11 +1049,23 @@ " NaN\n", " \n", " \n", - " AVOCAs0001\n", + " AVOCAs0007\n", + " 65.3283\n", + " 44.2821\n", + " 21.4544\n", + " 32.8409\n", + " NaN\n", " NaN\n", " NaN\n", " NaN\n", " NaN\n", + " \n", + " \n", + " AVOCAs0008\n", + " 52.3933\n", + " 45.2243\n", + " 7.1728\n", + " 13.6904\n", " NaN\n", " NaN\n", " NaN\n", @@ -1185,11 +1073,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0002\n", - " 97.3864\n", - " 26.6619\n", - " 70.7246\n", - " 72.6226\n", + " BILG0001\n", + " 20.3405\n", + " 7.6207\n", + " 12.7198\n", + " 62.5344\n", " NaN\n", " NaN\n", " NaN\n", @@ -1197,11 +1085,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0003\n", - " 70.7401\n", - " 40.0608\n", - " 30.7919\n", - " 43.5282\n", + " BILG0002\n", + " 156.4205\n", + " 98.1716\n", + " 58.1659\n", + " 37.1856\n", " NaN\n", " NaN\n", " NaN\n", @@ -1209,11 +1097,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0004\n", - " 97.8389\n", - " 45.4845\n", - " 52.2157\n", - " 53.3691\n", + " BOAT0001\n", + " 23.8361\n", + " 23.6865\n", + " -0.0926\n", + " -0.3885\n", " NaN\n", " NaN\n", " NaN\n", @@ -1221,11 +1109,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0005\n", - " 95.1711\n", - " 54.9722\n", - " 40.1706\n", - " 42.2088\n", + " BOAT0002\n", + " 38.8398\n", + " 14.0819\n", + " 24.7579\n", + " 63.7436\n", " NaN\n", " NaN\n", " NaN\n", @@ -1233,11 +1121,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0006\n", - " 112.5818\n", - " 67.8718\n", - " 44.8252\n", - " 39.8157\n", + " BOAT0003\n", + " 73.6809\n", + " 17.8545\n", + " 55.8264\n", + " 75.7678\n", " NaN\n", " NaN\n", " NaN\n", @@ -1245,11 +1133,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0007\n", - " 65.3531\n", - " 44.1537\n", - " 21.5228\n", - " 32.9331\n", + " BOAT0004\n", + " 73.1954\n", + " 23.1583\n", + " 50.0372\n", + " 68.3610\n", " NaN\n", " NaN\n", " NaN\n", @@ -1257,11 +1145,11 @@ " NaN\n", " \n", " \n", - " AVOCAs0008\n", - " 52.3940\n", - " 44.9152\n", - " 7.4803\n", - " 14.2770\n", + " BOAT0005\n", + " 53.5122\n", + " 22.4537\n", + " 31.0585\n", + " 58.0400\n", " NaN\n", " NaN\n", " NaN\n", @@ -1269,11 +1157,11 @@ " NaN\n", " \n", " \n", - " BILG0001\n", - " 19.4177\n", - " 7.5746\n", - " 11.8431\n", - " 60.9913\n", + " BOOM0001\n", + " 236.4540\n", + " 218.4918\n", + " 20.8725\n", + " 8.8273\n", " NaN\n", " NaN\n", " NaN\n", @@ -1281,11 +1169,23 @@ " NaN\n", " \n", " \n", - " BILG0002\n", - " 155.6998\n", - " 98.1693\n", - " 57.4340\n", - " 36.8876\n", + " CATHIE0024\n", + " 63.6452\n", + " 38.4261\n", + " 25.2191\n", + " 39.6245\n", + " 0.0\n", + " 0.0\n", + " 0.0\n", + " NaN\n", + " NaN\n", + " \n", + " \n", + " CATHIE0026\n", + " 75.1334\n", + " 43.7179\n", + " 31.0940\n", + " 41.3851\n", " NaN\n", " NaN\n", " NaN\n", @@ -1293,119 +1193,119 @@ " NaN\n", " \n", " \n", - " BILG0003\n", - " 83.5219\n", - " 52.8059\n", - " 30.5537\n", - " 36.5817\n", - " 41.1469\n", - " 40.6081\n", - " 0.0000\n", - " 0.0000\n", - " swash\n", + " CRESn0069\n", + " 37.5896\n", + " 8.3495\n", + " 29.2401\n", + " 77.7877\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BILG0004\n", - " 158.6283\n", - " 116.6189\n", - " 42.1178\n", - " 26.5512\n", - " 11.2211\n", - " 11.0892\n", - " -0.0132\n", - " -0.1179\n", - " swash\n", + " DEEWHYn0008\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BILG0005\n", - " 212.8478\n", - " 164.0044\n", - " 48.4312\n", - " 22.7539\n", - " 613.6156\n", - " 606.2766\n", - " 5.7738\n", - " 0.9410\n", - " collision\n", + " DEEWHYn0009\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BLUEYS0001\n", - " 65.4628\n", - " 19.2938\n", - " 46.1690\n", - " 70.5270\n", - " 130.7447\n", - " 120.5446\n", - " 9.5601\n", - " 7.3121\n", - " collision\n", + " DEEWHYs0005\n", + " 62.3514\n", + " 24.9797\n", + " 37.3716\n", + " 59.9372\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BLUEYS0002\n", - " 50.2084\n", - " 10.3009\n", - " 39.9074\n", - " 79.4836\n", - " 512.0154\n", - " 477.1774\n", - " 33.2825\n", - " 6.5003\n", - " collision\n", + " DEEWHYs0008\n", + " 1.0688\n", + " 1.3640\n", + " 0.0000\n", + " 0.0000\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BLUEYS0003\n", - " 50.6308\n", - " 11.1682\n", - " 39.4625\n", - " 77.9418\n", - " 443.0853\n", - " 414.3901\n", - " 24.8870\n", - " 5.6167\n", - " collision\n", + " DIAMONDn0023\n", + " 67.9416\n", + " 21.1812\n", + " 46.7603\n", + " 68.8244\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BLUEYS0004\n", - " 95.1608\n", - " 31.3330\n", - " 63.8279\n", - " 67.0737\n", - " 287.5805\n", - " 272.4267\n", - " 12.9641\n", - " 4.5080\n", - " collision\n", + " DIAMONDs0006\n", + " 74.9357\n", + " 42.4382\n", + " 32.2536\n", + " 43.0416\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BLUEYS0005\n", - " 141.0643\n", - " 58.2545\n", - " 82.8098\n", - " 58.7036\n", - " 539.3864\n", - " 520.0732\n", - " 12.0470\n", - " 2.2335\n", - " collision\n", + " DIAMONDs0007\n", + " 153.7639\n", + " 127.9469\n", + " 26.1595\n", + " 17.0128\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BLUEYS0006\n", - " 88.4207\n", - " 51.6205\n", - " 36.2553\n", - " 41.0032\n", - " 271.6036\n", - " 267.1954\n", - " 3.6045\n", - " 1.3271\n", - " collision\n", + " DUNBn0031\n", + " 36.5301\n", + " 6.3289\n", + " 30.2012\n", + " 82.6748\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " BOAT0001\n", - " 23.8510\n", - " 23.5660\n", - " -0.0264\n", - " -0.1108\n", + " DUNBn0055\n", + " 189.5283\n", + " 134.2760\n", + " 56.7139\n", + " 29.9237\n", " NaN\n", " NaN\n", " NaN\n", @@ -1413,11 +1313,11 @@ " NaN\n", " \n", " \n", - " BOAT0002\n", - " 37.6524\n", - " 14.0209\n", - " 23.6316\n", - " 62.7624\n", + " ELIZA0002\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " NaN\n", " NaN\n", " NaN\n", @@ -1437,215 +1337,215 @@ " ...\n", " \n", " \n", - " TREACH0014\n", - " 97.5323\n", - " 46.2994\n", - " 51.0816\n", - " 52.3740\n", - " 508.7400\n", - " 505.7877\n", - " 0.4254\n", - " 0.0836\n", - " swash\n", - " \n", - " \n", - " TREACH0015\n", - " 96.8327\n", - " 45.1962\n", - " 51.6364\n", - " 53.3254\n", - " 690.8275\n", - " 683.4458\n", - " 1.5086\n", - " 0.2184\n", + " STOCNs0170\n", + " 198.2785\n", + " 216.6368\n", + " -18.6067\n", + " -9.3841\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " NaN\n", " \n", " \n", - " TREACH0016\n", - " 106.9083\n", - " 66.0567\n", - " 40.3629\n", - " 37.7547\n", - " 508.0014\n", - " 499.6315\n", - " 0.3386\n", - " 0.0667\n", - " swash\n", + " STOCNs0175\n", + " 136.2126\n", + " 114.1715\n", + " 22.4984\n", + " 16.5171\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0001\n", - " 132.3413\n", - " 143.4459\n", - " -9.7255\n", - " -7.3488\n", - " 665.9898\n", - " 667.5923\n", - " 0.0410\n", - " 0.0062\n", - " swash\n", + " STOCNs0179\n", + " 67.7795\n", + " 45.3981\n", + " 22.3815\n", + " 33.0210\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0002\n", - " 151.1833\n", - " 126.9844\n", - " 23.9548\n", - " 15.8449\n", - " 385.8467\n", - " 386.7284\n", - " -0.0449\n", - " -0.0116\n", - " swash\n", + " STOCNs0180\n", + " 166.0813\n", + " 149.5195\n", + " 15.3441\n", + " 9.2389\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0003\n", - " 154.1788\n", - " 117.9441\n", - " 36.2425\n", - " 23.5068\n", - " 694.2226\n", - " 700.5105\n", - " -4.2136\n", - " -0.6070\n", - " swash\n", + " STOCNs0181\n", + " 90.1147\n", + " 98.9808\n", + " -9.1107\n", + " -10.1102\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0004\n", - " 137.8449\n", - " 76.6007\n", - " 61.2725\n", - " 44.4503\n", - " 559.5485\n", - " 569.8591\n", - " -4.4590\n", - " -0.7969\n", - " swash\n", + " STOCNs0182\n", + " 67.8622\n", + " 86.0118\n", + " -18.1671\n", + " -26.7705\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0005\n", + " STOCNs0183\n", + " 125.9085\n", + " 137.7342\n", + " -12.6233\n", + " -10.0257\n", " NaN\n", " NaN\n", " NaN\n", " NaN\n", " NaN\n", + " \n", + " \n", + " STOCNs0184\n", + " 146.6586\n", + " 123.5371\n", + " 23.2603\n", + " 15.8602\n", + " NaN\n", " NaN\n", " NaN\n", " NaN\n", " NaN\n", " \n", " \n", - " WAMBE0006\n", - " 205.8453\n", - " 186.0784\n", - " 22.5892\n", - " 10.9739\n", - " 55.0898\n", - " 55.8919\n", - " 0.0000\n", - " 0.0000\n", - " swash\n", + " STOCNs0185\n", + " 141.5421\n", + " 142.7279\n", + " -1.2619\n", + " -0.8915\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0007\n", - " 80.4674\n", - " 35.4614\n", - " 45.0059\n", - " 55.9307\n", - " 178.1005\n", - " 178.5439\n", - " 0.4727\n", - " 0.2654\n", - " swash\n", + " STOCNs0186\n", + " 115.9148\n", + " 123.1507\n", + " -7.8392\n", + " -6.7629\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0008\n", - " 88.4574\n", - " 40.3200\n", - " 48.1374\n", - " 54.4187\n", - " 258.7513\n", - " 258.3849\n", - " -1.2073\n", - " -0.4666\n", - " swash\n", + " STOCNs0187\n", + " 126.5519\n", + " 147.4371\n", + " -22.4452\n", + " -17.7359\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0009\n", - " 70.9159\n", - " 26.1742\n", - " 44.7418\n", - " 63.0913\n", - " 267.3725\n", - " 258.3720\n", - " 9.8041\n", - " 3.6668\n", - " collision\n", + " STOCNs0188\n", + " 345.5234\n", + " 353.8766\n", + " -11.3322\n", + " -3.2797\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0010\n", - " 58.6604\n", - " 18.0418\n", - " 40.6186\n", - " 69.2437\n", - " 187.5259\n", - " 161.9748\n", - " 25.3087\n", - " 13.4961\n", - " collision\n", + " STOCNs0189\n", + " 171.6354\n", + " 134.9192\n", + " 35.8697\n", + " 20.8988\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0011\n", - " 59.2415\n", - " 16.3165\n", - " 42.9250\n", - " 72.4577\n", - " 197.0129\n", - " 175.2512\n", - " 21.9882\n", - " 11.1608\n", - " collision\n", + " STOCNs0190\n", + " 151.4113\n", + " 116.6381\n", + " 35.0161\n", + " 23.1264\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0012\n", - " 74.4189\n", - " 23.0232\n", - " 51.3957\n", - " 69.0627\n", - " 178.4783\n", - " 168.3475\n", - " 10.0386\n", - " 5.6246\n", - " collision\n", + " STOCS0014\n", + " 98.4991\n", + " 57.3495\n", + " 40.9006\n", + " 41.5238\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0013\n", - " 70.4964\n", - " 22.7546\n", - " 47.7419\n", - " 67.7224\n", - " 231.1513\n", - " 195.2581\n", - " 35.8072\n", - " 15.4908\n", - " collision\n", + " STOCS0043\n", + " 36.4256\n", + " 11.7208\n", + " 24.7048\n", + " 67.8225\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", - " WAMBE0014\n", - " 68.0896\n", - " 24.1853\n", - " 43.9043\n", - " 64.4802\n", - " 82.4268\n", - " 61.2601\n", - " 21.1718\n", - " 25.6856\n", - " collision\n", + " WAMBE0005\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", + " NaN\n", " \n", " \n", " WAMBE0015\n", - " 55.0780\n", - " 16.0119\n", - " 39.0660\n", - " 70.9286\n", + " 56.2724\n", + " 16.0428\n", + " 40.2296\n", + " 71.4908\n", " NaN\n", " NaN\n", " NaN\n", @@ -1654,10 +1554,10 @@ " \n", " \n", " WAMBE0016\n", - " 96.7687\n", - " 39.8224\n", - " 56.9463\n", - " 58.8479\n", + " 97.8849\n", + " 39.8432\n", + " 58.0417\n", + " 59.2958\n", " NaN\n", " NaN\n", " NaN\n", @@ -1666,10 +1566,10 @@ " \n", " \n", " WAMBE0017\n", - " 35.2987\n", - " 8.5140\n", - " 26.7847\n", - " 75.8801\n", + " 36.5683\n", + " 8.5380\n", + " 28.0303\n", + " 76.6520\n", " NaN\n", " NaN\n", " NaN\n", @@ -1678,10 +1578,10 @@ " \n", " \n", " WAMBE0018\n", - " 40.9407\n", - " 10.5147\n", - " 30.4260\n", - " 74.3173\n", + " 42.1423\n", + " 10.5498\n", + " 31.5925\n", + " 74.9662\n", " NaN\n", " NaN\n", " NaN\n", @@ -1690,10 +1590,10 @@ " \n", " \n", " WAMBE0019\n", - " 38.2838\n", - " 9.2156\n", - " 29.0682\n", - " 75.9282\n", + " 39.6097\n", + " 9.2404\n", + " 30.3693\n", + " 76.6714\n", " NaN\n", " NaN\n", " NaN\n", @@ -1726,10 +1626,10 @@ " \n", " \n", " WAMBE0022\n", - " 0.5516\n", - " 0.2840\n", - " 0.2675\n", - " 48.5063\n", + " 1.1034\n", + " 0.4478\n", + " 0.6556\n", + " 59.4166\n", " NaN\n", " NaN\n", " NaN\n", @@ -1738,10 +1638,10 @@ " \n", " \n", " WAMBE0023\n", - " 3.3761\n", - " 0.3020\n", - " 3.0741\n", - " 91.0554\n", + " 4.4796\n", + " 0.3356\n", + " 4.1440\n", + " 92.5081\n", " NaN\n", " NaN\n", " NaN\n", @@ -1750,10 +1650,10 @@ " \n", " \n", " WAMBE0024\n", - " 60.8648\n", - " 31.2794\n", - " 29.5854\n", - " 48.6084\n", + " 61.8478\n", + " 31.3007\n", + " 30.5470\n", + " 49.3907\n", " NaN\n", " NaN\n", " NaN\n", @@ -1762,10 +1662,10 @@ " \n", " \n", " WAMBE0025\n", - " 45.1055\n", - " 14.6028\n", - " 30.5028\n", - " 67.6253\n", + " 45.9707\n", + " 14.6125\n", + " 31.3582\n", + " 68.2134\n", " NaN\n", " NaN\n", " NaN\n", @@ -1774,10 +1674,10 @@ " \n", " \n", " WAMBE0026\n", - " 32.1502\n", - " 12.9335\n", - " 19.2167\n", - " 59.7716\n", + " 32.8591\n", + " 12.9479\n", + " 19.9112\n", + " 60.5957\n", " NaN\n", " NaN\n", " NaN\n", @@ -1786,10 +1686,10 @@ " \n", " \n", " WAMBE0027\n", - " 26.2310\n", - " 18.6828\n", - " 7.5482\n", - " 28.7759\n", + " 26.4132\n", + " 18.7142\n", + " 7.6990\n", + " 29.1484\n", " NaN\n", " NaN\n", " NaN\n", @@ -1798,211 +1698,568 @@ " \n", " \n", "\n", - "

1768 rows × 9 columns

\n", + "

242 rows × 9 columns

\n", "" ], "text/plain": [ - " prestorm_swash_vol poststorm_swash_vol swash_vol_change \\\n", - "site_id \n", - "AVOCAn0001 113.9091 77.9830 35.6103 \n", - "AVOCAn0002 106.8959 67.0913 39.6373 \n", - "AVOCAn0003 99.0484 53.6563 45.1621 \n", - "AVOCAn0004 74.7543 29.3280 45.4262 \n", - "AVOCAn0005 70.5968 24.1071 46.4897 \n", - "AVOCAn0006 68.7583 23.9665 44.7918 \n", - "AVOCAn0007 75.7895 27.2715 48.5180 \n", - "AVOCAn0008 93.3107 42.4968 50.8139 \n", - "AVOCAn0009 3.6955 0.1038 3.5917 \n", - "AVOCAs0001 NaN NaN NaN \n", - "AVOCAs0002 97.3864 26.6619 70.7246 \n", - "AVOCAs0003 70.7401 40.0608 30.7919 \n", - "AVOCAs0004 97.8389 45.4845 52.2157 \n", - "AVOCAs0005 95.1711 54.9722 40.1706 \n", - "AVOCAs0006 112.5818 67.8718 44.8252 \n", - "AVOCAs0007 65.3531 44.1537 21.5228 \n", - "AVOCAs0008 52.3940 44.9152 7.4803 \n", - "BILG0001 19.4177 7.5746 11.8431 \n", - "BILG0002 155.6998 98.1693 57.4340 \n", - "BILG0003 83.5219 52.8059 30.5537 \n", - "BILG0004 158.6283 116.6189 42.1178 \n", - "BILG0005 212.8478 164.0044 48.4312 \n", - "BLUEYS0001 65.4628 19.2938 46.1690 \n", - "BLUEYS0002 50.2084 10.3009 39.9074 \n", - "BLUEYS0003 50.6308 11.1682 39.4625 \n", - "BLUEYS0004 95.1608 31.3330 63.8279 \n", - "BLUEYS0005 141.0643 58.2545 82.8098 \n", - "BLUEYS0006 88.4207 51.6205 36.2553 \n", - "BOAT0001 23.8510 23.5660 -0.0264 \n", - "BOAT0002 37.6524 14.0209 23.6316 \n", - "... ... ... ... \n", - "TREACH0014 97.5323 46.2994 51.0816 \n", - "TREACH0015 96.8327 45.1962 51.6364 \n", - "TREACH0016 106.9083 66.0567 40.3629 \n", - "WAMBE0001 132.3413 143.4459 -9.7255 \n", - "WAMBE0002 151.1833 126.9844 23.9548 \n", - "WAMBE0003 154.1788 117.9441 36.2425 \n", - "WAMBE0004 137.8449 76.6007 61.2725 \n", - "WAMBE0005 NaN NaN NaN \n", - "WAMBE0006 205.8453 186.0784 22.5892 \n", - "WAMBE0007 80.4674 35.4614 45.0059 \n", - "WAMBE0008 88.4574 40.3200 48.1374 \n", - "WAMBE0009 70.9159 26.1742 44.7418 \n", - "WAMBE0010 58.6604 18.0418 40.6186 \n", - "WAMBE0011 59.2415 16.3165 42.9250 \n", - "WAMBE0012 74.4189 23.0232 51.3957 \n", - "WAMBE0013 70.4964 22.7546 47.7419 \n", - "WAMBE0014 68.0896 24.1853 43.9043 \n", - "WAMBE0015 55.0780 16.0119 39.0660 \n", - "WAMBE0016 96.7687 39.8224 56.9463 \n", - "WAMBE0017 35.2987 8.5140 26.7847 \n", - "WAMBE0018 40.9407 10.5147 30.4260 \n", - "WAMBE0019 38.2838 9.2156 29.0682 \n", - "WAMBE0020 NaN NaN NaN \n", - "WAMBE0021 NaN NaN NaN \n", - "WAMBE0022 0.5516 0.2840 0.2675 \n", - "WAMBE0023 3.3761 0.3020 3.0741 \n", - "WAMBE0024 60.8648 31.2794 29.5854 \n", - "WAMBE0025 45.1055 14.6028 30.5028 \n", - "WAMBE0026 32.1502 12.9335 19.2167 \n", - "WAMBE0027 26.2310 18.6828 7.5482 \n", + " prestorm_swash_vol poststorm_swash_vol swash_vol_change \\\n", + "site_id \n", + "AVOCAn0009 4.5783 0.1110 4.4673 \n", + "AVOCAs0001 NaN NaN NaN \n", + "AVOCAs0002 97.9463 26.6638 71.2825 \n", + "AVOCAs0003 70.7306 40.2020 30.7232 \n", + "AVOCAs0004 98.2859 45.4986 52.6330 \n", + "AVOCAs0005 95.5841 54.9753 40.5733 \n", + "AVOCAs0006 113.0441 67.8912 45.2582 \n", + "AVOCAs0007 65.3283 44.2821 21.4544 \n", + "AVOCAs0008 52.3933 45.2243 7.1728 \n", + "BILG0001 20.3405 7.6207 12.7198 \n", + "BILG0002 156.4205 98.1716 58.1659 \n", + "BOAT0001 23.8361 23.6865 -0.0926 \n", + "BOAT0002 38.8398 14.0819 24.7579 \n", + "BOAT0003 73.6809 17.8545 55.8264 \n", + "BOAT0004 73.1954 23.1583 50.0372 \n", + "BOAT0005 53.5122 22.4537 31.0585 \n", + "BOOM0001 236.4540 218.4918 20.8725 \n", + "CATHIE0024 63.6452 38.4261 25.2191 \n", + "CATHIE0026 75.1334 43.7179 31.0940 \n", + "CRESn0069 37.5896 8.3495 29.2401 \n", + "DEEWHYn0008 NaN NaN NaN \n", + "DEEWHYn0009 NaN NaN NaN \n", + "DEEWHYs0005 62.3514 24.9797 37.3716 \n", + "DEEWHYs0008 1.0688 1.3640 0.0000 \n", + "DIAMONDn0023 67.9416 21.1812 46.7603 \n", + "DIAMONDs0006 74.9357 42.4382 32.2536 \n", + "DIAMONDs0007 153.7639 127.9469 26.1595 \n", + "DUNBn0031 36.5301 6.3289 30.2012 \n", + "DUNBn0055 189.5283 134.2760 56.7139 \n", + "ELIZA0002 NaN NaN NaN \n", + "... ... ... ... \n", + "STOCNs0170 198.2785 216.6368 -18.6067 \n", + "STOCNs0175 136.2126 114.1715 22.4984 \n", + "STOCNs0179 67.7795 45.3981 22.3815 \n", + "STOCNs0180 166.0813 149.5195 15.3441 \n", + "STOCNs0181 90.1147 98.9808 -9.1107 \n", + "STOCNs0182 67.8622 86.0118 -18.1671 \n", + "STOCNs0183 125.9085 137.7342 -12.6233 \n", + "STOCNs0184 146.6586 123.5371 23.2603 \n", + "STOCNs0185 141.5421 142.7279 -1.2619 \n", + "STOCNs0186 115.9148 123.1507 -7.8392 \n", + "STOCNs0187 126.5519 147.4371 -22.4452 \n", + "STOCNs0188 345.5234 353.8766 -11.3322 \n", + "STOCNs0189 171.6354 134.9192 35.8697 \n", + "STOCNs0190 151.4113 116.6381 35.0161 \n", + "STOCS0014 98.4991 57.3495 40.9006 \n", + "STOCS0043 36.4256 11.7208 24.7048 \n", + "WAMBE0005 NaN NaN NaN \n", + "WAMBE0015 56.2724 16.0428 40.2296 \n", + "WAMBE0016 97.8849 39.8432 58.0417 \n", + "WAMBE0017 36.5683 8.5380 28.0303 \n", + "WAMBE0018 42.1423 10.5498 31.5925 \n", + "WAMBE0019 39.6097 9.2404 30.3693 \n", + "WAMBE0020 NaN NaN NaN \n", + "WAMBE0021 NaN NaN NaN \n", + "WAMBE0022 1.1034 0.4478 0.6556 \n", + "WAMBE0023 4.4796 0.3356 4.1440 \n", + "WAMBE0024 61.8478 31.3007 30.5470 \n", + "WAMBE0025 45.9707 14.6125 31.3582 \n", + "WAMBE0026 32.8591 12.9479 19.9112 \n", + "WAMBE0027 26.4132 18.7142 7.6990 \n", "\n", - " swash_pct_change prestorm_dune_face_vol poststorm_dune_face_vol \\\n", - "site_id \n", - "AVOCAn0001 31.2620 165.4760 166.7296 \n", - "AVOCAn0002 37.0803 256.4137 258.9174 \n", - "AVOCAn0003 45.5960 372.7031 373.9198 \n", - "AVOCAn0004 60.7674 275.1689 276.0476 \n", - "AVOCAn0005 65.8524 268.5194 263.4262 \n", - "AVOCAn0006 65.1438 202.6770 198.2397 \n", - "AVOCAn0007 64.0168 149.8479 143.1312 \n", - "AVOCAn0008 54.4567 187.9201 187.3459 \n", - "AVOCAn0009 97.1908 NaN NaN \n", - "AVOCAs0001 NaN NaN NaN \n", - "AVOCAs0002 72.6226 NaN NaN \n", - "AVOCAs0003 43.5282 NaN NaN \n", - "AVOCAs0004 53.3691 NaN NaN \n", - "AVOCAs0005 42.2088 NaN NaN \n", - "AVOCAs0006 39.8157 NaN NaN \n", - "AVOCAs0007 32.9331 NaN NaN \n", - "AVOCAs0008 14.2770 NaN NaN \n", - "BILG0001 60.9913 NaN NaN \n", - "BILG0002 36.8876 NaN NaN \n", - "BILG0003 36.5817 41.1469 40.6081 \n", - "BILG0004 26.5512 11.2211 11.0892 \n", - "BILG0005 22.7539 613.6156 606.2766 \n", - "BLUEYS0001 70.5270 130.7447 120.5446 \n", - "BLUEYS0002 79.4836 512.0154 477.1774 \n", - "BLUEYS0003 77.9418 443.0853 414.3901 \n", - "BLUEYS0004 67.0737 287.5805 272.4267 \n", - "BLUEYS0005 58.7036 539.3864 520.0732 \n", - "BLUEYS0006 41.0032 271.6036 267.1954 \n", - "BOAT0001 -0.1108 NaN NaN \n", - "BOAT0002 62.7624 NaN NaN \n", - "... ... ... ... \n", - "TREACH0014 52.3740 508.7400 505.7877 \n", - "TREACH0015 53.3254 690.8275 683.4458 \n", - "TREACH0016 37.7547 508.0014 499.6315 \n", - "WAMBE0001 -7.3488 665.9898 667.5923 \n", - "WAMBE0002 15.8449 385.8467 386.7284 \n", - "WAMBE0003 23.5068 694.2226 700.5105 \n", - "WAMBE0004 44.4503 559.5485 569.8591 \n", - "WAMBE0005 NaN NaN NaN \n", - "WAMBE0006 10.9739 55.0898 55.8919 \n", - "WAMBE0007 55.9307 178.1005 178.5439 \n", - "WAMBE0008 54.4187 258.7513 258.3849 \n", - "WAMBE0009 63.0913 267.3725 258.3720 \n", - "WAMBE0010 69.2437 187.5259 161.9748 \n", - "WAMBE0011 72.4577 197.0129 175.2512 \n", - "WAMBE0012 69.0627 178.4783 168.3475 \n", - "WAMBE0013 67.7224 231.1513 195.2581 \n", - "WAMBE0014 64.4802 82.4268 61.2601 \n", - "WAMBE0015 70.9286 NaN NaN \n", - "WAMBE0016 58.8479 NaN NaN \n", - "WAMBE0017 75.8801 NaN NaN \n", - "WAMBE0018 74.3173 NaN NaN \n", - "WAMBE0019 75.9282 NaN NaN \n", - "WAMBE0020 NaN NaN NaN \n", - "WAMBE0021 NaN NaN NaN \n", - "WAMBE0022 48.5063 NaN NaN \n", - "WAMBE0023 91.0554 NaN NaN \n", - "WAMBE0024 48.6084 NaN NaN \n", - "WAMBE0025 67.6253 NaN NaN \n", - "WAMBE0026 59.7716 NaN NaN \n", - "WAMBE0027 28.7759 NaN NaN \n", + " swash_pct_change prestorm_dune_face_vol \\\n", + "site_id \n", + "AVOCAn0009 97.5750 NaN \n", + "AVOCAs0001 NaN NaN \n", + "AVOCAs0002 72.7771 NaN \n", + "AVOCAs0003 43.4369 NaN \n", + "AVOCAs0004 53.5509 NaN \n", + "AVOCAs0005 42.4478 NaN \n", + "AVOCAs0006 40.0359 NaN \n", + "AVOCAs0007 32.8409 NaN \n", + "AVOCAs0008 13.6904 NaN \n", + "BILG0001 62.5344 NaN \n", + "BILG0002 37.1856 NaN \n", + "BOAT0001 -0.3885 NaN \n", + "BOAT0002 63.7436 NaN \n", + "BOAT0003 75.7678 NaN \n", + "BOAT0004 68.3610 NaN \n", + "BOAT0005 58.0400 NaN \n", + "BOOM0001 8.8273 NaN \n", + "CATHIE0024 39.6245 0.0 \n", + "CATHIE0026 41.3851 NaN \n", + "CRESn0069 77.7877 NaN \n", + "DEEWHYn0008 NaN NaN \n", + "DEEWHYn0009 NaN NaN \n", + "DEEWHYs0005 59.9372 NaN \n", + "DEEWHYs0008 0.0000 NaN \n", + "DIAMONDn0023 68.8244 NaN \n", + "DIAMONDs0006 43.0416 NaN \n", + "DIAMONDs0007 17.0128 NaN \n", + "DUNBn0031 82.6748 NaN \n", + "DUNBn0055 29.9237 NaN \n", + "ELIZA0002 NaN NaN \n", + "... ... ... \n", + "STOCNs0170 -9.3841 NaN \n", + "STOCNs0175 16.5171 NaN \n", + "STOCNs0179 33.0210 NaN \n", + "STOCNs0180 9.2389 NaN \n", + "STOCNs0181 -10.1102 NaN \n", + "STOCNs0182 -26.7705 NaN \n", + "STOCNs0183 -10.0257 NaN \n", + "STOCNs0184 15.8602 NaN \n", + "STOCNs0185 -0.8915 NaN \n", + "STOCNs0186 -6.7629 NaN \n", + "STOCNs0187 -17.7359 NaN \n", + "STOCNs0188 -3.2797 NaN \n", + "STOCNs0189 20.8988 NaN \n", + "STOCNs0190 23.1264 NaN \n", + "STOCS0014 41.5238 NaN \n", + "STOCS0043 67.8225 NaN \n", + "WAMBE0005 NaN NaN \n", + "WAMBE0015 71.4908 NaN \n", + "WAMBE0016 59.2958 NaN \n", + "WAMBE0017 76.6520 NaN \n", + "WAMBE0018 74.9662 NaN \n", + "WAMBE0019 76.6714 NaN \n", + "WAMBE0020 NaN NaN \n", + "WAMBE0021 NaN NaN \n", + "WAMBE0022 59.4166 NaN \n", + "WAMBE0023 92.5081 NaN \n", + "WAMBE0024 49.3907 NaN \n", + "WAMBE0025 68.2134 NaN \n", + "WAMBE0026 60.5957 NaN \n", + "WAMBE0027 29.1484 NaN \n", "\n", - " dune_face_vol_change dune_face_pct_change storm_regime \n", - "site_id \n", - "AVOCAn0001 0.0000 0.0000 swash \n", - "AVOCAn0002 0.0000 0.0000 swash \n", - "AVOCAn0003 -0.3147 -0.0844 swash \n", - "AVOCAn0004 0.4104 0.1492 swash \n", - "AVOCAn0005 7.4196 2.7631 collision \n", - "AVOCAn0006 4.7944 2.3655 collision \n", - "AVOCAn0007 7.1323 4.7597 collision \n", - "AVOCAn0008 2.8297 1.5058 collision \n", - "AVOCAn0009 NaN NaN NaN \n", - "AVOCAs0001 NaN NaN NaN \n", - "AVOCAs0002 NaN NaN NaN \n", - "AVOCAs0003 NaN NaN NaN \n", - "AVOCAs0004 NaN NaN NaN \n", - "AVOCAs0005 NaN NaN NaN \n", - "AVOCAs0006 NaN NaN NaN \n", - "AVOCAs0007 NaN NaN NaN \n", - "AVOCAs0008 NaN NaN NaN \n", - "BILG0001 NaN NaN NaN \n", - "BILG0002 NaN NaN NaN \n", - "BILG0003 0.0000 0.0000 swash \n", - "BILG0004 -0.0132 -0.1179 swash \n", - "BILG0005 5.7738 0.9410 collision \n", - "BLUEYS0001 9.5601 7.3121 collision \n", - "BLUEYS0002 33.2825 6.5003 collision \n", - "BLUEYS0003 24.8870 5.6167 collision \n", - "BLUEYS0004 12.9641 4.5080 collision \n", - "BLUEYS0005 12.0470 2.2335 collision \n", - "BLUEYS0006 3.6045 1.3271 collision \n", - "BOAT0001 NaN NaN NaN \n", - "BOAT0002 NaN NaN NaN \n", - "... ... ... ... \n", - "TREACH0014 0.4254 0.0836 swash \n", - "TREACH0015 1.5086 0.2184 NaN \n", - "TREACH0016 0.3386 0.0667 swash \n", - "WAMBE0001 0.0410 0.0062 swash \n", - "WAMBE0002 -0.0449 -0.0116 swash \n", - "WAMBE0003 -4.2136 -0.6070 swash \n", - "WAMBE0004 -4.4590 -0.7969 swash \n", - "WAMBE0005 NaN NaN NaN \n", - "WAMBE0006 0.0000 0.0000 swash \n", - "WAMBE0007 0.4727 0.2654 swash \n", - "WAMBE0008 -1.2073 -0.4666 swash \n", - "WAMBE0009 9.8041 3.6668 collision \n", - "WAMBE0010 25.3087 13.4961 collision \n", - "WAMBE0011 21.9882 11.1608 collision \n", - "WAMBE0012 10.0386 5.6246 collision \n", - "WAMBE0013 35.8072 15.4908 collision \n", - "WAMBE0014 21.1718 25.6856 collision \n", - "WAMBE0015 NaN NaN NaN \n", - "WAMBE0016 NaN NaN NaN \n", - "WAMBE0017 NaN NaN NaN \n", - "WAMBE0018 NaN NaN NaN \n", - "WAMBE0019 NaN NaN NaN \n", - "WAMBE0020 NaN NaN NaN \n", - "WAMBE0021 NaN NaN NaN \n", - "WAMBE0022 NaN NaN NaN \n", - "WAMBE0023 NaN NaN NaN \n", - "WAMBE0024 NaN NaN NaN \n", - "WAMBE0025 NaN NaN NaN \n", - "WAMBE0026 NaN NaN NaN \n", - "WAMBE0027 NaN NaN NaN \n", + " poststorm_dune_face_vol dune_face_vol_change \\\n", + "site_id \n", + "AVOCAn0009 NaN NaN \n", + "AVOCAs0001 NaN NaN \n", + "AVOCAs0002 NaN NaN \n", + "AVOCAs0003 NaN NaN \n", + "AVOCAs0004 NaN NaN \n", + "AVOCAs0005 NaN NaN \n", + "AVOCAs0006 NaN NaN \n", + "AVOCAs0007 NaN NaN \n", + "AVOCAs0008 NaN NaN \n", + "BILG0001 NaN NaN \n", + "BILG0002 NaN NaN \n", + "BOAT0001 NaN NaN \n", + "BOAT0002 NaN NaN \n", + "BOAT0003 NaN NaN \n", + "BOAT0004 NaN NaN \n", + "BOAT0005 NaN NaN \n", + "BOOM0001 NaN NaN \n", + "CATHIE0024 0.0 0.0 \n", + "CATHIE0026 NaN NaN \n", + "CRESn0069 NaN NaN \n", + "DEEWHYn0008 NaN NaN \n", + "DEEWHYn0009 NaN NaN \n", + "DEEWHYs0005 NaN NaN \n", + "DEEWHYs0008 NaN NaN \n", + "DIAMONDn0023 NaN NaN \n", + "DIAMONDs0006 NaN NaN \n", + "DIAMONDs0007 NaN NaN \n", + "DUNBn0031 NaN NaN \n", + "DUNBn0055 NaN NaN \n", + "ELIZA0002 NaN NaN \n", + "... ... ... \n", + "STOCNs0170 NaN NaN \n", + "STOCNs0175 NaN NaN \n", + "STOCNs0179 NaN NaN \n", + "STOCNs0180 NaN NaN \n", + "STOCNs0181 NaN NaN \n", + "STOCNs0182 NaN NaN \n", + "STOCNs0183 NaN NaN \n", + "STOCNs0184 NaN NaN \n", + "STOCNs0185 NaN NaN \n", + "STOCNs0186 NaN NaN \n", + "STOCNs0187 NaN NaN \n", + "STOCNs0188 NaN NaN \n", + "STOCNs0189 NaN NaN \n", + "STOCNs0190 NaN NaN \n", + "STOCS0014 NaN NaN \n", + "STOCS0043 NaN NaN \n", + "WAMBE0005 NaN NaN \n", + "WAMBE0015 NaN NaN \n", + "WAMBE0016 NaN NaN \n", + "WAMBE0017 NaN NaN \n", + "WAMBE0018 NaN NaN \n", + "WAMBE0019 NaN NaN \n", + "WAMBE0020 NaN NaN \n", + "WAMBE0021 NaN NaN \n", + "WAMBE0022 NaN NaN \n", + "WAMBE0023 NaN NaN \n", + "WAMBE0024 NaN NaN \n", + "WAMBE0025 NaN NaN \n", + "WAMBE0026 NaN NaN \n", + "WAMBE0027 NaN NaN \n", "\n", - "[1768 rows x 9 columns]" + " dune_face_pct_change storm_regime \n", + "site_id \n", + "AVOCAn0009 NaN NaN \n", + "AVOCAs0001 NaN NaN \n", + "AVOCAs0002 NaN NaN \n", + "AVOCAs0003 NaN NaN \n", + "AVOCAs0004 NaN NaN \n", + "AVOCAs0005 NaN NaN \n", + "AVOCAs0006 NaN NaN \n", + "AVOCAs0007 NaN NaN \n", + "AVOCAs0008 NaN NaN \n", + "BILG0001 NaN NaN \n", + "BILG0002 NaN NaN \n", + "BOAT0001 NaN NaN \n", + "BOAT0002 NaN NaN \n", + "BOAT0003 NaN NaN \n", + "BOAT0004 NaN NaN \n", + "BOAT0005 NaN NaN \n", + "BOOM0001 NaN NaN \n", + "CATHIE0024 NaN NaN \n", + "CATHIE0026 NaN NaN \n", + "CRESn0069 NaN NaN \n", + "DEEWHYn0008 NaN NaN \n", + "DEEWHYn0009 NaN NaN \n", + "DEEWHYs0005 NaN NaN \n", + "DEEWHYs0008 NaN NaN \n", + "DIAMONDn0023 NaN NaN \n", + "DIAMONDs0006 NaN NaN \n", + "DIAMONDs0007 NaN NaN \n", + "DUNBn0031 NaN NaN \n", + "DUNBn0055 NaN NaN \n", + "ELIZA0002 NaN NaN \n", + "... ... ... \n", + "STOCNs0170 NaN NaN \n", + "STOCNs0175 NaN NaN \n", + "STOCNs0179 NaN NaN \n", + "STOCNs0180 NaN NaN \n", + "STOCNs0181 NaN NaN \n", + "STOCNs0182 NaN NaN \n", + "STOCNs0183 NaN NaN \n", + "STOCNs0184 NaN NaN \n", + "STOCNs0185 NaN NaN \n", + "STOCNs0186 NaN NaN \n", + "STOCNs0187 NaN NaN \n", + "STOCNs0188 NaN NaN \n", + "STOCNs0189 NaN NaN \n", + "STOCNs0190 NaN NaN \n", + "STOCS0014 NaN NaN \n", + "STOCS0043 NaN NaN \n", + "WAMBE0005 NaN NaN \n", + "WAMBE0015 NaN NaN \n", + "WAMBE0016 NaN NaN \n", + "WAMBE0017 NaN NaN \n", + "WAMBE0018 NaN NaN \n", + "WAMBE0019 NaN NaN \n", + "WAMBE0020 NaN NaN \n", + "WAMBE0021 NaN NaN \n", + "WAMBE0022 NaN NaN \n", + "WAMBE0023 NaN NaN \n", + "WAMBE0024 NaN NaN \n", + "WAMBE0025 NaN NaN \n", + "WAMBE0026 NaN NaN \n", + "WAMBE0027 NaN NaN \n", + "\n", + "[242 rows x 9 columns]" ] }, - "execution_count": 30, + "execution_count": 58, "metadata": {}, "output_type": "execute_result" } ], + "source": [ + "df_impacts = impacts['observed']\n", + "df_no_obs_impacts = df_impacts[df_impacts.storm_regime.isnull()]\n", + "no_obs_impacts_sites = df_no_obs_impacts.index\n", + "\n", + "df_no_obs_impacts\n" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "It looks like the problem comes from if we cannot identify the prestorm and post storm swash and berm volume changes." + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### Check pre and post storm profiles\n", + "It looks like, for some reason, there are no post storm profiles (`STUART0011`) for some of the profiles. Not sure if this is a processing error, or it hasn't been included in the dataset." + ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], "source": [] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "### Check prestorm dune crests\n", + "If there are no prestorm dune crests defined, we cannot define a mean slope and hence we won't be able to get observed impacts" + ] + }, + { + "cell_type": "code", + "execution_count": 74, + "metadata": { + "ExecuteTime": { + "end_time": "2018-12-13T02:37:07.225965Z", + "start_time": "2018-12-13T02:37:07.213921Z" + } + }, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "231 sites have no dune crests:\n", + "AVOCAn0009\n", + "AVOCAs0001\n", + "AVOCAs0002\n", + "AVOCAs0003\n", + "AVOCAs0004\n", + "AVOCAs0005\n", + "AVOCAs0006\n", + "AVOCAs0007\n", + "AVOCAs0008\n", + "BILG0001\n", + "BILG0002\n", + "BOAT0001\n", + "BOAT0002\n", + "BOAT0003\n", + "BOAT0004\n", + "BOAT0005\n", + "BOOM0001\n", + "CATHIE0026\n", + "CRESn0069\n", + "DEEWHYn0008\n", + "DEEWHYn0009\n", + "DEEWHYs0005\n", + "DEEWHYs0008\n", + "DIAMONDs0006\n", + "DIAMONDs0007\n", + "ENTRA0005\n", + "ENTRA0006\n", + "ENTRA0077\n", + "ENTRA0078\n", + "ENTRA0079\n", + "FOST0003\n", + "GRANTSn0004\n", + "GRANTSn0005\n", + "GRANTSn0006\n", + "GRANTSn0007\n", + "GRANTSn0008\n", + "GRANTSn0009\n", + "GRANTSn0021\n", + "GRANTSs0014\n", + "HARGs0003\n", + "HARGs0004\n", + "HARGs0005\n", + "HARR0056\n", + "LHOUSE0001\n", + "LHOUSE0002\n", + "LHOUSE0003\n", + "LHOUSE0004\n", + "LHOUSE0012\n", + "LHOUSE0013\n", + "LHOUSEs0015\n", + "MACM0008\n", + "MACM0012\n", + "MACM0013\n", + "MACM0014\n", + "MACM0015\n", + "MACM0016\n", + "MANNING0001\n", + "MANNING0002\n", + "MANNING0003\n", + "MANNING0004\n", + "MANNING0005\n", + "MANNING0101\n", + "MANNING0102\n", + "MANNING0103\n", + "MANNING0104\n", + "MANNING0105\n", + "MANNING0106\n", + "MANNING0107\n", + "MANNING0108\n", + "MANNING0109\n", + "MONA0001\n", + "MONA0002\n", + "MONA0003\n", + "MONA0014\n", + "MONA0015\n", + "MONA0016\n", + "MONA0017\n", + "MONA0018\n", + "MONA0019\n", + "MONA0020\n", + "MONA0021\n", + "NAMB0027\n", + "NAMB0041\n", + "NARRA0001\n", + "NARRA0028\n", + "NARRA0035\n", + "NINEMn0050\n", + "OLDBAR0035\n", + "PEARLn0001\n", + "PEARLn0002\n", + "PEARLn0003\n", + "PEARLn0004\n", + "PEARLs0003\n", + "PEARLs0004\n", + "PEARLs0005\n", + "STOCNn0012\n", + "STOCNn0013\n", + "STOCNn0014\n", + "STOCNn0015\n", + "STOCNn0016\n", + "STOCNn0017\n", + "STOCNn0018\n", + "STOCNn0019\n", + "STOCNn0020\n", + "STOCNn0021\n", + "STOCNn0022\n", + "STOCNn0023\n", + "STOCNn0024\n", + "STOCNn0025\n", + "STOCNn0026\n", + "STOCNn0027\n", + "STOCNn0028\n", + "STOCNn0029\n", + "STOCNn0030\n", + "STOCNn0031\n", + "STOCNn0032\n", + "STOCNn0033\n", + "STOCNn0034\n", + "STOCNn0035\n", + "STOCNn0036\n", + "STOCNn0037\n", + "STOCNn0038\n", + "STOCNn0039\n", + "STOCNn0044\n", + "STOCNn0059\n", + "STOCNn0062\n", + "STOCNn0063\n", + "STOCNn0064\n", + "STOCNn0065\n", + "STOCNs0022\n", + "STOCNs0025\n", + "STOCNs0026\n", + "STOCNs0031\n", + "STOCNs0045\n", + "STOCNs0048\n", + "STOCNs0049\n", + "STOCNs0053\n", + "STOCNs0055\n", + "STOCNs0056\n", + "STOCNs0057\n", + "STOCNs0058\n", + "STOCNs0059\n", + "STOCNs0060\n", + "STOCNs0061\n", + "STOCNs0062\n", + "STOCNs0073\n", + "STOCNs0079\n", + "STOCNs0088\n", + "STOCNs0089\n", + "STOCNs0090\n", + "STOCNs0091\n", + "STOCNs0092\n", + "STOCNs0093\n", + "STOCNs0094\n", + "STOCNs0095\n", + "STOCNs0096\n", + "STOCNs0097\n", + "STOCNs0098\n", + "STOCNs0099\n", + "STOCNs0100\n", + "STOCNs0101\n", + "STOCNs0102\n", + "STOCNs0103\n", + "STOCNs0104\n", + "STOCNs0105\n", + "STOCNs0106\n", + "STOCNs0107\n", + "STOCNs0108\n", + "STOCNs0109\n", + "STOCNs0110\n", + "STOCNs0111\n", + "STOCNs0112\n", + "STOCNs0113\n", + "STOCNs0114\n", + "STOCNs0115\n", + "STOCNs0116\n", + "STOCNs0117\n", + "STOCNs0118\n", + "STOCNs0119\n", + "STOCNs0120\n", + "STOCNs0121\n", + "STOCNs0122\n", + "STOCNs0123\n", + "STOCNs0124\n", + "STOCNs0125\n", + "STOCNs0126\n", + "STOCNs0127\n", + "STOCNs0128\n", + "STOCNs0137\n", + "STOCNs0141\n", + "STOCNs0144\n", + "STOCNs0150\n", + "STOCNs0155\n", + "STOCNs0156\n", + "STOCNs0157\n", + "STOCNs0158\n", + "STOCNs0159\n", + "STOCNs0160\n", + "STOCNs0167\n", + "STOCNs0168\n", + "STOCNs0169\n", + "STOCNs0170\n", + "STOCNs0175\n", + "STOCNs0179\n", + "STOCNs0180\n", + "STOCNs0181\n", + "STOCNs0182\n", + "STOCNs0183\n", + "STOCNs0184\n", + "STOCNs0185\n", + "STOCNs0186\n", + "STOCNs0187\n", + "STOCNs0188\n", + "STOCNs0189\n", + "STOCNs0190\n", + "STOCS0014\n", + "STOCS0043\n", + "WAMBE0005\n", + "WAMBE0015\n", + "WAMBE0016\n", + "WAMBE0017\n", + "WAMBE0018\n", + "WAMBE0019\n", + "WAMBE0020\n", + "WAMBE0021\n", + "WAMBE0022\n", + "WAMBE0023\n", + "WAMBE0024\n", + "WAMBE0025\n", + "WAMBE0026\n", + "WAMBE0027\n" + ] + } + ], + "source": [ + "df_no_crests = df_profile_features.query('profile_type==\"prestorm\" & (dune_crest_x != dune_crest_x)')\n", + "print('{} sites have no dune crests:'.format(len(df_no_crests)))\n", + "\n", + "for site_id in df_no_crests.index.tolist():\n", + " print(site_id)" + ] } ], "metadata": { @@ -2037,10 +2294,10 @@ "title_sidebar": "Contents", "toc_cell": false, "toc_position": { - "height": "calc(100% - 180px)", - "left": "10px", - "top": "150px", - "width": "275.797px" + "height": "656px", + "left": "508px", + "top": "90px", + "width": "218.797px" }, "toc_section_display": true, "toc_window_display": true