20       subroutine orbital  subroutine
30       integer i,j,k,m,ii  local
31       integer mi,mj,mk
32       integer iorb,jorb,korb
33       integer nlist,next
34       integer, allocatable :: list(:)
35       character*20 keyword
36       character*240 record
37       character*240 string
83    10       continue
247 c     ##  subroutine piplane  --  plane perpendicular to orbital  ##  argument