extract_grid_point Subroutine

  1. grid.f90
  2. grid
  3. extract_grid_point

private subroutine extract_grid_point(line, i, j, k)

Extract a grid point from a line of the grid file.

Arguments

Type IntentOptional AttributesName
character(len=STRING_BUFFER_LENGTH), intent(in) :: line
integer, intent(in) :: i
integer, intent(in) :: j
integer, intent(in) :: k

Calls

proc~~extract_grid_point~~CallsGraph proc~extract_grid_point extract_grid_point proc~dmsg dmsg proc~extract_grid_point->proc~dmsg

Called by

proc~~extract_grid_point~~CalledByGraph proc~extract_grid_point extract_grid_point proc~populate_grid_points populate_grid_points proc~populate_grid_points->proc~extract_grid_point proc~setup_grid setup_grid proc~setup_grid->proc~populate_grid_points proc~setup_solver setup_solver proc~setup_solver->proc~setup_grid proc~setup_all setup_all proc~setup_all->proc~setup_grid proc~start_run start_run proc~start_run->proc~setup_solver program~main main program~main->proc~start_run

Contents

Source Code

extract_grid_point

Source Code

        subroutine extract_grid_point(line, i, j, k)
            !< Extract a grid point from a line of the grid file. 
            !-----------------------------------------------------------

            implicit none
            character(len=STRING_BUFFER_LENGTH), intent(in) :: line
            integer, intent(in) :: i, j, k

            call dmsg(0, 'grid', 'extract_grid_point')

            if (kmx > 1) then
                read(line, *) grid_x(i, j, k), grid_y(i, j, k), grid_z(i, j, k)
            else    
                if (jmx > 1) then
                    read(line, *) grid_x(i, j, k), grid_y(i, j, k)
                else
                    read(line, *) grid_x(i, j, k)
                    grid_y(i, j, k) = 0.
                end if
                grid_z(i, j, k) = 0.
            end if
        end subroutine extract_grid_point