high_res_cube/eddy_flux/unmangleJPL1.m

function m = unmangleJPL1(A, ne, tx)

% Function m = unmangleJPL1(A, ne, tx)
% 
%
% INPUTS
% A        matrix in 1 layer of the rotten JPL format
% ne       number of values alone each edge of each face
% tx       number of values along each tile edge
%
% OUTPUTS
% m        output array of dimension ne,ne,6

if nargin ~= 3
  disp('Error: wrong number of arguments')
  exit
end

cx = ne;
cy = ne;
nt = ne*ne*6/(tx*tx);

m = zeros(ne,ne,6);
ite  = 0;
jte  = 1;
ntte = 1;
for i = 1:nt
  il = (i - 1)*tx + 1;
  ih = il + tx - 1;
  tepart = A(il:ih,:);
  ite = ite + 1;
  itel = (ite - 1)*tx + 1;
  iteh = itel + tx - 1;
  jtel = (jte - 1)*tx + 1;
  jteh = jtel + tx - 1;
  m(itel:iteh,jtel:jteh,ntte) = tepart;
  if iteh == cx
    ite = 0;
    jte = jte + 1;
    if jteh == cy
      jte = 1;
      ntte = ntte + 1;
    end
  end
end
1	function m = unmangleJPL1(A, ne, tx)
2
3	% Function m = unmangleJPL1(A, ne, tx)
4	%
5	%
6	% INPUTS
7	% A matrix in 1 layer of the rotten JPL format
8	% ne number of values alone each edge of each face
9	% tx number of values along each tile edge
10	%
11	% OUTPUTS
12	% m output array of dimension ne,ne,6
13
14	if nargin ~= 3
15	disp('Error: wrong number of arguments')
16	exit
17	end
18
19	cx = ne;
20	cy = ne;
21	nt = nene6/(tx*tx);
22
23	m = zeros(ne,ne,6);
24	ite = 0;
25	jte = 1;
26	ntte = 1;
27	for i = 1:nt
28	il = (i - 1)*tx + 1;
29	ih = il + tx - 1;
30	tepart = A(il:ih,:);
31	ite = ite + 1;
32	itel = (ite - 1)*tx + 1;
33	iteh = itel + tx - 1;
34	jtel = (jte - 1)*tx + 1;
35	jteh = jtel + tx - 1;
36	m(itel:iteh,jtel:jteh,ntte) = tepart;
37	if iteh == cx
38	ite = 0;
39	jte = jte + 1;
40	if jteh == cy
41	jte = 1;
42	ntte = ntte + 1;
43	end
44	end
45	end