write_vtu_unstructured Subroutine

   subroutine write_vtu_unstructured(params, flow, mesh, fields, species, energy, transport, step, time)
      type(case_params_t), intent(in) :: params
      type(flow_mpi_t), intent(inout) :: flow
      type(mesh_t), intent(in) :: mesh
      type(flow_fields_t), intent(in) :: fields
      type(species_fields_t), intent(in) :: species
      type(energy_fields_t), intent(in) :: energy
      type(transport_properties_t), intent(in) :: transport
      integer, intent(in) :: step
      real(rk), intent(in) :: time

      if (.not. params%write_vtu) return

      if (params%vtu_format == "binary") then
         call write_vtu_binary(params, flow, mesh, fields, species, energy, transport, step, time)
      else
         call write_vtu_ascii(params, flow, mesh, fields, species, energy, transport, step, time)
      end if

      ! Write global diagnostics and parallel XML master files
      call write_species_energy_conservation_diagnostics(mesh, flow, params, fields, species, energy, transport, &
                                                   step, time)
      call write_enthalpy_energy_budget_diagnostics(mesh, flow, params, fields, energy, transport, &
                                                   step, time)

      if (flow%rank == 0) then
         call write_pvtu_master(params, flow, step)
         call update_pvd_collection(params, flow, step, time)
      end if
   end subroutine write_vtu_unstructured