high_res_cube/eddy_flux/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	edhill	1.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