function malujdens(dfile_L,dfile_A,nn,x0,x1,y0,y1);
% Matlab function for drawing zoom of the grid
  
  ncfile_L = netcdf(dfile_L,'nowrite');
  ncfile_A = netcdf(dfile_A,'nowrite');
  rho_L = ncfile_L{'rho'};
  rho_A = ncfile_A{'rho'};
  x_L = ncfile_L{'x'};
  y_L = ncfile_L{'y'};
  xc_L = ncfile_L{'xc'};
  yc_L = ncfile_L{'yc'};
  x_A = ncfile_A{'x'};
  y_A = ncfile_A{'y'};
  xc_A = ncfile_A{'xc'};
  yc_A = ncfile_A{'yc'};
  
  rmin = min(min(min(min(rho_L(:),rho_A(:)))));
  rmax = max(max(max(max(rho_L(:),rho_A(:)))));
  
  prho_L  = rho_L(nn,:,:)';
  prho_A  = rho_A(nn,:,:)';
  nic = size(prho_A,1);
  njc = size(prho_A,2);
  nin = nic + 1;
  njn = njc + 1;
% move the problem to [zero,zero] position  
  px_L   = (x_L(nn,:,:)' - 1.05*ones(nin,njn))*10000;
  py_L   = (y_L(nn,:,:)' - 1.05*ones(nin,njn))*10000;
  pxc_L  = (xc_L(nn,:,:)' - 1.05*ones(nic,njc))*10000;
  pyc_L  = (yc_L(nn,:,:)' - 1.05*ones(nic,njc))*10000;
  px_A   = (x_A(nn,:,:)' - 1.05*ones(nin,njn))*10000;
  py_A   = (y_A(nn,:,:)' - 1.05*ones(nin,njn))*10000;
  pxc_A  = (xc_A(nn,:,:)' - 1.05*ones(nic,njc))*10000;
  pyc_A  = (yc_A(nn,:,:)' - 1.05*ones(nic,njc))*10000;

  
  subplot(1,2,1);
%  hf = figure(1);
%  aa = axes;
%  set(aa,'FontSize',16);
%  set(hf,'Position',[100,100,1000,1000]);
%  pp = get(hf,'PaperPosition');
%  pp(2) = pp(2) - 1.4;
%  pp(4) = pp(4) - 3.5;
%  set(hf,'Position',[10,10,1500,1500*(y1-y0)/(x1-x0)],'PaperPosition',pp);
  hold on;
  plot(px_L,py_L,'-k','LineWidth',1);
  plot(px_L',py_L','-k','LineWidth',1);
  axis equal;
%  axis off;
  axis([x0 x1 y0 y1]); 
  if (1==0)
  pcolor(pxc_L,pyc_L,prho_L);
  shading interp;
% this doesnot work very well, draws in dual cells
% instead, do it as polygons; slower, but exctely what we want
  else
  for i=1:nic
    for j=1:njc
      n1x = px_L(i,j);
      n1y = py_L(i,j);
      n2x = px_L(i+1,j);
      n2y = py_L(i+1,j);
      n3x = px_L(i+1,j+1);
      n3y = py_L(i+1,j+1);
      n4x = px_L(i,j+1);
      n4y = py_L(i,j+1);
  
      fill([n1x n2x n3x n4x],[n1y n2y n3y n4y],prho_L(i,j));
%      [i j e1(i,j) e2(i,j) e3(i,j) e(i,j)]
    end
  end
  end
  colormap(jet);
  caxis([rmin rmax]);

  
  subplot(1,2,2);
%  hf = figure(1);
%  aa = axes;
%  set(aa,'FontSize',16);
%  set(hf,'Position',[100,100,1000,1000]);
%  pp = get(hf,'PaperPosition');
%  pp(2) = pp(2) - 1.4;
%  pp(4) = pp(4) - 3.5;
%  set(hf,'Position',[10,10,1500,1500*(y1-y0)/(x1-x0)],'PaperPosition',pp);
  hold on;
  plot(px_A,py_A,'-k','LineWidth',1);
  plot(px_A',py_A','-k','LineWidth',1);
  axis equal;
%  axis off;
  axis([x0 x1 y0 y1]); 
  if (1==0)
  pcolor(pxc_A,pyc_A,prho_A);
  shading interp;
% this doesnot work very well, draws in dual cells
% instead, do it as polygons; slower, but exctely what we want
  else
  for i=1:nic
    for j=1:njc
      n1x = px_A(i,j);
      n1y = py_A(i,j);
      n2x = px_A(i+1,j);
      n2y = py_A(i+1,j);
      n3x = px_A(i+1,j+1);
      n3y = py_A(i+1,j+1);
      n4x = px_A(i,j+1);
      n4y = py_A(i,j+1);
  
      fill([n1x n2x n3x n4x],[n1y n2y n3y n4y],prho_A(i,j));
%      [i j e1(i,j) e2(i,j) e3(i,j) e(i,j)]
    end
  end
  end
  
  colormap(jet);
  caxis([rmin rmax]);
  
%  cb = colorbar;
%  cb_values = [0 1/3 2/3 0.99];
%  cb_labels = {'0'; '1.e10'; '1.e11'; '8.e11'};
%  set(cb,'YTick',cb_values);
%  set(cb,'YTickLabel',cb_labels);
%  cbp = get(cb,'Position');
%  cbp(2) = cbp(2) + 0.01;
%  cbp(4) = cbp(4) - 0.02;
%  set(cb,'Position',cbp,'FontSize',16);
%  set(cb,'FontSize',16);
  
%  xlabel('X [\mum]','FontSize',16);
%  ylabel('Y [\mum]','FontSize',16);
  