hdiff output

r33146/potential.f 2017-08-09 09:30:15.580969176 +0100 r33145/potential.f 2017-08-09 09:30:15.892973324 +0100
 51:          LOGICAL PTEST, BOXTEST, file_exists 51:          LOGICAL PTEST, BOXTEST, file_exists
 52:  52: 
 53:          INTEGER J1, J2, J3, NN, MM, IPOT, NELEMENTS, NTYPE(105), NCOUNT, ISTART, NDUM, NPCALL, ECALL, FCALL, SCALL, 53:          INTEGER J1, J2, J3, NN, MM, IPOT, NELEMENTS, NTYPE(105), NCOUNT, ISTART, NDUM, NPCALL, ECALL, FCALL, SCALL,
 54:      1   J4, J5, J6, JSTART, ISTAT, NDUMMY,IDUM1,IDUM2, NATOMSSAVE 54:      1   J4, J5, J6, JSTART, ISTAT, NDUMMY,IDUM1,IDUM2, NATOMSSAVE
 55:          DOUBLE PRECISION P2, P3, BA, XLAMBDA, AMAT(3,3), AINV(3,3), TEMPX, BOXLX, 55:          DOUBLE PRECISION P2, P3, BA, XLAMBDA, AMAT(3,3), AINV(3,3), TEMPX, BOXLX,
 56:      1   TEMPXX, TEMPH(3,3), TEMPV(3), ENERGY1, ENERGY2,DOUBLEDUM1, 56:      1   TEMPXX, TEMPH(3,3), TEMPV(3), ENERGY1, ENERGY2,DOUBLEDUM1,
 57:      1   EDOUBLE, XTEMP, YTEMP, ZTEMP, GAMESR(3,3), GAMEST(3), 57:      1   EDOUBLE, XTEMP, YTEMP, ZTEMP, GAMESR(3,3), GAMEST(3),
 58:      2   TEMPA, TEMP(6), COORDSO(3*NATOMS), GRADO(3*NATOMS), GEMAX, 58:      2   TEMPA, TEMP(6), COORDSO(3*NATOMS), GRADO(3*NATOMS), GEMAX,
 59:      4   ETIME, FTIME, STIME, TIME, TIME0, DUMMY1, DUMMY2, DUMMY3, EPLUS, EMINUS, DUMMY, DOTOPT, DIST, DIST1, DIST2, RMAT(3,3) 59:      4   ETIME, FTIME, STIME, TIME, TIME0, DUMMY1, DUMMY2, DUMMY3, EPLUS, EMINUS, DUMMY, DOTOPT, DIST, DIST1, DIST2, RMAT(3,3)
 60:          INTEGER NSTART, NFINISH, NSTEP 60:          INTEGER NSTART, NFINISH, NSTEP
 61:          INTEGER LUNIT, LUNIT2, GETUNIT 61:          INTEGER LUNIT, GETUNIT
 62:          DOUBLE PRECISION,ALLOCATABLE,DIMENSION(:,:) ::  HGAUSS 62:          DOUBLE PRECISION,ALLOCATABLE,DIMENSION(:,:) ::  HGAUSS
 63:          CHARACTER(LEN=87) ESTRING 63:          CHARACTER(LEN=87) ESTRING
 64:          CHARACTER(LEN=80) GPSTRING, NSTRING, FSTRING, FNAME, FNAME2, GSTRING,FROMFILE 64:          CHARACTER(LEN=80) GPSTRING, NSTRING, FSTRING, FNAME, FNAME2, GSTRING,FROMFILE
 65:          CHARACTER(LEN=132) STRING 65:          CHARACTER(LEN=132) STRING
 66:          COMMON /STRINGS/ ESTRING, GPSTRING, NSTRING, FSTRING 66:          COMMON /STRINGS/ ESTRING, GPSTRING, NSTRING, FSTRING
 67:          DOUBLE PRECISION C1, C2, C3, IZ, ROT, EDISP, EIND 67:          DOUBLE PRECISION C1, C2, C3, IZ, ROT, EDISP, EIND
 68:          LOGICAL ETEST, SSTEST, YESNO, AMIDEFAIL,GAUSSIANTEST 68:          LOGICAL ETEST, SSTEST, YESNO, AMIDEFAIL,GAUSSIANTEST
 69:          COMMON /CAS/ AMAT, AINV, NELEMENTS, NTYPE 69:          COMMON /CAS/ AMAT, AINV, NELEMENTS, NTYPE
 70:          COMMON /PCALL/ NPCALL, ECALL, FCALL, SCALL, ETIME, FTIME, STIME 70:          COMMON /PCALL/ NPCALL, ECALL, FCALL, SCALL, ETIME, FTIME, STIME
 71:          LOGICAL KNOWE, KNOWG, KNOWH 71:          LOGICAL KNOWE, KNOWG, KNOWH
3216:             ENDDO3216:             ENDDO
3217:             DO J1=1,3*NATOMS3217:             DO J1=1,3*NATOMS
3218:                DO J2=1,3*NATOMS3218:                DO J2=1,3*NATOMS
3219:                   HESS(J2,J1)=HGAUSS(J2,J1)3219:                   HESS(J2,J1)=HGAUSS(J2,J1)
3220:                ENDDO3220:                ENDDO
3221:             ENDDO3221:             ENDDO
3222:             CLOSE(15)3222:             CLOSE(15)
3223:          ELSE IF (GAUSSIAN03) THEN3223:          ELSE IF (GAUSSIAN03) THEN
3224:             WRITE(*,'(A)') ' potential> Calling GAUSSIAN03'3224:             WRITE(*,'(A)') ' potential> Calling GAUSSIAN03'
3225:             CALL SYSTEM('cp gaussian.inp.pre gaussian.inp' )3225:             CALL SYSTEM('cp gaussian.inp.pre gaussian.inp' )
3226:             LUNIT=GETUNIT()3226:             OPEN(15,file="gaussian.inp",STATUS="OLD",POSITION="APPEND")
3227:             OPEN(LUNIT,file="gaussian.inp",STATUS="OLD",POSITION="APPEND") 
3228:             INQUIRE(FILE='gaussian.chk',EXIST=GAUSSIANTEST) 3227:             INQUIRE(FILE='gaussian.chk',EXIST=GAUSSIANTEST) 
3229:             IF(GAUSSIANTEST)THEN3228:             IF(GAUSSIANTEST)THEN
3230:                WRITE(LUNIT,"(A10)") "GUESS=READ"3229:                WRITE(15,"(A10)") "GUESS=READ"
3231:             ENDIF3230:             ENDIF
3232:             WRITE(LUNIT,*) ""3231:             WRITE(15,*) ""
3233:             WRITE(LUNIT,*) "OPTIM INTERFACE"3232:             WRITE(15,*) "OPTIM INTERFACE"
3234:             WRITE(LUNIT,*) ""3233:             WRITE(15,*) ""
3235:             WRITE(LUNIT,"(2I4)") GAUSSIANCHARGE,GAUSSIANMULTI3234:             WRITE(15,"(2I4)") GAUSSIANCHARGE,GAUSSIANMULTI
3236:             LUNIT2=GETUNIT()3235:             OPEN(7,file="gaussian.name")
3237:             OPEN(LUNIT2,file="gaussian.name") 
3238:             DO J1=1,NATOMS3236:             DO J1=1,NATOMS
3239:                READ(LUNIT2,*) STRING3237:                READ(7,*) STRING
3240:                WRITE(LUNIT,"(A2,3F20.10)") ADJUSTL(TRIM(STRING)),COORDS(3*J1-2:3*J1)3238:                WRITE(15,"(A2,3F20.10)") ADJUSTL(TRIM(STRING)),COORDS(3*J1-2:3*J1)
3241:             ENDDO3239:             ENDDO
3242:             CLOSE(LUNIT2)3240:             CLOSE(7)
3243:             CLOSE(LUNIT)3241:             CLOSE(15)
3244:             INQUIRE(FILE='gaussian.inp.after',EXIST=GAUSSIANTEST)3242:             INQUIRE(FILE='gaussian.inp.after',EXIST=GAUSSIANTEST)
3245:             IF(.NOT.GAUSSIANTEST) THEN3243:             IF(.NOT.GAUSSIANTEST) THEN
3246:                CALL SYSTEM('echo "" > gaussian.inp.after' )3244:                CALL SYSTEM('echo "" > gaussian.inp.after' )
3247:                CALL SYSTEM('echo "" >> gaussian.inp.after' )3245:                CALL SYSTEM('echo "" >> gaussian.inp.after' )
3248:                CALL SYSTEM('echo "" >> gaussian.inp.after' )3246:                CALL SYSTEM('echo "" >> gaussian.inp.after' )
3249:             ENDIF3247:             ENDIF
3250:             CALL SYSTEM('cat gaussian.inp.after >> gaussian.inp' )3248:             CALL SYSTEM('cat gaussian.inp.after >> gaussian.inp' )
3251:             IF (SSTEST) THEN3249:             IF (SSTEST) THEN
3252:                CALL SYSTEM('sed "s/Force/Freq=Noraman/Ig" gaussian.inp >gaussian.inp.temp')3250:                CALL SYSTEM('sed "s/Force/Freq=Noraman/Ig" gaussian.inp >gaussian.inp.temp')
3253:                CALL SYSTEM('mv gaussian.inp.temp gaussian.inp')3251:                CALL SYSTEM('mv gaussian.inp.temp gaussian.inp')


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0