Lines Matching refs:e
128 integer e local
172 do e = 1, npro
174 nuInv = rho(e)/rmu(e)
175 k = abs(kay(e))
176 if (k.lt.1.e-32) k=0
177 eps = abs(epsilon(e))
178 y = dwi(e)
193 if ( k .gt.1.e-32 ) then
208 if ( eps .gt.1.e-32) then
214 ss = gradV(e,1,1) ** 2
215 & + gradV(e,2,2) ** 2
216 & + gradV(e,3,3) ** 2
218 & * ( (gradV(e,2,3) + gradV(e,3,2)) ** 2
219 & + (gradV(e,3,1) + gradV(e,1,3)) ** 2
220 & + (gradV(e,1,2) + gradV(e,2,1)) ** 2 )
247 if(fmuKe.gt.1e-32) then
262 if(ff1.gt.1.0e-32) ff1Inv=1./ff1
263 if(f2.gt.1.0e-32) f2Inv =1./f2
286 if(Tscale.gt.1.e-32) Tinv=1./Tscale
298 if(Tscale.gt.1.e-32) Tinv=1./Tscale
320 mut = rho(e)*CmuKE*fmuKE* k *Tscale
322 mut_k = rho(e)*CmuKE*(
328 mut_eps = rho(e) * CmuKE *k*(
337 pk2 = rho(e)*two3rdq*k*sqrt(ss)
351 src(e,1) = pk - rho(e) * eps
352 src(e,2) = ff1 * Ce1 * pk * Tinv
353 & -rho(e)* Ce2 * f2 * eps * Tinv
355 srcRat(e,1) = -kInv*( pk - rho(e) * eps )
356 srcRat(e,2) = -epsInv*(
358 & -rho(e)* Ce2 * f2 * eps * Tinv
361 srcJac(e,1) = -(pk_k) ! jacobian for k PDE alone
362 srcJac(e,3) = -(pk_eps-rho(e)) ! d(Fsrck)/d(epsilon)
363 srcJac(e,2) = -(
367 & -rho(e) * Ce2 * f2_k * eps * Tinv
368 & +rho(e) * Ce2 * f2 * eps * Tscale_k * Tinv**2 ! do not touch
370 srcJac(e,4) = -(
374 & -rho(e) * Ce2 * f2_eps * eps * Tinv
375 & -rho(e) * Ce2 * f2 * Tinv
376 & +rho(e) * Ce2 * f2* eps * Tscale_eps * Tinv**2
386 do e = 1, npro
389 srcJac(e,1) = max( srcJac(e,1), srcRat(e,1), 0.d0 )
390 srcJac(e,4) = max( srcJac(e,4), srcRat(e,2), 0.d0 )
395 tmp1 = min( srcJac(e,1) * srcJac(e,4),
396 & (srcJac(e,1)-srcRat(e,1)) *
397 & (srcJac(e,4)-srcRat(e,2)) )
398 tmp2 = srcJac(e,2) * srcJac(e,3)
403 srcJac(e,2) = tmp2 * srcJac(e,2)
404 srcJac(e,3) = tmp2 * srcJac(e,3)