Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- Q22 = Q22 + trd*pot*rho(l1,m1)*ILU(l1, m1, 0, 0, 0, 0, 0, 2, 1, 0)*r*it2 -&
- frd*pot*rho(l1,m1)*ILU(l1, m1, 0, 0, 0, 0, 1, 1, 2, 0)*r*it2 + &
- trd*pot*rho(l1,m1)*ILU(l1, m1, 0, 0, 0, 0, 2, 0, 0, 2)*r*it2 - &
- trd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 0, 2, 0, 0)*vr(l1,m2)*vr(l2,m2) - &
- trd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 0, 2, 0, 2)*r**2*vth(l2,m2)*vth(l3,m3) + &
- frd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 0, 2, 2, 0)*r**2*vth(l2,m2)*vth(l3,m3) + &
- frd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 1, 1, 0, 0)*r*vr(l2,m2)*vth(l3,m3) - &
- frd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 1, 1, 0, 2)*r*vr(l2,m2)*vth(l3,m3) + &
- 4.d0*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 1, 1, 1, 1)*r**2*vph(l2,m2)*vth(l3,m3) + &
- erd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 1, 1, 2, 0)*r*vr(l2,m2)*vth(l3,m3) -&
- trd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 2, 0, 0, 0)*r**2*vth(l2,m2)*vth(l3,m3) +&
- frd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 2, 0, 0, 2)*r**2*vph(l2,m2)*vph(l3,m3) - &
- trd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 2, 0, 0, 2)*vr(l2,m2)*vr(l3,m3) + &
- 4.d0*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 2, 0, 1, 1)*r*vph(l2,m2)*vr(l3,m3) - &
- trd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 2, 0, 2, 0)*r**2*vph(l2,m2)*vph(l3,m3) + &
- frd*rho(l1,m1)*ILU(l1, m1, l2, m2, l3, m3, 2, 0, 2, 0)*vr(l2,m2)*vr(l3,m3)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement