| 381 |
|
else |
| 382 |
|
g = ul1(1)*ul2(1) + ul1(2)*ul2(2) + ul1(3)*ul2(3) |
| 383 |
|
endif |
| 384 |
+ |
|
| 385 |
|
au = a / r |
| 386 |
|
bu = b / r |
| 386 |
– |
|
| 387 |
– |
au2 = au*au |
| 388 |
– |
bu2 = bu*bu |
| 389 |
– |
g2 = g*g |
| 387 |
|
|
| 388 |
+ |
au2 = au * au |
| 389 |
+ |
bu2 = bu * bu |
| 390 |
+ |
g2 = g * g |
| 391 |
|
|
| 392 |
|
H = (xa2 * au2 + xai2 * bu2 - 2.0_dp*x2*au*bu*g) / (1.0_dp - x2*g2) |
| 393 |
|
Hp = (xpap2*au2 + xpapi2*bu2 - 2.0_dp*xp2*au*bu*g) / (1.0_dp - xp2*g2) |
| 394 |
+ |
|
| 395 |
|
sigma = sigma0 / sqrt(1.0_dp - H) |
| 396 |
|
e1 = 1.0_dp / sqrt(1.0_dp - x2*g2) |
| 397 |
|
e2 = 1.0_dp - Hp |
| 411 |
|
|
| 412 |
|
pref1 = - 8.0_dp * eps * mu * (R12 - R6) / (e2 * r) |
| 413 |
|
|
| 413 |
– |
|
| 414 |
|
pref2 = 8.0_dp * eps * s3 * (6.0_dp*R13 - 3.0_dp*R7) / (dw*r*s03) |
| 415 |
|
|
| 416 |
|
dUdr = - (pref1 * Hp + pref2 * (sigma0*sigma0*r/s3 + H)) |
