[Nek5000-users] Boundary condition loop

[nek5000-users at lists.mcs.anl.gov]

Neks,

I need to loop over all element faces that are boundary faces (doesn't matter which type of boundary condition).  I know one way to do this would be

do e=1,nelv
  do f=1,2*ndim
    if (cbc(f,e,1).ne.'E  '.and.cbc(f,e,1).ne.'P  ') then
        ! Do computations
    endif
  enddo
enddo

Is there a more elegant solution for this, maybe some array that points to element/face combinations that are boundary conditions?  

Thanks,

Josh