Curvi/input/plot_grid.m

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

fin  = '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

fac = pi/180.0;
[ xg, yg, zg ] = sph2cart( grid.XG*fac, grid.YG*fac, 1 );

plot3(xg,yg,zg, 'b-')
hold on
plot3(xg',yg',zg', 'b-')
hold off

center = [ sum(sum(xg)) sum(sum(yg)) sum(sum(zg)) ];
center = center./norm(center);
hold on
plot3(center(1),center(2),center(3),'rx')
hold off

load coast
near = zeros([ length(lat) 3 ]);
[near(:,1) near(:,2) near(:,3)] = ...
    sph2cart(long*(pi/180),lat*(pi/180),1);
for i = 1:length(near)
  if norm(near(i,:) - center) > 0.4
    near(i,:) = NaN;
  end
end
hold on
plot3(near(:,1),near(:,2),near(:,3),'g-')
hold off

axis equal
view(32,36)


nr = 5;
xr = xg(1:nr:end,1:nr:end);
yr = yg(1:nr:end,1:nr:end);
zr = zg(1:nr:end,1:nr:end);
figure(2)
plot3(xr,yr,zr, 'b-')
hold on
plot3(xr',yr',zr', 'b-')
plot3(near(:,1),near(:,2),near(:,3),'g-')
hold off
axis equal
view(32,36)

print -depsc lab_sea_001.eps
1	%
2	% Ed Hill
3	%
4	% Quick script to double the lengths/areas around the face edges.
5
6	fin = 'b.face001.bin';
7
8	sNx = 102;
9	sNy = 41;
10
11	vnall = { 'XC','YC','dxF','dyF','rA','XG','YG','dxV', ...
12	'dyU','rAz','dxC','dyC','rAw','rAs','dxG','dyG' };
13
14	gvars = { 'XC','YC','DXF','DYF','RA','XG','YG','DXV', ...
15	'DYU','RAZ','DXC','DYC','RAW','RAS','DXG','DYG' };
16
17	nslab = (sNx + 1)*(sNy + 1);
18
19	grid = {};
20	for ivar = 1:length(gvars)
21
22	disp([ ' reading : ' gvars{ivar} ]);
23
24	gid = fopen(fin, 'r', 'ieee-be');
25	fseek(gid,(ivar-1)nslab8,'bof');
26	tmp = reshape(fread(gid,nslab,'real*8',0,'ieee-be'),[sNx+1,sNy+1]);
27	fclose(gid);
28	grid.(gvars{ivar}) = tmp;
29
30	end
31
32	fac = pi/180.0;
33	[ xg, yg, zg ] = sph2cart( grid.XGfac, grid.YGfac, 1 );
34
35	plot3(xg,yg,zg, 'b-')
36	hold on
37	plot3(xg',yg',zg', 'b-')
38	hold off
39
40	center = [ sum(sum(xg)) sum(sum(yg)) sum(sum(zg)) ];
41	center = center./norm(center);
42	hold on
43	plot3(center(1),center(2),center(3),'rx')
44	hold off
45
46	load coast
47	near = zeros([ length(lat) 3 ]);
48	[near(:,1) near(:,2) near(:,3)] = ...
49	sph2cart(long(pi/180),lat(pi/180),1);
50	for i = 1:length(near)
51	if norm(near(i,:) - center) > 0.4
52	near(i,:) = NaN;
53	end
54	end
55	hold on
56	plot3(near(:,1),near(:,2),near(:,3),'g-')
57	hold off
58
59	axis equal
60	view(32,36)
61
62
63	nr = 5;
64	xr = xg(1:nr:end,1:nr:end);
65	yr = yg(1:nr:end,1:nr:end);
66	zr = zg(1:nr:end,1:nr:end);
67	figure(2)
68	plot3(xr,yr,zr, 'b-')
69	hold on
70	plot3(xr',yr',zr', 'b-')
71	plot3(near(:,1),near(:,2),near(:,3),'g-')
72	hold off
73	axis equal
74	view(32,36)
75
76	print -depsc lab_sea_001.eps