Calculation of the local irregular solutions
VJLR
, VNL
, VHLR
, which appear in the integrands TAU(K,IPAN)
is used instead of TAU(K,IPAN)**2
, which directly gives RLL(r)
and SLL(r)
multiplied with r
. TAU
is the radial mesh.YR
, ZR
, YI
, ZI
because of the conventions used
by Gonzalez et al, Journal of Computational Physics 134, 134-149 (1997)
a factor is included in the source terms this factor is removed by
the definition of ZSLC1SUM
given below
i.e. prepare terms , appearing in 5.11, 5.12.Then determine the local solutions solve the equations and and and i.e., solve system , see eq. 5.68.
C a l c u l a t i o n
o f
t h e
l o c a l
i r r e g u l a r
s o l u t i o n s
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=dp) | :: | vll(lmsize*nvec,lmsize*nvec,nrmax) | ||||
real(kind=dp) | :: | tau(0:ncheb) | ||||
real(kind=dp) | :: | drpan2 | ||||
real(kind=dp) | :: | csrc1(0:ncheb,0:ncheb) | ||||
real(kind=dp) | :: | slc1sum(0:ncheb) | ||||
complex(kind=dp) | :: | mihvy(lmsize,lmsize) | ||||
complex(kind=dp) | :: | mihvz(lmsize,lmsize) | ||||
complex(kind=dp) | :: | mijvy(lmsize,lmsize) | ||||
complex(kind=dp) | :: | mijvz(lmsize,lmsize) | ||||
complex(kind=dp) | :: | yif(lmsize2,lmsize,0:ncheb) | ||||
complex(kind=dp) | :: | zif(lmsize2,lmsize,0:ncheb) | ||||
integer | :: | ncheb | ||||
integer | :: | ipan | ||||
integer | :: | lmsize | ||||
integer | :: | lmsize2 | ||||
integer | :: | nrmax | ||||
integer | :: | nvec | ||||
integer | :: | jlk_index(2*lmsize) | ||||
complex(kind=dp) | :: | hlk(lbessel,nrmax) | ||||
complex(kind=dp) | :: | jlk(lbessel,nrmax) | ||||
complex(kind=dp) | :: | hlk2(lbessel,nrmax) | ||||
complex(kind=dp) | :: | jlk2(lbessel,nrmax) | ||||
complex(kind=dp) | :: | gmatprefactor | ||||
character(len=1) | :: | cmodesll | ||||
integer | :: | lbessel | ||||
integer | :: | use_sratrick1 |
T h i s
s u b r o u t i n e
f a c i l i t a t e s
c o m p i l e
o p t i m i z a t i o n
b y
w o r k i n g
w i t t
o n l y
t w o - d i m e n s i o n a l
a r r a y s
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
complex(kind=dp) | :: | srv1(0:ncheb,lmsize) | ||||
complex(kind=dp) | :: | jlmkmn(0:ncheb,lmsize) | ||||
complex(kind=dp) | :: | hlmkmn(0:ncheb,lmsize) | ||||
complex(kind=dp) | :: | vhli(lmsize) | ||||
complex(kind=dp) | :: | vjli(lmsize) | ||||
integer | :: | ncheb | ||||
integer | :: | lmsize |