matlab_class/gcmfaces_convert/convert2vector.m

function [a]=convert2vector(b,method);
%[a]=CONVERT2VECTOR(b,method); 
%   converts gcmfaces object b to vector format a -- and vice versa 
%   when b is instead a vector obtained ealier using convert2vector.
%
%   With the 'new' method: each 2D field becomes a column vector, 
%   while the other dimensions remain the same; convert2gcmfaces is used.
%   With the 'old' method (default): convert2array is used and the third+
%   dimensions of b get conflated in the column vector length.
%
%   The 'old' method will get removed after updating routines that use convert2vector.

global mygrid;

if isa(b,'gcmfaces'); do_gcmfaces2vector=1; else; do_gcmfaces2vector=0; end;
if isempty(whos('method')); method='old'; end;

if strcmp(method,'new');
if do_gcmfaces2vector;
  bb=convert2gcmfaces(b);
  siz=size(bb); if length(siz)==2; siz=[siz 1]; end;
  a=reshape(bb,[prod(siz(1:2)) siz(3:end)]);
else;
  bb=convert2gcmfaces(mygrid.XC);
  siz=size(b);
  bb=reshape(b,[size(bb) siz(2:end)]);
  a=convert2gcmfaces(bb);
end;
return;
end;

if do_gcmfaces2vector;
  bb=convert2array(b);
  a=bb(:);
else;
  bb=convert2array(mygrid.XC);
  if mod(length(b(:)),length(bb(:)))~=0;
      error('vector length is inconsistent with gcmfaces objects');
  else;
      n3=length(b(:))/length(bb(:));
  end;
  b=reshape(b,[size(bb) n3]);
  a=convert2array(b);
end;


1	gforget	1.2	function [a]=convert2vector(b,method);
2			%[a]=CONVERT2VECTOR(b,method);
3			% converts gcmfaces object b to vector format a -- and vice versa
4			% when b is instead a vector obtained ealier using convert2vector.
5	gforget	1.1	%
6	gforget	1.2	% With the 'new' method: each 2D field becomes a column vector,
7			% while the other dimensions remain the same; convert2gcmfaces is used.
8			% With the 'old' method (default): convert2array is used and the third+
9			% dimensions of b get conflated in the column vector length.
10			%
11			% The 'old' method will get removed after updating routines that use convert2vector.
12	gforget	1.1
13			global mygrid;
14
15			if isa(b,'gcmfaces'); do_gcmfaces2vector=1; else; do_gcmfaces2vector=0; end;
16	gforget	1.2	if isempty(whos('method')); method='old'; end;
17
18			if strcmp(method,'new');
19			if do_gcmfaces2vector;
20			bb=convert2gcmfaces(b);
21			siz=size(bb); if length(siz)==2; siz=[siz 1]; end;
22			a=reshape(bb,[prod(siz(1:2)) siz(3:end)]);
23			else;
24			bb=convert2gcmfaces(mygrid.XC);
25			siz=size(b);
26			bb=reshape(b,[size(bb) siz(2:end)]);
27			a=convert2gcmfaces(bb);
28			end;
29			return;
30			end;
31	gforget	1.1
32			if do_gcmfaces2vector;
33			bb=convert2array(b);
34			a=bb(:);
35			else;
36			bb=convert2array(mygrid.XC);
37			if mod(length(b(:)),length(bb(:)))~=0;
38			error('vector length is inconsistent with gcmfaces objects');
39			else;
40			n3=length(b(:))/length(bb(:));
41			end;
42			b=reshape(b,[size(bb) n3]);
43			a=convert2array(b);
44			end;
45
46
47