hdiff output

r33409/optim_quip_wrapper.f90 2017-10-23 10:30:13.417206308 +0100 r33408/optim_quip_wrapper.f90 2017-10-23 10:30:13.637209214 +0100
  1: SUBROUTINE OPTIM_QUIP_WRAPPER(NAT,CRDS,GRAD,ENRG)  1: SUBROUTINE OPTIM_QUIP_WRAPPER(NAT,CRDS,GRAD,ENRG,BOXX,BOXY,BOXZ)
  2:   2: 
  3: USE COMMONS, ONLY : NATOMS  3: USE COMMONS, ONLY : NATOMS
  4: USE KEY, ONLY : QUIPZ, QUIPLATT 
  5:   4: 
  6: IMPLICIT NONE  5: IMPLICIT NONE
  7:   6: 
  8: INTEGER, INTENT(IN) :: NAT  7: INTEGER, INTENT(IN) :: NAT
  9: DOUBLE PRECISION, DIMENSION(3*NAT), INTENT(IN) :: CRDS  8: DOUBLE PRECISION, DIMENSION(3*NAT), INTENT(IN) :: CRDS
 10: DOUBLE PRECISION, DIMENSION(3*NAT), INTENT(OUT) :: GRAD  9: DOUBLE PRECISION, DIMENSION(3*NAT), INTENT(OUT) :: GRAD
 11: DOUBLE PRECISION, INTENT(OUT) :: ENRG 10: DOUBLE PRECISION, INTENT(OUT) :: ENRG
  11: DOUBLE PRECISION, INTENT(IN) :: BOXX, BOXY, BOXZ
 12:  12: 
 13: INTEGER :: I 13: INTEGER :: I
 14: DOUBLE PRECISION, DIMENSION(3,3) :: LATT, VIR 14: DOUBLE PRECISION, DIMENSION(3,3) :: LATT, VIR
 15: INTEGER ATOMICNUMBER(NATOMS) 15: INTEGER ATOMICNUMBER(NATOMS)
 16: DOUBLE PRECISION, DIMENSION(3,NATOMS) :: NEW_COORDS, FORCES 16: DOUBLE PRECISION, DIMENSION(3,NATOMS) :: NEW_COORDS, FORCES
 17: CHARACTER(LEN=3), DIMENSION(NATOMS) :: SYMB 17: CHARACTER(LEN=3), DIMENSION(NATOMS) :: SYMB
 18:  18: 
 19: ! LATT = RESHAPE((/ BOXX, 0.0D0, 0.0D0, 0.0D0, BOXY, 0.0D0, 0.0D0, 0.0D0, BOXZ /), SHAPE(LATT)) 19: LATT = RESHAPE((/ BOXX, 0.0D0, 0.0D0, 0.0D0, BOXY, 0.0D0, 0.0D0, 0.0D0, BOXZ /), SHAPE(LATT))
 20: LATT(1:3,1:3)=QUIPLATT(1:3,1:3) 
 21:  20: 
 22: ATOMICNUMBER(1:NATOMS)=QUIPZ 21: ATOMICNUMBER(1:NATOMS)=14
 23:  22: 
 24: NEW_COORDS = RESHAPE(CRDS, (/ 3, NAT /)) 23: NEW_COORDS = RESHAPE(CRDS, (/ 3, NAT /))
 25:  24: 
  25: !CALL QUIP_WRAPPER(NAT,LATT,SYMB,NEW_COORDS,QARGSTRING,ENRG,FORCES,VIR)
  26: 
 26: ! PRINT*,'before CALL QUIP_WRAPPER_SIMPLE' 27: ! PRINT*,'before CALL QUIP_WRAPPER_SIMPLE'
 27: CALL QUIP_WRAPPER_SIMPLE(NAT,LATT,ATOMICNUMBER,NEW_COORDS,ENRG,FORCES,VIR) 28: CALL QUIP_WRAPPER_SIMPLE(NAT,LATT,ATOMICNUMBER,NEW_COORDS,ENRG,FORCES,VIR)
 28: ! PRINT*,'after CALL QUIP_WRAPPER_SIMPLE' 29: ! PRINT*,'after CALL QUIP_WRAPPER_SIMPLE'
 29: ! PRINT*,'NEW_COORDS=',NEW_COORDS 30: ! PRINT*,'NEW_COORDS=',NEW_COORDS
 30: ! PRINT*,'ENRG=',ENRG 31: ! PRINT*,'ENRG=',ENRG
 31: ! PRINT*,'FORCES=',FORCES 32: ! PRINT*,'FORCES=',FORCES
 32:  33: 
 33: !Virial is not passed back to OPTIM 34: !Virial is not passed back to OPTIM
 34:  35: 
 35: GRAD = RESHAPE(FORCES, SHAPE(GRAD)) 36: GRAD = RESHAPE(FORCES, SHAPE(GRAD))


legend
Lines Added 
Lines changed
 Lines Removed

hdiff - version: 2.1.0