Skip to content

Commit 9fae993

Browse files
Fix bond calc for atoms starting with L, M or R.
Revert bug from March where L, M, R are ignored in bonding calc (but also accidentally ignored bonds from Li, Mg, Mn, Ru, Re, etc). Effect was sometimes missed because reverse direction bonds were allowed, so effect depended on atom order).
1 parent a9bbe83 commit 9fae993

File tree

1 file changed

+7
-7
lines changed

1 file changed

+7
-7
lines changed

crystals/distangl.F

Lines changed: 7 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -7107,10 +7107,10 @@ SUBROUTINE XBCALC(INTERN)
71077107
END IF
71087108
NOTFND = NOTFND + 1
71097109
WRITE(CATTYP,'(A4)')ISTORE(M5)
7110-
IF (( CATTYP(1:1).NE.'Q' ) .AND.
7111-
1 ( CATTYP(1:1).NE.'L' ) .AND.
7112-
1 ( CATTYP(1:1).NE.'M' ) .AND.
7113-
1 ( CATTYP(1:1).NE.'R' )) THEN
7110+
IF (( CATTYP(1:2).NE.'Q ' ) .AND.
7111+
1 ( CATTYP(1:2).NE.'L ' ) .AND.
7112+
1 ( CATTYP(1:2).NE.'M ' ) .AND.
7113+
1 ( CATTYP(1:2).NE.'R ' )) THEN
71147114
WRITE(CMON,'(3A)')
71157115
1 'FYI: Element not in L40 or L29: ',ISTORE(M5),ICOL
71167116
CALL XPRVDU(NCVDU,1,0)
@@ -7164,9 +7164,9 @@ SUBROUTINE XBCALC(INTERN)
71647164
DO I5 = 0,N5-1
71657165
M5 = L5 + I5*MD5
71667166
WRITE(CATTYP,'(A4)')ISTORE(M5)
7167-
IF (( CATTYP(1:1).EQ.'L' ) .OR.
7168-
1 ( CATTYP(1:1).EQ.'M' ) .OR.
7169-
2 ( CATTYP(1:1).EQ.'R' )) CYCLE
7167+
IF (( CATTYP(1:2).EQ.'L ' ) .OR.
7168+
1 ( CATTYP(1:2).EQ.'M ' ) .OR.
7169+
2 ( CATTYP(1:2).EQ.'R ' )) CYCLE
71707170

71717171
NEXTLC = NFL
71727172
NFOUND = MAKE41( M5, NEXTLC )

0 commit comments

Comments
 (0)