Lines Matching refs:BCtmp
1 subroutine genBC1 (BCtmp, iBC, BC) argument
26 dimension BCtmp(nshg,ndof+7), iBC(nshg),
34 where (btest(iBC,5+isclr)) BC(:,6+isclr) = BCtmp(:,12+isclr)
39 where (btest(iBC,0)) BC(:,1) = BCtmp(:,1) ! density
40 where (btest(iBC,1)) BC(:,2) = BCtmp(:,2) ! temperature
41 where (btest(iBC,2)) BC(:,1) = BCtmp(:,3) ! pressure
46 tmp = BCtmp(:,4)**2 + BCtmp(:,5)**2 + BCtmp(:,6)**2
47 BC(:,3) = tmp * BCtmp(:,7) / BCtmp(:,4)
48 BC(:,4) = BCtmp(:,5) / BCtmp(:,4)
49 BC(:,5) = BCtmp(:,6) / BCtmp(:,4)
55 tmp = BCtmp(:,4)**2 + BCtmp(:,5)**2 + BCtmp(:,6)**2
56 BC(:,3) = tmp * BCtmp(:,7) / BCtmp(:,5)
57 BC(:,4) = BCtmp(:,4) / BCtmp(:,5)
58 BC(:,5) = BCtmp(:,6) / BCtmp(:,5)
69 & (BCtmp(i,4).eq.0 .or. BCtmp(i,9).eq.0)) then !flip them
70 tmpbc(1:4)=BCtmp(i,4:7)
71 BCtmp(i,4:7)=BCtmp(i,8:11)
72 BCtmp(i,8:11)=tmpbc(1:4)
76 tmp = sqrt (BCtmp(:, 4)**2 + BCtmp(:, 5)**2
77 & + BCtmp(:, 6)**2)
78 BCtmp(:, 4) = BCtmp(:, 4) / tmp
79 BCtmp(:, 5) = BCtmp(:, 5) / tmp
80 BCtmp(:, 6) = BCtmp(:, 6) / tmp
81 BCtmp(:, 7) = BCtmp(:, 7) * tmp
83 tmp = sqrt (BCtmp(:, 8)**2 + BCtmp(:, 9)**2
84 & + BCtmp(:,10)**2)
85 BCtmp(:, 8) = BCtmp(:, 8) / tmp
86 BCtmp(:, 9) = BCtmp(:, 9) / tmp
87 BCtmp(:,10) = BCtmp(:,10) / tmp
88 BCtmp(:,11) = BCtmp(:,11) * tmp
90 BCtmp(:, 4) = BCtmp(:, 9) * BCtmp(:, 4)
91 & - BCtmp(:, 5) * BCtmp(:, 8)
92 BCtmp(:, 6) = BCtmp(:, 9) * BCtmp(:, 6)
93 & - BCtmp(:, 5) * BCtmp(:,10)
94 BCtmp(:, 7) = BCtmp(:, 9) * BCtmp(:, 7)
95 & - BCtmp(:, 5) * BCtmp(:,11)
96 BC(:,3) = BCtmp(:, 7) / BCtmp(:, 4)
97 BC(:,4) = BCtmp(:, 6) / BCtmp(:, 4)
99 BCtmp(:, 9) = BCtmp(:, 4) * BCtmp(:, 9)
100 BCtmp(:,10) = BCtmp(:, 4) * BCtmp(:,10)
101 & - BCtmp(:, 8) * BCtmp(:, 6)
102 BCtmp(:,11) = BCtmp(:, 4) * BCtmp(:,11)
103 & - BCtmp(:, 8) * BCtmp(:, 7)
104 BC(:,5) = BCtmp(:,11) / BCtmp(:, 9)
105 BC(:,6) = BCtmp(:,10) / BCtmp(:, 9)
112 tmp = BCtmp(:,4)**2 + BCtmp(:,5)**2 + BCtmp(:,6)**2
113 BC(:,3) = tmp * BCtmp(:,7) / BCtmp(:,6)
114 BC(:,4) = BCtmp(:,4) / BCtmp(:,6)
115 BC(:,5) = BCtmp(:,5) / BCtmp(:,6)
128 & (BCtmp(i,4).eq.0 .or. BCtmp(i,10).eq.0)) then !flip them
129 tmpbc(1:4)=BCtmp(i,4:7)
130 BCtmp(i,4:7)=BCtmp(i,8:11)
131 BCtmp(i,8:11)=tmpbc(1:4)
135 tmp = sqrt (BCtmp(:, 4)**2 + BCtmp(:, 5)**2
136 & + BCtmp(:, 6)**2)
137 BCtmp(:, 4) = BCtmp(:, 4) / tmp
138 BCtmp(:, 5) = BCtmp(:, 5) / tmp
139 BCtmp(:, 6) = BCtmp(:, 6) / tmp
140 BCtmp(:, 7) = BCtmp(:, 7) * tmp
142 tmp = sqrt (BCtmp(:, 8)**2 + BCtmp(:, 9)**2
143 & + BCtmp(:,10)**2)
144 BCtmp(:, 8) = BCtmp(:, 8) / tmp
145 BCtmp(:, 9) = BCtmp(:, 9) / tmp
146 BCtmp(:,10) = BCtmp(:,10) / tmp
147 BCtmp(:,11) = BCtmp(:,11) * tmp
149 BCtmp(:, 4) = BCtmp(:,10) * BCtmp(:, 4)
150 & - BCtmp(:, 6) * BCtmp(:, 8)
151 BCtmp(:, 5) = BCtmp(:,10) * BCtmp(:, 5)
152 & - BCtmp(:, 6) * BCtmp(:, 9)
153 BCtmp(:, 7) = BCtmp(:,10) * BCtmp(:, 7)
154 & - BCtmp(:, 6) * BCtmp(:,11)
155 BC(:,3) = BCtmp(:, 7) / BCtmp(:, 4)
156 BC(:,4) = BCtmp(:, 5) / BCtmp(:, 4)
158 BCtmp(:, 9) = BCtmp(:, 4) * BCtmp(:, 9)
159 & - BCtmp(:, 8) * BCtmp(:, 5)
160 BCtmp(:,10) = BCtmp(:, 4) * BCtmp(:,10)
161 BCtmp(:,11) = BCtmp(:, 4) * BCtmp(:,11)
162 & - BCtmp(:, 8) * BCtmp(:, 7)
163 BC(:,5) = BCtmp(:,11) / BCtmp(:,10)
164 BC(:,6) = BCtmp(:, 9) / BCtmp(:,10)
177 & BCtmp(i,5).eq.0 .or. BCtmp(i,10).eq.0)) then !flip them
178 tmpbc(1:4)=BCtmp(i,4:7)
179 BCtmp(i,4:7)=BCtmp(i,8:11)
180 BCtmp(i,8:11)=tmpbc(1:4)
184 tmp = sqrt (BCtmp(:, 4)**2 + BCtmp(:, 5)**2
185 & + BCtmp(:, 6)**2)
186 BCtmp(:, 4) = BCtmp(:, 4) / tmp
187 BCtmp(:, 5) = BCtmp(:, 5) / tmp
188 BCtmp(:, 6) = BCtmp(:, 6) / tmp
189 BCtmp(:, 7) = BCtmp(:, 7) * tmp
191 tmp = sqrt (BCtmp(:, 8)**2 + BCtmp(:, 9)**2
192 & + BCtmp(:,10)**2)
193 BCtmp(:, 8) = BCtmp(:, 8) / tmp
194 BCtmp(:, 9) = BCtmp(:, 9) / tmp
195 BCtmp(:,10) = BCtmp(:,10) / tmp
196 BCtmp(:,11) = BCtmp(:,11) * tmp
198 BCtmp(:, 4) = BCtmp(:,10) * BCtmp(:, 4)
199 & - BCtmp(:, 6) * BCtmp(:, 8)
200 BCtmp(:, 5) = BCtmp(:,10) * BCtmp(:, 5)
201 & - BCtmp(:, 6) * BCtmp(:, 9)
202 BCtmp(:, 7) = BCtmp(:,10) * BCtmp(:, 7)
203 & - BCtmp(:, 6) * BCtmp(:,11)
204 BC(:,3) = BCtmp(:, 7) / BCtmp(:, 5)
205 BC(:,4) = BCtmp(:, 4) / BCtmp(:, 5)
207 BCtmp(:, 8) = BCtmp(:, 5) * BCtmp(:, 8)
208 & - BCtmp(:, 9) * BCtmp(:, 4)
209 BCtmp(:,10) = BCtmp(:, 5) * BCtmp(:,10)
210 BCtmp(:,11) = BCtmp(:, 5) * BCtmp(:,11)
211 & - BCtmp(:, 9) * BCtmp(:, 7)
212 BC(:,5) = BCtmp(:,11) / BCtmp(:,10)
213 BC(:,6) = BCtmp(:, 8) / BCtmp(:,10)
221 BC(:,3) = BCtmp(:,7) * BCtmp(:,4)
222 BC(:,4) = BCtmp(:,7) * BCtmp(:,5)
223 BC(:,5) = BCtmp(:,7) * BCtmp(:,6)