/*LimitSh sets upper limit to subshell integrations */
long LimitSh(long int ion, 
  long int nshell, 
  long int nelem);