hdiff output

r25567/amberinterface.f 2017-03-30 13:31:23.869228199 +0100 r25566/amberinterface.f 2017-03-30 13:31:24.345234494 +0100
 65:  65: 
 66: ! sf344> some basic IO 66: ! sf344> some basic IO
 67:         facc = 'A'      ! try to append 67:         facc = 'A'      ! try to append
 68:         IF (ncalls == 3) THEN  68:         IF (ncalls == 3) THEN 
 69:                 mdin = 'min_md.in' 69:                 mdin = 'min_md.in'
 70:         ELSE 70:         ELSE
 71:                 mdin = 'min.in' 71:                 mdin = 'min.in'
 72:         END IF 72:         END IF
 73:         mdout = 'min.out'      ! not used yet 73:         mdout = 'min.out'      ! not used yet
 74:         parm = 'coords.prmtop' 74:         parm = 'coords.prmtop'
 75:         prmtop = 'coords.prmtop' ! khs26> to pass to modamber9 
 76:         inpcrd = trim(adjustl(inpcrd1))     ! inpcrd file name defined in keywords.f  75:         inpcrd = trim(adjustl(inpcrd1))     ! inpcrd file name defined in keywords.f 
 77:         mdcrd = 'coords.mdcrd' 76:         mdcrd = 'coords.mdcrd'
 78:         restrt = 'coords.rst'    77:         restrt = 'coords.rst'   
 79:         mdinfo = 'mdinfo' 78:         mdinfo = 'mdinfo'
 80:         inpdip = 'inpdip' 79:         inpdip = 'inpdip'
 81:         rstdip = 'rstdip' 80:         rstdip = 'rstdip'
 82:         refc = 'refc' 81:         refc = 'refc'
 83:         mdvel = 'mdvel' 82:         mdvel = 'mdvel'
 84:         mden = 'mden' 83:         mden = 'mden'
 85:  84: 
318:         master = .false.317:         master = .false.
319:         natoms = natom318:         natoms = natom
320: 319: 
321:         IF(.NOT.allocated(coords1)) ALLOCATE(COORDS1(3*natom))320:         IF(.NOT.allocated(coords1)) ALLOCATE(COORDS1(3*natom))
322:         IF(.NOT.allocated(atmass1)) ALLOCATE(atmass1(natom))321:         IF(.NOT.allocated(atmass1)) ALLOCATE(atmass1(natom))
323: 322: 
324:         coords1(1:3*natom) = x(lcrd:lcrd+3*natom-1)323:         coords1(1:3*natom) = x(lcrd:lcrd+3*natom-1)
325: !  write(MYUNITNEW,*) 'sf344> check ', COORDS1(1:100)324: !  write(MYUNITNEW,*) 'sf344> check ', COORDS1(1:100)
326: ! Atomic masses:325: ! Atomic masses:
327:         atmass1(1:natom) = x(lmass:lmass+natom-1)326:         atmass1(1:natom) = x(lmass:lmass+natom-1)
328: 327: !        write(MYUNITNEW,*) 'atomic masses in amberinterface', atmass1(1:10)
329: !        IF (DEBUG) THEN 
330: !            write(MYUNITNEW,*) 'atomic masses in amberinterface', atmass1(1:10) 
331: !        END IF 
332: !close(5)328: !close(5)
333: !open (5,file='odata',status='old')329: !open (5,file='odata',status='old')
334: !write(MYUNITNEW,*) ih(1:lasth)330: !write(MYUNITNEW,*) ih(1:lasth)
335: !write(MYUNITNEW,*) 'sf344> exiting amber interface routine...',natoms,natom331: !write(MYUNITNEW,*) 'sf344> exiting amber interface routine...',natoms,natom
336: END SUBROUTINE AMBERINTERFACE332: END SUBROUTINE AMBERINTERFACE
337: 333: 
338: 334: 
339: SUBROUTINE AMBERENERGIES(y,grad,ereal,gradt,stest)335: SUBROUTINE AMBERENERGIES(y,grad,ereal,gradt,stest)
340: !  use commons336: !  use commons
341:   use modamber9337:   use modamber9


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0