femwind/swap23.m

   1 function u=swap23(v)
   2 % swap indices to/from WRF ikj ordering
   3 % arrays of rank up to 4
   4 [n1,n2,n3,n4]=size(v);
   5 u=zeros(n1,n3,n2,n4);
   6 for l=1:n4
   7     for k=1:n3
   8         for j=1:n2
   9             for i=1:n1
  10                  u(i,k,j,l)=v(i,j,k,l);
  11             end
  12         end
  13     end
  14 end
  15 end
  16