diff --git a/Plots_Droplet_FRAP.ipynb b/Plots_Droplet_FRAP.ipynb
index a325c4cb4889dfab47f76e09651f9fc337a0c8d4..c09a66792e880b44a93e8ee994e5c71b0fd4630d 100644
--- a/Plots_Droplet_FRAP.ipynb
+++ b/Plots_Droplet_FRAP.ipynb
@@ -1,12 +1,5 @@
 {
  "cells": [
-  {
-   "cell_type": "markdown",
-   "metadata": {},
-   "source": [
-    "### FRAP geometries"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -22,11 +15,26 @@
     "import pandas as pd\n",
     "import seaborn as sns\n",
     "fol = '/Users/hubatsch/Desktop/DropletFRAP/Latex/Figures/'\n",
-    "sns.set_style(\"ticks\")\n",
+    "# sns.set_style(\"ticks\")\n",
     "rcParams['axes.linewidth'] = 0.75\n",
     "rcParams['xtick.major.width'] = 0.75\n",
     "rcParams['ytick.major.width'] = 0.75\n",
-    "rcParams['text.usetex']=True"
+    "# rcParams['text.usetex']=True"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Define colors\n",
+    "sns.set_palette(sns.color_palette(\"Set2\"))\n",
+    "grey = (0.6, 0.6, 0.6)\n",
+    "dark_grey = (0.2, 0.2, 0.2)\n",
+    "green = pa[0]\n",
+    "blue = pa[2]\n",
+    "red = pa[1]"
    ]
   },
   {
@@ -55,6 +63,25 @@
     "                transparent=True)"
    ]
   },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "### FRAP geometries"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# For radial average: define angles and radial spacing\n",
+    "alphas = np.linspace(0,2*np.pi, 20)\n",
+    "ns = np.c_[np.cos(alphas), np.sin(alphas), np.zeros(len(alphas))]\n",
+    "eps = np.linspace(0, 0.23, 100)"
+   ]
+  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -91,9 +118,6 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "alphas = np.linspace(0,2*np.pi, 20)\n",
-    "ns = np.c_[np.cos(alphas), np.sin(alphas), np.zeros(len(alphas))]\n",
-    "eps = np.linspace(0, 0.23, 100)\n",
     "profs = []\n",
     "for i in range(len(f_i)):\n",
     "#     if i>2:\n",
@@ -150,7 +174,7 @@
     "plt.plot([np.mean(x)/f_i[2].phi_tot_int for x in profs[2]],\n",
     "         lw=2, label='d=1.5', ls=':')\n",
     "plt.plot(range(0, 100), np.ones(100), linestyle='--', color='k')\n",
-    "plt.title('$\\Phi_{out}=0.01}$', size=12)\n",
+    "plt.title('$P=99}$', size=12)\n",
     "plt.gca().get_yaxis().set_visible(False)\n",
     "save_nice_fig(fol+'Fig3/tot_recov_neighbours_bad.pdf')"
    ]
@@ -169,7 +193,7 @@
     "plt.plot([np.mean(x)/f_i[5].phi_tot_int for x in profs[5]],\n",
     "         lw=2, label='d=1.5', ls=':')\n",
     "plt.plot(range(0, 100), np.ones(100), linestyle='--', color='k')\n",
-    "plt.title('$\\Phi_{out}=0.1}$', size=12)\n",
+    "# plt.title('$P=9$', size=12)\n",
     "plt.legend(prop={'size': 9}, frameon=False)\n",
     "save_nice_fig(fol+'Fig3/tot_recov_neighbours_good.pdf')"
    ]
@@ -255,7 +279,7 @@
     "    plt.plot([np.mean(x)/f.phi_tot_int for x in profs_cs[i]],\n",
     "             label='d='+str(z[i]), ls=ls[i], lw=2)\n",
     "plt.plot(range(0, 100), np.ones(100), linestyle='--', color='k')\n",
-    "plt.title('$\\Phi_{out}=0.01}$', size=12)\n",
+    "plt.title('$P=99$', size=12)\n",
     "plt.gca().get_yaxis().set_visible(False)\n",
     "save_nice_fig(fol+'Fig3/tot_recov_cs_bad.pdf')"
    ]
@@ -272,7 +296,7 @@
     "    plt.plot([np.mean(x)/f.phi_tot_int for x in profs_cs[i+3]],\n",
     "             label='h='+str(z[i]), lw=2, ls=ls[i])\n",
     "plt.plot(range(0, 100), np.ones(100), linestyle='--', color='k')\n",
-    "plt.title('$\\Phi_{out}=0.1}$', size=12)\n",
+    "plt.title('$P=9$', size=12)\n",
     "plt.legend(prop={'size': 9}, frameon=False)\n",
     "save_nice_fig(fol+'Fig3/tot_recov_cs_good.pdf')"
    ]