MODULE Schem_MODULE
!  scheme   1 for molecular problems.
!           2 for coulombic problems.
   INTEGER :: scheme=1
END MODULE Schem_MODULE