From 58545ff0900753e2066d96978a92b60d1035a923 Mon Sep 17 00:00:00 2001
From: Joern Behrens <joern.behrens@uni-hamburg.de>
Date: Wed, 16 Feb 2022 19:30:09 +0100
Subject: [PATCH] fixed a bug in vtu velo plotting
---
.../src/options-sphere/IO_vtuplot_stacked.F90 | 20 +++++++++----------
1 file changed, 10 insertions(+), 10 deletions(-)
diff --git a/flash2d/src/options-sphere/IO_vtuplot_stacked.F90 b/flash2d/src/options-sphere/IO_vtuplot_stacked.F90
index 52dec7a..0f0cb1c 100644
--- a/flash2d/src/options-sphere/IO_vtuplot_stacked.F90
+++ b/flash2d/src/options-sphere/IO_vtuplot_stacked.F90
@@ -151,19 +151,19 @@ MODULE IO_vtuplot
!---------- encode velocity fields for layers
DO i_lay=1,i_stcnumlayers
- i_vcnt= (i_lay-1)*GRID_dimension+1+i_nodnum
+ i_vcnt= (i_lay-1)*GRID_dimension+1
write(c_tmp,*) 'u-wnd',i_lay
- nodedata(i_lay+5)%c_name = c_tmp
- nodedata(i_lay+5)%i_size = 1
- nodedata(i_lay+5)%p_vdata => r_velo(i_vcnt:i_vcnt,:)
+ nodedata(i_vcnt+i_nodnum)%c_name = c_tmp
+ nodedata(i_vcnt+i_nodnum)%i_size = 1
+ nodedata(i_vcnt+i_nodnum)%p_vdata => r_velo(i_vcnt:i_vcnt,:)
write(c_tmp,*) 'v-wnd',i_lay
- nodedata(i_lay+5)%c_name = c_tmp
- nodedata(i_lay+5)%i_size = 1
- nodedata(i_lay+5)%p_vdata => r_velo(i_vcnt+1:i_vcnt+1,:)
+ nodedata(i_vcnt+i_nodnum+1)%c_name = c_tmp
+ nodedata(i_vcnt+i_nodnum+1)%i_size = 1
+ nodedata(i_vcnt+i_nodnum+1)%p_vdata => r_velo(i_vcnt+1:i_vcnt+1,:)
write(c_tmp,*) 'w-wnd',i_lay
- nodedata(i_lay+5)%c_name = c_tmp
- nodedata(i_lay+5)%i_size = 1
- nodedata(i_lay+5)%p_vdata => r_velo(i_vcnt+2:i_vcnt+2,:)
+ nodedata(i_vcnt+i_nodnum+2)%c_name = c_tmp
+ nodedata(i_vcnt+i_nodnum+2)%i_size = 1
+ nodedata(i_vcnt+i_nodnum+2)%p_vdata => r_velo(i_vcnt+2:i_vcnt+2,:)
END DO
i_nodnum= i_nodnum+ GRID_dimension*i_stcnumlayers
--
GitLab