/[MITgcm]/MITgcm_contrib/high_res_cube/eddy_flux/unmangleJPL1.m
ViewVC logotype

Contents of /MITgcm_contrib/high_res_cube/eddy_flux/unmangleJPL1.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Tue Aug 3 15:54:26 2004 UTC (20 years, 11 months ago) by edhill
Branch: MAIN
CVS Tags: HEAD
 o sanitizing and checking in parts of the eddy flux diagnostics

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 = ne*ne*6/(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

  ViewVC Help
Powered by ViewVC 1.1.22