diff --git a/flash2d/src/options-sphere/IO_vtu_stacked.F90 b/flash2d/src/options-sphere/IO_vtu_stacked.F90
index 264780769828f6714db8f27917254d52588f5e62..cb83bae98ec4aa005b8def94771aacbe9f2ae088 100644
--- a/flash2d/src/options-sphere/IO_vtu_stacked.F90
+++ b/flash2d/src/options-sphere/IO_vtu_stacked.F90
@@ -127,7 +127,7 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
   INTEGER( KIND = GRID_SI)                       :: i_numberofvtupoints
   INTEGER( KIND = GRID_SI)                       :: i_numberofvtufaces
   INTEGER( KIND = GRID_SI)                       :: i_numberofvtucells
-  REAL ( KIND = GRID_SR )                        :: r_heightfactor
+  REAL ( KIND = GRID_SR )                        :: r_heightfactor, r_aux
   TYPE(t_vtu_data), DIMENSION(6)                 :: p_stcnodedata
   REAL (KIND = GRID_SR), DIMENSION(:,:), POINTER :: r_stcval
   INTEGER( KIND = GRID_SI)                       :: i_stcbeg
@@ -394,7 +394,7 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
 	END IF
     DO i_cnt = 1, i_nodedata
       IF (p_nodedata(i_cnt)%c_name(1:5) == 'layer') THEN
-        READ(i_hcnt,*) p_nodedata(i_cnt)%c_name(6:32)
+        READ(i_hcnt,'(I8)') p_nodedata(i_cnt)%c_name(6:32)
         i_stcbeg = p_mesh%i_nnumber*(i_hcnt-1)+1
         i_stcend = p_mesh%i_nnumber*i_hcnt
         r_stcval(1:1,i_stcbeg:i_stcend) = p_nodedata(i_cnt)%p_vdata(1:1,:)
@@ -443,7 +443,8 @@ SUBROUTINE plot_vtu(p_mesh, c_filename, &
       DO i_hcnt=1,i_numberofcells
         i_stcbeg = (i_hcnt-1)*(i_lay-1)+1
         i_stcend = i_hcnt*(i_lay-1)
-        r_stcval(1:1,i_stcbeg:i_stcend) = p_celldata(i_cnt)%p_vdata(1:1,i_hcnt)
+        r_aux    = p_celldata(i_cnt)%p_vdata(1,i_hcnt)
+        r_stcval(1,i_stcbeg:i_stcend) = r_aux
       END DO
       p_stcnodedata(1)%c_name = p_celldata(i_cnt)%c_name
       p_stcnodedata(1)%i_size = p_celldata(i_cnt)%i_size