comsol_evol_plot_sim_data.m

function p = comsol_evol_plot_sim_data(model,model_data,sim_data,params)
%function p = comsol_evol_plot_sim_data(model,model_data,sim_data,params)
%
% function plotting simulation results for evolution problems.
% This function needs a plot group being defined in the comsol_model!
% it uses COMSOL plot routines to plot the results!
%
% Optional fileds of params:
%   plot_title: Title of the plot
%   plot_output_title: title of the output plot
%   axis_equal: axis equal
%
% Oliver Zeeb, 2013/01/14



 sim_data.comsol_model=model_data.comsol_model;
% %get the current solution stored in the comsol model
% PVals_original = sim_data.comsol_model.sol('sol1').getPVals;
% U_original = zeros((sim_data.comsol_model.sol('sol1').getSize)');
% for i = 1:length(PVals_original)
%     U_original(:,i) = sim_data.comsol_model.sol('sol1').getU(i);
% end

%Write this solution into the model
PVals = 0:model.T/model.nt:model.T;
%set new solution
for i = 1:length(PVals)
    sim_data.comsol_model.sol('sol1').setU(i,sim_data.U(:,i));
end
%set new time list
sim_data.comsol_model.sol('sol1').setPVals(PVals);
%create new solution vector:
sim_data.comsol_model.sol('sol1').createSolution;


p = comsol_plot_sequence(sim_data,params);

% plot output functional
if model.compute_output_functional
  p2 = figure;
  plot(sim_data.s);
  title('output functional sim\_data.s')
  p = [p,p2];
end;



%reset original solution
% for i = 1:length(PVals)
%     sim_data.comsol_model.sol('sol1').setU(i,U_original(:,i));
% end
% %set new time list
% sim_data.comsol_model.sol('sol1').setPVals(PVals_original);
% %create new solution vector:
% sim_data.comsol_model.sol('sol1').createSolution;

%%%