Hi again,

I have tried to solve the problem that I described last Friday but without success. I was looking more carefully the log file (that was submitted in my first post) and there is something really strange. Maybe some of the more expert users could help me with this. When the program starts to run it generate a file pkci (that is the name of my sequence) that contains coordinates for the model. I assume this is the ini file (I am not sure why the ini extension is omitted). The number of atoms in this file is 890 (including 2 Zn atoms). I used this file to extract the ID for the residues coordinated with the Zn, information that is input in the ali file. The special restraints are read correctly as is evident from the log file:

Dynamically allocated memory at amaxseq [B,kB,MB]: 5097041 4977.579 4.861

TOP_________> 216 163 CALL ROUTINE = 'special_restraints'

TOP_________> 217 756 ADD_RESTRAINT ATOM_IDS = 'SG:14' 'ZN:55', RESTRAINT_PARAM;

ETERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 218 757 ADD_RESTRAINT ATOM_IDS = 'SG:17' 'ZN:55', RESTRAINT_PARAM;

ETERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 219 758 ADD_RESTRAINT ATOM_IDS = 'NE2:39' 'ZN:55', RESTRAINT_PARA;

METERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 220 759 ADD_RESTRAINT ATOM_IDS = 'SG:42' 'ZN:55', RESTRAINT_PARAM;

ETERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 221 760 ADD_RESTRAINT ATOM_IDS = 'SG:31' 'ZN:54', RESTRAINT_PARAM;

ETERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 222 761 ADD_RESTRAINT ATOM_IDS = 'SG:34' 'ZN:54', RESTRAINT_PARAM;

ETERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 223 762 ADD_RESTRAINT ATOM_IDS = 'SG:50' 'ZN:54', RESTRAINT_PARAM;

ETERS = 2 1 1 22 2 2 0 3.5 0.1

Dynamically allocated memory at amprmcns [B,kB,MB]: 5097191 4977.726 4.861

TOP_________> 224 763 ADD_RESTRAINT ATOM_IDS = 'NE2:53' 'ZN:54', RESTRAINT_PARA;

METERS = 2 1 1 22 2 2 0 3.5 0.1

After that the model file is read again:

openf5__224_> Open 18 UNKNOWN SEQUENTIAL pkci.D00000001

TOP_________> 243 218 READ_MODEL FILE = MODEL

openf5__224_> Open 11 OLD SEQUENTIAL pkci.ini

openf5__224_> Open 11 OLD SEQUENTIAL pkci.ini

rdatm___297_> Segments, residues, atoms: 2 55 438

rdatm___298_> Segment: 1 1 53 436

rdatm___298_> Segment: 2 54 55 2

TOP_________> 244 219 CALL ROUTINE = 'select_atoms'

But the strange thing is that the number of atoms in the file is 438 not 890 as in the pkci file generated initially

So after a few lines the log file gives an out of bounds message:

csrrng__299E> Atom index is out of bounds: 441 438

Any ideas of what is going on here? By the way, I am running in windows.

Thanks

Alfredo Cardenas