src/coordinates/coordinates.cpp
//----------------------------------------------------------------------------------------
// Function for determining if index corresponds to a polar boundary
// Inputs:
//   j: x2-index
// Outputs:
//   returned value: true if face indexed with j is on a pole; false otherwise
Function for determining if index corresponds to a polar boundary
Inputs:
j: x2-index
Outputs:
returned value: true if face indexed with j is on a pole; false otherwise
Function for determining if index corresponds to a polar boundary Inputs : j : x2 - index Outputs : returned value : true if face indexed with j is on a pole ; false otherwise
bool Coordinates::IsPole(int j)
{
  if ((pmy_block->block_bcs[INNER_X2] == POLAR_BNDRY
      or pmy_block->block_bcs[INNER_X2] == POLAR_BNDRY_WEDGE) and j == pmy_block->js) {
    return true;
  }
  if ((pmy_block->block_bcs[OUTER_X2] == POLAR_BNDRY
      or pmy_block->block_bcs[OUTER_X2] == POLAR_BNDRY_WEDGE) and j == pmy_block->je+1) {
    return true;
  }
  return false;
}
bool Coordinates :: IsPole ( int j ) { if ( ( pmy_block -> block_bcs [ INNER_X2 ] == POLAR_BNDRY or pmy_block -> block_bcs [ INNER_X2 ] == POLAR_BNDRY_WEDGE ) and j == pmy_block -> js ) { return true ; } if ( ( pmy_block -> block_bcs [ OUTER_X2 ] == POLAR_BNDRY or pmy_block -> block_bcs [ OUTER_X2 ] == POLAR_BNDRY_WEDGE ) and j == pmy_block -> je + 1 ) { return true ; } return false ; }
Function for determining if index corresponds to a polar boundary
Inputs:
j: x2-index
Outputs:
returned value: true if face indexed with j is on a pole; false otherwise
Function for determining if index corresponds to a polar boundary Inputs : j : x2 - index Outputs : returned value : true if face indexed with j is on a pole ; false otherwise