Purpose: Enables 2D plotting of the wavefunction or PES.
 Usage:   showsys [-i -D -f -p -o -pes -nopes -ddraw -n -skip -step
                     -rst -nw -nofs -u -pop2 -pop2all -plgrd -DB
                     -shepord -ddlocal -ver -h -?]

Options : 
 -i DIR  : data is stored in directory DIR
 -f FILE : the wavefunction is read from FILE rather than from ./psi
           The string FILE may be a relative or a full path-name.
 -n num  : only the first num WFs are read from psi-file.
 -skip m : the first m WFs are skipped.
 -p FILE : the PES is read from FILE rather than from ./oper
           The string FILE may be a relative or a full path-name.
 -o FILE : The ouput is written to file FILE.pl rather than to ./surface.pl
           The string FILE may be a relative or a full path-name.
 -D DIR  : The output is written to directory DIR. This is useful,
           if one has no right to write to the name directory.
           The directory DIR must exist.
           The string DIR may be a relative or a full path-name.
 -pes    : only the OPER file is read,  i.e. no wavefunction plotting.
 -nopes  : the OPER file is not read, no potential plotting, no dia->ad trafo.
 -ddraw    : the raw adiabatic PES from a DD,  calculation is plotted.
 -rst    : the restart file is read, rather than the psi file.
 -n num  : only num WFs are processed.
 -skip m : the first m WFs are skipped.
 -step step : only every step's WF will be processed.
 -nw     : No weights are employed. Plot the "naked" DVR populations. (WF only).
 -nofs   : No transformation to fs. Use when 'time-not-fs' was set in mctdh.
 -u UNIT : The energy unit UNIT is applied. (default for pes: eV).
 -pop2   : the basis set occupations are plotted rather than grid
           populations. For FFT this implies momentum space representation.
           Note: -pop2 sets -nw.
 -pop2all: Similar to -pop2, but now all DOFs are transformed to second population.
           This is only useful for cuts. Note: -pop2all sets -pop2 and -nw.
 -plgrd MAXPLGRID : The maximum number of points for a grid.
 -DB PATH : The path to a direct dynamics DB.
 -shepord I : Order I will be used for the Shephard interpolation.
 -ddlocal : The local rather than interpolated energy is plotted.
 -ver    : Version information about the program.
 -h      : Print this help text. 
 -?      : Print this help text. 

  The program is menu driven. Just type "showsys" or "showsys -pop2"
  or "showsys -pes" or "showsys -rst" and then follow the menu.
  If there are more WFs on psi file than you want to inspect, use the
  options -n, -skip, or -step. Example: "showsys -skip 10 -n 2" will plot
  the WFs no. 11 and 12 (i.e. t=10 and 11 if tpsi=1fs),
  and "showsys -step 10 -n 2" will plot the WFs no. 1 and 11.
  If one wants to plot a potential one should use "showsys -pes",
  because otherwise the WF is read, which consumes a large amount of memory.
  NB1: The option -step may also be set by the menu-point 285.
       The WFs to be processed may also be selected by setting time bounds
       (menu-point 30).
  NB2: The option -pop2 and -pop2all are ignored for phifbr,sphfbr,k,external bases.
--------------------------------------------------------------------------------