94
95
96
97#include "implicit_f.inc"
98
99
100
101#include "scr05_c.inc"
102
103
104
105 INTEGER NB, LEN, NSVSI(*), IXLINS(2,*)
107 . bufr(len,*), alphak(3,*)
108 DOUBLE PRECISION DAANC6(3,6,*)
109
110
111
112 INTEGER I, N, NL1, NL2, INC ,ISIGN
113 DOUBLE PRECISION DAANC6L(3,6)
114
115
116
117
118 DO i = 1, nb
119 n = nsvsi(i)
120 nl1 = ixlins(1,n)
121 nl2 = ixlins(2,n)
122 isign = 1
123 IF(alphak(2,nl1) < zero .or. bufr(1,i) < zero)isign = -1
124 alphak(2,nl1) = isign*
min(abs(alphak(2,nl1)),abs(bufr(1,i)))
125 isign = 1
126 IF(alphak(2,nl2) < zero .or. bufr(2,i) < zero)isign = -1
127 alphak(2,nl2) = isign*
min(abs(alphak(2,nl2)),abs(bufr(2,i)))
129 daanc6(1,1,nl1) = daanc6(1,1,nl1)+daanc6l(1,1)
130 daanc6(1,2,nl1) = daanc6(1,2,nl1)+daanc6l(1,2)
131 daanc6(1,3,nl1) = daanc6(1,3,nl1)+daanc6l(1,3)
132 daanc6(1,4,nl1) = daanc6(1,4,nl1)+daanc6l(1,4)
133 daanc6(1,5,nl1) = daanc6(1,5,nl1)+daanc6l(1,5)
134 daanc6(1,6,nl1) = daanc6(1,6,nl1)+daanc6l(1,6)
135 daanc6(2,1,nl1) = daanc6(2,1,nl1)+daanc6l(2,1)
136 daanc6(2,2,nl1) = daanc6(2,2,nl1)+daanc6l(2,2)
137 daanc6(2,3,nl1) = daanc6(2,3,nl1)+daanc6l(2,3)
138 daanc6(2,4,nl1) = daanc6(2,4,nl1)+daanc6l(2,4)
139 daanc6(2,5,nl1) = daanc6(2,5,nl1)+daanc6l(2,5)
140 daanc6(2,6,nl1) = daanc6(2,6,nl1)+daanc6l(2,6)
141 daanc6(3,1,nl1) = daanc6(3,1,nl1)+daanc6l(3,1)
142 daanc6(3,2,nl1) = daanc6(3,2,nl1)+daanc6l(3,2)
143 daanc6(3,3,nl1) = daanc6(3,3,nl1)+daanc6l(3,3)
144 daanc6(3,4,nl1) = daanc6(3,4,nl1)+daanc6l(3,4)
145 daanc6(3,5,nl1) = daanc6(3,5,nl1)+daanc6l(3,5)
146 daanc6(3,6,nl1) = daanc6(3,6,nl1)+daanc6l(3,6)
147 CALL getdpdaanc(daanc6l,bufr(3+inc,i),iresp,inc)
148 daanc6(1,1,nl2) = daanc6(1,1,nl2)+daanc6l(1,1)
149 daanc6(1,2,nl2) = daanc6(1,2,nl2)+daanc6l(1,2)
150 daanc6(1,3,nl2) = daanc6(1,3,nl2)+daanc6l(1,3)
151 daanc6(1,4,nl2) = daanc6(1,4,nl2)+daanc6l(1,4)
152 daanc6(1,5,nl2) = daanc6(1,5,nl2)+daanc6l(1,5)
153 daanc6(1,6,nl2) = daanc6(1,6,nl2)+daanc6l(1,6)
154 daanc6(2,1,nl2) = daanc6(2,1,nl2)+daanc6l(2,1)
155 daanc6(2,2,nl2) = daanc6(2,2,nl2)+daanc6l(2,2)
156 daanc6(2,3,nl2) = daanc6(2,3,nl2)+daanc6l(2,3)
157 daanc6(2,4,nl2) = daanc6(2,4,nl2)+daanc6l(2,4)
158 daanc6(2,5,nl2) = daanc6(2,5,nl2)+daanc6l(2,5)
159 daanc6(2,6,nl2) = daanc6(2,6,nl2)+daanc6l(2,6)
160 daanc6(3,1,nl2) = daanc6(3,1,nl2)+daanc6l(3,1)
161 daanc6(3,2,nl2) = daanc6(3,2,nl2)+daanc6l(3,2)
162 daanc6(3,3,nl2) = daanc6(3,3,nl2)+daanc6l(3,3)
163 daanc6(3,4,nl2) = daanc6(3,4,nl2)+daanc6l(3,4)
164 daanc6(3,5,nl2) = daanc6(3,5,nl2)+daanc6l(3,5)
165 daanc6(3,6,nl2) = daanc6(3,6,nl2)+daanc6l(3,6)
166 END DO
167
168 RETURN