Hi,
I would not use the wildcard option specifying g1tetra_prod_mdx* and simply
just specify the second trajectory and put startframe= 1000, endframe= 3000
or just simply make one big trajectory of both simulations and put
startframe=6000, endframe=8000 so that should work fine. With the way you
tried to run MMPBSA.py it simply wrote:
g1tetra_prod_mdx1 6000 8000 1
g1tetra_prod_mdx2 6000 8000 1
hence is the reason for the error
*WARNING in ptrajSetupIO(): trajectory start is > stop; no
configurations will be processed*
On Wed, Mar 23, 2011 at 12:44 PM, Chris Chris <alpharecept.yahoo.com> wrote:
> Please help me understand my issue here. I have run two simulations in
> series,
> the first for 5ns and the second for 5 ns. The second 5ns simulation is
> merely a
> continuation of the first. My prod.in file for the production simulation
> is as
> follows:
>
> &cntrl
> imin=0,irest=1,ntx=5,
> nstlim=2500000,dt=0.002,
> ntc=2,ntf=2,
> cut=8.0, ntb=2, ntp=1, taup=2.0,
> ntpr=500, ntwx=500,
> ntt=3, gamma_ln=2.0,
> temp0=300.0,
> /
> thus, I am generatring a frame per ps, correct? If so, the 2 simulations
> must
> comprise 10,000 frames- is this assumption correct?
>
> I am interested in reading frames 6,000-8,000 into mmpbsa.py and thus
> obtain
> binding data. I have the submit file for mmpbsa.py as follows:
>
> mpirun -np 4 /gpfs1/apps/chemistry/amber/amber11/bin/MMPBSA.py -O -i
> mmpbsa_g1tetra.in -o G1TETRA_RESULTS_MMPBSA.dat -sp g1tetra.top -cp
> g1tetra_nosolv.top -rp g1tetra_dimer1_no
> solv.top -lp g1tetra_dimer2_nosolv.top -y g1tetra_prod_mdx*
>
> My thinking is that the program will read in both trajectories and choose
> frames
> 6,000-8,000 and process them. However, this does not work as I get the
> following
> error message:
>
>
>
> FYI ptrajSetupIO(): trajin stop value (8000) is greater than the number of
> sets
> read in.
> Setting stop to the maximum value (5000)
> WARNING in ptrajSetupIO(): trajectory start is > stop; no
> configurations will be processed
>
> Thus, only 5,000 frames are read. However- 10,000 should be present.
>
> Can anyone help me with this?
>
> Thanks
>
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber
>
--
T. Dwight McGee Jr.
Quantum Theory Project
University of Florida
Graduate Student
dwight.mcgee.gmail.com
"Problems cannot be solved at the same level of awareness that created
them."
Albert Einstein
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Wed Mar 23 2011 - 10:30:02 PDT
