Commit a3620fcc authored by Mickael Delcey's avatar Mickael Delcey

Fix SUBS (old REDL) keyword for SOC

parent 10abab18
Pipeline #66020554 passed with stages
in 283 minutes and 30 seconds
......@@ -482,20 +482,8 @@ C printing threshold
! At the moment memory is not reduced
!
IF(REDUCELOOP) THEN
EX=ENSOR(1)
L=1
LD=1
DO ISO = 2, NSS
If (ABS(ENSOR(ISO)-EX).gt.1.0D-8) Then
LD = LD + 1
EX = ENSOR(ISO)
Else
L = L + 1
End If
If (LD.gt.LOOPDIVIDE) Exit
End Do
IEND = L
JSTART = L+1
IEND = LOOPDIVIDE
JSTART = LOOPDIVIDE+1
ELSE
IEND = NSS
JSTART = 1
......
......@@ -194,7 +194,7 @@ C Nr of spin states and division of loops:
IF(ISTATE.GT.LOOPDIVIDE) CYCLE
LOOPDIVIDE_TEMP = LOOPDIVIDE_TEMP + MPLET
END DO
* LOOPDIVIDE = LOOPDIVIDE_TEMP
LOOPDIVIDE = LOOPDIVIDE_TEMP
CALL GETMEM('UTOTR','ALLO','REAL',LUTOTR,NSS**2)
CALL GETMEM('UTOTI','ALLO','REAL',LUTOTI,NSS**2)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment