source: trunk/src/script_readnc.m

Last change on this file was 1119, checked in by sommeria, 20 months ago
File size: 557 bytes
Line 
1DataFolder='.fsnet/project/meige/2019/19TORE';
2fileinput={'Tore_a8b40_T22s_N047_A1cm.png.civ.mproj.tfilter_1'}
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.