source: trunk/src/script_readnc.m @ 975

Last change on this file since 975 was 890, checked in by sommeria, 10 years ago

bugs solved for reading general netcdf files

File size: 610 bytes
Line 
1DataFolder=''; %TODO: put the actual path to start browser
2DataFolder='C:\Users\sommeria\Documents\MATLAB_WORK\UVMAT_DEMO_SOURCES\UVMAT_DEMO03_PIVchallenge_2005C';
3fileinput=uigetfile_uvmat('pick an input file',DataFolder);
4Data=nc2struct(fileinput);
5figure
6npy=numel(Data.coord_y);
7npx=numel(Data.coord_x);
8X=ones(npy,1)*Data.coord_x';
9Sum=sum(Data.VMean,2);
10XMean=sum(X.*Data.VMean,2)./Sum;
11XMean=XMean*ones(1,npx);
12Xrms=sum(((X-XMean).*(X-XMean).*Data.VMean),2)./Sum;
13Xrm=sqrt(Xrms);
14plot(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.