Merge pull request #246 from kostahoraites/read_fsgrid_variable

read_interpolated_fsgrid_variable() bug fix

pyVlsv/vlsvreader.py

@@ -1143,6 +1143,8 @@ def read_interpolated_fsgrid_variable(self, name, coordinates, operator="pass",p
       .. seealso:: :func:`read` :func:`read_variable_info`
       '''
 
+      warnings.warn("read_interpolated_fsgrid_variable: face- vs. edge- centered variables not accounted for!")
+
       if name[0:3] != 'fg_':
          raise ValueError("Interpolation of FsGrid called on non-FsGrid data; exiting.")
 
@@ -1152,7 +1154,6 @@ def read_interpolated_fsgrid_variable(self, name, coordinates, operator="pass",p
       #First off let's fetch the data and some meta
       fg_data=self.read_fsgrid_variable( name,operator=operator)
       fg_size=self.get_fsgrid_mesh_size()
-      fg_data=np.reshape(fg_data,fg_size)
       nx,ny,nz=fg_size
       extents=self.get_fsgrid_mesh_extent()
       xmin,ymin,zmin,xmax,ymax,zmax=extents