function [xu,yu,zu,xv,yv,zv,xw,yw,zw,xp,yp,zp,...
  nxu,nyu,nzu,nxv,nyv,nzv,nxw,nyw,nzw,nxp,nyp,nzp] = read_grid_legacy(file,a,b,c,d,e,f,x_periodic,y_periodic,z_periodic)
  % [xu,yu,zu,xv,yv,zv,xw,yw,zw,xp,yp,zp,...
  %  nxu,nyu,nzu,nxv,nyv,nzv,nxw,nyw,nzw,nxp,nyp,nzp] = read_grid_legacy(file,a,b,c,d,e,f,x_periodic,y_periodic,z_periodic)
  % Reconstructs grid from 'current_points.tmp'.
  % Input
  %   file            path to 'current_points.tmp'
  %   a,b,c,d,e,f     domain boundaries
  %   x/y/z_periodic  periodicity
  % Output
  % xu,yu,zu,...      grid vectors
  % nxu,nyu,nzu,...   number of points
  
  nn=load(file);
  nxp=nn(1);
  nyp=nn(2);
  nzp=nn(3);
  [nxu,nyu,nzu,nxv,nyv,nzv,nxw,nyw,nzw,...
       xu,yu,zu,xv,yv,zv,xw,yw,zw,xp,yp,zp]=...
              generate_grid(a,b,c,d,e,f,nxp,nyp,nzp,x_periodic,y_periodic,z_periodic);
end