plot_sim can now work with two species.

@Ternary_model/plot_sim.m

@@ -18,26 +18,37 @@ if strcmp(mode, 'mov')
         axis([0, T.system_size, 0, max(T.sol(:))]);
         ax = gca; ax.FontSize = 12;
         xlabel('position'); ylabel('probability [not normalized]');
-        plot(T.x, T.phi_tot(T.x, T.a, T.b, T.e, T.u0, T.t(i)*T.v), 'LineWidth', 2, ...
-            'LineStyle', '--');
-        plot(T.x, T.sol(i, :), 'LineWidth', 2);
+        plot(T.x, T.phi_tot(T.x, T.a, T.b, T.e, T.u0, T.t(i)*T.v),...
+            'LineWidth', 2, 'LineStyle', '--');
+        if length(size(T.sol))==2
+            plot(T.x, T.sol(i, :), 'LineWidth', 2);
+        elseif length(size(T.sol))==3
+            plot(T.x, T.sol(i, :, 1), 'LineWidth', 2, 'Color', color);
+            plot(T.x, T.sol(i, :, 2), 'LineWidth', 2, 'Color', 'k');
+        end
         text(abs(-T.a+2), 0.7, num2str(T.t(i)));
         shg; pause();
 %         print([num2str(i),'.png'],'-dpng')
     end
 elseif strcmp(mode, 'plot')
-    figure(20);
-%     cla;
-    hold on;
+    figure(20); hold on;
     xlim([-inf, -2*T.a]); %ylim([-inf, max(T.sol(:))]);
     ax = gca;
     ax.FontSize = 12;
     xlabel('x [\mum]'); ylabel('volume fraction'); 
 %     N = max(max(T.sol([1, 2:nth:si(1)], :)));
     N = 1;
-    plot(T.x*x_scale, T.sol([1, start:nth:si(1)], :)/N, 'LineWidth', 1, 'Color', color);
     plot(T.x*x_scale, 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]);
+    if length(size(T.sol))==2
+        plot(T.x*x_scale, T.sol([1, start:nth:si(1)], :)/N, 'LineWidth', 1,...
+            'Color', color);
+    elseif length(size(T.sol))==3
+        plot(T.x*x_scale, T.sol([1, start:nth:si(1)], :, 1)/N, 'LineWidth', 1,...
+            'Color', color);
+        plot(T.x*x_scale, T.sol([1, start:nth:si(1)], :, 2)/N, 'LineWidth', 1,...
+            'Color', 'k');
+    end
 end
 end