working on plot functions

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
 
 #----------------------------------------------------------------#

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)