Lines Matching refs:iel
139 do iel = 1, npro
144 iface = abs(iBCB(iel,2))
154 flxID(1,iface) = flxID(1,iface) + WdetJb(iel)! measure area too
155 flxID(2,iface) = flxID(2,iface) - WdetJb(iel) * unm(iel)
157 & - ( tau1n(iel) - bnorm(iel,1)*pres(iel))
158 & * WdetJb(iel)
160 & - ( tau2n(iel) - bnorm(iel,2)*pres(iel))
161 & * WdetJb(iel)
163 & - ( tau3n(iel) - bnorm(iel,3)*pres(iel))
164 & * WdetJb(iel)
172 if (btest(iBCB(iel,1),0)) then
173 unm(iel) = zero
176 unm(iel) = unm(iel)
177 & + shape(iel,nodlcl) * BCB(iel,n,1)
183 if (btest(iBCB(iel,1),1)) then
184 pres(iel) = zero
187 pres(iel) = pres(iel)
188 & + shape(iel,nodlcl) * BCB(iel,n,2)
194 if (btest(iBCB(iel,1),2)) then
195 tau1n(iel) = zero
196 tau2n(iel) = zero
197 tau3n(iel) = zero
200 tau1n(iel) = tau1n(iel)
201 & + shape(iel,nodlcl)*BCB(iel,n,3)
202 tau2n(iel) = tau2n(iel)
203 & + shape(iel,nodlcl)*BCB(iel,n,4)
204 tau3n(iel) = tau3n(iel)
205 & + shape(iel,nodlcl)*BCB(iel,n,5)
212 if (btest(iBCB(iel,1),4)) then
213 rlKwall(iel,:,:) = rlKwall(iel,:,:) / ngaussb ! divide by number of gauss points
214 pres(iel) = zero ! to avoid the gauss point loop
215 tau1n(iel) = zero ! and make the traction contribution
216 tau2n(iel) = zero ! zero
217 tau3n(iel) = zero
219 rlKwall(iel,:,:) = zero ! this is not a deformable element
220 vdot(iel,:) = zero
221 xKebe(iel,:,:,:) = zero
222 rKwall_glob(iel,:,:,:) = zero ! no stiffness: not a wall element
235 do iel = 1, npro
236 if (nsrflist(iBCB(iel,2)) .ne. 0) then
237 unm(iel) = zero
238 tau1n(iel) = zero
239 tau2n(iel) = zero
240 tau3n(iel) = zero
255 vn=dot_product(vdot(iel,:),bnorm(iel,:))
256 vdot(iel,:)=vn*bnorm(iel,:)
257 walln1=dot_product(rlkwall(iel,1,:),bnorm(iel,:))
258 walln2=dot_product(rlkwall(iel,2,:),bnorm(iel,:))
259 walln3=dot_product(rlkwall(iel,3,:),bnorm(iel,:))
260 rlkwall(iel,1,:)=walln1*bnorm(iel,:)
261 rlkwall(iel,2,:)=walln2*bnorm(iel,:)
262 rlkwall(iel,3,:)=walln3*bnorm(iel,:)
429 do iel=1, npro
434 if (nsrflist(iBCB(iel,2)).ne.0 .and. ires.ne.2) then
435 iface = abs(iBCB(iel,2))
437 & - WdetJb(iel) * flux(iel)
440 if (btest(iBCB(iel,1),ib-1)) then
441 flux(iel) = zero
444 flux(iel) = flux(iel)
445 & + shape(iel,nodlcl) * BCB(iel,n,ibb)