Nbody6
Main Page
Files
File List
File Members
Nbody6
Files
File List
Nbody6
Chain
Nchain
Ncode
adjust.f
assess.f
bindat.f
binev.f
binout.f
binpop.f
block.f
bodies.f
brake.f
brake2.f
brake3.f
brake4.f
bsetid.f
chaos.f
chaos0.f
check.f
checkl.f
chrect.f
circ.f
clint.f
cloud.f
cloud0.f
cmbody.f
cmcorr.f
cmfirr.f
cmfreg.f
coal.f
comenv.f
common6.h
core.f
corerd.f
cputim.f
data.f
decide.f
define.f
deform.f
deform2.f
degen.f
delay.f
derqp3.f
derqp4.f
dgcore.f
difsy3.f
difsy4.f
dtchck.f
eccmod.f
ecirc.f
edot.f
efac2.f
efac3.f
endreg.f
energy.f
erel3.f
erel4.f
escape.f
events.f
evolve.f
expand.f
expand2.f
expel.f
extend.f
fbulge.f
fclose.f
fcloud.f
fcorr.f
fdisk.f
fhalo.f
ficorr.f
findj.f
findm.f
flyby.f
fnuc.f
fpcorr.f
fpert.f
fpoly1.f
fpoly2.f
freeze.f
gcinit.f
gcint.f
giant.f
giant3.f
gntage.f
grrad.f
hcorr.f
hiarch.f
hicirc.f
hidat.f
higrow.f
himax.f
himax2.f
himod.f
hipop.f
hirect.f
histab.f
hivel.f
hmdot.f
hmdot2.f
hotsys.f
hrdiag.f
hrplot.f
hut.f
hut2.f
iblock.f
ichain.f
imf.f
imf2.f
imfbd.f
impact.f
induce.f
inext.f
input.f
insert.f
instar.f
intgrt.f
intide.f
jacobi.f
kepler.f
kick.f
kick2.f
kozai.f
ksapo.f
kscorr.f
ksin2.f
ksinit.f
ksint.f
ksintp.f
kslist.f
ksmod.f
ksperi.f
kspert.f
kspoly.f
kspred.f
ksrect.f
ksreg.f
ksres.f
ksres2.f
ksrot.f
ksterm.f
kstide.f
lagr.f
lagr2.f
levels.f
magbrk.f
matrix.f
mdot.f
merge.f
merge2.f
mix.f
mloss.f
mlwind.f
modify.f
mrenv.f
mtrace.f
mydump.f
nbhist.f
nbint.f
nblist.f
nbody6.f
nbpot.f
nbrem.f
nbrest.f
nbsort.f
nbtide.f
newhut.f
newreg.f
newsys.f
newtev.f
newton.f
nstab.f
ntint.f
offset.f
orbit.f
output.f
params.h
peri.f
permit.f
pfac.f
plpot1.f
pnpert2.f
poarams.h
poti.f
premsf.f
proto_star.f
qpmod3.f
qpmod4.f
qtides.f
quad.f
ran2.f
rchain.f
reflct.f
regint.f
remove.f
rename.f
reset.f
reset2.f
resolv.f
rkint.f
rl.f
roche.f
rpmax.f
rpmax2.f
rpmin.f
rsort.f
scale.f
search.f
setup.f
setup2.f
shrink.f
sort1.f
spinup.f
spiral.f
stability.f
stabl3.f
stabl4.f
stablz.f
star.f
start.f
start3.f
start4.f
status.f
stepi.f
stepk.f
steps.f
stumpf.f
subint.f
subsys.f
swap.f
sweep.f
synch.f
tail0.f
tcirc.f
test.f
tides.f
tides2.f
tides3.f
touch.f
tperi.f
tpert.f
trans3.f
trans4.f
trdot.f
trdot2.f
trflow.f
triple.f
tstab.f
tstep.f
units.f
unpert.f
update.f
verify.f
xtrnl0.f
xtrnld.f
xtrnlf.f
xtrnlp.f
xtrnlt.f
xtrnlv.f
xvpred.f
zare.f
zcnsts.f
zdata.h
zero.f
zfuncs.f
File Members
•
All
Files
Functions
Variables
pfac.f
Go to the documentation of this file.
1
REAL*8
FUNCTION
pfac
(A,Z)
2
3
* Precession factor for hierarchy.
4
* --------------------------------
5
*
6
* Maths by Douglas Heggie; function by Rosemary Mardling (11/96).
7
* @@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
8
*
9
REAL*8
m,m1,a,z,c,c1,c2,c3,d1,d2,d3,z1
10
11
*The function is singular for m1=0, but the following holds:
12
*For 0 < m1 <=1 (m1 never > 1), \infty > PFAC(m1) > 1.6, with P(0.01)=3.6.
13
*So taking it as 1 isn't too bad for an order of magnitude estimate.
14
15
16
c1 = 0.5*(z + a**2)
17
z1 = 5.0 - z + 4*a**2
18
c2 = (z1 + sqrt(z1))/6.0
19
c3 = (z1 - sqrt(z1))/6.0
20
IF
(c1.GE.c2)
THEN
21
d1 = c1
22
d2 = c2
23
d3 = c3
24
ELSE
IF
(c2.GE.c1.AND.c1.GE.c3)
THEN
25
d1 = c2
26
d2 = c1
27
d3 = c3
28
ELSE
29
d1 = c2
30
d2 = c3
31
d3 = c1
32
END IF
33
*
34
c = 1.0/sqrt(d1 - d3)
35
m = (d2 - d3)/(d1 - d3)
36
*
37
* Evaluate elliptic integral by approximate expression.
38
a0=1.3862944
39
a1=0.1119723
40
a2=0.0725296
41
b0=0.5
42
b1=0.1213478
43
b2=0.0288729
44
m1 = 1.0 - m
45
m1 = max(m1,1.0d-10)
46
47
pfac
= (a0+a1*m1+a2*m1**2)+(b0+b1*m1+b2*m1**2)*log(1.0/m1)
48
*
49
* Include the square root of (D1 - D3).
50
pfac
= c*
pfac
51
52
END
Nbody6
Ncode
pfac.f
Generated on Sat Nov 30 2013 19:04:37 for Nbody6 by
1.8.2