ecco-godae/input/rotate_grid.m

%
%  Rotate the grid by -37 degrees longitude
%

rlon = -37.0;

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

i = 1;
for i = [ 1:5 ]
  
  fbname = sprintf('llc_p90_%d',i);
  eval([ '!cp ' fbname '.nc ' fbname '_rot.nc' ]);
  % eval([ '!cp ' fbname '.mitgrid ' fbname '_rot.mitgrid' ]);
  
  nc = netcdf([ fbname '_rot.nc' ],'write');
  nxp1 = length( nc('Xp1') );  nx = nxp1 - 1;
  nyp1 = length( nc('Yp1') );  ny = nyp1 - 1;
  XCo = [];
  XGo = [];
  XCo = nc{'XC'}(:);
  XGo = nc{'XG'}(:);
  XC = XCo + rlon;
  XG = XGo + rlon;
  icp = find( XC(:) < -180.0 );  XC(icp) = XC(icp) + 360.0;
  icm = find( XC(:) >  180.0 );  XC(icm) = XC(icm) - 360.0;
  igp = find( XG(:) < -180.0 );  XG(icp) = XG(icp) + 360.0;
  igm = find( XG(:) >  180.0 );  XG(icm) = XG(icm) - 360.0;
  nc{'XC'}(:) = XC(:);
  nc{'XG'}(:) = XG(:);

  gid = fopen([ fbname '.mitgrid' ], 'r', 'ieee-be');
  tmp = reshape(fread(gid,inf,'real*8',0,'ieee-be'),[nxp1*nyp1,16]);
  XC = tmp(:,1) + rlon;
  XG = tmp(:,6) + rlon;
  icp = find( XC(:) < -180.0 );  XC(icp) = XC(icp) + 360.0;
  icm = find( XC(:) >  180.0 );  XC(icm) = XC(icm) - 360.0;
  igp = find( XG(:) < -180.0 );  XG(icp) = XG(icp) + 360.0;
  igm = find( XG(:) >  180.0 );  XG(icm) = XG(icm) - 360.0;
  tmp(:,1) = XC;
  tmp(:,6) = XG;
  fclose( gid );

  gid = fopen([ fbname '_rot.mitgrid' ], 'w', 'ieee-be');
  fwrite(gid,tmp,'real*8');
  fclose( gid );
  
end

1	edhill	1.1	%
2			% Rotate the grid by -37 degrees longitude
3			%
4
5			rlon = -37.0;
6
7			gvars = { 'XC','YC','DXF','DYF','RA','XG','YG','DXV', ...
8			'DYU','RAZ','DXC','DYC','RAW','RAS','DXG','DYG' };
9
10			i = 1;
11			for i = [ 1:5 ]
12
13			fbname = sprintf('llc_p90_%d',i);
14			eval([ '!cp ' fbname '.nc ' fbname '_rot.nc' ]);
15			% eval([ '!cp ' fbname '.mitgrid ' fbname '_rot.mitgrid' ]);
16
17			nc = netcdf([ fbname '_rot.nc' ],'write');
18			nxp1 = length( nc('Xp1') ); nx = nxp1 - 1;
19			nyp1 = length( nc('Yp1') ); ny = nyp1 - 1;
20			XCo = [];
21			XGo = [];
22			XCo = nc{'XC'}(:);
23			XGo = nc{'XG'}(:);
24			XC = XCo + rlon;
25			XG = XGo + rlon;
26			icp = find( XC(:) < -180.0 ); XC(icp) = XC(icp) + 360.0;
27			icm = find( XC(:) > 180.0 ); XC(icm) = XC(icm) - 360.0;
28			igp = find( XG(:) < -180.0 ); XG(icp) = XG(icp) + 360.0;
29			igm = find( XG(:) > 180.0 ); XG(icm) = XG(icm) - 360.0;
30			nc{'XC'}(:) = XC(:);
31			nc{'XG'}(:) = XG(:);
32
33			gid = fopen([ fbname '.mitgrid' ], 'r', 'ieee-be');
34			tmp = reshape(fread(gid,inf,'real8',0,'ieee-be'),[nxp1nyp1,16]);
35			XC = tmp(:,1) + rlon;
36			XG = tmp(:,6) + rlon;
37			icp = find( XC(:) < -180.0 ); XC(icp) = XC(icp) + 360.0;
38			icm = find( XC(:) > 180.0 ); XC(icm) = XC(icm) - 360.0;
39			igp = find( XG(:) < -180.0 ); XG(icp) = XG(icp) + 360.0;
40			igm = find( XG(:) > 180.0 ); XG(icm) = XG(icm) - 360.0;
41			tmp(:,1) = XC;
42			tmp(:,6) = XG;
43			fclose( gid );
44
45			gid = fopen([ fbname '_rot.mitgrid' ], 'w', 'ieee-be');
46			fwrite(gid,tmp,'real*8');
47			fclose( gid );
48
49			end
50