Lines Matching refs:n

86        do n = 1, nshl 
87 pres = pres + shpfun(:,n) * yl(:,n,1)
88 u1 = u1 + shpfun(:,n) * yl(:,n,2)
89 u2 = u2 + shpfun(:,n) * yl(:,n,3)
90 u3 = u3 + shpfun(:,n) * yl(:,n,4)
94 do n = 1, nshl
95 Temp = Temp + shpfun(:,n) * yl(:,n,5)
101 do n = 1,nenl
102 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
103 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
104 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3)
110 do n = 1, nenl
111 dist2w = dist2w + shpfun(:,n) * dwl(:,n)
118 do n = 1, nshl
120 rlsli(:,1) = rlsli(:,1) + shpfun(:,n) * rlsl(:,n,1)
121 rlsli(:,2) = rlsli(:,2) + shpfun(:,n) * rlsl(:,n,2)
122 rlsli(:,3) = rlsli(:,3) + shpfun(:,n) * rlsl(:,n,3)
123 rlsli(:,4) = rlsli(:,4) + shpfun(:,n) * rlsl(:,n,4)
124 rlsli(:,5) = rlsli(:,5) + shpfun(:,n) * rlsl(:,n,5)
125 rlsli(:,6) = rlsli(:,6) + shpfun(:,n) * rlsl(:,n,6)
135 do n = 1, nshl
136 aci(:,1) = aci(:,1) + shpfun(:,n) * acl(:,n,2)
137 aci(:,2) = aci(:,2) + shpfun(:,n) * acl(:,n,3)
138 aci(:,3) = aci(:,3) + shpfun(:,n) * acl(:,n,4)
152 do n = 1, nshl
153 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1)
154 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
155 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
156 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4)
158 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * yl(:,n,1)
159 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2)
160 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3)
161 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4)
163 g3yi(:,1) = g3yi(:,1) + shg(:,n,3) * yl(:,n,1)
164 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2)
165 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3)
166 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4)
176 do n=1, nshl
177 divqi(:,1) = divqi(:,1) + shg(:,n,1)*ql(:,n,1 )
178 & + shg(:,n,2)*ql(:,n,4 )
179 & + shg(:,n,3)*ql(:,n,7 )
181 divqi(:,2) = divqi(:,2) + shg(:,n,1)*ql(:,n,2 )
182 & + shg(:,n,2)*ql(:,n,5 )
183 & + shg(:,n,3)*ql(:,n,8)
185 divqi(:,3) = divqi(:,3) + shg(:,n,1)*ql(:,n,3 )
186 & + shg(:,n,2)*ql(:,n,6 )
187 & + shg(:,n,3)*ql(:,n,9 )
195 do n=1, nshl
197 & + shg(:,n,1)*ql(:,n,idflx-2)
198 & + shg(:,n,2)*ql(:,n,idflx-1)
199 & + shg(:,n,3)*ql(:,n,idflx)
207 do n = 1, nshl
208 Sclr(i) = Sclr(i) + shpfun(i,n) * yl(i,n,6) !scalar
212 gyti(i,1) = gyti(i,1) + shg(i,n,1) * yl(i,n,6)
213 gyti(i,2) = gyti(i,2) + shg(i,n,2) * yl(i,n,6)
214 gyti(i,3) = gyti(i,3) + shg(i,n,3) * yl(i,n,6)
291 do n = 1, nenl
292 dkei = dkei + shpfun(:,n) * (1.0-xl(:,n,2)**2) !u_ex^~ (in FEM space)
360 do n = 1, nshl
361 u1 = u1 + shpfun(:,n) * yl(:,n,2)
362 u2 = u2 + shpfun(:,n) * yl(:,n,3)
363 u3 = u3 + shpfun(:,n) * yl(:,n,4)
364 Sclr = Sclr + shpfun(:,n) * yl(:,n,id)
371 do n = 1, nshl
372 Sdot = Sdot + shpfun(:,n) * acl(:,n,id)
386 do n = 1, nshl
387 gradS(:,1) = gradS(:,1) + shg(:,n,1) * yl(:,n,id)
388 gradS(:,2) = gradS(:,2) + shg(:,n,2) * yl(:,n,id)
389 gradS(:,3) = gradS(:,3) + shg(:,n,3) * yl(:,n,id)
397 do n=1, nshl
398 divS(:) = divS(:) + shg(:,n,1)*ql(:,n,1 )
399 & + shg(:,n,2)*ql(:,n,2 )
400 & + shg(:,n,3)*ql(:,n,3 )
410 do n = 1, nshl
411 Temp = Temp + shpfun(:,n) * yl(:,n,5)
417 do n = 1,nenl
418 xx(:,1) = xx(:,1) + shpfun(:,n) * xl(:,n,1)
419 xx(:,2) = xx(:,2) + shpfun(:,n) * xl(:,n,2)
420 xx(:,3) = xx(:,3) + shpfun(:,n) * xl(:,n,3)
424 do n = 1, nshl
425 aci(:,1) = aci(:,1) + shpfun(:,n) * acl(:,n,2)
426 aci(:,2) = aci(:,2) + shpfun(:,n) * acl(:,n,3)
427 aci(:,3) = aci(:,3) + shpfun(:,n) * acl(:,n,4)
432 do n = 1, nshl
433 g1yi(:,1) = g1yi(:,1) + shg(:,n,1) * yl(:,n,1)
434 g1yi(:,2) = g1yi(:,2) + shg(:,n,1) * yl(:,n,2)
435 g1yi(:,3) = g1yi(:,3) + shg(:,n,1) * yl(:,n,3)
436 g1yi(:,4) = g1yi(:,4) + shg(:,n,1) * yl(:,n,4)
438 g2yi(:,1) = g2yi(:,1) + shg(:,n,2) * yl(:,n,1)
439 g2yi(:,2) = g2yi(:,2) + shg(:,n,2) * yl(:,n,2)
440 g2yi(:,3) = g2yi(:,3) + shg(:,n,2) * yl(:,n,3)
441 g2yi(:,4) = g2yi(:,4) + shg(:,n,2) * yl(:,n,4)
443 g3yi(:,1) = g3yi(:,1) + shg(:,n,3) * yl(:,n,1)
444 g3yi(:,2) = g3yi(:,2) + shg(:,n,3) * yl(:,n,2)
445 g3yi(:,3) = g3yi(:,3) + shg(:,n,3) * yl(:,n,3)
446 g3yi(:,4) = g3yi(:,4) + shg(:,n,3) * yl(:,n,4)