add map

Author: ocefpaf

jupyterbook/content/code_gallery/data_access_notebooks/2022-11-25-plan_surf_trip_EDS.ipynb

@@ -209,7 +209,7 @@
     "    print(f\"Getting data from {urls[key]}\")\n",
     "\n",
     "    ds = xr.open_dataset(urls[key])\n",
-    "    \n",
+    "\n",
     "    data = (\n",
     "        ds.filter_by_attrs(abbreviation=\"HTSGW\")\n",
     "        .sel(lat=lat, lon=lon360, method=\"nearest\")\n",
@@ -592,20 +592,19 @@
     "import matplotlib.dates as dates\n",
     "import pytz\n",
     "\n",
-    "\n",
     "kw = {\n",
     "    \"marker\": \"o\",\n",
     "    \"linestyle\": \"none\",\n",
     "}\n",
-    "        \n",
+    "\n",
+    "\n",
     "def station_scatter(df, station, ax):\n",
     "    df.plot(\n",
     "        ax=ax[0],\n",
     "        x=\"time (UTC)\",\n",
     "        y=\"sea_surface_wave_significant_height (m)\",\n",
     "        label=station,\n",
     "        **kw,\n",
-    "\n",
     "    )\n",
     "\n",
     "    df.plot(\n",
@@ -658,15 +657,25 @@
    ],
    "source": [
     "import datetime\n",
-    "import pytz\n",
     "\n",
+    "import pytz\n",
     "\n",
     "dfs[\"gdfs\"] = df_out[\n",
-    "    [\"time\", \"Significant_height_of_combined_wind_waves_and_swell_surface_gfs\", \"lat_gfs\", \"lon_gfs\"]\n",
+    "    [\n",
+    "        \"time\",\n",
+    "        \"Significant_height_of_combined_wind_waves_and_swell_surface_gfs\",\n",
+    "        \"lat_gfs\",\n",
+    "        \"lon_gfs\",\n",
+    "    ]\n",
     "].copy()\n",
     "\n",
     "dfs[\"nwps\"] = df_out[\n",
-    "    [\"time\", \"Significant_height_of_combined_wind_waves_and_swell_surface_nwps\", 'lat_nwps', 'lon_nwps']\n",
+    "    [\n",
+    "        \"time\",\n",
+    "        \"Significant_height_of_combined_wind_waves_and_swell_surface_nwps\",\n",
+    "        \"lat_nwps\",\n",
+    "        \"lon_nwps\",\n",
+    "    ]\n",
     "].copy()\n",
     "\n",
     "dfs[\"gdfs\"].rename(\n",
@@ -733,6 +742,7 @@
     "tools = \"pan,box_zoom,reset\"\n",
     "width, height = 750, 250\n",
     "\n",
+    "\n",
     "def make_plot(series, station):\n",
     "    p = figure(\n",
     "        toolbar_location=\"above\",\n",
@@ -751,6 +761,7 @@
     "    )\n",
     "    return p, line\n",
     "\n",
+    "\n",
     "def make_marker(p, location, fname):\n",
     "    html = file_html(p, CDN, fname)\n",
     "    iframe = IFrame(html, width=width + 45, height=height + 80)\n",
@@ -1077,7 +1088,6 @@
     "import folium\n",
     "import folium.plugins\n",
     "\n",
-    "\n",
     "m = folium.Map()\n",
     "folium.plugins.Fullscreen().add_to(m)\n",
     "\n",
@@ -1086,9 +1096,9 @@
     "    lat = df[\"latitude (degrees_north)\"].dropna().unique().squeeze()\n",
     "    cols = [\"time (UTC)\", \"sea_surface_wave_significant_height (m)\"]\n",
     "    df = df[cols].set_index(\"time (UTC)\").squeeze()\n",
-    "    \n",
+    "\n",
     "    df.index = df.index.tz_localize(None)\n",
-    "    \n",
+    "\n",
     "    p, _ = make_plot(df, station=station)\n",
     "    marker = make_marker(p, location=(lat, lon), fname=station)\n",
     "    marker.add_to(m)\n",