| 417 |
|
gpi = 1.0d0 / gp |
| 418 |
|
gmum = gmu*gpi |
| 419 |
|
|
| 420 |
< |
curlyE = 1.0d0/dsqrt(1.0d0 - Chi*Chi*u1dotu2*u1dotu2) |
| 420 |
> |
curlyE = 1.0d0/sqrt(1.0d0 - Chi*Chi*u1dotu2*u1dotu2) |
| 421 |
|
dcE = (curlyE**3)*Chi*Chi*u1dotu2 |
| 422 |
|
|
| 423 |
|
dcEdu1x = dcE*ul2(1) |
| 638 |
|
l2 = (gb_sigma*gb_l2b_ratio)**2 |
| 639 |
|
d2 = gb_sigma**2 |
| 640 |
|
lj2 = ljsigma**2 |
| 641 |
< |
s0 = dsqrt(d2 + lj2) |
| 641 |
> |
s0 = sqrt(d2 + lj2) |
| 642 |
|
|
| 643 |
|
chioalpha2 = (l2 - d2)/(l2 + lj2) |
| 644 |
|
|
| 645 |
< |
eE = dsqrt(gb_eps*gb_eps_ratio*ljeps) |
| 646 |
< |
eS = dsqrt(gb_eps*ljeps) |
| 645 |
> |
eE = sqrt(gb_eps*gb_eps_ratio*ljeps) |
| 646 |
> |
eS = sqrt(gb_eps*ljeps) |
| 647 |
|
moom = 1.0d0 / gb_mu |
| 648 |
|
mum1 = gb_mu-1 |
| 649 |
|
chipoalphap2 = 1 - (eE/eS)**moom |
| 651 |
|
!! mess matches cleaver (eq 20) |
| 652 |
|
|
| 653 |
|
mess = 1-rdotu*rdotu*chioalpha2 |
| 654 |
< |
sab = 1.0d0/dsqrt(mess) |
| 654 |
> |
sab = 1.0d0/sqrt(mess) |
| 655 |
|
|
| 656 |
|
dsabdct = s0*sab*sab*sab*rdotu*chioalpha2 |
| 657 |
|
|
