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' ]);
  
  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
16			nc = netcdf([ fbname '_rot.nc' ],'write');
17			nxp1 = length( nc('Xp1') ); nx = nxp1 - 1;
18			nyp1 = length( nc('Yp1') ); ny = nyp1 - 1;
19			XCo = [];
20			XGo = [];
21			XCo = nc{'XC'}(:);
22			XGo = nc{'XG'}(:);
23			XC = XCo + rlon;
24			XG = XGo + rlon;
25			icp = find( XC(:) < -180.0 ); XC(icp) = XC(icp) + 360.0;
26			icm = find( XC(:) > 180.0 ); XC(icm) = XC(icm) - 360.0;
27			igp = find( XG(:) < -180.0 ); XG(icp) = XG(icp) + 360.0;
28			igm = find( XG(:) > 180.0 ); XG(icm) = XG(icm) - 360.0;
29			nc{'XC'}(:) = XC(:);
30			nc{'XG'}(:) = XG(:);
31
32			gid = fopen([ fbname '.mitgrid' ], 'r', 'ieee-be');
33			tmp = reshape(fread(gid,inf,'real8',0,'ieee-be'),[nxp1nyp1,16]);
34			XC = tmp(:,1) + rlon;
35			XG = tmp(:,6) + rlon;
36			icp = find( XC(:) < -180.0 ); XC(icp) = XC(icp) + 360.0;
37			icm = find( XC(:) > 180.0 ); XC(icm) = XC(icm) - 360.0;
38			igp = find( XG(:) < -180.0 ); XG(icp) = XG(icp) + 360.0;
39			igm = find( XG(:) > 180.0 ); XG(icm) = XG(icm) - 360.0;
40			tmp(:,1) = XC;
41			tmp(:,6) = XG;
42			fclose( gid );
43
44			gid = fopen([ fbname '_rot.mitgrid' ], 'w', 'ieee-be');
45			fwrite(gid,tmp,'real*8');
46			fclose( gid );
47
48			end
49