Writing statistical averages for postprocessing in .vtu files.

nikola-m · Aug 17, 2023 · 8dc9dda · 8dc9dda
1 parent 65c46a8
commit 8dc9dda
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/src/io/writefiles.f90 b/src/io/writefiles.f90
@@ -127,6 +127,19 @@ subroutine writefiles
   ! call vtu_write_XML_scalar_field_node ( output_unit, 'mueff', uN )
   ! !#--
 
+  !#--
+  ! Write ensemble averaged velocity components values in nodes.
+  ! Interpolate to nodes using interpolation weights generated at init stage.
+  if(ltransient .and. lturb) then
+
+  call interpolate_to_nodes(u_aver, uN)
+  call interpolate_to_nodes(v_aver, vN)
+  call interpolate_to_nodes(w_aver, wN)
+  call vtu_write_XML_vector_field_node( output_unit, 'Uavg', uN, vN, wN )
+
+  endif
+  ! !#--
+
   deallocate( uN,vN, wN )
 
   do i=1,2! Lines where we close PointData tag and open CellData