plot_sim_data.m

function p = plot_sim_data(dmodel, model_data, sim_data, plot_params)
%function plot_sim_data(descr, model_data, sim_data, plot_params)
% plots the simulation data as returned by detailed_simulation()
%
% Parameters:
%   sim_data:    simulation data structure as returned by detailed_simulation()
%   plot_params: structure which controls the plot output
%
% Optional fields of plot_params:
%  plot:            function pointer to a function plotting a single snapshot
%                   for a specific time instant. (default = 'descr.plot')
%
% Return values:
%   p         : Matlab GUI handle to figure window
%
% function plotting simulation results

  if isfield(plot_params,'plot')
    pplot = plot_params.plot;
  else
    if isfield(dmodel.descr, 'plot')
      pplot = dmodel.descr.plot;
    else
      pplot = @fv_plot;
    end
  end

  UP_data.P = sim_data.P;
  UP_data.U = sim_data.U;
  UP_data.gEI = model_data.gEI;

  data = { sim_data.S, sim_data.S, sim_data.P, UP_data };
  plot_params.index = { 1, 2, 3, 4 };
  plot_params.plot_type = { 'contour', 'patch', 'patch', 'contour' };
  plot_params.plot = { pplot, pplot, pplot, @pu_plot };
  plot_params.sp_x = 2;
  plot_params.sp_y = 2;

  p = plot_subplot_sequence(data, model_data.grid, plot_params);

function p = pu_plot(grid, sim_data, plot_params)

  plot_params.plot_type = 'contour';
  p = fv_plot(grid, sim_data.P(:, sim_data.t), plot_params);
  hold on;
  NBI_pos_dir = grid.NBI > repmat((1:grid.nelements)', 1, 4) | grid.NBI < 0;
  UU = sim_data.U(:, sim_data.t);
  UU = UU(sim_data.gEI(NBI_pos_dir));
  ECX = grid.ECX(NBI_pos_dir);
  ECY = grid.ECY(NBI_pos_dir);
  NX  = grid.NX (NBI_pos_dir);
  NY  = grid.NY (NBI_pos_dir);
  quiver(ECX(:), ECY(:), UU.*NX(:), UU.*NY(:));
  hold off;