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

bugfix: nxl/nyl/nzl not correctly written

This commit is contained in:
Michael Stumpf (ifhcluster) 2018-12-10 10:41:04 +00:00
parent e05dcdc34b
commit fa769ec4b4
1 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
matlab/write_scal_complete_ucf.m
View File

@ -39,12 +39,16 @@ function [] = write_scal_complete_ucf(file,s,ibegp,iendp,jbegp,jendp,kbegp,kendp
fname = sprintf('%s.%05d',fbasepath,iproc); fname = sprintf('%s.%05d',fbasepath,iproc);
% Extract local indices % Extract local indices
ighost = 0;
ib = ibegp(ixproc+1); ib = ibegp(ixproc+1);
ie = iendp(ixproc+1); ie = iendp(ixproc+1);
jb = jbegp(iyproc+1); jb = jbegp(iyproc+1);
je = jendp(iyproc+1); je = jendp(iyproc+1);
kb = kbegp(izproc+1); kb = kbegp(izproc+1);
ke = kendp(izproc+1); ke = kendp(izproc+1);
nxl = ie-ib+1;
nyl = je-jb+1;
nzl = ke-kb+1;
% Create a new chunk file % Create a new chunk file
obj = ucf('verbosity',par.Results.verbosity,'debug',par.Results.debug); obj = ucf('verbosity',par.Results.verbosity,'debug',par.Results.debug);
@ -53,10 +57,6 @@ function [] = write_scal_complete_ucf(file,s,ibegp,iendp,jbegp,jendp,kbegp,kendp
% Write data field by field % Write data field by field
for iscal=1:nscal for iscal=1:nscal
ighost = 0;
nxl = size(s,1)-2*ighost;
nyl = size(s,2)-2*ighost;
nzl = size(s,3)-2*ighost;
params = int64([ighost,ib,jb,kb,nxl,nyl,nzl,nx,ny,nz]); params = int64([ighost,ib,jb,kb,nxl,nyl,nzl,nx,ny,nz]);
data = s(ib:ie,jb:je,kb:ke,iscal); data = s(ib:ie,jb:je,kb:ke,iscal);
obj.appendSet(data,params); obj.appendSet(data,params);