Ignore:
Timestamp:
Jul 19, 2012, 11:50:39 AM (12 years ago)
Author:
sommeria
Message:

various bugs corrected after testing in Windows OS. Introduction
of filter tps

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/transform_field/phys.m

    r433 r494  
    159159        DataOut.V=(YOut_2-YOut_1)/Dt;
    160160    end
    161     if ~strcmp(Calib.CalibrationType,'rescale') && isfield(Data,'X_tps') && isfield(Data,'Y_tps')
    162         [DataOut.X_tps,DataOut.Y_tps]=phys_XYZ(Calib,Data.X,Data.Y,ZIndex);
    163     end
    164 end
     161%     if ~strcmp(Calib.CalibrationType,'rescale') && isfield(Data,'X_tps') && isfield(Data,'Y_tps')
     162%         [DataOut.X_tps,DataOut.Y_tps]=phys_XYZ(Calib,Data.X,Data.Y,ZIndex);
     163%     end
     164end
     165
     166%% suppress tps
     167list_tps={'Coord_tps'  'U_tps'  'V_tps'  'SubRange'  'NbSites'};
     168ind_remove=[];
     169for ilist=1:numel(list_tps)
     170    ind_tps=find(strcmp(list_tps{ilist},Data.ListVarName));
     171    if ~isempty(ind_tps)
     172        ind_remove=[ind_remove ind_tps];
     173        DataOut=rmfield(DataOut,list_tps{ilist});
     174    end
     175end
     176DataOut.ListVarName(ind_remove)=[];
     177DataOut.VarDimName(ind_remove)=[];
     178DataOut.VarAttribute(ind_remove)=[];
     179   
     180   
    165181
    166182%% transform of spatial derivatives: TODO check the case with plane angles
Note: See TracChangeset for help on using the changeset viewer.