function TL_ovlay(Odat,Mdat,freq_lims,tsize,lcols,scols,clim,iinewfig)
%
% Call TL_ovlay(Odat,Mdat,freq_lims,tsize,lcols,scols,clim,iinewfig);
%
% Odat      = Structure containing ORCA TL at frequencies desired;
% Mdat      = Structure containing measured TL or, if a string, .mat file root;
% freq_lims = Frequency limits [fmin,fmax];
% tsize     = Plot size in [width,height];
% lcols     = Line colors for modeled TL;
% scols     = Symbol colors for measured TL;
% clim      = Vertical axis limits [TLmin,TLmax];
% iinewfig  = Flag to make new figure;
%
%
% COPYRIGHT (c) 2008
% 
% Applied Research Laboratories, The University of Texas at Austin
% (ARL:UT)
% All rights reserved.
% 
% This material may be reproduced by or for the U.S. Government
% only. This notice must appear in all copies of this file and its
% derivatives. Further distribution outside the U.S. Government is
% prohibited without the explicit written authorization from the
% Director of the Environmental Sciences Laboratory, ARL:UT.
% 
% THIS SOFTWARE IS PROVIDED "AS IS" WITHOUT WARRANTY OF ANY KIND,
% EITHER EXPRESSED OR IMPLIED, INCLUDING, BUT NOT LIMITED TO THE
% IMPLIED WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR
% PURPOSE, OR NONINFRINGEMENT.
% 
% IN NO EVENT WILL ARL:UT BE LIABLE FOR ANY DAMAGES WHATSOEVER ARISING
% OUT OF THE USE OR INABILITY TO USE THIS SOFTWARE, EVEN IF ARL:UT HAS
% BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGES.
if nargin<3, freq_lims=[]; end
if nargin<4 | isempty(tsize), 
    tsize=[3.5,1.75]; tsize=[5,3]; 
end
if nargin<5 | isempty(lcols), lcols='rbgmck'; end
if nargin<6 | isempty(scols), scols='rbgmck'; end
if nargin<7 | isempty(clim), clim=[]; end
if nargin<8 | isempty(iinewfig), iinewfig=1; end
%
if ischar(Mdat),
    str=['load ',Mdat],
    eval(str); Mdat=DATA;
end
jzr=1; jzs=1;
ncol=length(lcols);
if iinewfig, figure; else; clf; end
[nzr,nrng,nzs,nf]=size(Odat.data);
rmin_data=min(abs(Mdat.scales{1}));
rmax_data=max(abs(Mdat.scales{1}));
iir=find(Odat.scales{2}>=rmin_data & Odat.scales{2}<=rmax_data);
%
iif=ekw_var_spec(freq_lims,Odat.scales{4}); niif=length(iif);
if niif==1, lcols='r'; scols='b'; end
%
for jff=1:niif
    jf=iif(jff);
    jcol=mod(jff-1,length(lcols)) + 1;
    plot_cfill(Odat.scales{2}(iir),-Odat.data(jzr,iir,jzs,jf),'none',lcols(jcol),5,'-',1.75); hold on;
end
for jff=1:niif
    jf=iif(jff);
    freq=Odat.scales{4}(jf);
    [df,jf2]=min(abs(freq - Mdat.scales{2}));
    jcol=mod(jff-1,length(scols)) + 1;
    plot_cfill(abs(Mdat.scales{1}),Mdat.data(:,jf2),'o',scols(jcol),5);
end
axis('ij'); grid on; zoom on;
plot_size(gca,'in',tsize(1),tsize(2)); plot_set3(gca,1.5,'Arial',14,1);
%
ax=axis; ax(2)=rmax_data; 
if ~isempty(clim), 
    if any(clim<0), clim=-clim(2:-1:1); end 
    ax(3:4)=clim; 
end; 
axis(ax);
%
xlabel(Odat.labels{2}); ylabel(Odat.dlab);
str=[Odat.labels{1},' = ',num2str(Odat.scales{1}(jzr)),'; ',Odat.labels{3},' = ',num2str(Odat.scales{3}(jzs)), ...
    '; ',Odat.labels{4},' = ',num2str(Odat.scales{4}(iif))];
h=title(str); set(h,'FontSize',12);
