28 |
C. ****************************************************************** |
C. ****************************************************************** |
29 |
C. |
C. |
30 |
IMPLICIT DOUBLE PRECISION(A-H,O-Z) |
IMPLICIT DOUBLE PRECISION(A-H,O-Z) |
31 |
COMMON/DELTAB/DELTA0,DELTA1 |
COMMON/DELTAB/DELTA0,DELTA1,DLT |
32 |
* |
* |
33 |
REAL VVV(3),FFF(3) |
REAL VVV(3),FFF(3) |
34 |
REAL*8 CHARGE, STEP, VECT(*), VOUT(*), F(4) |
REAL*8 CHARGE, STEP, VECT(*), VOUT(*), F(4) |
38 |
+ (XT,XYZT(1)),(YT,XYZT(2)),(ZT,XYZT(3)) |
+ (XT,XYZT(1)),(YT,XYZT(2)),(ZT,XYZT(3)) |
39 |
* |
* |
40 |
PARAMETER (MAXIT = 1992, MAXCUT = 11) |
PARAMETER (MAXIT = 1992, MAXCUT = 11) |
41 |
PARAMETER (EC=2.9979251D-4,DLT=1D-4,DLT32=DLT/32) |
cPP PARAMETER (EC=2.9979251D-4,DLT=1D-4,DLT32=DLT/32) |
42 |
PARAMETER (ZERO=0, ONE=1, TWO=2, THREE=3) |
PARAMETER (EC=2.99792458D-4) |
43 |
|
cPP PARAMETER (ZERO=0, ONE=1, TWO=2, THREE=3) |
44 |
|
PARAMETER (ZERO=0.D0, ONE=1.D0, TWO=2.D0, THREE=3.D0) |
45 |
PARAMETER (THIRD=ONE/THREE, HALF=ONE/TWO) |
PARAMETER (THIRD=ONE/THREE, HALF=ONE/TWO) |
46 |
PARAMETER (PISQUA=.986960440109D+01) |
PARAMETER (PISQUA=.986960440109D+01) |
47 |
PARAMETER (IX=1,IY=2,IZ=3,IPX=4,IPY=5,IPZ=6) |
PARAMETER (IX=1,IY=2,IZ=3,IPX=4,IPY=5,IPZ=6) |
48 |
|
|
49 |
REAL*8 DELTAB(3) |
REAL*8 DELTAB(3) |
50 |
|
REAL*8 DLT32 |
51 |
|
DLT32=DLT/32. |
52 |
|
|
53 |
*. |
*. |
54 |
*. ------------------------------------------------------------------ |
*. ------------------------------------------------------------------ |