Adjusting to new input parameters.

@Ternary_model/create_mesh.m

@@ -25,7 +25,7 @@ if strcmp(T.mode, 'Constituent')
 
     % Fine mesh close to x0 if ic = 'Gauss'
     [~, ind] = min(abs(T.x-T.x0));
-    if T.x(ind+20)-T.x(ind) > 0.1
+    if (T.x(ind+20)-T.x(ind) > 0.1) && strcmp(T.ic, 'Gauss')
         T.x = [T.x(1:ind-21), linspace(T.x(ind-20), T.x(ind+20),...
                                        ceil((T.x(ind+20)-T.x(ind-20))*200)),...
                T.x(ind+21:end)];

@Ternary_model/plot_sim.m

@@ -32,7 +32,7 @@ elseif strcmp(mode, 'plot')
     xlabel('x [\mum]'); ylabel('volume fraction'); 
     plot(T.x, T.sol([1, 2:nth:si(1)], :), 'LineWidth', 1, 'Color', color);
 %     plot(T.x, T.sol(nth, :), 'LineWidth', 1, 'Color', color);
-    plot(T.x, Ternary_model.phi_tot(T.x, T.a, T.b, T.e, T.u0),...
+    plot(T.x, Ternary_model.phi_tot(T.x, T.a, T.b, T.e, T.u0, 0),...
         'LineWidth', 1, 'LineStyle', '--', 'Color', [0.97, 0.55, 0.03]);
 end
 end