function xg = plotall(direc, saveplot_fmt, plotfun, xg, visible)

narginchk(1,5)

cwd = fileparts(mfilename('fullpath'));

addpath([cwd, '/plotfunctions'])

addpath([cwd, '/../script_utils'])

validateattributes(direc, {'char'}, {'vector'}, mfilename, 'path to data', 1)

direc = absolute_path(direc);

assert(isfolder(direc), [direc, ' is not a directory'])

if nargin<2, saveplot_fmt={}; end %e.g. {'png'} or {'png', 'eps'}

if nargin<3

plotfun=[];

else

validateattributes(plotfun, {'char', 'function_handle'}, {'nonempty'}, mfilename, 'plotting function',3)

end

if nargin<4

xg=[];

else

validateattributes(xg, {'struct'}, {'scalar'}, mfilename, 'grid structure', 4)

end

if nargin<5, visible = 'on'; end

validateattributes(visible, {'char'}, {'vector'}, mfilename, 'plot visibility: on/off', 5)

lxs = simsize(direc);

disp(['sim grid dimensions: ',num2str(lxs)])

%% NEED TO READ INPUT FILE TO GET DURATION OF SIMULATION AND START TIME

params = read_config([direc, '/inputs']);

%% CHECK WHETHER WE NEED TO RELOAD THE GRID (check if one is given because this can take a long time)

if isempty(xg)

disp('Reloading grid...')

xg = readgrid([direc, '/inputs']);

end

plotfun = grid2plotfun(plotfun, xg);

%% TIMES OF INTEREST

times=params.UTsec0:params.dtout:params.UTsec0+params.tdur;

Nt=numel(times);

%% MAIN FIGURE LOOP

% NOTE: keep figure() calls in case plotfcn misses a graphics handle, and

% for Octave...

ymd(1,:) = params.ymd;

UTsec(1) = params.UTsec0;

for i = 2:Nt

[ymd(i,:), UTsec(i)] = dateinc(params.dtout, ymd(i-1,:), UTsec(i-1)); %#ok<AGROW>

end

h = plotinit(xg, visible);

if ~isempty(saveplot_fmt) % plot and save as fast as possible.

if ~isoctave && ~isinteractive

% if isempty(gcp('nocreate'))

% parpool(2); % don't use too much RAM loading huge data

% end

% parfor i = 1:Nt

for i=1:Nt

plotframe(direc, ymd(i,:), UTsec(i), saveplot_fmt, plotfun, xg, h);

end

else

for i = 1:Nt

plotframe(direc, ymd(i,:), UTsec(i), saveplot_fmt, plotfun, xg, h);

end

end

else % displaying interactively, not saving

for i = 1:Nt

xg = plotframe(direc, ymd(i,:), UTsec(i), saveplot_fmt, plotfun, xg, h);

drawnow % need this here to ensure plots update (race condition)

disp(''), disp('** press any key to plot next time step, or Ctrl C to stop**')

pause

end

end % if saveplots

if is_folder([direc, '/aurmaps']) % glow sim

plotglow(direc, saveplot_fmt)

end

%% Don't print

if nargout==0, clear('xg'), end

end % function