Parent Directory
|
Revision Log
|
Revision Graph
|
Patch
revision 1.1 by jmc, Tue Jun 28 02:34:54 2005 UTC | revision 1.2 by jmc, Thu Aug 18 18:26:15 2005 UTC | |
---|---|---|
# | Line 12 function [uE,vN] = rotate_csAg_EN(u,v) | Line 12 function [uE,vN] = rotate_csAg_EN(u,v) |
12 | Rac='grid_files/'; | Rac='grid_files/'; |
13 | ||
14 | NN=size(u); | NN=size(u); |
15 | [nnx nc nz]=size(u); | if length(NN) < 3, nz=1; |
16 | nPg=nnx*nc; | else nz=prod(NN(3:end)); |
17 | end | |
18 | nnx=NN(1); nc=NN(2); nPg=nnx*nc; | |
19 | if nnx ~= 6*nc, | if nnx ~= 6*nc, |
20 | fprintf('Error in CS-dim: %i %i %i \n',NN); | fprintf('Error in CS-dim: %i %i %i \n',NN); |
21 | return | return |
# | Line 34 for k=1:nz; | Line 36 for k=1:nz; |
36 | vN(:,k)=uvEN(:,2).*u(:,k)+uvEN(:,1).*v(:,k); | vN(:,k)=uvEN(:,2).*u(:,k)+uvEN(:,1).*v(:,k); |
37 | end | end |
38 | ||
39 | uE=reshape(uE,[nnx NN(2:end)]); | uE=reshape(uE,NN); |
40 | vN=reshape(vN,[nnx NN(2:end)]); | vN=reshape(vN,NN); |
41 | ||
42 | return | return |
|
ViewVC Help | |
Powered by ViewVC 1.1.22 |