1 | function DataOut=calc_tps(DataIn) |
---|
2 | DataOut=DataIn;%default |
---|
3 | SubDomain=1000; %default, estimated nbre of vectors in a subdomain used for tps |
---|
4 | if isfield(DataIn,'SubDomain') |
---|
5 | SubDomain=DataIn.SubDomain;% |
---|
6 | end |
---|
7 | [DataOut.SubRange,DataOut.NbSites,DataOut.Coord_tps,DataOut.U_tps,DataOut.V_tps] =... |
---|
8 | filter_tps([DataIn.X(DataIn.FF==0) DataIn.Y(DataIn.FF==0)],DataIn.U(DataIn.FF==0),DataIn.V(DataIn.FF==0),[],SubDomain,0); |
---|
9 | nbvar=numel(DataIn.ListVarName); |
---|
10 | DataOut.ListVarName=[DataIn.ListVarName {'SubRange','NbSites','Coord_tps','U_tps','V_tps'}]; |
---|
11 | DataOut.VarDimName=[DataIn.VarDimName {{'nb_coord','nb_bounds','nb_subdomain'},{'nb_subdomain'},... |
---|
12 | {'nb_tps','nb_coord','nb_subdomain'},{'nb_tps','nb_subdomain'},{'nb_tps','nb_subdomain'}}]; |
---|
13 | DataOut.VarAttribute{nbvar+3}.Role='coord_tps'; |
---|
14 | DataOut.VarAttribute{nbvar+4}.Role='vector_x'; |
---|
15 | DataOut.VarAttribute{nbvar+5}.Role='vector_y'; |
---|
16 | if isfield(DataOut,'ListDimName')%cleaning |
---|
17 | DataOut=rmfield(DataOut,'ListDimName'); |
---|
18 | end |
---|
19 | if isfield(DataOut,'DimValue')%cleaning |
---|
20 | DataOut=rmfield(DataOut,'DimValue'); |
---|
21 | end |
---|