function [ierr] = plotELCOMinputFile(filename,datelim,bcsetlist,pltTheme)

% function [data,data_type,bcset_no,n_data] = plotELCOMinputFile(filename)
%
% funcion to plot ELCOM input file
%
% Inputs:
%		filename  : filename string
%		datelim	  : (optional) array of min and max dates to plot in CWR format
%		bcsetlist : (optional) structure giving bcset_no and bc_name
%						eg:  bcsetlist.bcset_no = [1 2]
%						eg:  bcsetlist.bc_name = [{'River A'},{'Dam Wall Outflow'} ]
%
% Outputs
%		ierr      : error flag
%
% Uses:
%		plotELCOMprofileBCFile.m
%		plotELCOMicFile.m
%		plotELCOMbcFile.m
%		plotELCOMupdateFile.m
%		readELCOMinputFile.m
%		readELCOMinputFileHeader.m
%		line2strs.m
%
% Written by C. Dallimore 27 July 04

ierr = 0;
close all

% Find pts within date range
if (nargin == 1 )
	datelim = [1970001 2020001];
	bcsetlist = [];
	pltTheme = 'default';
elseif (nargin == 2)
	bcsetlist = [];
	pltTheme = 'default';
elseif (nargin == 3)
	pltTheme = 'default';
end

if isempty(datelim)
	datelim = [1970001 2020001];
end

% Read file
[data,data_type,bcset_no,n_data,file_type,ierr] = readELCOMinputFile(filename);

% Check read was successfull
if (isempty(data))
	ierr = 1;
	return
end


% Read data
if (strcmp(file_type,'PROFILE'))
	disp(['Plotting profile bc file ',filename]);
	ierr = plotELCOMprofileBCFile(filename,data,data_type,bcset_no,n_data,datelim,bcsetlist,['ELCOMprofileBCFile',pltTheme]);
elseif (strcmp(file_type,'INITIAL'))
	disp(['Plotting initial profile file ',filename]);
	ierr = plotELCOMicFile(filename,data,data_type,bcset_no,n_data,['ELCOMicFile',pltTheme]);
elseif (strcmp(file_type,'BC'))
	disp(['Plotting bc file ',filename])
	ierr = plotELCOMbcFile(filename,data,data_type,bcset_no,n_data,datelim,bcsetlist,['ELCOMbcFile',pltTheme]);
elseif (strcmp(file_type,'UPDATE'))
	disp(['Plotting update file ',filename])
	ierr = plotELCOMupdateFile(filename,data,data_type,bcset_no,n_data,datelim,bcsetlist,['ELCOMupdateFile',pltTheme]);
elseif (strcmp(file_type,'INIT_HORIZ'))
	% We can't plot horizontal initial condition files because we don't have dx and dy
	disp(['Plotting of horizontal initial files not currently enabled: ',filename])
else
	disp(['ERROR in plotELCOMinputFile unknown file type: ',file_type,' for file ',filename]);
	ierr = 1;
	return
end