Purpose: Calculates the diabatic and adiabatic population.
          And also one or two dimensional adiabatic densities.
 Usage:   adpop [-i -f -d -p -o -r -n -skip -step -tcut -w 
            -q -mc -mcd -mcs -mcsob -noq -qtol -mctol -nogwp 
            -loose -tight -saddle -saddle0 -saddle1 
            -cdvr -tddvr -DB -dblocal -batch -ver -h -?] [Degrees of Freedom]

Options  : 
 -i DIR  : data is stored in directory DIR
 -f FILE : the PSI is read from FILE rather than from ./psi
           The string FILE may be a relative or a full path-name.
 -d FILE : the DVR is read from FILE rather than from ./dvr
           The string FILE may be a relative or a full path-name.
 -p FILE : the PES is read from FILE rather than from ./pes
           The string FILE may be a relative or a full path-name.
 -o FILE : The ouput is written to file FILE rather than to ./adp
           The string FILE may be a relative or a full path-name.
           If FILE=no, i.e. -o no, no output file will be written.
 -r      : The PSI is read from restart file rather than from ./psi
 -n num  : only the first num WFs are read  from psi-file.
 -skip m : the first m WFs are skipped.
 -step step : Compute the population only every step-th WF.
 -tcut tcut : Compute the population only at time tcut.
 -w      : Overwrite the output-file adp.
 -q      : The quick algorithm is used.
 -mc     : Monte-Carlo integration is used.
 -nosob  : Turn off Sobol sampling.
 -noq    : Quick algorithm not used.
 -qtol qtol   : The quick algorithm cutoff set to qtol.
 -mctol mctol : The Monte-Carlo convergence set to mctol.
 -nogwp : GWP integration turned off
      Defaults: qtol=1.0d-4, mctol=1.0d-2
 -sleazy : qtol and mctol relaxed   to:  qtol=1.0d-3, mctol=1.0d-1
 -tight  : qtol and mctol tightened to:  qtol=1.0d-8, mctol=1.0d-3
 -DB PATH : The path to a direct dynamics DB.
 -dblocal : Take nearest point in DB - no interpolation.
 -saddle : Use saddle point approximation for vMCG.
 -saddle0 : Uses RHS saddle point approximation for vMCG.
 -saddle1 sadtol : Uses saddle point approximation for vMCG
          with MC for points with energy gap less than sadgap (in eV).
 -cdvr : Uses CDVR integration. Only works for single-set WFs
 -tddvr : Uses TDDVR integration. Only works for single-set WFs
 -batch : Do not output anything to the screen
 -ver    : Version information about the program 
 -h      : Print this help text. 
 -?      : Print this help text. 

 If an exact WF is used then the pes file must be in exact representation too!

 The arguments "Degrees of Freedom" have to be replaced by the names of the
 degrees of freedom for which densities are to be calculated.
 To calculate 2-dim densities separate the two DOF with a comma.
 EXAMPLES: adpop dof1 
           adpop dof3,dof4
           adpop dof1 dof2 dof3,dof4
 The reduced densities can be plotted with pladpop.
------------------------------------------------------------------------------