From 949eea162400939b5701b71b3685c4ac114b5a7c Mon Sep 17 00:00:00 2001
From: Martin Diehl <mail@martin-diehl.net>
Date: Wed, 1 Sep 2021 11:07:39 +0200
Subject: [PATCH] bugfix: did not work with VTK 8

---
 PRIVATE                          | 2 +-
 src/grid/discretization_grid.f90 | 3 ++-
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/PRIVATE b/PRIVATE
index 8a6ab2d16..65c453c46 160000
--- a/PRIVATE
+++ b/PRIVATE
@@ -1 +1 @@
-Subproject commit 8a6ab2d16b01ca6fa4cd8cef62f68cd1b4f0cadd
+Subproject commit 65c453c46d157f6448ab48829d059b1d641fde47
diff --git a/src/grid/discretization_grid.f90 b/src/grid/discretization_grid.f90
index 1c88e5817..1758fd855 100644
--- a/src/grid/discretization_grid.f90
+++ b/src/grid/discretization_grid.f90
@@ -266,7 +266,8 @@ subroutine readVTI(grid,geomSize,origin,material, &
     integer :: i
 
 
-    if (getXMLValue(header,'Direction') /= '1 0 0 0 1 0 0 0 1') &
+    temp = getXMLValue(header,'Direction')
+    if (temp /= '1 0 0 0 1 0 0 0 1' .and. temp /= '') &                                             ! https://discourse.vtk.org/t/vti-specification/6526
       call IO_error(error_ID = 844, ext_msg = 'coordinate order')
 
     temp = getXMLValue(header,'WholeExtent')