function [H,bathStruct] = plotELCOMbathy(filename)
% function [H, bathStruct] = plotELCOMbathy(filename)
%
% Inputs:
%		filename   : filename of bathymetry file
% Outputs
%       H          : Handle to graphic object
%		bathStruct : a matlab structure that contains all the data in the file
%
% Uses:
%		readELCOMbathy.m
%		pcolor_CD.m
%
% Written by C. Dallimore 1 Sept 04
%

% Get Bathymetry from file
bathStruct = readELCOMbathy(filename)
data = (bathStruct.bathData);

% Sort out x and y
if ~isempty(bathStruct.dx)
	dx1 = bathStruct.dx(1);
	x(1) = bathStruct.dx(1)/2.0;
	for i =2:bathStruct.x_rows
		x(i) = x(i-1)+bathStruct.dx(i-1)/2.0+bathStruct.dx(i)/2.0;
	end
else
	dx1 = bathStruct.x_grid;
	x(1) = bathStruct.x_grid/2.0;
	for i =2:bathStruct.x_rows
		x(i) = x(i-1)+bathStruct.x_grid;
	end
end
if ~isempty(bathStruct.dy)
	dy1 = bathStruct.dy(1);
	y(1) = bathStruct.dy(1)/2.0;
	for i =2:bathStruct.y_columns
		y(i) = y(i-1)+bathStruct.dy(i-1)/2.0+bathStruct.dy(i)/2.0;
	end
else
	dy1 = bathStruct.y_grid;
	y(1) = bathStruct.y_grid/2.0;
	for i =2:bathStruct.y_columns
		y(i) = y(i-1)+bathStruct.y_grid;
	end
end

% Remove land and open bc cells
pts = find(data == bathStruct.land_value | ...
		   data == bathStruct.open_value);
data(pts) = NaN;

H = pcolor_CD(y/1000,x/1000,data,dy1/1000,dx1/1000);
%H = pcolor(data);
set (gca,'ydir','reverse')
axis equal