Change of parameters. Write solutions to file.

interface_resistance.jl

 include("eff_Drop.jl")
 
-P = 5; Δt = 0.001; D_i = 0.2; D_o = 0.2; BP=0.5; x0=0.5; 
-T=2; perc=0.5; L=0; R=2; k="none";
+P = 10; Δt = 0.001; D_i = 1; D_o = 50; BP=1; x0=0.5; 
+T=20; perc=0.5; L=0; R=10; k="none";
 
 ## Solve diffusion equation in phase-separated media
 A, u, x, b, Nb, t_perc, αi, αo= @time solve_D("3D","FRAP","Y","N",zeros(9),
 P,D_i,D_o,0,Δt,L,R,BP,x0,T,perc,zeros(9),zeros(9), k)
+u = u[:, vcat(1:Nb-1, Nb+2:size(u, 2))]
+x = vcat(collect(x[1:Nb-1]), collect(x[Nb+2:end]))
+writedlm(string("interface_ref_none.csv"), hcat(x, u'), ',')
 
 ##
-k = 0.001
+k = 0.05
 A, u1, x, b, Nb, t_perc, αi, αo= @time solve_D("3D","FRAP","Y","N",zeros(9),
 P,D_i,D_o,0,Δt,L,R,BP,x0,T,perc,zeros(9),zeros(9), k)
-##
-# Write results to file
-# writedlm( "x_solution.csv",  vcat(x', u), ',')
+u1 = u1[:, vcat(1:Nb-1, Nb+2:size(u1, 2))]
+x = vcat(collect(x[1:Nb-1]), collect(x[Nb+2:end]))
+writedlm(string("interface_ref_k_", k, ".csv"),hcat(x, u1'), ',')
 
 plot_model(x, u, Nb, T, Δt)
 plot_model(x, u1, Nb, T, Δt)