From 9637edbee5a4f41df7449b3524dac29723c5305d Mon Sep 17 00:00:00 2001
From: Joern Behrens <joern.behrens@uni-hamburg.de>
Date: Thu, 27 Jan 2022 08:26:54 +0100
Subject: [PATCH] working on plot functions
---
flash2d/compile/Makefiles/Makefile.common | 4 +++-
flash2d/src/options-sphere/IO_vtuplot_stacked.F90 | 5 +++--
2 files changed, 6 insertions(+), 3 deletions(-)
diff --git a/flash2d/compile/Makefiles/Makefile.common b/flash2d/compile/Makefiles/Makefile.common
index 49b85db..eb4f9a6 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 aebd0d4..c1ea6cc 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)
--
GitLab