mesh_points.m

hodge decomposition for 2D vector fields.

function xy=mesh_points(mv_x,mv_y) 
%% generate the Node Table 
%% mv_x: horizontal coordinates of the grids 
%% mv_y: vertical coordinates of the grids 
%% xy: the coordinates of all nodes 
sz=size(mv_x); 
MN=prod(sz); 
xx=reshape(mv_x,MN,1); 
yy=reshape(mv_y,MN,1); 
xy=zeros(MN,2);  
for i=1:MN 
xy(i,1)=mv_y(i); 
xy(i,2)=mv_x(i); 
end; 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%