MOLPRO

4 GENERAL PROGRAM STRUCTURE 16
in dump record 2100.2, the second SCF in record 2101.2, the first MCSCF in 2140.2, the second
MCSCF in 2141.2 and so on. Note that these numbers refer to the occurrence in the input and
not on the order in which the calculations are performed in the actual run. If an input or part of
it is repeated using DO loops, this ensures that each calculation will start with the orbitals from
the corresponding orbitals from the previous cycle, as long as the order of the commands in the
input remains unchanged. If for instance the first SCF would be skipped in the second cycle
using some IF / ENDIF structure, the second SCF would still use record 2101.2. Thus, under
most circumstances the program defaults are appropriate, and the user does not have to specify
the records.
After a restart this logic will still work correctly if the number and sequence of SCF and MCSCF
commands is kept unchanged. Thus, if you want to skip certain parts of the input after a restart,
it is recommended to use IF / ENDIF structures or the GOTO command rather than to delete or
comment certain commands. If for some reason this is not possible, the START and ORBITAL
directives can be used to specify explicitely the records to be used.
In general we recommend the use of program defaults whenever possible, since this minimizes
the probability of input errors and frustration!
After completion of each program step, MOLPRO prints a summary of the records on each file.
4.4 Restart
Information from the permanent files is automatically recovered in subsequent calculations.
This can be controlled using the RESTART directive.
4.5 Data set manipulation
It is possible to truncate files and rename or copy records using the DATA command. Several
standard matrix operations can be performed with MATROP, e.g., printing records, linearly
combining or multiplying matrices, or forming the trace of a product of two matrices.
4.6 Memory allocation
MOLPRO allocates memory dynamically as required by the user on the MEMORY card. Thus it is
not necessary to maintain different versions of the program with different memory sizes. If the
MEMORY command is omitted, the program will use a default memory size, which depends on
the hardware used and how the program was installed. Note that, on Unix machines, the default
memory can be set on the molpro command line using the flag -m.
4.7 Multiple passes through the input
It is possible to perform loops over parts of the input using DO loops, very much as in FORTRAN
programs. DO loops may be nested to any reasonable depth. This can be conveniently used, for
instance, to compute automatically whole potential energy surfaces.
4.8 Symmetry
MOLPRO can use Abelian point group symmetry only. For molecules with degenerate symmetry,
an Abelian subgroup must be used — e.g., C 2v or D 2h for linear molecules. The symmetry