include("eff_Drop.jl")

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.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)
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)