------------------------------------------------------- Amber 9 SANDER 2006 ------------------------------------------------------- | Run on 03/19/2009 at 13:35:33 [-O]verwriting output File Assignments: | MDIN: dna_wat_mdfiinal.in | MDOUT: dna_wat_mdfiinal.out |INPCRD: dna_wat_4min5.rst | PARM: dna_wat.prmtop |RESTRT: dna_wat_final.rst | REFC: refc | MDVEL: mdvel | MDEN: mden | MDCRD: dna_wat_mdfinal.mdcrd |MDINFO: mdinfo |INPDIP: inpdip |RSTDIP: rstdip |INPTRA: inptraj | Here is the input file: dna2: 10ps MD &cntrl imin = 0, irest = 1, ntx = 7, ntb = 2, pres0 = 1.0, ntp = 1, taup = 2.0, cut = 10.0, ntr = 0, ntc = 2, ntf = 2, tempi = 100.0, temp0 = 300.0, ntt = 3, gamma_ln = 1.0, nstlim = 5000, dt = 0.002, ntpr = 100, ntwx = 100, ntwr = 1000 / -------------------------------------------------------------------------------- 1. RESOURCE USE: -------------------------------------------------------------------------------- | Flags: getting new box info from bottom of inpcrd | INFO: Old style inpcrd file read | peek_ewald_inpcrd: Box info found |Largest sphere to fit in unit cell has radius = 24.756 | New format PARM file being parsed. | Version = 1.000 Date = 03/18/09 Time = 12:30:17 NATOM = 14545 NTYPES = 16 NBONH = 14029 MBONA = 550 NTHETH = 636 MTHETA = 844 NPHIH = 1354 MPHIA = 1480 NHPARM = 0 NPARM = 0 NNB = 22586 NRES = 4635 NBONA = 550 NTHETA = 844 NPHIA = 1480 NUMBND = 37 NUMANG = 73 NPTRA = 41 NATYP = 26 NPHB = 1 IFBOX = 2 NMXRS = 33 IFCAP = 0 NEXTRA = 0 NCOPY = 0 | Memory Use Allocated | Real 782098 | Hollerith 91907 | Integer 391395 | Max Pairs 8377920 | nblistReal 174540 | nblist Int 487552 | Total 43992 kbytes | Duplicated 0 dihedrals | Duplicated 0 dihedrals BOX TYPE: TRUNCATED OCTAHEDRON -------------------------------------------------------------------------------- 2. CONTROL DATA FOR THE RUN -------------------------------------------------------------------------------- General flags: imin = 0, nmropt = 0 Nature and format of input: ntx = 7, irest = 1, ntrx = 1 Nature and format of output: ntxo = 1, ntpr = 100, ntrx = 1, ntwr = 1000 iwrap = 0, ntwx = 100, ntwv = 0, ntwe = 0 ioutfm = 0, ntwprt = 0, idecomp = 0, rbornstat= 0 Potential function: ntf = 2, ntb = 2, igb = 0, nsnb = 25 ipol = 0, gbsa = 0, iesp = 0 dielc = 1.00000, cut = 10.00000, intdiel = 1.00000 scnb = 2.00000, scee = 1.20000 Frozen or restrained atoms: ibelly = 0, ntr = 0 Molecular dynamics: nstlim = 5000, nscm = 1000, nrespa = 1 t = 0.00000, dt = 0.00200, vlimit = 20.00000 Langevin dynamics temperature regulation: ig = 71277 temp0 = 300.00000, tempi = 100.00000, gamma_ln= 1.00000 Pressure regulation: ntp = 1 pres0 = 1.00000, comp = 44.60000, taup = 2.00000 SHAKE: ntc = 2, jfastw = 0 tol = 0.00001 Ewald parameters: verbose = 0, ew_type = 0, nbflag = 1, use_pme = 1 vdwmeth = 1, eedmeth = 1, netfrc = 1 Box X = 60.640 Box Y = 60.640 Box Z = 60.640 Alpha = 109.471 Beta = 109.471 Gamma = 109.471 NFFT1 = 60 NFFT2 = 60 NFFT3 = 60 Cutoff= 10.000 Tol =0.100E-04 Ewald Coefficient = 0.27511 Interpolation order = 4 -------------------------------------------------------------------------------- 3. ATOMIC COORDINATES AND VELOCITIES -------------------------------------------------------------------------------- FATAL: Could not read velocities from dna_wat_4min5.rst