/[MITgcm]/MITgcm/utils/cs_grid/rotate_csAg_EN.m
ViewVC logotype

Diff of /MITgcm/utils/cs_grid/rotate_csAg_EN.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch 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

Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2

  ViewVC Help
Powered by ViewVC 1.1.22