osse/utils/test.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

% load data
U=rdmds('U');
V=rdmds('V');
W=rdmds('W');
T=rdmds('T');
% set null-space areas to NaN (makes interpolation happier)
U(find(U==0))=NaN;
V(find(V==0))=NaN;
W(find(W==0))=NaN;
T(find(T==0))=NaN;

% 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;

% convert theta/rho vector fields to u/v
[u v]=cyl2cartuv(U,V,-30:step:30,-30:step:30);

% convert w (effectively a scalar)
 w=cyl2cart(W,-30:step:30,-30:step:30);
% Tc=cyl2cart(T,-30:step:30,-30:step:30);

% interpolate theta and rho into annulus
 azim=cyl2cart(U,-30:step:30,-30:step:30);
 radial=cyl2cart(V,-30:step:30,-30:step:30);
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	% load data
6	U=rdmds('U');
7	V=rdmds('V');
8	W=rdmds('W');
9	T=rdmds('T');
10	% set null-space areas to NaN (makes interpolation happier)
11	U(find(U==0))=NaN;
12	V(find(V==0))=NaN;
13	W(find(W==0))=NaN;
14	T(find(T==0))=NaN;
15
16	% this determines resolution of interpolated fields
17	% 0.5 is perhaps a bit fine, but it looks pretty
18	step=0.5;
19	% good for quiver plots
20	% step=1;
21
22	% convert theta/rho vector fields to u/v
23	[u v]=cyl2cartuv(U,V,-30:step:30,-30:step:30);
24
25	% convert w (effectively a scalar)
26	w=cyl2cart(W,-30:step:30,-30:step:30);
27	% Tc=cyl2cart(T,-30:step:30,-30:step:30);
28
29	% interpolate theta and rho into annulus
30	azim=cyl2cart(U,-30:step:30,-30:step:30);
31	radial=cyl2cart(V,-30:step:30,-30:step:30);