Line data Source code
1 : MODULE m_crtail
2 : c
3 : c initial point for income regular solution of dirac eq.
4 : c order kap1=-L-1, kap2=L
5 : c
6 : CONTAINS
7 332 : SUBROUTINE crtail(
8 332 : > mrad,e,rc,nsol,nzero,csq,
9 332 : X gc,fc)
10 :
11 : IMPLICIT NONE
12 : C ..
13 : C .. Scalar Arguments ..
14 : INTEGER, INTENT (IN) :: mrad,nsol,nzero
15 : REAL , INTENT (IN) :: e,csq
16 : C ..
17 : C .. Array Arguments ..
18 : REAL , INTENT (IN) :: rc(mrad)
19 : REAL , INTENT (INOUT) :: fc(2,2,mrad),gc(2,2,mrad)
20 : C ..
21 : C .. Local Scalars ..
22 : REAL beta,rr
23 : INTEGER i,ir,j
24 : C ..
25 : C .. Intrinsic Functions ..
26 : INTRINSIC exp,sqrt
27 : c
28 : C----->infinity point solutions construction
29 : c
30 332 : beta = sqrt(-e-e*e/csq)
31 834 : DO j = 1,nsol
32 1676 : DO i = 1,nsol
33 78256 : DO ir = nzero + 1,mrad
34 76912 : rr = rc(ir) - rc(nzero)
35 76912 : gc(i,j,ir) = gc(i,j,nzero)*exp(-beta*rr)
36 77754 : fc(i,j,ir) = fc(i,j,nzero)*exp(-beta*rr)
37 : END DO
38 : END DO
39 : END DO
40 :
41 332 : END SUBROUTINE crtail
42 : END MODULE m_crtail
|