8 common/galaxy/ gmg,rg(3),vg(3),fg(3),fgd(3),tg,
9 & omega,disk,a,b,v02,rl2,gmb,ar,gam,zdum(7)
14 IF (gmg.GT.0.0d0)
THEN
15 rin2 = 1.0/(rg(1)**2 + rg(2)**2 + rg(3)**2)
16 rin3 = rin2*sqrt(rin2)
17 rgvg = 3.0*(rg(1)*vg(1) + rg(2)*vg(2) + rg(3)*vg(3))*rin2
20 fg(k) = -gmg*rg(k)*rin3
21 fgd(k) = -gmg*(vg(k) - rgvg*rg(k))*rin3
31 IF (gmb.GT.0.0d0)
THEN
35 fgd(k) = fgd(k) + fdd(k)
40 IF (disk.GT.0.0d0)
THEN
41 CALL
fdisk(rg,vg,fd,fdd)
44 fgd(k) = fgd(k) + fdd(k)
49 IF (v02.GT.0.0d0)
THEN
50 CALL
fhalo(rg,vg,fd,fdd)
53 fgd(k) = fgd(k) + fdd(k)