!cc
      indlsod(n,l) = jbasod + (n-l)/2 + 1
!cc