Changeset 1068 for trunk/src/series/merge_proj_polar.m
 Timestamp:
 Jul 9, 2019, 10:10:18 AM (5 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/src/series/merge_proj_polar.m
r1061 r1068 63 63 ParamOut.AllowInputSort='on';% allow alphabetic sorting of the list of input file SubDir (options 'off'/'on', 'off' by default) 64 64 ParamOut.WholeIndexRange='off';% prescribes the file index ranges from min to max (options 'off'/'on', 'off' by default) 65 ParamOut.NbSlice='o n'; %nbre of slices ('off' by default)65 ParamOut.NbSlice='off'; %nbre of slices ('off' by default) 66 66 ParamOut.VelType='one';% menu for selecting the velocity type (options 'off'/'one'/'two', 'off' by default) 67 67 ParamOut.FieldName='off';% menu for selecting the field (s) in the input file(options 'off'/'one'/'two', 'off' by default) … … 268 268 CheckOverwrite=Param.CheckOverwrite; 269 269 end 270 270 NbField 271 271 for index=1:NbField 272 disp(['index=' num2str(index)])273 disp(['ellapsed time ' num2str(toc(tstart)/60,4) ' minutes'])274 272 update_waitbar(WaitbarHandle,index/NbField) 275 273 if ~isempty(RUNHandle) && ~strcmp(get(RUNHandle,'BusyAction'),'queue') … … 352 350 %% calculate tps coefficients 353 351 Data{iview}=tps_coeff_field(Data{iview},1); 354 352 'tps_coeff done' 355 353 %% projection on the polar grid 356 354 [DataOut,VarAttribute,errormsg]=calc_field_tps(Data{iview}.Coord_tps,Data{iview}.NbCentre,Data{iview}.SubRange,... 357 355 cat(3,Data{iview}.U_tps,Data{iview}.V_tps),FieldNames,cat(3,XI,YI)); 356 if ~isempty(errormsg) 357 disp(errormsg) 358 end 359 358 360 % set to NaN interpolation points which are too far from any initial data (more than 2 CoordMesh) 359 361 Coord=permute(Data{iview}.Coord_tps,[1 3 2]); … … 366 368 G=TriScatteredInterp(Coord,Coord(:,2),'nearest'); 367 369 end 370 'Interp done' 368 371 Distx=F(XI,YI)XI;% diff of x coordinates with the nearest measurement point 369 372 Disty=G(XI,YI)YI;% diff of y coordinates with the nearest measurement point … … 390 393 391 394 %% merge the NbView fields 395 ProjData 392 396 [MergeData,errormsg]=merge_field(ProjData); 393 397 if ~isempty(errormsg) … … 446 450 TimeData.div=MergeData.div; 447 451 448 [error,ncid]=struct2nc(OutputFile,TimeData);%save result file452 error=struct2nc(OutputFile,TimeData);%save result file 449 453 if isempty(error) 450 454 disp(['output file ' OutputFile ' written']) … … 453 457 end 454 458 ellapsed_time=toc(tstart); 455 disp(['total ellapsed time ' num2str(ellapsed_time/60,2) ' minutes']) 456 end 457 458 ellapsed_time=toc(tstart); 459 disp(['total ellapsed time ' num2str(ellapsed_time/60,2) ' minutes']) 459 disp(['ellapsed time since start ' num2str(ellapsed_time/60,2) ' minutes']) 460 end 461 460 462 disp([ num2str(ellapsed_time/(60*NbField),3) ' minutes per iteration']) 461 463
Note: See TracChangeset
for help on using the changeset viewer.