Purpose: check the sampling points
         of a mccpd/mcpotfit run

  Usage:  chksampling [-h|-?]  [-ver -d -e -oi -oe
                                -rem -C -hist] 

   vers:  Program version (e.g. 86 for version 8.6)
    d/D:  Indicates a debug version

 optn: -h          : Print this help-text.
       -?          : Print this help-text.
       -ver        : Version information.

 -d       : read dvr from  rather
                     than from ./dvr
 -e  : File containing energies matching the
                     . Required for -rem and -rem2.
 -oi     : write cleaned index to 
                     rather than to _chk
 -oe  : write cleaned energies to 
                     rather than to _chk

 -rem     : remove points if at least 
                     DOF are at most  points away from 
                     the grid boundary (i=0: first or last
                     point ) AND the energy is less then .

 -rem2    : remove points if at least 
                     DOF are at most  points away from 
                     the grid boundary (i=0: first or last
                     point) AND the energy is larger then .

 -C             : If -rem or rem2 is set, also remove  lines
                     of context before and after a removed point.

 -hist             : Write histogram files of the distribution
                     of sampling points. If -rem or -rem2 is set,
                     histograms will be made before and after the
                     trajectory was cleaned. The histgram files
                     are named:
                     _hist_     (before -rem)
                     _hist__chk (after -rem) 

 -hist2            : Write 2D histogram files of the distribution
                     of sampling points. If -rem or -rem2 is set,
                     histograms will be made before and after the
                     trajectory was cleaned. The histgram files
                     are named:
                     _hist2d__     (before -rem)
                     _hist2d___chk (after -rem) 

 -cov           : Write the normalized covariance matrix of the
                     sampling points to file.  is the exponent
                     of the variables, i.e. cov(q_i^n, q_j^n) will
                     be calculated. If -rem or rem2 is is set
                     the matrix is calculated before and after
                     the trajectory was cleaned. The files are named:
                     _covariance     (before -rem)
                     _covariance_chk (after -rem).
                     Also creates a sorted list of covariances.
                     and a file with means and variances