From fa769ec4b4c55f40234ec2d2702bf7e06e57ba30 Mon Sep 17 00:00:00 2001
From: "Michael Stumpf (ifhcluster)" <michael.stumpf@student.kit.edu>
Date: Mon, 10 Dec 2018 10:41:04 +0000
Subject: [PATCH] bugfix: nxl/nyl/nzl not correctly written

---
 matlab/write_scal_complete_ucf.m | 8 ++++----
 1 file changed, 4 insertions(+), 4 deletions(-)

diff --git a/matlab/write_scal_complete_ucf.m b/matlab/write_scal_complete_ucf.m
index 90fdf41..713e683 100644
--- a/matlab/write_scal_complete_ucf.m
+++ b/matlab/write_scal_complete_ucf.m
@@ -39,12 +39,16 @@ function [] = write_scal_complete_ucf(file,s,ibegp,iendp,jbegp,jendp,kbegp,kendp
           fname = sprintf('%s.%05d',fbasepath,iproc);
 
           % Extract local indices
+          ighost = 0;
           ib = ibegp(ixproc+1);
           ie = iendp(ixproc+1);
           jb = jbegp(iyproc+1);
           je = jendp(iyproc+1);
           kb = kbegp(izproc+1);
           ke = kendp(izproc+1);
+          nxl = ie-ib+1;
+          nyl = je-jb+1;
+          nzl = ke-kb+1;
 
           % Create a new chunk file
           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
           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]);
             data = s(ib:ie,jb:je,kb:ke,iscal);
             obj.appendSet(data,params);