Rev  Line  

[1119]  1  DataFolder='.fsnet/project/meige/2019/19TORE';


 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(((XXMean).*(XXMean).*Data.VMean),2)./Sum;


 13  Xrm=sqrt(Xrms);


 14  plot(Data.coord_y(3:end1),Xrms(3:end1))


 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.