New mobility possible: D=D0*(1-a*phi). Not tested.

flory_hugg_pde.m

@@ -22,5 +22,7 @@ elseif strcmp(mode, 'Client')
 elseif strcmp(mode, 'Const_flux')
     g0 = Ternary_model.gamma0(x, a+t*v, b, e_g0, u_g0, v*t, e, u0);
     f = g0.*(1-pt)./pt.*(pt.*dudx-u.*gra_a)-j*u/pt;
+elseif strcmp(mode, 'Mobility_ratio')
+    f = u_g0.*(1-e_g0*pt).*(dudx-u.*gra_a./pt);
 end
 end
\ No newline at end of file

prob_laplace.m

@@ -274,7 +274,7 @@ for j = 1:length(bp)
 end
 
 %% Now for d = 0.5 across the entire domain bp=0, lb=0.02
-params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 'const', 0.5, 10, 7, 0,...
+params = {-5, b(7.7/3, 10^-6), 0.5, e(7.7/3), 'const', 0.5, 10, 7, 0, 0,...
           'Constituent', 0};
 t = [0, 0.05, 0.1];
 direc = 1;

run_jump_lengths.m

@@ -3,7 +3,7 @@ function run_jump_lengths(a, ad, chi, direc, l_max, name, nu, t_ind, u0)
 % ad ... prefactor of quadratic term in Stefano's mobility ansatz
 b = @(chi, nu) nu^(1/3)*sqrt(chi/(chi-2));
 e = @(chi) sqrt(3/8*(chi-2));
-params = {a, b(chi, nu), u0, e(chi), ad, 0, 10, 7, 0, 'Const_mob', 0};
+params = {a, b(chi, nu), u0, e(chi), ad, 0, 10, 7, 0, 0, 'Const_mob', 0};
 t = [0, 0.05, 0.1, 1];
 x0 = sort(5-direc*(0:0.001:4.01));
 %% Run simulations for 'delta' IC across outside