
%this script asks you to enter the simrad mat file processed by 
%simwaymthen.pl and simthenplotter.m 
%you need to change the yearday for mingmt and maxgmt
% 16 sept 2002 smg

clear all
close all

file=input('What simrad mat file do you want to plot?   ','s');

load(file) %load the simrad mat file of choice

%you must set the appropriate yearday begin and end
mingmt=257
maxgmt=258

%Define color map same as for BiomaperII
	colormap(cmap(1,100));         % make color map for -100 to -55 dB when
	colmap=colormap;              % datrange = 1e-4 to 1e-8
	%colmap=[colmap;ones(10,3)];   % add white for range -55 to -50dB
	m=32;
	g=[12:m-1]'/max(m-1,1);
	g=1-g;
	g=[g g g];
	colmap=[colmap;g];            % add gray-scale for range -55 to -45dB
	colmap=[colmap;zeros(20,3)];  % add black for range -45 to -40dB
	%colmap=[colmap;zeros(19,3)];  % add black for range -45 to -41.25dB
	%colmap=[colmap;ones(1,3)];  % add white for -40 dB
	colormap(colmap);
	% datrange= max and min values for 10log(sv) data (Albatross IV data series)
	%datrange =[10*log10(1e-4) 10*log10(1e-10)]; %this sets range of color bar.
	%clegend5(datrange,'sv (dB)',.87)                %see below same values applied
%End Define color map

depth=([1:250]);

figure (1)
%PLOT 38
subplot(3,1,1)
%GET RIGHT DEPTHS DEPENDING ON DAY
if DAY<= 227,depth38=depth;
else depth38=depth.*2;
end;   
colormap(colmap);
imagesc(gmt(1:length(echo38),1),depth38,echo38,[-100 -40])
colorbar
title(['NBP0204 Simrad EK500 38 kHz on Day:  ', num2str(DAY(q,:)) ]);
axis([mingmt maxgmt 0 max(depth38)]);
ylabel('Depth (m)')
drawnow
orient landscape

%PLOT 120
subplot(3,1,2)
colormap(colmap);
imagesc(gmt(1:length(echo120),1),depth,echo120,[-100 -40])
colorbar
title(['NBP0204 Simrad EK500 120 kHz on Day:  ', num2str(DAY(q,:)) ]);
axis([mingmt maxgmt 0 250]);
ylabel('Depth (m)')
drawnow
orient landscape

%PLOT 200
subplot(3,1,3)
colormap(colmap);
imagesc(gmt(1:length(echo200),1),depth,echo200,[-100 -40])
colorbar
title(['NBP0204 Simrad EK500 200 kHz on Day:  ', num2str(DAY(q,:)) ]);
axis([mingmt maxgmt 0 250]);
xlabel('Year Day (GMT)')
ylabel('Depth (m)')
drawnow
orient landscape

