Can calculate important experimental quantities.

@Ternary_model/Ternary_model.m

@@ -23,6 +23,7 @@ classdef Ternary_model < handle
         u_g0 = 0.2;
         system_size = 300;
         x0 = 1; % Center of Gauss initial condition
+        real_params;
     end
     methods
         function T = Ternary_model(geometry, ic, params, t, prec)
@@ -44,11 +45,13 @@ classdef Ternary_model < handle
             end
             T.create_mesh();
             T.phi_t = Ternary_model.phi_tot(T.x, T.a, T.b, T.e, T.u0);
+            T.calc_real_params;
         end
         create_mesh(T)
         plot_sim(T, mode, nth, color)
         solve_tern_frap(T)
         calc_time_scales(T);
+        calc_real_params(T);
     end
     methods(Static)
         pt = phi_tot(x, a, b, e, u0)

@Ternary_model/calc_real_params.m

+function calc_real_params(T)
+% Calculate important real-life parameters, such as D_in, D_out,
+% partitioning etc..
+T.real_params.phi_in = Ternary_model.phi_tot(-1000, T.a, T.b, T.e, T.u0);
+T.real_params.phi_out = Ternary_model.phi_tot(1000, T.a, T.b, T.e, T.u0);
+T.real_params.Gamma_in = Ternary_model.gamma0(-1000, T.a, T.b, T.e_g0,...
+                                               T.u_g0);
+T.real_params.Gamma_out = Ternary_model.gamma0(1000, T.a, T.b, T.e_g0,...
+                                               T.u_g0);
+T.real_params.D_in = (1-T.real_params.phi_in)*T.real_params.Gamma_in;
+T.real_params.D_out = (1-T.real_params.phi_out)*T.real_params.Gamma_out;
+T.real_params.partitioning = T.real_params.phi_in/T.real_params.phi_out;
+end
\ No newline at end of file