fixed a bug in vtu velo plotting

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