eddy_flux/from_Dimitris/foo.m

% convert phi(nt*tx,ty,nz) to te(cx,cy,6,nz)

te=zeros(cx,cy,6,nz);
ite=0;jte=1;ntte=1;
for i=1:nt
 il=(i-1)*tx+1;
 ih=il+tx-1;
 tepart=phi(il:ih,:,:);
 ite=ite+1;
 itel=(ite-1)*tx+1;
 iteh=itel+tx-1;
 jtel=(jte-1)*ty+1;
 jteh=jtel+ty-1;
 te(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	% convert phi(nt*tx,ty,nz) to te(cx,cy,6,nz)
2
3	te=zeros(cx,cy,6,nz);
4	ite=0;jte=1;ntte=1;
5	for i=1:nt
6	il=(i-1)*tx+1;
7	ih=il+tx-1;
8	tepart=phi(il:ih,:,:);
9	ite=ite+1;
10	itel=(ite-1)*tx+1;
11	iteh=itel+tx-1;
12	jtel=(jte-1)*ty+1;
13	jteh=jtel+ty-1;
14	te(itel:iteh,jtel:jteh,ntte,:)=tepart;
15	if iteh == cx
16	ite = 0;
17	jte = jte+1;
18	if jteh == cy
19	jte=1;
20	ntte=ntte+1;
21	end
22	end
23	end