[molpro-user] Overriding geometry (and basis) reordering in molpro?

[Dr. Seth Olsen]

Hi Molpro-Users,

Is there any way to override MolPro's default tendency of reordering the 
atoms in a molecule to group all members of a given element together?  
I've been trying to use MERGE to splice together orbitals of benzene 
dimer from monomer calculations, and everything gets mangled because as 
soon as the dimer geometry is read in, the basis is reordered and the 
merged orbitals become useless.  I have turned all symmetry and extra 
symmetry mechanisms off, and given the carbons and hydrogens from each 
monomer their own basis group number, but it doesn't seem to turn off 
the default of reordering so that all C's are in one block and all H's 
in another.  I've attached the output, which also displays the output of 
MERGE before and after orthonormalization (I've gzipped the file).

I think that I can get this done with MERGE or MATROP by the end of the 
day, now that I know what the problem is.  However, this issue also 
comes up regularly when (for example) MOLEKEL or GABEDIT is used to 
visualize MOLDEN output created by MolPro.  Neither MOLEKEL nor GABEDIT 
seem capable of mapping basis to atoms if the ordering of the atoms is 
changed relative to the ordering of the basis (as in MolPro-output 
MOLDEN files, where the atoms are rearranged but the basis is not - 
although there is a number attached to each function to identify it's 
atom, it seems that neither MOLEKEL nor GABEDIT pay attention). 

Can this behavior be turned off?

Cheers,

Seth

-------------- next part --------------
A non-text attachment was scrubbed...
Name: dimer-frags-rohf-631p.out.gz
Type: application/x-gzip
Size: 82061 bytes
Desc: not available
URL: <http://www.molpro.net/pipermail/molpro-user/attachments/20060324/272a6808/attachment.gz>