ph=menu('Which phase?','build file','actual plot');
if ph==1;
  %load /data5/jmanning/ship/oc303/cnv/head.dat
  load /home3/ocn/jmanning/gbs95/barb/bongo97.mtl
  d=bongo97;

% for k=1:length(d);
% [la(k),lo(k)]=ddmm2dec(d(k,7),d(k,8));
% end;
  la=d(:,1);
  lo=d(:,2);

  load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_333.dat
  b=bongm97_333;


  y=[97*ones(length(lo),1) d(:,3) lo la b(:,9) b(:,10)];
% fploopy(y)

% sizesym(y(:,3),y(:,4),-1*y(:,6),1,100*y(:,6),'g');
end;

if ph==2;
       wnet=questdlg('Which net?','333 ','505 ','both');
       if wnet=='333 ';
         wsize=menu('Which size catagory?','1','2','3','combined');
         if wsize<4;
           % get lat/lon from one file & data from another (build "b")
           load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_505.dat
           b(:,1:4)=bongm97_505(:,1:4);
           load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_333_bysize.dat
           b(:,5:6)=bongm97_333_bysize(wsize:3:length(b)*3,10:11);
         else;
           load /home3/ocn/jmanning/gbs95/bongo/92_95/bongall.dat
           b=bongall;
           id=find(b(:,1)==97&b(:,2)>=96&b(:,2)<=145);
           b=b(id,:);
         end;
       elseif wnet=='505 ';
         wsize=menu('Which size catagory?','1','2','3','combined');
         if wsize<4;
           % get lat/lon from one file & data from another (build "b")
           load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_505.dat
           b(:,1:4)=bongm97_505(:,1:4);
           load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_505_bysize.dat
           b(:,5:6)=bongm97_505_bysize(wsize:3:length(b)*3,10:11);
         else;
           load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_505.dat
           b=bongm97_505;
         end;
       elseif wnet=='both';
         load /home3/ocn/jmanning/gbs95/bongo/92_95/bongall.dat
         b=bongall;
         id=find(b(:,1)==97);
         b=b(id,:);
         load /home3/ocn/jmanning/gbs95/bongo/92_95/bongm97_505.dat
         b(:,5)=(b(:,5)+bongm97_505(:,5))/2;
         b(:,6)=(b(:,6)+bongm97_505(:,6))/2;
       end;
       x=b(:,3);
       y=b(:,4);
       wfish=menu('Which fish?',' cod',' had ');
       wfish=wfish+4;
       z=b(:,wfish);
       ptype=menu('What type of plot?','Contour','sizeplot');
       if ptype==1;
         hsec       %(with two sections)
       elseif ptype==2;
         hold on;
         % add a value to standardize the symbol size
         x=[x;-70];y=[y;43];z=[z;5];
         sizesym(x,y,-1*z,3,1,'g');
       end;
end;

if wfish==5;fish='Cod';elseif wfish==6;fish='Haddock';end;
title('May 1997 Tidal Front Entrainment','fontsize',18)
global lw
lw=1;
sumplt
plotvel
text(-68.3,41.4,'CTD-derived sigma-t @30m (line contours)')
text(-68.3,41.35,'CTD-derived vertically-averaged velocity (arrows)')
text(-68.3,41.3,'333 Bongo-net larval abundance (shades)')
text(-68.3,41.25,'Stations positions (+)')
text(-67.95,40.77,'drifters','color','k')
axis([-68.75 -67.2 40.333 41.4])
xlabel('Longitude','fontsize',14)
ylabel('Latitude','fontsize',14)
colormap(jet)