
%load data
 load ../data/processed/ref/surfmom.mat
 EEref=EE;
 load ../data/processed/noed/surfmom.mat
 EEnoe=EE;

%more ini
 xlim = [12.5 26] %[12 max(x)];
 ylim = [54 61] %min(y) max(y)];
 [X,Y]=meshgrid(x,y); clear x y

%define contour interval
 in = [50:10:150];

for k=[7] % 11]
 figure
  drawnow
        m_proj('Mercator','long', ...
                xlim, ...
               'lat',ylim);
         hold on
         [xxtt,yytt]=m_ll2xy(X,Y);
         flat_pcolor(xxtt,yytt,squee(EEref(:,:,k))./squee(EEnoe(:,:,k)).*100), shading flat;

         amcolormap(jet,in); amcaxis(in);

         [cc]=amcolorbar('horiz',in);
         set(cc,'linewidth',2,'tickdir','out','fontweight','bold', ...
                 'xtick',in(1:2:end),'xticklabel',in(1:2:end),'position', ...
                 [0.130+0.15 0.1100+0.075 0.7750-0.3 0.0611/3], ...
                 'fontsize',14)
%         m_gshhs_f('line','color','k')

         m_grid('box','fancy','tickdir','in','linestyle','none', ...
                'fontsize',8);

 brighten(0.5)
 if k==7
  print -r400 -depsc ../ps/DEKEnoed_jul.eps
 elseif k==11
  print -r400 -depsc ../ps/DEKEnoed_nov.eps
 end
end
