IUBio GIL .. BIOSCI/Bionet News .. Biosequences .. Software .. FTP

%GENRES-ERR: exceeded MAXA (PSF) parameter --> recompile program

Eugene Leitl ui22204 at mail.lrz-muenchen.de
Wed May 6 13:17:23 EST 1998


Dear X-PLORERS,

I have a fairly large system:


ATOM  45084  OH2 TIP3 1563      28.614  -4.056 -55.687  0.00  0.00      W101
ATOM  45085  H1  TIP3 1563      28.619  -4.272 -54.754  0.00  0.00      W101
ATOM  45086  H2  TIP3 1563      29.497  -3.736 -55.872  0.00  0.00      W101
END

which I'd like to create a structure file from. Using script:


!******************************************************
!* Make structure file for water shell                *
!* Created by SOLVATE 1.0 (C) 1996 Helmut Grubmueller *
!* http://www.imo.physik.uni-muenchen.de/~grubi       *
!* Helmut.Grubmueller at Physik.uni-muenchen.de          *
!******************************************************


struct @dry-lipidbox.psf end

topo @toph19.sol end
topo @toph19.nacl end
param @param19.sol end

param @param19.nacl end

segment
  name=W100
  molecule
    name=TIP3
    number=9999
  end
end
vector do (resid = encode(decode(resid) + 0))   (segid W100)

segment
  name=W101
  molecule
    name=TIP3
    number=1563
  end
end
vector do (resid = encode(decode(resid) + 0))   (segid W101)

coord @wetted-dry-lipidbox.pdb

vector do (q = 0.0) (resname TIP3)
vector do (b = 0.0) (resname TIP3)

vector do (b = 0.0) (resname INA)

vector do (b = 0.0) (resname ICL)

write struct output=wetted-dry-lipidbox.psf end
write coors  output=wetted-dry-lipidbox.pdb end

stop

the system aborts with:

MOLECULE>    name=TIP3 
 MOLECULE>    number=9999 
 MOLECULE>  end 
 %GENRES-ERR: exceeded MAXA (PSF) parameter --> recompile program
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 BOMLEV=    0 reached.  Program execution will be terminated.
 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
 Subroutine DIE called . Terminating

Having set MAXA and MAXA2 to 300000 and  X-PLOR(online) Version 3.851
recompiled with fort77 generated a binary which fails in the same way.

[md at liposome cleanbox]$ /usr/home/md/XPLOR/object_library/xplor-fort77.exe
 X-PLOR(online) Version 3.851

 Author: Axel T. Brunger
 X-PLOR(online): Copyright 1996 by Yale University
 X-PLOR: Copyright 1987 by Harvard U., 1988-1996 by Yale U.

 User: md           on: linux           at:  6-May-98 22:13:22
X-PLOR>?
 XPLOR: current counts (number in parenthesis is maximum)
 NATOM=          0(MAXA=       24000)  NBOND=          0(MAXB=       24000)
 NTHETA=         0(MAXT=       40000)  NGRP=           0(MAXGRP=     24000)
 NPHI=           0(MAXP=       50000)  NIMPHI=         0(MAXIMP=     24000)  
 NDON=           0(MAXPAD=      4000)  NACC=           0(MAXPAD=      4000)
 NNB=            0(MAXNB=       3000) 
 X-PLOR: list of open files:
         INPUT read             
         OUTPUT write            

Is 24000 really the end of the line or am I stupid?
Please help. I'd settle for an ELF binary for Linux with the constraint 
relaxed.

Thanks in advance,

Regards,
Eugene Leitl, MSU Polymer Sci. Dept.



More information about the X-plor mailing list

Send comments to us at archive@iubioarchive.bio.net