Re: [AMBER] error while running mmpbsa.py in amber12

From: Alessandro Contini <alessandro.contini.unimi.it>
Date: Wed, 30 Jan 2013 15:51:43 +0100

Hi Jason,
here is the output when "print output" is placed after outputlines =
output.split('\n')"



Reading command-line arguments and input files...
Loading and checking parameter files for compatibility...
sander found! Using /usr/local/amber12/bin/sander
mmpbsa_py_energy found! Using /usr/local/amber12/bin/mmpbsa_py_energy
cpptraj found! Using /usr/local/amber12/bin/cpptraj
Preparing trajectories for simulation...

CPPTRAJ: Trajectory Analysis. V12.5
    ___ ___ ___ ___
     | \/ | \/ | \/ |
    _|_/\_|_/\_|_/\_|_

INPUT: Reading Input from STDIN, type "go" to run, "quit" to exit:
>
[trajin /archive/acontini/tubuline/md/tetramer/cp1/from_empty_equil12/csearch/tetramer_cp1_csearch.prod4.mdcrd.bz2]
Warning: /archive/acontini/tubuline/md/tetramer/cp1/from_empty_equil12/csearch/tetramer_cp1_csearch.prod4.mdcrd.bz2: Uncompressed size of trajectory could not be determined.
         (This is normal for bzipped files)
         Number of frames could not be calculated.
         Frames will be read until EOF.
        [tetramer_cp1_csearch.prod4.mdcrd.bz2] contains an unknown number of
frames.
>
INPUT TRAJECTORIES:
  [tetramer_cp1_csearch.prod4.mdcrd.bz2] is an AMBER trajectory, Parm 0
(with box info), unknown #frames, start=0 offset=1
  Coordinate processing will occur on an unknown number of frames.

PARAMETER FILES:

0: /archive/acontini/tubuline/md/tetramer/cp1/from_empty_equil12/csearch/tetramer_cp1_csearch.top, 254878 atoms, 77696 res, ortho. box, 75972 mol, 75920 solvent mol

REFERENCE COORDS:
  No reference coordinates.
  No frames defined.

OUTPUT TRAJECTORIES:
  No files.

ACTIONS: Initializing 0 actions:
BEGIN TRAJECTORY PROCESSING:
Progress: '+' = 200 frames processed.
.....................................................
PARM [tetramer_cp1_csearch.top]: Setting up 0 actions.
----- [tetramer_cp1_csearch.prod4.mdcrd.bz2] (1-EOF, 1) -----
         0 +++++Error: FileIO_Bzip2::Read: BZ2_bzRead error:
[BZ_SEQUENCE_ERROR]
                        size=1 count=6193561

Read 1000 frames and processed 1000 frames.

ACTION OUTPUT:

DATASETS:
  There are no data sets set up for analysis.


  File "/usr/local/amber12/bin/MMPBSA.py", line 516, in <module>
    str(external_progs['cpptraj']))
  File "/usr/local/amber12/bin/MMPBSA_mods/make_trajs.py", line 57, in
make_trajectories
    traj = Trajectory(FILES.solvated_prmtop, FILES.mdcrd, cpptraj)
  File "/usr/local/amber12/bin/MMPBSA_mods/make_trajs.py", line 453, in
__init__
    self.Query()
  File "/usr/local/amber12/bin/MMPBSA_mods/make_trajs.py", line 592, in
Query
    num_frames = int(line.split()[5])
ValueError: invalid literal for int() with base 10: 'an'
Exiting. All files have been retained.


Might be an issue with bzip2? In the meanwhile I try to decompress
before running mmpbsa.py

Thanks again

Alessandro



Il giorno mer, 30/01/2013 alle 08.35 -0500, Jason Swails ha scritto:
> At line 587, can you add
>
> print output
>
> either after or before "outputlines = output.split('\n')
>
> in $AMBERHOME/bin/MMPBSA_mods/make_trajs.py then re-run it again and send
> me the output?
>
> Thanks!
> Jason
>
> On Wed, Jan 30, 2013 at 7:54 AM, Alessandro Contini <
> alessandro.contini.unimi.it> wrote:
>
> > Dear Jason, thanks for your answer.
> > I've rerun the job as suggested and here is the output:
> >
> > Reading command-line arguments and input files...
> > Loading and checking parameter files for compatibility...
> > sander found! Using /usr/local/amber12/bin/sander
> > mmpbsa_py_energy found! Using /usr/local/amber12/bin/mmpbsa_py_energy
> > cpptraj found! Using /usr/local/amber12/bin/cpptraj
> > Preparing trajectories for simulation...
> > File "/usr/local/amber12/bin/MMPBSA.py", line 516, in <module>
> > str(external_progs['cpptraj']))
> > File "/usr/local/amber12/bin/MMPBSA_mods/make_trajs.py", line 57, in
> > make_trajectories
> > traj = Trajectory(FILES.solvated_prmtop, FILES.mdcrd, cpptraj)
> > File "/usr/local/amber12/bin/MMPBSA_mods/make_trajs.py", line 453, in
> > __init__
> > self.Query()
> > File "/usr/local/amber12/bin/MMPBSA_mods/make_trajs.py", line 591, in
> > Query
> > num_frames = int(line.split()[5])
> > ValueError: invalid literal for int() with base 10: 'an'
> > Exiting. All files have been retained.
> >
> >
> > Thanks,
> >
> > Alessandro
> >
> > Il giorno mer, 30/01/2013 alle 07.29 -0500, Jason Swails ha scritto:
> >
> > > On Wed, Jan 30, 2013 at 6:47 AM, Alessandro Contini <
> > > alessandro.contini.unimi.it> wrote:
> > >
> > > > Dear community,
> > > > I've recently installed amber12 (composer XE 2013.1.117) and I'm
> > > > experiencing an error while running mmpbsa.py on a system where I've
> > > > previously run the very same calculation with amber11.
> > > >
> > > > here follows the input:
> > > >
> > > > GB/PBSA analysis
> > > > &general
> > > > startframe=1, endframe=1000, interval=20, verbose=2, keep_files=0,
> > > > strip_mask=:WAT,MG2,Na+,GTP,GDP,:1-438
> > > > /
> > > > &gb
> > > > igb=2, saltcon=0.150,
> > > > /
> > > > &pb
> > > > istrng=0.15, fillratio=4.0
> > > > /
> > > >
> > > >
> > > > I've launched mmpbsa.py trough the following command:
> > > >
> > > > /usr/local/amber12/bin/MMPBSA.py -O -i mmpbsa_py.in -o
> > MMPBSA_PB-GB.txt
> > > > -sp solvated_complex.top -cp complex.top -rp rec.top -lp lig.top -y
> > > > traj.mdcrd.bz2 >& mmpbsa_py.out
> > > >
> > > > and this is the output:
> > > >
> > > > Reading command-line arguments and input files...
> > > > Loading and checking parameter files for compatibility...
> > > > mmpbsa_py_energy found! Using /usr/local/amber12/bin/mmpbsa_py_energy
> > > > cpptraj found! Using /usr/local/amber12/bin/cpptraj
> > > > Preparing trajectories for simulation...
> > > > ValueError: invalid literal for int() with base 10: 'an'
> > > >
> > >
> > > This is an uncaught Python exception. I'm not sure where it's coming
> > from.
> > > Can you put "debug_printlevel=2" in the &general section of your input
> > > file and run it again?
> > >
> > > Thanks,
> > > Jason
> > >
> > _______________________________________________
> > AMBER mailing list
> > AMBER.ambermd.org
> > http://lists.ambermd.org/mailman/listinfo/amber
> >
>
>
>


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Jan 30 2013 - 07:00:03 PST
Custom Search