corner_pair_for_side Subroutine

  1. SELF_Mesh_2D_t.f90
  2. SELF_Mesh_2D_t
  3. corner_pair_for_side

public subroutine corner_pair_for_side(corners, iSide, cN1, cN2)

Returns the two corner-node IDs delimiting a given local side using SELF's 2D side convention: Side 1 South = [CN1, CN2] Side 2 East = [CN2, CN3] Side 3 North = [CN4, CN3] Side 4 West = [CN1, CN4]

Arguments

TypeIntentOptionalAttributesName
integer, intent(in) :: corners(1:4)
integer, intent(in) :: iSide
integer, intent(out) :: cN1
integer, intent(out) :: cN2

Called by

proc~~corner_pair_for_side~~CalledByGraph proc~corner_pair_for_side corner_pair_for_side proc~read_hohqmesh_mesh2d_t Read_HOHQMesh_Mesh2D_t proc~read_hohqmesh_mesh2d_t->proc~corner_pair_for_side

Contents

Source Code

corner_pair_for_side

Source Code

  subroutine corner_pair_for_side(corners,iSide,cN1,cN2)
    !! Returns the two corner-node IDs delimiting a given local side
    !! using SELF's 2D side convention:
    !!   Side 1 South = [CN1, CN2]
    !!   Side 2 East  = [CN2, CN3]
    !!   Side 3 North = [CN4, CN3]
    !!   Side 4 West  = [CN1, CN4]
    implicit none
    integer,intent(in) :: corners(1:4)
    integer,intent(in) :: iSide
    integer,intent(out) :: cN1,cN2
    select case(iSide)
    case(1); cN1 = corners(1); cN2 = corners(2)
    case(2); cN1 = corners(2); cN2 = corners(3)
    case(3); cN1 = corners(4); cN2 = corners(3)
    case(4); cN1 = corners(1); cN2 = corners(4)
    case default; cN1 = 0; cN2 = 0
    endselect
  endsubroutine corner_pair_for_side