564 |
|
vpair = vpair + vterm |
565 |
|
epot = epot + sw * vterm |
566 |
|
|
567 |
< |
dudx = dudx - 5.0_dp*sw*vterm*riji*xhat - pref * sw * ri4 * ( & |
567 |
> |
dudx = dudx - 5.0_dp*sw*vterm*riji*xhat + pref * sw * ri4 * ( & |
568 |
|
qxx_j*(6.0_dp*cx_j*ux_j(1) - 2.0_dp*xhat) + & |
569 |
|
qyy_j*(6.0_dp*cy_j*uy_j(1) - 2.0_dp*xhat) + & |
570 |
|
qzz_j*(6.0_dp*cz_j*uz_j(1) - 2.0_dp*xhat) ) |
571 |
< |
dudy = dudy - 5.0_dp*sw*vterm*riji*yhat - pref * sw * ri4 * ( & |
571 |
> |
dudy = dudy - 5.0_dp*sw*vterm*riji*yhat + pref * sw * ri4 * ( & |
572 |
|
qxx_j*(6.0_dp*cx_j*ux_j(2) - 2.0_dp*yhat) + & |
573 |
|
qyy_j*(6.0_dp*cy_j*uy_j(2) - 2.0_dp*yhat) + & |
574 |
|
qzz_j*(6.0_dp*cz_j*uz_j(2) - 2.0_dp*yhat) ) |
575 |
< |
dudz = dudz - 5.0_dp*sw*vterm*riji*zhat - pref * sw * ri4 * ( & |
575 |
> |
dudz = dudz - 5.0_dp*sw*vterm*riji*zhat + pref * sw * ri4 * ( & |
576 |
|
qxx_j*(6.0_dp*cx_j*ux_j(3) - 2.0_dp*zhat) + & |
577 |
|
qyy_j*(6.0_dp*cy_j*uy_j(3) - 2.0_dp*zhat) + & |
578 |
|
qzz_j*(6.0_dp*cz_j*uz_j(3) - 2.0_dp*zhat) ) |
690 |
|
vpair = vpair + vterm |
691 |
|
epot = epot + sw * vterm |
692 |
|
|
693 |
< |
dudx = dudx - 5.0_dp*sw*vterm*riji*xhat - pref * sw * ri4 * ( & |
693 |
> |
dudx = dudx - 5.0_dp*sw*vterm*riji*xhat + pref * sw * ri4 * ( & |
694 |
|
qxx_i*(6.0_dp*cx_i*ux_i(1) - 2.0_dp*xhat) + & |
695 |
|
qyy_i*(6.0_dp*cy_i*uy_i(1) - 2.0_dp*xhat) + & |
696 |
|
qzz_i*(6.0_dp*cz_i*uz_i(1) - 2.0_dp*xhat) ) |
697 |
< |
dudy = dudy - 5.0_dp*sw*vterm*riji*yhat - pref * sw * ri4 * ( & |
697 |
> |
dudy = dudy - 5.0_dp*sw*vterm*riji*yhat + pref * sw * ri4 * ( & |
698 |
|
qxx_i*(6.0_dp*cx_i*ux_i(2) - 2.0_dp*yhat) + & |
699 |
|
qyy_i*(6.0_dp*cy_i*uy_i(2) - 2.0_dp*yhat) + & |
700 |
|
qzz_i*(6.0_dp*cz_i*uz_i(2) - 2.0_dp*yhat) ) |
701 |
< |
dudz = dudz - 5.0_dp*sw*vterm*riji*zhat - pref * sw * ri4 * ( & |
701 |
> |
dudz = dudz - 5.0_dp*sw*vterm*riji*zhat + pref * sw * ri4 * ( & |
702 |
|
qxx_i*(6.0_dp*cx_i*ux_i(3) - 2.0_dp*zhat) + & |
703 |
|
qyy_i*(6.0_dp*cy_i*uy_i(3) - 2.0_dp*zhat) + & |
704 |
|
qzz_i*(6.0_dp*cz_i*uz_i(3) - 2.0_dp*zhat) ) |