Changeset 491 for trunk/src/calc_field.m


Ignore:
Timestamp:
Jul 11, 2012, 2:56:40 PM (9 years ago)
Author:
sommeria
Message:

fix the the way to deal with filter fields using tps
fix the main projection plane in uvmat

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/calc_field.m

    r412 r491  
    6767for ilist=1:length(FieldList)
    6868    switch FieldList{ilist}
    69         case {'u','v'}
     69        case {'u','v','velocity','norm_vel','ima_cor'}
    7070            check_grid=1;% needs a regular grid
    7171        case{'vort','div','strain'}% needs spatial derivatives spatial derivatives
    7272            check_der=1;
    73         case {'velocity','norm_vel','ima_cor'};
     73%         case {'velocity','norm_vel','ima_cor'};
    7474        otherwise
    7575            check_calc(ilist)=0;
     
    119119        coord_x=XMin:DataIn.Mesh:XMax;
    120120        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);
    123123        DataOut.coord_x=[XMin XMax];
    124124        DataOut.coord_y=[YMin YMax];
    125125        [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);
    131134    %initialise output
    132135    nb_sites=size(Coord_interp,1);
Note: See TracChangeset for help on using the changeset viewer.