diff --git a/flash2d/src/options-sphere/IO_vtu_stacked.F90 b/flash2d/src/options-sphere/IO_vtu_stacked.F90
index 7ec16c980f220b1c756b72b1515a7832f48e67d9..30612ad55fc489ab520652d2dcdaff0bff508d82 100644
--- a/flash2d/src/options-sphere/IO_vtu_stacked.F90
+++ b/flash2d/src/options-sphere/IO_vtu_stacked.F90
@@ -26,20 +26,20 @@ PUBLIC :: plot_vtu, t_vtu_data
   INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodesperface = GRID_edgenodes
 
   ! This defines a triangle in VTK format
-  INTEGER(KIND=GRID_SI), PARAMETER         :: i_vtucelltype = 5
+  INTEGER(KIND=GRID_SI), PARAMETER         :: i_vtucelltype  = 5
 #elif defined VTU_OUTPUT3D
   INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodespercell = GRID_tetranodes
   INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodesperface = GRID_elementnodes
 
   ! This defines a tetrahedron in VTK format
-  INTEGER(KIND=GRID_SI), PARAMETER         :: i_vtucelltype = 10
+  INTEGER(KIND=GRID_SI), PARAMETER         :: i_vtucelltype  = 10
 #elif defined VTU_OUTPUTSTC
-  INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodespercell = GRID_elementnodes
-  INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodesperface = GRID_edgenodes
+  INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodespercell    = GRID_elementnodes
+  INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodesperface    = GRID_edgenodes
 
   ! This defines a wedge/prism in VTK format
   INTEGER(KIND=GRID_SI), PARAMETER         :: i_nodespervtucell = GRID_elementnodes*2
-  INTEGER(KIND=GRID_SI), PARAMETER         :: i_vtucelltype = 13
+  INTEGER(KIND=GRID_SI), PARAMETER         :: i_vtucelltype     = 13
 #endif
 
 
@@ -355,14 +355,26 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
 
   ! write the cell type. Tetrahedra are represented by 10 and triangles by 5
   WRITE(i_fhandle, *) '<DataArray type="UInt8" Name="types" format="ascii">'
+#if defined VTU_OUTPUT2D || defined VTU_OUTPUTSPH || defined VTU_OUTPUT2D
   DO i_cnt = 1, i_numberofcells
     WRITE(i_fhandle, *) i_vtucelltype
   END DO
+#elif defined VTU_OUTPUTSTC
+  DO i_cnt = 1, i_numberofvtucells
+    WRITE(i_fhandle, *) i_vtucelltype
+  END DO
+#endif
   WRITE(i_fhandle, *) '</DataArray>'
   WRITE(i_fhandle, *) '<DataArray type="Int32" Name="offsets" format="ascii">'
+#if defined VTU_OUTPUT2D || defined VTU_OUTPUTSPH || defined VTU_OUTPUT2D
   DO i_cnt = 1, i_numberofcells
     WRITE(i_fhandle, *) i_cnt * i_nodespercell
   END DO
+#elif defined VTU_OUTPUTSTC
+  DO i_cnt = 1, i_numberofvtucells
+    WRITE(i_fhandle, *) i_cnt * i_nodespervtucell
+  END DO
+#endif
   WRITE(i_fhandle, *) '</DataArray>'
   WRITE(i_fhandle, *) '</Cells>'