%
%  Ed Hill
%
%  Quick script to double the lengths/areas around the face edges.

fin  = 'a.face001.bin';
fout = 'b.face001.bin';

sNx   = 102;
sNy   =  41;

vnall = { 'XC','YC','dxF','dyF','rA','XG','YG','dxV', ...
          'dyU','rAz','dxC','dyC','rAw','rAs','dxG','dyG' };

gvars = { 'XC','YC','DXF','DYF','RA','XG','YG','DXV', ...
          'DYU','RAZ','DXC','DYC','RAW','RAS','DXG','DYG' };

nslab = (sNx + 1)*(sNy + 1);

grid = {};
for ivar = 1:length(gvars)
  
  disp([ '  reading   :  ' gvars{ivar} ]);
  
  gid = fopen(fin, 'r', 'ieee-be');
  fseek(gid,(ivar-1)*nslab*8,'bof');
  tmp = reshape(fread(gid,nslab,'real*8',0,'ieee-be'),[sNx+1,sNy+1]);
  fclose(gid);
  grid.(gvars{ivar}) = tmp;
  
end

vn_fonly = { 'XC' 'YC' 'dxF' 'dyF' 'rA' 'XG' 'YG' 'dxG' 'dyG' };
vn_fedge = { 'rAw' 'rAs' 'rAz' 'dxV' 'dyU' 'dxC' 'dyC' };

gnew = grid;

%  RAW
gnew.RAW(1  ,:) = 2*grid.RAW(1  ,:);  % Bot
gnew.RAW(end,:) = 2*grid.RAW(end,:);  % Top
subplot(1,2,1), surf(grid.RAW)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.RAW)
title('new'), shading interp, caxis(cax);
%  RAS
gnew.RAS(:,1  ) = 2*grid.RAS(:,1  );  % L
gnew.RAS(:,end) = 2*grid.RAS(:,end);  % R
subplot(1,2,1), surf(grid.RAS)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.RAS)
title('new'), shading interp, caxis(cax);
%  RAZ
gnew.RAZ(1  ,:) = 2*gnew.RAZ(1  ,:);  % Bot
gnew.RAZ(end,:) = 2*gnew.RAZ(end,:);  % Top
gnew.RAZ(:,1  ) = 2*gnew.RAZ(:,1  );  % L
gnew.RAZ(:,end) = 2*gnew.RAZ(:,end);  % R
subplot(1,2,1), surf(grid.RAZ)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.RAZ)
title('new'), shading interp, caxis(cax);
%  DXV
gnew.DXV(1  ,:) = 2*grid.DXV(1  ,:);  % Bot
gnew.DXV(end,:) = 2*grid.DXV(end,:);  % Top
subplot(1,2,1), surf(grid.DXV)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.DXV)
title('new'), shading interp, caxis(cax);
%  DYU
gnew.DYU(:,1  ) = 2*grid.DYU(:,1  );  % L
gnew.DYU(:,end) = 2*grid.DYU(:,end);  % R
subplot(1,2,1), surf(grid.DYU)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.DYU)
title('new'), shading interp, caxis(cax);
%  DXC
gnew.DXC(1  ,:) = 2*grid.DXC(1  ,:);  % Bot
gnew.DXC(end,:) = 2*grid.DXC(end,:);  % Top
subplot(1,2,1), surf(grid.DXC)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.DXC)
title('new'), shading interp, caxis(cax);
%  DYC
gnew.DYC(:,1  ) = 2*grid.DYC(:,1  );  % L
gnew.DYC(:,end) = 2*grid.DYC(:,end);  % R
subplot(1,2,1), surf(grid.DYC)
title('old'), shading interp, cax = caxis;
subplot(1,2,2), surf(gnew.DYC)
title('new'), shading interp, caxis(cax);

fid = fopen(fout, 'w', 'ieee-be');
for i = 1:length(gvars)
  fwrite(fid,gnew.(gvars{i}),'real*8',0,'ieee-be');
end
fclose(fid);