[876] | 1 | %'diff_vel': calculate the difference of two input velocity fields.
|
---|
[870] | 2 | %
|
---|
[876] | 3 | % the second velocity field is linearly interpolated
|
---|
| 4 | % (after elimination of the vectors marked with an error flag) to the positions of
|
---|
| 5 | % the first one before subtraction. The ancilary data of the first field
|
---|
| 6 | % are preserved while those of the second one are lost.
|
---|
| 7 |
|
---|
[870] | 8 | %-----------------------------------------------------------------------
|
---|
[876] | 9 | % function SubData=diff_vel(Field,XmlData,Field_1)
|
---|
[870] | 10 | %
|
---|
| 11 | % OUPUT:
|
---|
| 12 | % SubData: structure representing the resulting field
|
---|
| 13 | %
|
---|
| 14 | % INPUT:
|
---|
| 15 | % Field: matlab structure representing the first field
|
---|
[876] | 16 | % XmlData: not used, needed for consistency with the call of transform fct.
|
---|
[870] | 17 | % Field_1:matlab structure representing the second field
|
---|
| 18 |
|
---|
| 19 | %=======================================================================
|
---|
| 20 | % Copyright 2008-2014, LEGI UMR 5519 / CNRS UJF G-INP, Grenoble, France
|
---|
| 21 | % http://www.legi.grenoble-inp.fr
|
---|
| 22 | % Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
|
---|
| 23 | %
|
---|
| 24 | % This file is part of the toolbox UVMAT.
|
---|
| 25 | %
|
---|
| 26 | % UVMAT is free software; you can redistribute it and/or modify
|
---|
| 27 | % it under the terms of the GNU General Public License as published
|
---|
| 28 | % by the Free Software Foundation; either version 2 of the license,
|
---|
| 29 | % or (at your option) any later version.
|
---|
| 30 | %
|
---|
| 31 | % UVMAT is distributed in the hope that it will be useful,
|
---|
| 32 | % but WITHOUT ANY WARRANTY; without even the implied warranty of
|
---|
| 33 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
---|
| 34 | % GNU General Public License (see LICENSE.txt) for more details.
|
---|
| 35 | %=======================================================================
|
---|
| 36 |
|
---|
| 37 | function SubData=diff_vel(Field,XmlData,Field_1)
|
---|
| 38 |
|
---|
| 39 | SubData=Field;
|
---|
| 40 | if exist('Field_1','var')
|
---|
| 41 | F.U=scatteredInterpolant(Field_1.X,Field_1.Y,Field_1.U,'linear');
|
---|
| 42 | SubData.U=Field.U-F.U(Field.X,Field.Y);%substract the interpolated ref to U
|
---|
| 43 | F.V=scatteredInterpolant(Field_1.X,Field_1.Y,Field_1.V,'linear');
|
---|
| 44 | SubData.V=Field.V-F.V(Field.X,Field.Y);%substract the interpolated ref to V
|
---|
| 45 | end
|
---|
| 46 | |
---|