double precision qlonshellcutoff
      parameter(qlonshellcutoff=1d-10)