Stylistic tutorial edits in 09_DataCorrection_Becker_Adusumilli

Author: mayakbecker

09_DataCorrection_Becker_Adusumilli/data_correction_rendered.ipynb

@@ -99,7 +99,7 @@
     "                        DS = '/gt%d%s/%s/%s' % (pair, beam, group, dataset)\n",
     "                        try:\n",
     "                            temp[dataset] = np.array(h5f[DS])                            \n",
-    "                            if '_FillValue' in h5f[DS].attrs: # Converting Fortran NaNs to Python NaNs\n",
+    "                            if '_FillValue' in h5f[DS].attrs: # Convert Fortran NaNs to Python NaNs\n",
     "                                fill_value = h5f[DS].attrs['_FillValue']\n",
     "                                temp[dataset][temp[dataset] == fill_value] = np.NaN\n",
     "                        except KeyError:\n",
@@ -137,7 +137,8 @@
    ],
    "source": [
     "%%bash\n",
-    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/byrd_moa.mat /home/jovyan/external_data/ "
+    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/byrd_moa.mat \\\n",
+    "    /home/jovyan/external_data/ "
    ]
   },
   {
@@ -199,16 +200,17 @@
      "output_type": "stream",
      "text": [
       "download: s3://pangeo-data-upload-oregon/icesat2/data-correction/IceShelf_Antarctica_v1.cpg to ../external_data/IceShelf_Antarctica_v1.cpg\n",
-      "download: s3://pangeo-data-upload-oregon/icesat2/data-correction/IceShelf_Antarctica_v1.shx to ../external_data/IceShelf_Antarctica_v1.shx\n",
       "download: s3://pangeo-data-upload-oregon/icesat2/data-correction/IceShelf_Antarctica_v1.dbf to ../external_data/IceShelf_Antarctica_v1.dbf\n",
+      "download: s3://pangeo-data-upload-oregon/icesat2/data-correction/IceShelf_Antarctica_v1.shx to ../external_data/IceShelf_Antarctica_v1.shx\n",
       "download: s3://pangeo-data-upload-oregon/icesat2/data-correction/IceShelf_Antarctica_v1.prj to ../external_data/IceShelf_Antarctica_v1.prj\n",
       "download: s3://pangeo-data-upload-oregon/icesat2/data-correction/IceShelf_Antarctica_v1.shp to ../external_data/IceShelf_Antarctica_v1.shp\n"
      ]
     }
    ],
    "source": [
     "%%bash\n",
-    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/ /home/jovyan/external_data/ --recursive --exclude \"*\" --include \"IceShelf_Antarctica_v1.*\""
+    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/ /home/jovyan/external_data/ \\\n",
+    "    --recursive --exclude \"*\" --include \"IceShelf_Antarctica_v1.*\""
    ]
   },
   {
@@ -278,7 +280,8 @@
     "coords = np.squeeze(np.dstack((xp,yp)))\n",
     "poly_mpl = path.Path(coords)\n",
     "\n",
-    "for Di in D6: # Add 'mask' dataset (containing the x-y coordinate pairs within the mask) to list D6 for each beam/pair combination\n",
+    "# Add 'mask' dataset (containing the x-y coordinate pairs within the mask) to list D6 for each beam/pair combination\n",
+    "for Di in D6:\n",
     "    Di['mask'] = poly_mpl.contains_points(np.transpose([Di['x'],Di['y']]))\n",
     "\n",
     "plt.figure(figsize = (10,8))\n",
@@ -346,7 +349,8 @@
     "plt.subplot(122)\n",
     "plot_moa()\n",
     "for Di in D6:\n",
-    "    plt.scatter(Di['x'][::20],Di['y'][::20],c = Di['tide_load'][::20],s = 0.1,vmin=-0.01,vmax=0.01)\n",
+    "    plt.scatter(Di['x'][::20],Di['y'][::20],c = Di['tide_load'][::20],s = 0.1,\\\n",
+    "                vmin=-0.01,vmax=0.01)\n",
     "\n",
     "cbar = plt.colorbar()\n",
     "cbar.set_label('ocean load tide height [m]', rotation=270, labelpad = 12)\n",
@@ -401,7 +405,8 @@
     "plt.figure(figsize = (10,8))\n",
     "plot_moa()\n",
     "for Di in D6:\n",
-    "    plt.scatter(Di['x'][Di['mask']][::20],Di['y'][Di['mask']][::20],c = Di['dac'][Di['mask']][::20],s = 0.1,vmin=0,vmax=0.5)\n",
+    "    plt.scatter(Di['x'][Di['mask']][::20],Di['y'][Di['mask']][::20],\\\n",
+    "                c = Di['dac'][Di['mask']][::20],s = 0.1,vmin=0,vmax=0.5)\n",
     "cbar = plt.colorbar()\n",
     "cbar.set_label('dynamic atmosphere correction [m]', rotation=270, labelpad = 12)\n",
     "\n",
@@ -474,7 +479,8 @@
    ],
    "source": [
     "%%bash\n",
-    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/ /home/jovyan/external_data/ --recursive --exclude \"*\" --include \"FDM_*.nc\""
+    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/ /home/jovyan/external_data/ \\\n",
+    "    --recursive --exclude \"*\" --include \"FDM_*.nc\""
    ]
   },
   {
@@ -496,8 +502,10 @@
     "FAC = Dataset(external_data + '/FDM_FirnAir_ANT27_1979-2016.nc')\n",
     "\n",
     "print(FAC.file_format)\n",
-    "print(FAC.dimensions.keys()) # In this case, dimensions represent the physical dimensions assigned to each variable\n",
-    "print(FAC.variables.keys()) # Variables are arrays of values where the bulk of the information is stored"
+    "print(FAC.dimensions.keys())\n",
+    "print(FAC.variables.keys())\n",
+    "\n",
+    "# In this case, dimensions represent the physical dimensions assigned to each variable. Variables are arrays of values where the bulk of the information is stored"
    ]
   },
   {
@@ -593,7 +601,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "inds = np.where((lon_fac > 147) & (lon_fac < 163) & (lat_fac > -82) & (lat_fac < -78)) # Subsetting area around Byrd\n",
+    "# Subset area around Byrd\n",
+    "inds = np.where((lon_fac > 147) & (lon_fac < 163) & (lat_fac > -82) & (lat_fac < -78))\n",
     "mean_fac_byrd = np.squeeze(meanFirnAir_16[inds]).ravel()\n",
     "xy = np.empty((len(x_fac[inds]),2))\n",
     "xy[:,0] = np.squeeze(x_fac[inds]).ravel()\n",
@@ -679,7 +688,8 @@
     "for Di in D6:\n",
     "    h_MSL = Di['h_li'] - MDT_corr - Di['geoid_h']\n",
     "    thickness = ((1/rho_ice - 1/rho_ocean)**(-1)) * ((h_MSL - Di['FAC_interp'])/rho_ice)\n",
-    "    plt.scatter(Di['x'][Di['mask']][::20],Di['y'][Di['mask']][::20],c = thickness[Di['mask']][::20],s = 0.1)\n",
+    "    plt.scatter(Di['x'][Di['mask']][::20],Di['y'][Di['mask']][::20],\\\n",
+    "                c = thickness[Di['mask']][::20],s = 0.1)\n",
     "    plt.clim([0,1000])\n",
     "cbar = plt.colorbar()\n",
     "cbar.set_label('thickness [m]', rotation=270, labelpad = 12)\n",
@@ -723,7 +733,8 @@
    ],
    "source": [
     "%%bash\n",
-    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/ /home/jovyan/external_data/ --recursive --exclude \"*\" --include \"SMB_*.nc\""
+    "aws s3 cp s3://pangeo-data-upload-oregon/icesat2/data-correction/ /home/jovyan/external_data/ \\\n",
+    "    --recursive --exclude \"*\" --include \"SMB_*.nc\""
    ]
   },
   {