osse/utils/testscalar.m

% converts U (rho), V (theta) and W in polar coords into 
% proper u,v, and w in cartesian coords, plus into theta and
% rho vector components in a more physical annulus layout

if (1) 
% load data
Wwhole=rdmds('W');
%Twhole=rdmds('T');
W=Wwhole(:,:,15);
% set null-space areas to NaN (makes interpolation happier)
W(find(W==0))=NaN;
%T(find(T==0))=NaN;
%W(120,:)=0.005;
[ydim xdim]=size(W);
else
xdim=31;
ydim=120;
[X Y ]=meshgrid(1:xdim,1:ydim);
W=X;
W(:,1:8)=nan;
W(:,31)=nan;
%W(:,9)=26;


end


% this determines resolution of interpolated fields
% 0.5 is perhaps a bit fine, but it looks pretty
%step=0.5;
% good for quiver plots
% step=1;
step=(xdim*2+1)/(xdim*2);


rhoi=1:31;
thetai=1:120;
z=1;


% convert w (effectively a scalar)
 w=cyl2cart(W,-xdim:step:xdim,-xdim:step:xdim);
 Wback=cart2cyl(w,thetai,rhoi);
% Tc=cyl2cart(T,-30:step:30,-30:step:30);
figure(1);imagesc(W(:,:,z));colorbar;
figure(2);imagesc(w(:,:,z));colorbar;
figure(3);imagesc(Wback(:,:,z));colorbar;
Werr=Wback-W;
figure(4);imagesc(Werr(:,:,z));colorbar;


1	afe	1.1	% converts U (rho), V (theta) and W in polar coords into
2			% proper u,v, and w in cartesian coords, plus into theta and
3			% rho vector components in a more physical annulus layout
4
5			if (1)
6			% load data
7			Wwhole=rdmds('W');
8			%Twhole=rdmds('T');
9			W=Wwhole(:,:,15);
10			% set null-space areas to NaN (makes interpolation happier)
11			W(find(W==0))=NaN;
12			%T(find(T==0))=NaN;
13			%W(120,:)=0.005;
14			[ydim xdim]=size(W);
15			else
16			xdim=31;
17			ydim=120;
18			[X Y ]=meshgrid(1:xdim,1:ydim);
19			W=X;
20			W(:,1:8)=nan;
21			W(:,31)=nan;
22			%W(:,9)=26;
23
24
25			end
26
27
28
29			% this determines resolution of interpolated fields
30			% 0.5 is perhaps a bit fine, but it looks pretty
31			%step=0.5;
32			% good for quiver plots
33			% step=1;
34			step=(xdim2+1)/(xdim2);
35
36
37			rhoi=1:31;
38			thetai=1:120;
39			z=1;
40
41
42
43
44			% convert w (effectively a scalar)
45			w=cyl2cart(W,-xdim:step:xdim,-xdim:step:xdim);
46			Wback=cart2cyl(w,thetai,rhoi);
47			% Tc=cyl2cart(T,-30:step:30,-30:step:30);
48			figure(1);imagesc(W(:,:,z));colorbar;
49			figure(2);imagesc(w(:,:,z));colorbar;
50			figure(3);imagesc(Wback(:,:,z));colorbar;
51			Werr=Wback-W;
52			figure(4);imagesc(Werr(:,:,z));colorbar;
53
54
55