write_energy_reconciliation_scalar Subroutine

   subroutine write_energy_reconciliation_scalar(unit_id, mesh, flow, energy, transport, name)
      integer, intent(in) :: unit_id
      type(mesh_t), intent(in) :: mesh
      type(flow_mpi_t), intent(in) :: flow
      type(energy_fields_t), intent(in) :: energy
      type(transport_properties_t), intent(in) :: transport
      character(len=*), intent(in) :: name

      integer :: c
      real(rk) :: value

      write(unit_id,'(a)') '        <DataArray type="Float64" Name="' // trim(name) // '" format="ascii">'
      do c = 1, mesh%ncells
         if (.not. flow%owned(c)) cycle
         value = energy_reconciliation_value(energy, transport, c, trim(name))
         write(unit_id,'(ES26.16E4)') value
      end do
      write(unit_id,'(a)') '        </DataArray>'

   end subroutine write_energy_reconciliation_scalar