%make as plot of "honey insertion"

%define file with grid information
 ifil = "ocean_physic_2000.nc"

%get grid information
 x=ncread(ifil,'xt_ocean');
 y=ncread(ifil,'yt_ocean');
 [X,Y]=meshgrid(x,y);
 d=ncread(ifil,'ht')'; I=find(d<-1e10); d(I)=NaN; d=double(d);

%choose lims
 xlims=360+[-70 -50];
 ylims=[-60 -54];

 xlims=360+[-70-20 -50+20];
 ylims=[-72 -54+20]; 

 ylims=[-69 -45];
 xlims=[280 310] 
% xlims=[min(x) max(x)];
% ylims=[min(y) max(y)];
 
    figure; drawnow
         in=[0:100:5200];
         in=[3200:100:4500];
         m_proj('Mercator','long',xlims,'lat',ylims);
          hold on
          [xxtt,yytt]=m_ll2xy(X,Y);
          flat_pcolor(xxtt,yytt,d), shading flat;
          amcolormap(flipud(ferretc('ocean_blue',length(in))),in); amcaxis(in);

          %laplacian horizontal viscosity
           %whole watercolumn
            %600m^2/s
             m_plot([0 360 360 0 0],[-10 -10 50 50 -10],'g','linewidth',2)
            %1200m^2/s
             m_plot([0 360 360 0 0],[50   50 60 60 50],'r--','linewidth',2)    
            %1800m^2/s
             m_plot([0 360 360 0 0],[60   60 max(y) max(y) 60],'y.','linewidth',2)

           %below level 45 (below ~3250m)
            %1000m^2/s

            m_plot([301.6 303.2 303.2 301.6 301.6],[-56 -56 -58.5 -58.5 -56],'r','linewidth',3.5)
      
          [cc]=amcolorbar('horiz',in);
          set(cc,'linewidth',2,'tickdir','out','fontweight','bold', ...
                  'xtick',in([1:10:end]),'xticklabel',in([1:10:end]),'position', ...
                  [0.130+0.185 0.1100+0.075 0.7750-0.4 0.0611/4], ...
                  'fontsize',10)
       %   m_grid('xtick',12,'tickdir','out','ytick',[-70:20:20],'linest','-', ...
       %          'color','k');

          m_grid('color','k')
          m_coast('patch',[0.7 0.7 0.7])
          %m_gshhs_h('patch',[0.7 0.7 0.7])
          legend('off')
          brighten(0.1)
          print -depsc drakeHoney.eps

