Skip to content

Commit

Permalink
Merge pull request #242 from ChemCryst/QLMR-atom-type-bonding-fix
Browse files Browse the repository at this point in the history
Fix bond calc for atoms starting with L, M or R.
  • Loading branch information
richardicooper authored Sep 6, 2021
2 parents a9bbe83 + 9fae993 commit 0e3dcff
Showing 1 changed file with 7 additions and 7 deletions.
14 changes: 7 additions & 7 deletions crystals/distangl.F
Original file line number Diff line number Diff line change
Expand Up @@ -7107,10 +7107,10 @@ SUBROUTINE XBCALC(INTERN)
END IF
NOTFND = NOTFND + 1
WRITE(CATTYP,'(A4)')ISTORE(M5)
IF (( CATTYP(1:1).NE.'Q' ) .AND.
1 ( CATTYP(1:1).NE.'L' ) .AND.
1 ( CATTYP(1:1).NE.'M' ) .AND.
1 ( CATTYP(1:1).NE.'R' )) THEN
IF (( CATTYP(1:2).NE.'Q ' ) .AND.
1 ( CATTYP(1:2).NE.'L ' ) .AND.
1 ( CATTYP(1:2).NE.'M ' ) .AND.
1 ( CATTYP(1:2).NE.'R ' )) THEN
WRITE(CMON,'(3A)')
1 'FYI: Element not in L40 or L29: ',ISTORE(M5),ICOL
CALL XPRVDU(NCVDU,1,0)
Expand Down Expand Up @@ -7164,9 +7164,9 @@ SUBROUTINE XBCALC(INTERN)
DO I5 = 0,N5-1
M5 = L5 + I5*MD5
WRITE(CATTYP,'(A4)')ISTORE(M5)
IF (( CATTYP(1:1).EQ.'L' ) .OR.
1 ( CATTYP(1:1).EQ.'M' ) .OR.
2 ( CATTYP(1:1).EQ.'R' )) CYCLE
IF (( CATTYP(1:2).EQ.'L ' ) .OR.
1 ( CATTYP(1:2).EQ.'M ' ) .OR.
2 ( CATTYP(1:2).EQ.'R ' )) CYCLE

NEXTLC = NFL
NFOUND = MAKE41( M5, NEXTLC )
Expand Down

0 comments on commit 0e3dcff

Please sign in to comment.