%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