read_velocity Subroutine

  1. read_output_vtk.f90
  2. read_output_vtk
  3. read_velocity

private subroutine read_velocity()

Read velocity vector from the vtk file

Arguments

None

Calls

proc~~read_velocity~~CallsGraph proc~read_velocity read_velocity debugcall debugcall proc~read_velocity->debugcall

Called by

proc~~read_velocity~~CalledByGraph proc~read_velocity read_velocity proc~read_file read_file proc~read_file->proc~read_velocity

Contents

Source Code

read_velocity

Source Code

    subroutine read_velocity()
      !< Read velocity vector from the vtk file
      implicit none

      DebugCall('read_output_vtk: read_velocity')
      read(IN_FILE_UNIT, *) !'CELL_DATA ', (imx-1)*(jmx-1)*(kmx-1)
      read(IN_FILE_UNIT, *) !'VECTORS Velocity FLOAT'
      do k = 1, kmx - 1
       do j = 1, jmx - 1
        do i = 1, imx - 1
          read(IN_FILE_UNIT, *) x_speed(i, j, k), y_speed(i, j, k), z_speed(i, j, k)
        end do
       end do
      end do
      read(IN_FILE_UNIT, *) 

    end subroutine read_velocity