fixing channels.ipynb

vcantarella · Sep 17, 2024 · 5eb09a8 · 5eb09a8
1 parent c47c52d
commit 5eb09a8
Showing 1 changed file with 19 additions and 19 deletions.
diff --git a/examples/channels/test_channels.ipynb b/examples/channels/test_channels.ipynb
@@ -549,26 +549,26 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "    tau = np.interp(500, s, curvature)\n",
-    "    Ro = k_1 *tau * W\n",
-    "    s_prime = np.where(s < 500, s, 0)\n",
-    "    s_prime = s_prime[s_prime != 0]\n",
-    "    curv_prime = curvature[:len(s_prime)]\n",
-    "    s_prime = np.concatenate((s_prime, [500]))\n",
-    "    curv_prime = np.concatenate((curv_prime, [tau]))\n",
-    "    sau = 500 - s_prime\n",
-    "    G_sau = np.exp(-2* Cf * sau/D)\n",
-    "    Ro_prime = k_1 * curv_prime * W\n",
-    "    index = np.where(G_sau>1e-13)\n",
-    "    index = index[0]\n",
+    "tau = np.interp(500, s, curvature)\n",
+    "Ro = k_1 *tau * W\n",
+    "s_prime = np.where(s < 500, s, 0)\n",
+    "s_prime = s_prime[s_prime != 0]\n",
+    "curv_prime = curvature[:len(s_prime)]\n",
+    "s_prime = np.concatenate((s_prime, [500]))\n",
+    "curv_prime = np.concatenate((curv_prime, [tau]))\n",
+    "sau = 500 - s_prime\n",
+    "G_sau = np.exp(-2* Cf * sau/D)\n",
+    "Ro_prime = k_1 * curv_prime * W\n",
+    "index = np.where(G_sau>1e-13)\n",
+    "index = index[0]\n",
     "\n",
-    "    G_sau = G_sau[index]\n",
-    "    Ro_prime = Ro_prime[index]\n",
-    "    sau = sau[index]\n",
-    "    G_sau = np.flip(G_sau)\n",
-    "    Ro_prime = np.flip(Ro_prime)\n",
-    "    sau = np.flip(sau)\n",
-    "    integration = np.trapz(Ro_prime * G_sau, sau)/(1/(2*Cf/D))"
+    "G_sau = G_sau[index]\n",
+    "Ro_prime = Ro_prime[index]\n",
+    "sau = sau[index]\n",
+    "G_sau = np.flip(G_sau)\n",
+    "Ro_prime = np.flip(Ro_prime)\n",
+    "sau = np.flip(sau)\n",
+    "integration = np.trapz(Ro_prime * G_sau, sau)/(1/(2*Cf/D))"
    ]
   },
   {