mwtkrayer
/
ucftools
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity

renamed

This commit is contained in:
Michael Stumpf (ifhcluster) 2019-03-26 12:25:50 +00:00
parent d8a650884a
commit 8064b589d2
1 changed files with 5 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/read_field_complete_ucftar.m → matlab/read_field_complete_hucf.m
View File

@ -1,8 +1,8 @@
function [q,x,y,z] = read_field_complete_ucftar(htar,field,varargin) function [q,x,y,z] = read_field_complete_hucf(hucf,field,varargin)
% [q,x,y,z] = read_field_complete_ucftar(file,field,varargin) % [q,x,y,z] = read_field_complete_ucftar(file,field,varargin)
% Reads a specific field from all processor chunks.(UCF tar version) % Reads a specific field from all processor chunks.(UCF tar version)
% Input % Input
% htar handle of tar file (ustar object) % hucf handle of UCF object (ustar,ucfmulti)
% field field to be read % field field to be read
% {'u','v','w','p','s1','s2',...} % {'u','v','w','p','s1','s2',...}
% ? step index of step to be read (default: 1) % ? step index of step to be read (default: 1)
@ -31,7 +31,7 @@ function [q,x,y,z] = read_field_complete_ucftar(htar,field,varargin)
end end
% Read parameters % Read parameters
params = read_parameters_ucf(htar,'tarmode',1); params = read_parameters_ucf(hucf);
% Construct an array with final mesh size % Construct an array with final mesh size
nx = params.mesh.(sprintf('nx%s',cmesh)); nx = params.mesh.(sprintf('nx%s',cmesh));
@ -40,7 +40,7 @@ function [q,x,y,z] = read_field_complete_ucftar(htar,field,varargin)
q = zeros(nx,ny,nz); q = zeros(nx,ny,nz);
% Read grid % Read grid
[x.u,y.u,z.u,x.v,y.v,z.v,x.w,y.w,z.w,x.p,y.p,z.p] = read_grid_ucf(htar,'tarmode',1); [x.u,y.u,z.u,x.v,y.v,z.v,x.w,y.w,z.w,x.p,y.p,z.p] = read_grid_ucf(hucf);
x = x.(cmesh); x = x.(cmesh);
y = y.(cmesh); y = y.(cmesh);
z = z.(cmesh); z = z.(cmesh);
@ -50,7 +50,7 @@ function [q,x,y,z] = read_field_complete_ucftar(htar,field,varargin)
% Now read chunk by chunk % Now read chunk by chunk
for iproc=0:nprocs-1 for iproc=0:nprocs-1
[data,ib,jb,kb,nxl,nyl,nzl] = read_field_chunk_ucf(htar,field,'tarmode',1,'rank',iproc,'ghost',0,'step',istep); [data,ib,jb,kb,nxl,nyl,nzl] = read_field_chunk_ucf(hucf,field,'rank',iproc,'ghost',0,'step',istep);
q(ib:ib+nxl-1,jb:jb+nyl-1,kb:kb+nzl-1) = data; q(ib:ib+nxl-1,jb:jb+nyl-1,kb:kb+nzl-1) = data;
end end
end end