From 085b1be9b2c02cf1716bd25fae9ab02cfbd71819 Mon Sep 17 00:00:00 2001
From: Lars Hubatsch <hubatsch@pks.mpg.de>
Date: Thu, 4 Feb 2021 20:03:30 +0100
Subject: [PATCH] Plot Anatol's PGL-3 partitioning data on top of mine.

---
 Plots_Droplet_FRAP.ipynb | 18 +++++++++++-------
 1 file changed, 11 insertions(+), 7 deletions(-)

diff --git a/Plots_Droplet_FRAP.ipynb b/Plots_Droplet_FRAP.ipynb
index 3bb41a4..ff1e379 100644
--- a/Plots_Droplet_FRAP.ipynb
+++ b/Plots_Droplet_FRAP.ipynb
@@ -963,26 +963,30 @@
    "metadata": {},
    "outputs": [],
    "source": [
+    "P_anatol = np.array([686, 533, 590, 474])/np.array([1.6, 5.5, 7.5, 19.3])\n",
+    "salt_anatol = np.array([100, 150, 175, 200])\n",
     "P = PGL_50['P']\n",
     "def i_P(df, Pa):\n",
     "    return interp1d(np.array(df.groupby('P').mean()).flatten(),\n",
     "                    P.unique(), bounds_error=False,\n",
     "                    fill_value=\"extrapolate\")(Pa)\n",
-    "Ps_15 = [i_P(PGL, 30) for PGL in [PGL_50, PGL_60, PGL_75, PGL_90,\n",
+    "Ps_15 = [i_P(PGL, 2) for PGL in [PGL_50, PGL_60, PGL_75, PGL_90,\n",
     "                           PGL_100, PGL_120, PGL_150, PGL_180]]\n",
-    "Ps_30 = [i_P(PGL, 60) for PGL in [PGL_50, PGL_60, PGL_75, PGL_90,\n",
+    "Ps_30 = [i_P(PGL, 2) for PGL in [PGL_50, PGL_60, PGL_75, PGL_90,\n",
     "                           PGL_100, PGL_120, PGL_150, PGL_180]]\n",
     "salts = [50, 60, 75, 90, 100, 120, 150, 180]\n",
     "nice_fig('$c_\\mathrm{salt} \\; [\\mathrm{mM}]$', 'Partition coefficient $P$',\n",
-    "         [50, 180], [0,20000], [2.3,2])\n",
+    "         [50, 180], [50,2000], [2.3,2])\n",
     "plt.plot(salts, Ps_15, color=red)\n",
     "plt.plot(salts, Ps_30, color=red)\n",
+    "plt.plot(salt_anatol, P_anatol)\n",
     "plt.gca().fill_between(salts, Ps_15, Ps_30, color=red, alpha=0.2)\n",
     "plt.yticks(rotation=45)\n",
-    "plt.text(110, 6200, '$D_\\mathrm{out}=30\\;\\mathrm{\\mu m^2/s}$',\n",
-    "         color=red, size=8)\n",
-    "plt.text(54, 700, '$D_\\mathrm{out}=15\\;\\mathrm{\\mu m^2/s}$',\n",
-    "         color=red, size=8)\n",
+    "plt.yscale('log')\n",
+    "# plt.text(110, 6200, '$D_\\mathrm{out}=30\\;\\mathrm{\\mu m^2/s}$',\n",
+    "#          color=red, size=8)\n",
+    "# plt.text(54, 700, '$D_\\mathrm{out}=15\\;\\mathrm{\\mu m^2/s}$',\n",
+    "#          color=red, size=8)\n",
     "# save_nice_fig(fol+'Fig4/PGL-3_part.pdf')"
    ]
   },
-- 
GitLab