C ... Control vector
      parameter( MXCNTL=100 )
      common /cntlvec/ icntl(MXCNTL), rcntl(MXCNTL)
      parameter( MAXMEM=100 000 )
      common /reala  / ra(MAXMEM)
      common /inta   / ia(MAXMEM)
      common /ptrs   / ip1,ip2,ip3,ip4,ip5,ip6,locrnd,locind