Lines Matching refs:b
49 integer aa, b local
79 do b = 1, nshl
80 t1(:,1) = uBar(:,1) * shg(:,b,1)
81 & + uBar(:,2) * shg(:,b,2)
82 & + uBar(:,3) * shg(:,b,3)
87 tmp1 = tsFct * shpfun(:,aa) * shpfun(:,b)
94 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp2
95 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp2
96 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp2
102 do b = 1, nshl
104 t1(:,1) = tauC * shg(:,b,1)
105 t1(:,2) = tauC * shg(:,b,2)
106 t1(:,3) = tauC * shg(:,b,3)
110 t2(:,1) = rmu * shg(:,b,1)
111 t2(:,2) = rmu * shg(:,b,2)
112 t2(:,3) = rmu * shg(:,b,3)
115 tmp1 = tauM * ( u1 * shg(:,b,1)
116 & + u2 * shg(:,b,2)
117 & + u3 * shg(:,b,3) )*rho
120 tmp2 = tauBar * ( rLui(:,1) * shg(:,b,1)
121 & + rLui(:,2) * shg(:,b,2)
122 & + rLui(:,3) * shg(:,b,3) )
135 aa = b
142 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp
145 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp
148 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp
154 xKebe(:,2,aa,b) = xKebe(:,2,aa,b) + tmp1
155 xKebe(:,4,b,aa) = xKebe(:,4,b,aa) + tmp1
159 xKebe(:,3,aa,b) = xKebe(:,3,aa,b) + tmp1
160 xKebe(:,7,b,aa) = xKebe(:,7,b,aa) + tmp1
164 xKebe(:,6,aa,b) = xKebe(:,6,aa,b) + tmp1
165 xKebe(:,8,b,aa) = xKebe(:,8,b,aa) + tmp1
169 do aa = b+1, nshl
177 xKebe(:,1,aa,b) = xKebe(:,1,aa,b) + tmp1
178 xKebe(:,1,b,aa) = xKebe(:,1,b,aa) + tmp1
183 xKebe(:,5,aa,b) = xKebe(:,5,aa,b) + tmp1
184 xKebe(:,5,b,aa) = xKebe(:,5,b,aa) + tmp1
189 xKebe(:,9,aa,b) = xKebe(:,9,aa,b) + tmp1
190 xKebe(:,9,b,aa) = xKebe(:,9,b,aa) + tmp1
196 xKebe(:,2,aa,b) = xKebe(:,2,aa,b) + tmp1
197 xKebe(:,4,b,aa) = xKebe(:,4,b,aa) + tmp1
201 xKebe(:,3,aa,b) = xKebe(:,3,aa,b) + tmp1
202 xKebe(:,7,b,aa) = xKebe(:,7,b,aa) + tmp1
206 xKebe(:,4,aa,b) = xKebe(:,4,aa,b) + tmp1
207 xKebe(:,2,b,aa) = xKebe(:,2,b,aa) + tmp1
211 xKebe(:,6,aa,b) = xKebe(:,6,aa,b) + tmp1
212 xKebe(:,8,b,aa) = xKebe(:,8,b,aa) + tmp1
216 xKebe(:,7,aa,b) = xKebe(:,7,aa,b) + tmp1
217 xKebe(:,3,b,aa) = xKebe(:,3,b,aa) + tmp1
221 xKebe(:,8,aa,b) = xKebe(:,8,aa,b) + tmp1
222 xKebe(:,6,b,aa) = xKebe(:,6,b,aa) + tmp1
230 do b = 1, nshl
231 t1(:,1) = tlW * shg(:,b,1)
232 t1(:,2) = tlW * shg(:,b,2)
233 t1(:,3) = tlW * shg(:,b,3)
235 xGoC(:,1,aa,b) = xGoC(:,1,aa,b) + t1(:,1) * shpfun(:,aa)
236 xGoC(:,2,aa,b) = xGoC(:,2,aa,b) + t1(:,2) * shpfun(:,aa)
237 xGoC(:,3,aa,b) = xGoC(:,3,aa,b) + t1(:,3) * shpfun(:,aa)
246 do b = 1, nshl
247 t1(:,1) = tauM * shg(:,b,1)
248 t1(:,2) = tauM * shg(:,b,2)
249 t1(:,3) = tauM * shg(:,b,3)
250 do aa = b, nshl
251 xGoC(:,4,aa,b) = xGoC(:,4,aa,b)
295 integer aa, b local
311 do b = 1, nshl
312 t1(:,1) = WdetJ * ( uMod(:,1) * shg(:,b,1)
313 & + uMod(:,2) * shg(:,b,2)
314 & + uMod(:,3) * shg(:,b,3) )
317 tmp1 = shpfun(:,aa) * shpfun(:,b)
319 xSebe(:,aa,b) = xSebe(:,aa,b) + tmp1 * (tsFct + srcL)
325 xSebe(:,b,aa) = xSebe(:,b,aa) + t1(:,2)*shpfun(:,aa)
331 do b = 1, nshl
333 & * ( uMod(:,1) * shg(:,b,1)
334 & + uMod(:,2) * shg(:,b,2)
335 & + uMod(:,3) * shg(:,b,3) )
337 t1(:,1) = kptmp * shg(:,b,1) + uMod(:,1) * tmp
338 t1(:,2) = kptmp * shg(:,b,2) + uMod(:,2) * tmp
339 t1(:,3) = kptmp * shg(:,b,3) + uMod(:,3) * tmp
353 t1(:,1) = t1(:,1) + giju(:,1) * shg(:,b,1)
354 2 + giju(:,4) * shg(:,b,2)
355 3 + giju(:,6) * shg(:,b,3)
356 t1(:,2) = t1(:,2) + giju(:,4) * shg(:,b,1)
357 2 + giju(:,2) * shg(:,b,2)
358 3 + giju(:,5) * shg(:,b,3)
359 t1(:,3) = t1(:,3) + giju(:,6) * shg(:,b,1)
360 2 + giju(:,5) * shg(:,b,2)
361 3 + giju(:,3) * shg(:,b,3)
367 aa = b
369 xSebe(:,aa,b) = xSebe(:,aa,b) + t1(:,1) * shg(:,aa,1)
376 do aa = b+1, nshl
381 xSebe(:,aa,b) = xSebe(:,aa,b) + tmp
382 xSebe(:,b,aa) = xSebe(:,b,aa) + tmp