diff --git a/flash2d/compile/Makefiles/Makefile.common b/flash2d/compile/Makefiles/Makefile.common
index 49b85db9c8d1580cc068323ecb240e5d1367ccd3..eb4f9a65101beca7c985e6726f32b08e65d6ab71 100644
--- a/flash2d/compile/Makefiles/Makefile.common
+++ b/flash2d/compile/Makefiles/Makefile.common
@@ -243,14 +243,16 @@ ShallowWater:
@cp $(OPTDIR)/SLM_initial.SWCosine.f90 SLM_initial.f90
@cp $(OPTDIR)/IO_vtuplot.SW.F90 IO_vtuplot.F90
-Stacked:
+StackedSpherical:
@make maincopy
@cp $(OPTDIR)/FLASH_parameters_stacked.f90 FLASH_parameters.f90
@cp $(OPTDIR)/ADV_rhs_stacked.f90 ADV_rhs.f90
@cp $(OPTDIR)/ADV_semilagrange_stacked.F90 ADV_semilagrange.F90
@cp $(OPTDIR)/ADV_wind_stacked.f90 ADV_wind.f90
+ @cp $(OPTDIR)/IO_vtuplot_stacked.F90 IO_vtuplot.F90
@cp $(OPTDIR)/SLM_advanced_stacked.f90 SLM_advanced.f90
@cp $(OPTDIR)/SLM_errorestimate_stacked.f90 SLM_errorestimate.f90
+ @cp $(OPTDIR)/SLM_initial_stacked.f90 SLM_initial.f90
@cp $(OPTDIR)/SLM_simple_stacked.f90 SLM_simple.f90
#----------------------------------------------------------------#
diff --git a/flash2d/src/options-sphere/IO_vtuplot_stacked.F90 b/flash2d/src/options-sphere/IO_vtuplot_stacked.F90
index aebd0d4b13537459da479355ab5f19b8f479e541..c1ea6ccf0a42c76944fa3577d2d7bbda794973cc 100644
--- a/flash2d/src/options-sphere/IO_vtuplot_stacked.F90
+++ b/flash2d/src/options-sphere/IO_vtuplot_stacked.F90
@@ -139,7 +139,7 @@ MODULE IO_vtuplot
nodedata(i_lay+5)%i_size = 1
nodedata(i_lay+5)%p_vdata => r_val(i_lay+5:i_lay+5,:)
END DO
- i_nodnum= i_stcnumlayers+ 4
+ i_nodnum= i_stcnumlayers+ 5
!---------- extract cells grid data
@@ -170,7 +170,8 @@ MODULE IO_vtuplot
! plot the vtu data
CALL plot_vtu(p_handle, c_mfile, i_nodedata = i_nodnum, p_nodedata = nodedata, &
- i_celldata = 2, p_celldata = celldata) !, i_zcoordinate = 4)
+ i_celldata = 2, p_celldata = celldata, i_layers = i_stcnumlayers, &
+ r_height = r_stclayerpressure) !, i_zcoordinate = 4)
DEALLOCATE(i_ids, i_sta, r_sta, r_lvl, &
r_val, r_velo)