Lines Matching refs:i
55 do i = 1, nelint
57 yintl(i,:,j+1) = y(ien2D(i,:),j)
59 yintl(i,:,1) = y(ien2D(i,:),4)
62 yintl(i,:,j) = y(ien2D(i,:),j)
72 do i=1,nfint
73 velbarR(i,:)=0
74 do j=1,imax(i)+1
75 call shptet(ipord,xsinfin(i,j,:),shapef(:),shgradl(:,:))
77 velbarR(i,:)=velbarR(i,:)
78 & + yintl(elcnfin(i,j),k,:)*shapef(k)
81 velbarR(i,:)=velbarR(i,:) / (imax(i)+1)
95 do i=2,nfint+1
96 v1(i)=velbarR(i-1,dir)-0.99*vel
97 if((v1(i).gt.0).and.(v1(i-1).le.0)) then
98 label=i-1
102 label=i-1
145 do i=1, nfint
146 if (xyn(i).ge.deltaint) then
147 freestream(:) = freestream(:) + velbarR(i,:)
157 do i=1, nfint
158 if (xyn(i).ge.deltaint) then
159 velbarR(i,:) = freestream(:)
167 do i=1,nelint
171 if ((x(ien2D(i,k),2).ge.xyn(j)) .and.
172 & (x(ien2D(i,k),2).le.(xyn(j+1)+0.000001))) then
173 tmp = (x(ien2D(i,k),2) - xyn(j)) /
176 velbarl(i,k,l) =
182 if ((xcyl(ien2D(i,k),1).ge.xcyl(nrint(j+1),1)) .and.
183 & (xcyl(ien2D(i,k),1).le.xcyl(nrint(j),1))) then
184 tmp = (xcyl(ien2D(i,k),1) - xcyl(nrint(j+1),1)) /
187 velbarl(i,k,l) =
219 do i=1,npin
226 xts1 = x(nen1(i),1) + plandist
227 if (xynin(i)*rdelta.gt.ymax) then
230 xts2 = xynin(i)*rdelta
232 xts3 = x(nen1(i),3)
239 if (xynin(i).le.0.00001) then
240 xts1 = (radcyl-xynin(i)*rdelta*sang-tolerence)
241 & *cos(xcyl(nen1(i),2))
242 xts2 = (radcyl-xynin(i)*rdelta*sang-tolerence)
243 & *sin(xcyl(nen1(i),2))
244 xts3 = (aR-(radcyl-xynin(i)*rdelta*sang-tolerence)
245 & * (xnrml*cos(xcyl(nen1(i),2))
246 & + ynrml*sin(xcyl(nen1(i),2))))/znrml
247 elseif (xynin(i)*rdelta.gt.ymax) then
249 & *cos(xcyl(nen1(i),2))
251 & *sin(xcyl(nen1(i),2))
253 & * (xnrml*cos(xcyl(nen1(i),2))
254 & + ynrml*sin(xcyl(nen1(i),2))))/znrml
256 xts1 = (radcyl-xynin(i)*rdelta*sang)
257 & *cos(xcyl(nen1(i),2))
258 xts2 = (radcyl-xynin(i)*rdelta*sang)
259 & *sin(xcyl(nen1(i),2))
260 xts3 = (aR-(radcyl-xynin(i)*rdelta*sang)
261 & * (xnrml*cos(xcyl(nen1(i),2))
262 & + ynrml*sin(xcyl(nen1(i),2))))/znrml
280 fluc2(i,:)= 0 !fluc2(i,:) + flucl(element,k,:)*shapef(k)
281 ubar2(i,:)=ubar2(i,:) + ubarintl(element,k,:)*shapef(k)
322 do i = 1, npin
323 if (xynin(i).ge.rbltin) then
324 ubar2(i,:) = freestream(:)
326 fluc2(i,:) = 0