Lines Matching refs:rGNa
56 & tmp(npro), rGNa(npro,nsd,nsd), rNa(npro,nsd), local
64 rGNa = zero
125 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp
126 rGNa(:,1,2) = tmp1
127 rGNa(:,1,3) = tmp3
128 rGNa(:,2,1) = tmp1
129 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp
130 rGNa(:,2,3) = tmp2
131 rGNa(:,3,1) = tmp3
132 rGNa(:,3,2) = tmp2
133 rGNa(:,3,3) = two * rmu * g3yi(:,4) + tmp
139 rGNa(:,1,1) = two * rmu * g1yi(:,2) + tmp - u1(:)*u1(:)*rho(:)
140 rGNa(:,1,2) = tmp1 - u1(:)*u2(:)*rho(:)
141 rGNa(:,1,3) = tmp3 - u1(:)*u3(:)*rho(:)
142 rGNa(:,2,1) = tmp1 - u1(:)*u2(:)*rho(:)
143 rGNa(:,2,2) = two * rmu * g2yi(:,3) + tmp - u2(:)*u2(:)*rho(:)
144 rGNa(:,2,3) = tmp2 - u3(:)*u2(:)*rho(:)
145 rGNa(:,3,1) = tmp3 - u1(:)*u3(:)*rho(:)
146 rGNa(:,3,2) = tmp2 - u3(:)*u2(:)*rho(:)
147 rGNa(:,3,3) = two * rmu * g3yi(:,4) + tmp - u3(:)*u3(:)*rho(:)
150 rGNa(:,1,1) = rGNa(:,1,1) - rlsli(:,1)*rho(:)
151 rGNa(:,1,2) = rGNa(:,1,2) - rlsli(:,4)*rho(:)
152 rGNa(:,1,3) = rGNa(:,1,3) - rlsli(:,5)*rho(:)
153 rGNa(:,2,1) = rGNa(:,2,1) - rlsli(:,4)*rho(:)
154 rGNa(:,2,2) = rGNa(:,2,2) - rlsli(:,2)*rho(:)
155 rGNa(:,2,3) = rGNa(:,2,3) - rlsli(:,6)*rho(:)
156 rGNa(:,3,1) = rGNa(:,3,1) - rlsli(:,5)*rho(:)
157 rGNa(:,3,2) = rGNa(:,3,2) - rlsli(:,6)*rho(:)
158 rGNa(:,3,3) = rGNa(:,3,3) - rlsli(:,3)*rho(:)
165 rGNa(:,1,1) = rGNa(:,1,1) + tmp1 * u1
166 rGNa(:,1,2) = rGNa(:,1,2) + tmp1 * u2
167 rGNa(:,1,3) = rGNa(:,1,3) + tmp1 * u3
168 rGNa(:,2,1) = rGNa(:,2,1) + tmp2 * u1
169 rGNa(:,2,2) = rGNa(:,2,2) + tmp2 * u2
170 rGNa(:,2,3) = rGNa(:,2,3) + tmp2 * u3
171 rGNa(:,3,1) = rGNa(:,3,1) + tmp3 * u1
172 rGNa(:,3,2) = rGNa(:,3,2) + tmp3 * u2
173 rGNa(:,3,3) = rGNa(:,3,3) + tmp3 * u3
180 rGNa(:,1,1) = rGNa(:,1,1) + tmp1 * u1
181 rGNa(:,1,2) = rGNa(:,1,2) + tmp2 * u1
182 rGNa(:,1,3) = rGNa(:,1,3) + tmp3 * u1
183 rGNa(:,2,1) = rGNa(:,2,1) + tmp1 * u2
184 rGNa(:,2,2) = rGNa(:,2,2) + tmp2 * u2
185 rGNa(:,2,3) = rGNa(:,2,3) + tmp3 * u2
186 rGNa(:,3,1) = rGNa(:,3,1) + tmp1 * u3
187 rGNa(:,3,2) = rGNa(:,3,2) + tmp2 * u3
188 rGNa(:,3,3) = rGNa(:,3,3) + tmp3 * u3
205 rGNa(:,1,1) = rGNa(:,1,1) + tmp1 * rLui(:,1)
206 rGNa(:,1,2) = rGNa(:,1,2) + tmp1 * rLui(:,2)
207 rGNa(:,1,3) = rGNa(:,1,3) + tmp1 * rLui(:,3)
208 rGNa(:,2,1) = rGNa(:,2,1) + tmp2 * rLui(:,1)
209 rGNa(:,2,2) = rGNa(:,2,2) + tmp2 * rLui(:,2)
210 rGNa(:,2,3) = rGNa(:,2,3) + tmp2 * rLui(:,3)
211 rGNa(:,3,1) = rGNa(:,3,1) + tmp3 * rLui(:,1)
212 rGNa(:,3,2) = rGNa(:,3,2) + tmp3 * rLui(:,2)
213 rGNa(:,3,3) = rGNa(:,3,3) + tmp3 * rLui(:,3)
240 & + shg(:,aa,1) * rGNa(:,1,1)
241 & + shg(:,aa,2) * rGNa(:,1,2)
242 & + shg(:,aa,3) * rGNa(:,1,3) )
245 & + shg(:,aa,1) * rGNa(:,2,1)
246 & + shg(:,aa,2) * rGNa(:,2,2)
247 & + shg(:,aa,3) * rGNa(:,2,3) )
250 & + shg(:,aa,1) * rGNa(:,3,1)
251 & + shg(:,aa,2) * rGNa(:,3,2)
252 & + shg(:,aa,3) * rGNa(:,3,3) )
287 real*8 rGNa(npro,nsd), rNa(npro), rcp(npro), tmp(npro) local
294 rGNa = zero
313 rGNa(:,1) = diffus * gradS(:,1) + uMod(:,1) * tmp
314 rGNa(:,2) = diffus * gradS(:,2) + uMod(:,2) * tmp
315 rGNa(:,3) = diffus * gradS(:,3) + uMod(:,3) * tmp
323 rGNa(:,1) = rGNa(:,1) + gGradS(:,1) ! gGradS is
324 rGNa(:,2) = rGNa(:,2) + gGradS(:,2) ! g^{ij}*Y_{j}*dcFct
325 rGNa(:,3) = rGNa(:,3) + gGradS(:,3) ! calculated in e3dc.f
336 & + shg(:,aa,1) * rGNa(:,1)
337 & + shg(:,aa,2) * rGNa(:,2)
338 & + shg(:,aa,3) * rGNa(:,3) )