Lines Matching refs:i
49 do i=1,numnp
50 if ((x(i,1).eq.0).and.(x(i,2).eq.0)) then
52 nula(j) = i
54 xcyl(i,2) = atan2(x(i,2),x(i,1))
66 do i=2,npin
67 if((x(nen1(i),1).eq.0).and.(x(nen1(i),2).eq.0)) then
70 if(thmin.gt.xcyl(nen1(i),2)) thmin = xcyl(nen1(i),2)
71 if(thmax.lt.xcyl(nen1(i),2)) thmax = xcyl(nen1(i),2)
76 do i=1,j
77 xcyl(nula(i),2) = thmin
86 do i=1,npin
87 error1 = abs(xcyl(nen1(i),2) - thmin)
90 nrin(j)=nen1(i)
95 do i=1,nfin
96 xrtmp(i) = xcyl(nrin(i),1)
106 do i=1,nfin
107 if (radinl.lt.xcyl(nrin(i),1)) radinl = xcyl(nrin(i),1)
129 do i=1,npin
130 xynin(i) = (radinl - xcyl(nen1(i),1))/sang
137 do i=1,nfint
138 xyn(i) = (radcyl - xcyl(nrint(i),1))/sang
160 do i=1,imax(jj)
163 & *cos(1.0*i/imax(jj)*(thmax-thmin)+thmin)
165 & *sin(1.0*i/imax(jj)*(thmax-thmin)+thmin)
169 & *cos(1.0*i/imax(jj)*(thmax-thmin)+thmin)
171 & *sin(1.0*i/imax(jj)*(thmax-thmin)+thmin)
176 xsinfin(jj,1+i,:) = xtmp(:)
177 elcnfin(jj,1+i) = etmp
225 do i=2,npin
226 if(x(nen1(i),3).lt.zmin) zmin = x(nen1(i),3)
227 if(x(nen1(i),3).gt.zmax) zmax = x(nen1(i),3)
230 do i=1,npin
231 if (x(nen1(i),3).eq.zmin) then
233 nyin(j) = nen1(i)
237 do i=1,nfin
238 yrtmp(i) = x(nyin(i),2)
247 do i=1,nfin
265 do i=1,npin
266 xynin(i) = x(nen1(i),2)
273 do i=1,nfint
274 xyn(i) = x(nrint(i),2)
296 do i=1,imax(jj)
299 xts3 = x(nrint(jj),3) + i*ds
302 xsinfin(jj,1+i,:) = xtmp(:)
303 elcnfin(jj,1+i) = etmp