pole Subroutine

private subroutine pole(face, bc, dims)

Boundary condition for the block face with zero area; turning into a pole

Arguments

Type IntentOptional AttributesName
character(len=*), intent(in) :: face

Name of the face at which boundary condition is called

type(boundarytype), intent(in) :: bc
type(extent), intent(in) :: dims

Calls

proc~~pole~~CallsGraph proc~pole pole proc~copy3 copy3 proc~pole->proc~copy3

Called by

proc~~pole~~CalledByGraph proc~pole pole proc~populate_ghost_primitive populate_ghost_primitive proc~populate_ghost_primitive->proc~pole proc~get_total_conservative_residue get_total_conservative_Residue proc~get_total_conservative_residue->proc~populate_ghost_primitive proc~get_next_solution get_next_solution proc~get_next_solution->proc~get_total_conservative_residue proc~iterate_one_more_time_step iterate_one_more_time_step proc~iterate_one_more_time_step->proc~get_next_solution program~main main program~main->proc~iterate_one_more_time_step

Contents

Source Code


Source Code

      subroutine pole(face, bc, dims)
        !< Boundary condition for the block face
        !< with zero area; turning into a pole
        implicit none
        type(extent), intent(in) :: dims
        type(boundarytype), intent(in) :: bc
        character(len=*), intent(in) :: face
        !< Name of the face at which boundary condition is called
        call copy3(density , "flat", face, bc, dims)
        call copy3(x_speed , "flat", face, bc, dims)
        call copy3(y_speed , "flat", face, bc, dims)
        call copy3(z_speed , "flat", face, bc, dims)
        call copy3(pressure, "flat", face, bc, dims)
        select case (turbulence)
          case('none')
            !do nothing
            continue
          case('sa', 'saBC')
            call copy3(tv, "flat", face, bc, dims) 
          case('sst', 'sst2003')
            call copy3(tk, "flat", face, bc, dims)
            call copy3(tw, "flat", face, bc, dims)
          case('kkl')
            call copy3(tk, "flat", face, bc, dims)
            call copy3(tkl, "flat", face, bc, dims)
          case DEFAULT
            Fatal_error
        end select
        select case(trim(transition))
          case('lctm2015')
            call copy3(tgm, "flat", face, bc, dims)
          case DEFAULT
            continue
        end select
      end subroutine pole