|
Last change
on this file since 1137 was
1131,
checked in by sommeria, 8 months ago
|
|
updated to read pivdata from fluidimage
|
|
File size:
578 bytes
|
| Rev | Line | |
|---|
| [1131] | 1 | DataFolder='.fsnet/project/coriolis/2024/24PLUME/VectorsEyemotion';
|
|---|
| [1119] | 2 | fileinput={'Tore_a8b40_T22s_N047_A1cm.png.civ.mproj.tfilter_1'}
|
|---|
| [890] | 3 | fileinput=uigetfile_uvmat('pick an input file',DataFolder);
|
|---|
| 4 | Data=nc2struct(fileinput);
|
|---|
| 5 | figure
|
|---|
| 6 | npy=numel(Data.coord_y);
|
|---|
| 7 | npx=numel(Data.coord_x);
|
|---|
| 8 | X=ones(npy,1)*Data.coord_x';
|
|---|
| 9 | Sum=sum(Data.VMean,2);
|
|---|
| 10 | XMean=sum(X.*Data.VMean,2)./Sum;
|
|---|
| 11 | XMean=XMean*ones(1,npx);
|
|---|
| 12 | Xrms=sum(((X-XMean).*(X-XMean).*Data.VMean),2)./Sum;
|
|---|
| 13 | Xrm=sqrt(Xrms);
|
|---|
| 14 | plot(Data.coord_y(3:end-1),Xrms(3:end-1))
|
|---|
| 15 | %0.024 y+0.14
|
|---|
| 16 | % spread= sqrt(2*log(2))*0.024=0.028
|
|---|
| 17 | %y0=-5.8;
|
|---|
| 18 |
|
|---|
Note: See
TracBrowser
for help on using the repository browser.
