Changeset 937 for trunk/src/series
- Timestamp:
- Mar 30, 2016, 11:04:43 PM (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/series/merge_proj.m
r924 r937 381 381 end 382 382 383 else 384 MergeData.ListGlobalAttribute={'Conventions','Project','InputFile_1','InputFile_end',' nb_coord','nb_dim'};383 else %output as netcdf files 384 MergeData.ListGlobalAttribute={'Conventions','Project','InputFile_1','InputFile_end','NbCoord','NbDim'}; 385 385 MergeData.Conventions='uvmat'; 386 MergeData.nb_coord=2; 387 MergeData.nb_dim=2; 388 dt=[]; 389 if isfield(Data{1},'dt')&& isnumeric(Data{1}.dt) 390 dt=Data{1}.dt; 386 MergeData.NbCoord=2; 387 MergeData.NbDim=2; 388 % time interval of PIV 389 Dt=[]; 390 if isfield(Data{1},'Dt')&& isnumeric(Data{1}.Dt) 391 Dt=Data{1}.Dt; 391 392 end 392 393 for iview =2:numel(Data) 393 if ~(isfield(Data{iview},' dt')&& isequal(Data{iview}.dt,dt))394 dt=[];%dt not the same for all fields394 if ~(isfield(Data{iview},'Dt')&& isequal(Data{iview}.Dt,Dt)) 395 Dt=[];%dt not the same for all fields 395 396 end 396 397 end … … 399 400 MergeData.Time=timeread; 400 401 end 401 if ~isempty(dt) 402 MergeData.ListGlobalAttribute=[MergeData.ListGlobalAttribute {'dt'}]; 403 MergeData.dt=dt; 402 % position of projection plane 403 if isfield(Data{1},'ProjObjectCoord')&& isfield(Data{1},'ProjObjectAngle') 404 'test' 405 ProjObjectCoord=Data{1}.ProjObjectCoord; 406 ProjObjectAngle=Data{1}.ProjObjectAngle; 407 for iview =2:numel(Data) 408 if ~(isfield(Data{iview},'ProjObjectCoord')&& isequal(Data{iview}.ProjObjectCoord,ProjObjectCoord))... 409 ||~(isfield(Data{iview},'ProjObjectAngle')&& isequal(Data{iview}.ProjObjectAngle,ProjObjectAngle)) 410 ProjObjectCoord=[];%dt not the same for all fields 411 end 412 end 413 if ~isempty(ProjObjectCoord) 414 MergeData.ListGlobalAttribute=[MergeData.ListGlobalAttribute {'ProjObjectCoord'} {'ProjObjectAngle'}]; 415 MergeData.ProjObjectCoord=ProjObjectCoord; 416 MergeData.ProjObjectAngle=ProjObjectAngle; 417 end 418 end 419 % coord unit 420 if isfield(Data{1},'CoordUnit') 421 CoordUnit=Data{1}.CoordUnit; 422 for iview =2:numel(Data) 423 if ~(isfield(Data{iview},'CoordUnit')&& isequal(Data{iview}.CoordUnit,CoordUnit)) 424 CoordUnit=[];%CoordUnit not the same for all fields 425 end 426 end 427 if ~isempty(CoordUnit) 428 MergeData.ListGlobalAttribute=[MergeData.ListGlobalAttribute {'CoordUnit'}]; 429 MergeData.CoordUnit=CoordUnit; 430 end 431 end 432 % time unit 433 if isfield(Data{1},'TimeUnit') 434 TimeUnit=Data{1}.TimeUnit; 435 for iview =2:numel(Data) 436 if ~(isfield(Data{iview},'TimeUnit')&& isequal(Data{iview}.TimeUnit,TimeUnit)) 437 TimeUnit=[];%TimeUnit not the same for all fields 438 end 439 end 440 if ~isempty(TimeUnit) 441 MergeData.ListGlobalAttribute=[MergeData.ListGlobalAttribute {'TimeUnit'}]; 442 MergeData.TimeUnit=TimeUnit; 443 end 404 444 end 405 445 error=struct2nc(OutputFile,MergeData);%save result file
Note: See TracChangeset
for help on using the changeset viewer.