hdiff output

r30747/newconnodata.f 2017-01-21 10:38:48.898330642 +0000 r30746/newconnodata.f 2017-01-21 10:38:49.182351762 +0000
 23: C  Make an odata tspath file for the minimum in LOCALPOINTS1 23: C  Make an odata tspath file for the minimum in LOCALPOINTS1
 24: C  Note that the last line of odata.tspath should be POINTS, like odata.connect. 24: C  Note that the last line of odata.tspath should be POINTS, like odata.connect.
 25: C  For historical reasons the POINTS line should not be present for the old style 25: C  For historical reasons the POINTS line should not be present for the old style
 26: C  odata.tsseach file! 26: C  odata.tsseach file!
 27: C 27: C
 28:       SUBROUTINE NEWCONNODATA(CONNID,LOCALPOINTS1) 28:       SUBROUTINE NEWCONNODATA(CONNID,LOCALPOINTS1)
 29:       USE COMMONS 29:       USE COMMONS
 30:       USE UTILS,ONLY : GETUNIT 30:       USE UTILS,ONLY : GETUNIT
 31:       IMPLICIT NONE 31:       IMPLICIT NONE
 32:       INTEGER J2, CONNID, STATUS, LUNIT 32:       INTEGER J2, CONNID, STATUS, LUNIT
 33:       DOUBLE PRECISION LOCALPOINTS1(3*NATOMS), SAVEPOINTS(3*NATOMS), RANDOM, DPRAND 33:       DOUBLE PRECISION LOCALPOINTS1(3*NATOMS), SAVEPOINTS(3*NATOMS)
 34:       CHARACTER(LEN=20) UNSTRING 34:       CHARACTER(LEN=20) UNSTRING
 35:       CHARACTER(LEN=10) CONNSTR 35:       CHARACTER(LEN=10) CONNSTR
 36:       CHARACTER(LEN=80) FPOO 36:       CHARACTER(LEN=80) FPOO
 37:  37: 
 38:       WRITE(CONNSTR,'(I10)') CONNID 38:       WRITE(CONNSTR,'(I10)') CONNID
 39:       CALL MYSYSTEM(STATUS,DEBUG,'cat odata.tspath > odata.'//TRIM(ADJUSTL(CONNSTR))) 39:       CALL MYSYSTEM(STATUS,DEBUG,'cat odata.tspath > odata.'//TRIM(ADJUSTL(CONNSTR)))
 40:  40: 
 41:       DO J2=1,3*NATOMS 
 42:          RANDOM=DPRAND() 
 43:          LOCALPOINTS1(J2)=LOCALPOINTS1(J2)+(RANDOM-0.5D0)*2.0D0*PERTVALUE 
 44:       ENDDO 
 45:  
 46:       IF (CHARMMT) THEN 41:       IF (CHARMMT) THEN
 47:          IF (MACHINE) THEN ! SAT 42:          IF (MACHINE) THEN ! SAT
 48:             DO J2=1,3*NATOMS 43:             DO J2=1,3*NATOMS
 49:                SAVEPOINTS(J2)=LOCALPOINTS1(J2) 44:                SAVEPOINTS(J2)=LOCALPOINTS1(J2)
 50:             ENDDO 45:             ENDDO
 51:             CALL CHARMMDUMP(SAVEPOINTS,'points1.inp.'//TRIM(ADJUSTL(CONNSTR))) 46:             CALL CHARMMDUMP(SAVEPOINTS,'points1.inp.'//TRIM(ADJUSTL(CONNSTR)))
 52:          ELSE 47:          ELSE
 53:             DO J2=1,3*NATOMS 48:             DO J2=1,3*NATOMS
 54:                SAVEPOINTS(J2)=LOCALPOINTS1(J2) 49:                SAVEPOINTS(J2)=LOCALPOINTS1(J2)
 55:             ENDDO 50:             ENDDO


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0