diff --git a/Plots_Droplet_FRAP.ipynb b/Plots_Droplet_FRAP.ipynb index ff1e3798783c52bd78064e5eba9c2bdeddbb4867..988644b8c2b3f73160756b763a54215d391a053b 100644 --- a/Plots_Droplet_FRAP.ipynb +++ b/Plots_Droplet_FRAP.ipynb @@ -100,12 +100,12 @@ "me = ['Meshes/multi_drop_gauss.xml', 'Meshes/multi_drop_gauss_med.xml',\n", " 'Meshes/multi_drop_gauss_far.xml', 'Meshes/multi_drop_gauss.xml',\n", " 'Meshes/multi_drop_gauss_med.xml', 'Meshes/multi_drop_gauss_far.xml']\n", - "point_lists = [[[4, 4.5, 0.5], [4, 3.5, 0.5], [3.5, 4, 0.5], [4.5, 4, 0.5]],\n", - " [[4, 5, 0.5], [4, 3, 0.5], [3, 4, 0.5], [5, 4, 0.5]],\n", - " [[4, 5.5, 0.5], [4, 2.5, 0.5], [2.5, 4, 0.5], [5.5, 4, 0.5]],\n", - " [[4, 4.5, 0.5], [4, 3.5, 0.5], [3.5, 4, 0.5], [4.5, 4, 0.5]],\n", - " [[4, 5, 0.5], [4, 3, 0.5], [3, 4, 0.5], [5, 4, 0.5]],\n", - " [[4, 5.5, 0.5], [4, 2.5, 0.5], [2.5, 4, 0.5], [5.5, 4, 0.5]]]\n", + "point_lists = [[[4, 4.5, 0.25], [4, 3.5, 0.25], [3.5, 4, 0.25], [4.5, 4, 0.25]],\n", + " [[4, 5, 0.25], [4, 3, 0.25], [3, 4, 0.25], [5, 4, 0.25]],\n", + " [[4, 5.5, 0.25], [4, 2.5, 0.25], [2.5, 4, 0.25], [5.5, 4, 0.25]],\n", + " [[4, 4.5, 0.25], [4, 3.5, 0.25], [3.5, 4, 0.25], [4.5, 4, 0.25]],\n", + " [[4, 5, 0.25], [4, 3, 0.25], [3, 4, 0.25], [5, 4, 0.25]],\n", + " [[4, 5.5, 0.25], [4, 2.5, 0.25], [2.5, 4, 0.25], [5.5, 4, 0.25]]]\n", "phi_tot_int = [.99, .99, .99, .9, .9, .9]\n", "phi_tot_ext = [.01, .01, .01, .1, .1, .1]\n", "G_in = [1, 1, 1, .1, .1, .1]\n", @@ -115,8 +115,8 @@ "\n", "for p, m, p_i, p_e, G_i, G_o in zip(point_lists, me, phi_tot_int,\n", " phi_tot_ext, G_in, G_out):\n", - " f = frap_solver([4, 4, 0.5], m, name='FRAP_multi_'+m[:-4]+str(G_i), point_list=p,\n", - " T=50, phi_tot_int=p_i, phi_tot_ext=p_e, G_in=G_i, G_out=G_o)\n", + " f = frap_solver([4, 4, 0.25], m, name='FRAP_multi/FRAP_multi_'+m[:-4]+str(G_i), point_list=p,\n", + " T=240, dt=0.02, phi_tot_int=p_i, phi_tot_ext=p_e, G_in=G_i, G_out=G_o)\n", " f.solve_frap()\n", " f_i.append(f)" ]