3 PARAMETER(MAXDIM
=4,MAXVEC
=4**3*8,MAXT
=20)
4 REAL(KIND
=8) B
,W1
,W2
,BNORM
,BINV
,WT
,W0
,C1
,C2
,R1
,R2
5 DIMENSION B
(MAXVEC
,0:3),W1
(MAXVEC
,0:3),W2
(MAXVEC
,0:3)
6 DIMENSION BNORM
(MAXVEC
),BINV
(MAXVEC
),WT
(MAXVEC
),W0
(MAXVEC
)
7 DIMENSION C1
(MAXVEC
),C2
(MAXVEC
),R1
(MAXVEC
),R2
(MAXVEC
)
11 CALL INTACT
(ILAT
,I1
,I1
,W1
)
15 WT
(I
)=ONE
+ C1
(I
)*LOG
(EPS
+R1
(I
))
16 IF( R2
(I
)**2 .LE
. (ONE
-WT
(I
)**2) )THEN