dyncore_ASP/csgrids/grid_shift_x.m


kwr=1;
nc=32; np=nc+1; n6x=6*nc; nPg=n6x*nc;
xShift=+180;

%inpName='tile.mitgrid';
%rDir='/home/jmc/exp/inp_CS_grid/';
inpName='dxC3_dXYa';
rDir='../grid_cs_input/';
%outpName='tile.newgrid';
outpName=[inpName,'_x0'];

for n=1:6,
%-- read :
%if strcmp(inpName(1:4),'tile'),
%if strcmp(inpName(5:end),'.mitgrid'),
 p=strfind(inpName,'.');
 if isempty(p),
   namF=sprintf([rDir,inpName,'.face%3.3i.bin'],n);
 else
   i=min(p);
   namF=sprintf([rDir,inpName(1:i-1),'%3.3i',inpName(i:end)],n);
 end
 fid=fopen(namF,'r','b');
 vv1=fread(fid,'real*8');
 fclose(fid);
 s=size(vv1,1); k1=s/np/np;
 fprintf(['read: ',namF,' : size: %i (%ix%ix%i)\n'],s,np,np,k1);
 vv1=reshape(vv1,[np np k1]);
%-- do not add angle:
 vv2=vv1; k2=k1;
%- XC YC DXF DYF RA XG YG DXV DYU RAZ DXC DYC RAW RAS DXG DYG:
%------------------------------------------------------------------
%- shift long:
vv2(:,:,1)=vv2(:,:,1)+xShift;
vv2(:,:,6)=vv2(:,:,6)+xShift;

%-- write:
 if kwr == 1,
% if ~strcmp(outpName(1:4),'tile'),
  p=strfind(outpName,'.');
  if isempty(p),
    namW=sprintf([outpName,'.face%3.3i.bin'],n);
  else
    namW=sprintf([outpName(1:p-1),'%3.3i',outpName(p:end)],n);
  end
  fid=fopen(namW,'w','b');
  fwrite(fid,vv2,'real*8');
  fclose(fid);
  fprintf([' write to file: ',namW,' %i 2D.var(%ix%i)\n'],k2,np,np);
 end
end
1
2	kwr=1;
3	nc=32; np=nc+1; n6x=6nc; nPg=n6xnc;
4	xShift=+180;
5
6	%inpName='tile.mitgrid';
7	%rDir='/home/jmc/exp/inp_CS_grid/';
8	inpName='dxC3_dXYa';
9	rDir='../grid_cs_input/';
10	%outpName='tile.newgrid';
11	outpName=[inpName,'_x0'];
12
13	for n=1:6,
14	%-- read :
15	%if strcmp(inpName(1:4),'tile'),
16	%if strcmp(inpName(5:end),'.mitgrid'),
17	p=strfind(inpName,'.');
18	if isempty(p),
19	namF=sprintf([rDir,inpName,'.face%3.3i.bin'],n);
20	else
21	i=min(p);
22	namF=sprintf([rDir,inpName(1:i-1),'%3.3i',inpName(i:end)],n);
23	end
24	fid=fopen(namF,'r','b');
25	vv1=fread(fid,'real*8');
26	fclose(fid);
27	s=size(vv1,1); k1=s/np/np;
28	fprintf(['read: ',namF,' : size: %i (%ix%ix%i)\n'],s,np,np,k1);
29	vv1=reshape(vv1,[np np k1]);
30	%-- do not add angle:
31	vv2=vv1; k2=k1;
32	%- XC YC DXF DYF RA XG YG DXV DYU RAZ DXC DYC RAW RAS DXG DYG:
33	%------------------------------------------------------------------
34	%- shift long:
35	vv2(:,:,1)=vv2(:,:,1)+xShift;
36	vv2(:,:,6)=vv2(:,:,6)+xShift;
37
38	%-- write:
39	if kwr == 1,
40	% if ~strcmp(outpName(1:4),'tile'),
41	p=strfind(outpName,'.');
42	if isempty(p),
43	namW=sprintf([outpName,'.face%3.3i.bin'],n);
44	else
45	namW=sprintf([outpName(1:p-1),'%3.3i',outpName(p:end)],n);
46	end
47	fid=fopen(namW,'w','b');
48	fwrite(fid,vv2,'real*8');
49	fclose(fid);
50	fprintf([' write to file: ',namW,' %i 2D.var(%ix%i)\n'],k2,np,np);
51	end
52	end