2 ! { dg
-options
"-O3 -std=legacy" }
4 SUBROUTINE EFTORD
(DM
,CHDINT
,L4
)
5 IMPLICIT DOUBLE PRECISION (A
-H
,O
-Z
)
6 PARAMETER (MXPT
=100, MXFRG
=50, MXFGPT
=MXPT*MXFRG
)
7 DIMENSION DM
(*),CHDINT
(L4
)
8 COMMON /FGRAD
/ DEF0
,DEFT0
,TORQ0
10 COMMON /CSSTV
/ CX
,CY
,CZ
11 * EFBTRM
(MXFGPT
),EFATRM2
(MXFGPT
),EFBTRM2
(MXFGPT
),
12 * EFDIP
(3,MXFGPT
),EFQAD
(6,MXFGPT
),
13 * EFOCT
(10,MXFGPT
),FRGNME
(MXFGPT
)
14 IF(NROOTS
.EQ
.5) CALL ROOT5
15 IF(NROOTS
.EQ
.6) CALL ROOT6
20 CHDINT
(ICC
)=CHDINT
(ICC
)-DUM*DUMY
23 CHDINT
(ICC
)=CHDINT
(ICC
)-DUM*DUMY
26 IF (LI
-LJ
) 920,940,940
30 960 NN
= (ID*
(ID
-1))/2+JD
32 ATORQ
(1,INF
)=ATORQ
(1,INF
)-DUM*
(CHDINT
(ICC
+1)*EFDIP
(3,IC
)
33 $
-CHDINT
(ICC
+2)*EFDIP
(2,IC
))