/[PAMELA software]/DarthVader/TrackerLevel2/src/F77/grkuta.f
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/F77/grkuta.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.5 by pam-fi, Fri Apr 27 11:19:51 2007 UTC revision 1.7 by pam-fi, Tue Jan 15 14:23:55 2008 UTC
# Line 28  C.    *                                 Line 28  C.    *                                
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,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)
# Line 37  C. Line 38  C.
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)
50          REAL*8 DLT32
51          DLT32=DLT/32.
52    
53  *.  *.
54  *.    ------------------------------------------------------------------  *.    ------------------------------------------------------------------
55  *.  *.
56  *             This constant is for units CM,GEV/C and KGAUSS  *             This constant is for units CM,GEV/C and KGAUSS
57  *  *
58    
59        ITER = 0        ITER = 0
60        NCUT = 0        NCUT = 0
61        DO 10 J=1,7        DO 10 J=1,7
# Line 69  C. Line 77  C.
77        DO I=1,3        DO I=1,3
78         F(I)=DBLE(FFF(I))         F(I)=DBLE(FFF(I))
79        ENDDO        ENDDO
80          DELTAB(2) = -F(2)*VECT(7)*CHARGE*(DELTA0+DELTA1*VVV(2))
81          F(2) = F(2)+DELTAB(2)
82    cPP   -----------------
83  *  *
84  *             Start of integration  *             Start of integration
85  *  *
# Line 106  C. Line 117  C.
117        CALL GUFLD(VVV,FFF)        CALL GUFLD(VVV,FFF)
118        DO I=1,3        DO I=1,3
119         F(I)=DBLE(FFF(I))         F(I)=DBLE(FFF(I))
120        ENDDO            ENDDO  
121          DELTAB(2) = -F(2)*VECT(7)*CHARGE*(DELTA0+DELTA1*VVV(2))
122          F(2) = F(2)+DELTAB(2)
123    cPP   -----------------
124  C      CALL GUFLD(XYZT,F)  C      CALL GUFLD(XYZT,F)
125        AT     = A + SECXS(1)        AT     = A + SECXS(1)
126        BT     = B + SECYS(1)        BT     = B + SECYS(1)
# Line 141  C      CALL GUFLD(XYZT,F) Line 155  C      CALL GUFLD(XYZT,F)
155        DO I=1,3        DO I=1,3
156         F(I)=DBLE(FFF(I))         F(I)=DBLE(FFF(I))
157        ENDDO        ENDDO
158          DELTAB(2) = -F(2)*VECT(7)*CHARGE*(DELTA0+DELTA1*VVV(2))
159          F(2) = F(2)+DELTAB(2)
160    cPP   -----------------
161  C      CALL GUFLD(XYZT,F)  C      CALL GUFLD(XYZT,F)
162  *  *
163        Z      = Z + (C + (SECZS(1) + SECZS(2) + SECZS(3)) * THIRD) * H        Z      = Z + (C + (SECZS(1) + SECZS(2) + SECZS(3)) * THIRD) * H
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.23