From 24868c99c17cd67082f7b04679aa89110bebda86 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Sat, 20 Feb 2021 10:15:29 +0100
Subject: [PATCH] Add supplementary plot for D~P and D~P^2.
---
Plots_Droplet_FRAP.ipynb | 75 +++++++++++++++++++++++++++++++++++++++-
1 file changed, 74 insertions(+), 1 deletion(-)
diff --git a/Plots_Droplet_FRAP.ipynb b/Plots_Droplet_FRAP.ipynb
index a0ea2eb..5f1b226 100644
--- a/Plots_Droplet_FRAP.ipynb
+++ b/Plots_Droplet_FRAP.ipynb
@@ -409,7 +409,8 @@
"D = lars.D[lars.conc==75].mean()\n",
"eta = 10.4\n",
"kBT = 4.11*10**-21\n",
- "r = kBT/(D*6*np.pi*eta)"
+ "r = kBT/(D*6*np.pi*eta)\n",
+ "print(r)"
]
},
{
@@ -1179,6 +1180,78 @@
"clb.ax.set_title('Cost')"
]
},
+ {
+ "cell_type": "markdown",
+ "metadata": {},
+ "source": [
+ "### Figure: supplement. Scaling at early and late times"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "scal1 = np.loadtxt('scaling_test_1.csv', delimiter=',')\n",
+ "scal2 = np.loadtxt('scaling_test_2.csv', delimiter=',')\n",
+ "scal3 = np.loadtxt('scaling_test_3.csv', delimiter=',')\n",
+ "s1 = plt.plot(scal1[:, 0], scal1[:, 14::2], '-', lw=1, c=dark_grey)\n",
+ "s2 = plt.plot(scal2[:, 0], scal2[:, 14::2], '-', lw=1, c=red, alpha=0.3)\n",
+ "s3 = plt.plot(scal3[:, 0], scal3[:, 14::2], '-', lw=1, c=green)\n",
+ "nice_fig('radial distance $r$ [$\\mathrm{\\mu m}$]', 'intensity [a.u]',\n",
+ " [0, 2], [0,1], [3,2])\n",
+ "plt.legend([s1[0], s2[0], s3[0]], ['Reference', '$D_\\mathrm{out}\\propto P^2$',\n",
+ " '$D_\\mathrm{out}\\propto P$'], frameon=False,\n",
+ " fontsize=9, handletextpad=0.4, handlelength=0.8, loc=(0.57, 0.5))\n",
+ "save_nice_fig(fol+'Supp/scal_late.pdf')\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "s1 = plt.plot(scal1[:, 0], scal1[:, 14::2], '-', lw=1, c=dark_grey)\n",
+ "s2 = plt.plot(scal2[:, 0], scal2[:, 14::2], '-', lw=1, c=red, alpha=0.3)\n",
+ "s3 = plt.plot(scal3[:, 0], scal3[:, 14::2], '-', lw=1, c=green)\n",
+ "nice_fig('radial distance $r$ [$\\mathrm{\\mu m}$]', 'intensity [a.u]',\n",
+ " [0, 1.2], [0.35,1], [3,2])\n",
+ "save_nice_fig(fol+'Supp/scal_late_zoom.pdf')"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "s1 = plt.plot(scal1[:, 0], scal1[:, 2:10:2], '-', lw=1, c=dark_grey)\n",
+ "s2 = plt.plot(scal2[:, 0], scal2[:, 2:10:2], '-', lw=1, c=red)\n",
+ "s3 = plt.plot(scal3[:, 0], scal3[:, 2:10:2], '-', lw=1, c=green, alpha=0.4)\n",
+ "nice_fig('radial distance $r$ [$\\mathrm{\\mu m}$]', 'intensity [a.u]',\n",
+ " [0, 2], [0,1], [3,2])\n",
+ "plt.legend([s1[0], s2[0], s3[0]], ['Reference', '$D_\\mathrm{out}\\propto P^2$',\n",
+ " '$D_\\mathrm{out}\\propto P$'], frameon=False,\n",
+ " fontsize=9, handletextpad=0.4, handlelength=0.8, loc=(0.57, 0.54))\n",
+ "save_nice_fig(fol+'Supp/scal_early.pdf')"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "s1 = plt.plot(scal1[:, 0], scal1[:, 4:11:2], '-', lw=1, c=dark_grey, label='Original')\n",
+ "s2 = plt.plot(scal2[:, 0], scal2[:, 4:11:2], '-', lw=1, c=red, label='Original')\n",
+ "s3 = plt.plot(scal3[:, 0], scal3[:, 4:11:2], '-', lw=1, c=green, label='Original', alpha=0.4)\n",
+ "nice_fig('radial distance $r$ [$\\mathrm{\\mu m}$]', 'intensity [a.u]',\n",
+ " [0.97,1.03], [0, 0.35], [3,2])\n",
+ "save_nice_fig(fol+'Supp/scal_early_zoom.pdf')"
+ ]
+ },
{
"cell_type": "code",
"execution_count": null,
--
GitLab