Merge pull request #649 from mkondratyev85/mkondratyev85_issue641

Fix export to vtk
gempy-project · Oct 29, 2021 · 2378578 · 2378578
2 parents 5958fed + ac0c2e0
commit 2378578
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/gempy/plot/visualization_3d.py b/gempy/plot/visualization_3d.py
@@ -1177,7 +1177,7 @@ def export_vtk_lith_block(geo_data, lith_block=None, path="./default"):
 
         if lith_block is None:
             lith_block = geo_data.solutions.lith_block
-        lith = lith_block.reshape((nx, ny, nz))
+        lith = lith_block.reshape((nx, ny, nz)).flatten(order='F')
         out = pv.RectilinearGrid(x,y,z)
         out['Lithology'] = lith
         out.save(path+'_lith_block.vtr')