AMBER Archive (2003)Subject: Re: AMBER: nmanal's error message
From: David Case (case_at_scripps.edu)
Date: Mon Aug 18 2003 - 18:24:11 CDT
On Mon, Aug 18, 2003, Qiang Lu wrote:
>
> I used new2oldparm to convert parm7 to parm6, then ran nmanal.
>
Try the following patch to amber7/src/nmode/rdparm2.f:
*** rdparm2.f 25 Jul 2001 23:51:39 -0000 7.2
--- rdparm2.f 18 Aug 2003 23:16:49 -0000
***************
*** 27,33 ****
READ(20,9108) (IHEAD1(I),I=1,20)
READ(20,9118) NATOM,NTYPES,NBONH,MBONA,NTHETH,MTHETA,NPHIH,
+ MPHIA,NHPARM,NPARM,NNB,NRES,NBONA,NTHETA,NPHIA,
! + NUMBND,NUMANG,NPTRA,NATYP,NPHB
C
C ----- CHECK FOR POSSIBLE OVERFLOWS OF ARRAY BOUNDS -----
C
--- 27,34 ----
READ(20,9108) (IHEAD1(I),I=1,20)
READ(20,9118) NATOM,NTYPES,NBONH,MBONA,NTHETH,MTHETA,NPHIH,
+ MPHIA,NHPARM,NPARM,NNB,NRES,NBONA,NTHETA,NPHIA,
! + NUMBND,NUMANG,NPTRA,NATYP,NPHB,ifpert,idum,idum,
! + idum,idum,idum,idum,idum,idum,idum,idum
C
C ----- CHECK FOR POSSIBLE OVERFLOWS OF ARRAY BOUNDS -----
C
Let me know this helps. [Note: if you have the amber6 version of nmanal,
that should also work. In amber6 there was a test to see if one had the
Amber 2 (!!) or Amber 3 and later prmtop files. This test got lost in the
conversion to Amber 7, resulting in code that was looking for Amber 2 (ca.
1985) formatted files....]
...thanks for pointing this out....dac
--
==================================================================
David A. Case | e-mail: case_at_scripps.edu
Dept. of Molecular Biology, TPC15 | fax: +1-858-784-8896
The Scripps Research Institute | phone: +1-858-784-9768
10550 N. Torrey Pines Rd. | home page:
La Jolla CA 92037 USA | http://www.scripps.edu/case
==================================================================
-----------------------------------------------------------------------
The AMBER Mail Reflector
To post, send mail to amber_at_scripps.edu
To unsubscribe, send "unsubscribe amber" to majordomo_at_scripps.edu
|