8 procedure,
nopass:: init
15 SUBROUTINE init (k,xi,wi)
36 integer,
parameter:: kmax = 10
38 real,
dimension(k):: wi, xi
40 real(kind=8),
dimension(55):: xn = (/ &
41 -1.00000000000000d0, &
42 -1.00000000000000d0, .33333333333333d0, &
43 -1.00000000000000d0,-.28989794855663d0, .68989794855663d0,&
44 -1.00000000000000d0,-.57531892352169d0, .18106627111853d0,&
46 -1.00000000000000d0,-.72048027131244d0,-.16718086473783d0,&
47 .44631397272376d0, .88579160777096d0, &
48 -1.00000000000000d0,-.80292982840235d0,-.39092854670727d0,&
49 .12405037950522d0, .60397316425279d0, .92038028589707d0,&
50 -1.00000000000000d0,-.85389134263949d0,-.53846772406011d0,&
51 -.11734303754309d0, .32603061943769d0, .70384280066303d0,&
53 -1.00000000000000d0,-.88747487892616d0,-.63951861652622d0,&
54 -.29475056577366d0, .09430725266111d0, .46842035443083d0,&
55 .77064189367819d0, .95504122712258d0, &
56 -1.00000000000000d0,-.91073208942006d0,-.71126748591571d0,&
57 -.42635048571114d0,-.09037336960685d0, .25613567083346d0,&
58 .57138304120873d0, .81735278420041d0, .96444016970528d0,&
59 -1.00000000000000d0,-.92748437423359d0,-.76384204242000d0,&
60 -.52564603037008d0,-.23623446939059d0, .07605919783798d0,&
61 .38066484014473d0, .64776668767401d0, .85122522058160d0,&
63 real(kind=8),
dimension(55):: wn = (/ &
65 0.50000000000000d0,1.50000000000005d0, &
66 0.22222222222222d0,1.02497165237682d0,0.75280612540100d0,&
67 0.12500000000000d0,0.65768863996010d0,0.77638693768637d0,&
69 0.08000000000000d0,0.44620780216715d0,0.62365304595145d0,&
70 0.56271203029887d0,0.28742712158247d0, &
71 0.05555555555556d0,0.31964075322051d0,0.48538718846895d0,&
72 0.52092678318961d0,0.41690133431186d0,0.20158838525329d0,&
73 0.04081632653061d0,0.23922748922532d0,0.38094987364422d0,&
74 0.44710982901453d0,0.42470377900601d0,0.31820423146727d0,&
76 0.03125000000000d0,0.18535815480298d0,0.30413062064680d0,&
77 0.37651754538912d0,0.39157216745253d0,0.34701479563444d0,&
78 0.24964790132996d0,0.11450881474417d0, &
79 0.02469135802469d0,0.14765401904631d0,0.24718937820459d0,&
80 0.31684377567046d0,0.34827300277294d0,0.33769396697589d0,&
81 0.28638669635730d0,0.20055329802460d0,0.09071450492309d0,&
82 0.02000000000000d0,0.12029667055749d0,0.20427013187899d0,&
83 0.26819483784117d0,0.30585928772442d0,0.31358245722693d0,&
84 0.29061016483286d0,0.23919343171435d0,0.16437601273700d0,&
86 integer,
dimension(kmax+1):: k0 = (/1,2,4,7,11,16,22,29,37,46,56/)
93 print *,
' WARNING(radaui): k=',k,
'< 1 does not make sense!' 94 print *,
' Setting k=1.' 98 print *,
' WARNING(radaui): k=',k,
'>',kmax,
' is not implemented.' 99 print *,
' Setting k=',kmax,
'.' 110 xi(i) = 1.0 - (c + d*xn(k0(k+1)-i))
111 wi(i) = d*wn(k0(k+1)-i)