Yawn... more typos

Author: ifenty

Tutorials_as_Jupyter_Notebooks/ECCO_v4_Accessing_and_Subsetting_Variables.ipynb

@@ -724,7 +724,7 @@
     "\n",
     "To demonstrate, we'll pull out the ``numpy`` array of *SSH* and then select the first time record and second tile\n",
     "\n",
-    "  > **Note:** *Note `numpy` array indexing starts with 0*"
+    "  > **Note:**  `numpy` array indexing starts with 0*"
    ]
   },
   {

Tutorials_as_Jupyter_Notebooks/ECCO_v4_Volume_budget_closure.ipynb

@@ -19,7 +19,9 @@
     "\n",
     "ECCOv4 uses the $z^*$ coordinate system in which the depth of the vertical coordinate, $z^*$ varies with time as:\n",
     "\n",
-    "$$ z^* = \\frac{z - \\eta(x,y,t)}{H(x,y) + \\eta(x,y,t)} H(x,y)$$\n",
+    "\\begin{equation}\n",
+    "z^* = \\frac{z - \\eta(x,y,t)}{H(x,y) + \\eta(x,y,t)} H(x,y)\n",
+    "\\end{equation}\n",
     "\n",
     "With $H$ being the model depth, $\\eta$ being the model sea level anomaly, and $z$ being depth.\n",
     "\n",
@@ -48,15 +50,15 @@
     "\n",
     "To make budget calculations easier we provide the scaled velocities quantities ``UVELMASS`` and ``VVELMASS``,\n",
     "\n",
-    "\\begin{equation}\n",
+    "\\begin{align}\n",
     "\\mathit{UVELMASS}(x,y,k) = \\mathit{UVEL}(x,y,k) \\times \\mathit{hFacW}(x,y,k) \\times s^{*}(x,y,k,t)\n",
-    "\\end{equation}\n",
+    "\\end{align}\n",
     "and \n",
-    "\\begin{equation}\n",
+    "\\begin{align}\n",
     "\\mathit{VVELMASS}(x,y,k) = \\mathit{VVEL}(x,y,k) \\times \\mathit{hFacS}(x,y,k) \\times s^{*}(x,y,k,t)\n",
-    "\\end{equation}\n",
+    "\\end{align}\n",
     "\n",
-    "> **Note:** The word **mass** in ``UVELMASS`` and ``VVELMASS`` is confusing since there is no mass involved here.  Think of these terms as simply being ``UVEL`` and ``VVEL`` multiplied by the scaled grid cell area fraction.\n",
+    "It is worth noting that the word **mass** in ``UVELMASS`` and ``VVELMASS`` is confusing since there is no mass involved here.  Think of these terms as simply being ``UVEL`` and ``VVEL`` multiplied by the fraction of the grid cell height that is open grid cell face across which the volume transport occurs.  Partial cell bathymetry can make this fraction (hFacW, hFacS) less than one, and the $s^*$ scaling factor further adjusts this fraction higher or lower through time.\n",
     "\n",
     "Fully closing the budget requires the vertical volume fluxes across the top and bottom 'w' faces of the grid cell and surface freshwater fluxes.  Regarding vertical volume fluxes, there are no $s^*$ or ``hFac`` equivalent scaling factors that modify our top and bottom grid cell areas.  Therefore, vertical volume fluxes through 'w' faces are simply:\n",
     "\n",
@@ -136,7 +138,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": 62,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -145,7 +147,8 @@
     "\n",
     "## needed to convert surface mass fluxes to volume fluxes\n",
     "\n",
-    "# lat/lon resolution in degrees to interpolate the model fields for the purposes of plotting\n",
+    "# lat/lon resolution in degrees to interpolate the model \n",
+    "# fields for the purposes of plotting\n",
     "map_dx = .2\n",
     "map_dy = .2"
    ]
@@ -235,8 +238,8 @@
     "\n",
     "# load one extra year worth of snapshots\n",
     "ecco_monthly_snaps = ecco.recursive_load_ecco_var_from_years_nc(data_dir, \\\n",
-    "                                           vars_to_load=['ETAN'],\\\n",
-    "                                           years_to_load=range(year_start, year_end+1)).load()\n",
+    "                         vars_to_load=['ETAN'],\\\n",
+    "                         years_to_load=range(year_start, year_end+1)).load()\n",
     "\n",
     "num_months = len(ecco_monthly_snaps.time.values)\n",
     "# drop the last 11 months so that we have one snapshot at the beginning and end of each month within the \n",
@@ -278,7 +281,8 @@
    ],
    "source": [
     "# find the record of the last ETAN snapshot\n",
-    "last_record_date = ecco.extract_yyyy_mm_dd_hh_mm_ss_from_datetime64(ecco_monthly_snaps.time[-1].values)\n",
+    "last_record_date = \n",
+    "    ecco.extract_yyyy_mm_dd_hh_mm_ss_from_datetime64(ecco_monthly_snaps.time[-1].values)\n",
     "print(last_record_date)\n",
     "last_record_year = last_record_date[0]"
    ]
@@ -768,6 +772,7 @@
     "                              cmin=-1, cmax=1, \\\n",
     "                              cmap=cmocean.cm.balance, user_lon_0=-67,\\\n",
     "                              dx=map_dx,dy=map_dy);\n",
+    "\n",
     "plt.title('Actual $\\Delta \\eta$ [m]', fontsize=20);"
    ]
   },
@@ -857,6 +862,7 @@
     "                              cmin=-1e-7, cmax=1e-7,\\\n",
     "                              cmap=cmocean.cm.balance, user_lon_0=-67,\\\n",
     "                              dx=map_dx,dy=map_dy);\n",
+    "\n",
     "plt.title('$\\partial \\eta / \\partial t$ [m/s] during ' + \n",
     "          str(tmp[0]) +'/' + str(tmp[1]), fontsize=20);"
    ]
@@ -1729,7 +1735,8 @@
     "plt.grid()\n",
     "\n",
     "plt.subplots_adjust(hspace = .5, wspace=.2)\n",
-    "plt.suptitle('Volume Budget for one Arctic Ocean point',fontsize=20);"
+    "plt.suptitle('Volume Budget for one Arctic Ocean point',\n",
+    "                 fontsize=20);"
    ]
   },
   {
@@ -1803,12 +1810,15 @@
     "area_masked = ecco_grid.rA.where(ecco_grid.maskC.isel(k=0) == 1)\n",
     "\n",
     "dETA_per_month_predicted_from_surf_fluxes = \\\n",
-    "    ((G_surf_fluxes * area_masked).sum(dim=('i','j','tile'))/area_masked.sum())*secs_per_month\n",
+    "    ((G_surf_fluxes * area_masked).sum(dim=('i','j','tile')) / \n",
+    "     area_masked.sum())*secs_per_month\n",
     "\n",
     "ETA_predicted_by_surf_fluxes = \\\n",
     "    np.cumsum(dETA_per_month_predicted_from_surf_fluxes.values)\n",
     "\n",
-    "ETA_from_ETAN = (ecco_monthly_snaps.ETAN * area_masked).sum(dim=('i','j','tile'))/ area_masked.sum()\n",
+    "ETA_from_ETAN = \n",
+    "    (ecco_monthly_snaps.ETAN * area_masked).sum(dim=('i','j','tile')) /\n",
+    "    area_masked.sum()\n",
     "\n",
     "# plotting\n",
     "plt.figure(figsize=(14,5));\n",
@@ -1819,14 +1829,15 @@
     "plt.grid()\n",
     "plt.ylabel('global mean $\\eta$');\n",
     "plt.legend(('predicted', 'actual'));\n",
-    "plt.title('$\\eta(t)$ as predicted from net surface fluxes and model ETAN [m]', fontsize=20);"
+    "plt.title('$\\eta(t)$ as predicted from net surface fluxes and model ETAN [m]', \n",
+    "          fontsize=20);"
    ]
   },
   {
    "cell_type": "markdown",
    "metadata": {},
    "source": [
-    "> **Note:** the first predicted $\\eta$ occurs at the end of the first month (one month time integral of $\\partial \\eta / \\partial t$.  The first *actual* $\\eta$ is set to be zero.\n",
+    "The first predicted $\\eta$ occurs at the end of the first month (one month time integral of $\\partial \\eta / \\partial t$.  The first *actual* $\\eta$ is set to be zero.\n",
     "\n",
     "The above plot is another way of confirming that in our Boussinesq model the only term that can change global mean model sea level anomaly $\\eta$ is net surface freshwater flux.  To account for changes in global mean density we must apply the Greatbatch correction, inverse-barometer correction, and a correction term to account for the fact that sea-ice does not 'float' on top of the ocean but in fact displaces seawater upwards.  All of these corrections are made for the term ``SSH``, dynamic sea surface height anomaly (not shown here)."
    ]
@@ -1881,7 +1892,8 @@
     "        1000*(annual_mean_GMSL_due_to_mass_fluxes.values[12:num_years]-.017), \\\n",
     "        color='k')\n",
     "plt.grid()\n",
-    "plt.xticks(np.arange(2005, annual_mean_GMSL_due_to_mass_fluxes.year.values[-1]+1,step=1));\n",
+    "plt.xticks(np.arange(2005, \n",
+    "                     annual_mean_GMSL_due_to_mass_fluxes.year.values[-1]+1,step=1));\n",
     "plt.title('Sea level (mm) caused by mass fluxes');"
    ]
   },