Changeset 491 for trunk/src/calc_field.m
- Timestamp:
- Jul 11, 2012, 2:56:40 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/calc_field.m
r412 r491 67 67 for ilist=1:length(FieldList) 68 68 switch FieldList{ilist} 69 case {'u','v' }69 case {'u','v','velocity','norm_vel','ima_cor'} 70 70 check_grid=1;% needs a regular grid 71 71 case{'vort','div','strain'}% needs spatial derivatives spatial derivatives 72 72 check_der=1; 73 case {'velocity','norm_vel','ima_cor'};73 % case {'velocity','norm_vel','ima_cor'}; 74 74 otherwise 75 75 check_calc(ilist)=0; … … 119 119 coord_x=XMin:DataIn.Mesh:XMax; 120 120 coord_y=YMin:DataIn.Mesh:YMax; 121 npx=length(coord_x);122 npy=length(coord_y);121 % npx=length(coord_x); 122 % npy=length(coord_y); 123 123 DataOut.coord_x=[XMin XMax]; 124 124 DataOut.coord_y=[YMin YMax]; 125 125 [XI,YI]=meshgrid(coord_x,coord_y); 126 XI=reshape(XI,[],1); 127 YI=reshape(YI,[],1); 128 Coord_interp=[XI YI]; 129 end 130 126 % XI=reshape(XI,[],1); 127 % YI=reshape(YI,[],1); 128 Coord_interp=cat(3,XI,YI);%[XI YI]; 129 end 130 npx=size(Coord_interp,2); 131 npy=size(Coord_interp,1); 132 Coord_interp=reshape(Coord_interp,npx*npy,size(Coord_interp,3)); 133 % npy=length(coord_y); 131 134 %initialise output 132 135 nb_sites=size(Coord_interp,1);
Note: See TracChangeset
for help on using the changeset viewer.