- focus on 3 channels!
- better simplifications!
- using occupation number representation
- in recalcop[]: m should contain conserved quantum numbers (q,s,p, etc),
  while the lists should contain those additional quantum numbers that
  drop out of the problem
  
- argument 'r' is totally useless (dirpdt)
