Lines Matching refs:y
1 subroutine itrBC (y,ac, iBC, BC, iper, ilwork) argument
25 dimension y(nshg,nflow), iBC(nshg),
40 locmax=maxloc(y(:,i))
41 locmin=minloc(y(:,i))
42 ymax=maxval(y(:,i))
43 ymin=minval(y(:,i))
47 ypc=y(j,i)
48 y(j,i)=min(ylimit(3,i),max(ylimit(2,i),y(j,i)))
49 if(ypc.ne.y(j,i) )limitcount(i)=limitcount(i)+one
63 y(:,5) = BC(:,2)
79 y(:,1) = BC(:,3) - BC(:,4) * y(:,2)
80 & - BC(:,5) * y(:,3)
86 y(:,2) = BC(:,3) - BC(:,4) * y(:,1)
87 & - BC(:,5) * y(:,3)
93 y(:,1) = BC(:,3) - BC(:,4) * y(:,3)
94 y(:,2) = BC(:,5) - BC(:,6) * y(:,3)
100 y(:,3) = BC(:,3) - BC(:,4) * y(:,1)
101 & - BC(:,5) * y(:,2)
107 y(:,1) = BC(:,3) - BC(:,4) * y(:,2)
108 y(:,3) = BC(:,5) - BC(:,6) * y(:,2)
114 y(:,2) = BC(:,3) - BC(:,4) * y(:,1)
115 y(:,3) = BC(:,5) - BC(:,6) * y(:,1)
121 y(:,1) = BC(:,3)
122 y(:,2) = BC(:,4)
123 y(:,3) = BC(:,5)
148 call getthm (y1, y(:,5), tmp,
155 y(:,1) = y1
168 y(:,4) = BC(:,1) ! pressure here
176 y(:,i) = y(iper(:),i)
183 call commu (y, ilwork, nflow, 'out')
190 call rotabc(y, iBC, 'out')
201 subroutine itrBCSclr (y, ac, iBC, BC, iper, ilwork) argument
211 dimension y(nshg,ndof), iBC(nshg),
226 & y(:,id)=min(ylimit(3,id),max(ylimit(2,id),y(:,id)))
232 y(:,id) = BC(:,ibb)
239 y(i,id) = y(iper(i),id)
247 T=y(:,id)
249 y(:,id)=T