Index: wflow-py/wflow/wf_netcdfio.py
===================================================================
diff -u -r5be5af755c609a0230420418bb558e9eeab0f6bc -ra31d19fc1aef843b047059c648ed259a801110f8
--- wflow-py/wflow/wf_netcdfio.py	(.../wf_netcdfio.py)	(revision 5be5af755c609a0230420418bb558e9eeab0f6bc)
+++ wflow-py/wflow/wf_netcdfio.py	(.../wf_netcdfio.py)	(revision a31d19fc1aef843b047059c648ed259a801110f8)
@@ -530,7 +530,7 @@
         for var in vars:
             try:
                 #self.alldat[var] = self.dataset.variables[var][self.fstep:self.maxsteps]
-                self.alldat[var] = reshape(self.dataset.variables[var],(self.maxlentime,len(y),len(x)))
+                self.alldat[var] = self.dataset.variables[var]
             except:
                 self.alldat.pop(var, None)
                 logging.warn("Variable " + var + " not found in netcdf file: " + netcdffile)
@@ -578,7 +578,10 @@
             if ncindex == self.lstep:  # Read new block of data in mem
                 logging.debug("reading new netcdf data block starting at: " + str(ncindex))
                 for vars in self.alldat:
-                    self.alldat[vars] = reshape(self.dataset.variables[vars],(self.maxlentime,len(self.y),len(self.x)))[ncindex:ncindex + self.maxsteps]
+                    if len(vars) == 3:
+                        self.alldat[vars] = self.dataset.variables[vars][ncindex:ncindex + self.maxsteps]
+                    elif len(vars) == 4:
+                        self.alldat[vars] = self.dataset.variables[vars][ncindex:ncindex + self.maxsteps,0]
 
                 self.fstep = ncindex
                 self.lstep = ncindex + self.maxsteps