real(kind=qp) qponshellcutoff
      parameter(qponshellcutoff=1.E-10_qp)