%  function to contour may95 bottemp
function bott(ld);
%  ld is the dist beweem label typical 150
%  of contour bot temp

% % % % read in the gridded results from surfer and contour
        fid=fopen('/home3/ocn/jmanning/gbs95/bongo/sj9507/bot.grd');
        fscanf(fid,'%s',1);
        nx=fscanf(fid,'%f',1);
        ny=fscanf(fid,'%f',1);
        nmilo=fscanf(fid,'%f',1);
        nmalo=fscanf(fid,'%f',1);
        nmila=fscanf(fid,'%f',1);
        nmala=fscanf(fid,'%f',1);
        mi=fscanf(fid,'%f',1);
        ma=fscanf(fid,'%f',1);
        loinc=(nmalo-nmilo)/nx;
        lainc=(nmala-nmila)/ny;

% % % % now we can form the xyz points
        for i=1:ny;
          for k=1:nx;
            t(i,k)=fscanf(fid,'%f',1);
            if t(i,k)>200000;
               t(i,k)=NaN;
            end;
            lo(i,k)=nmilo+k*loinc-loinc;
            la(i,k)=nmila+i*lainc-lainc;
          end;
        end;

%        econtour(lo(1,:),la(:,1),t,'fill','label',300);
%        ld=300;
        econtour(lo(1,:),la(:,1),t,'label',ld,'linewidth',2,'color','k');
        colormap('jet');
