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
efac2.f
Go to the documentation of this file.
1
FUNCTION
efac2
(ZETA,ICASE)
2
*
3
*
4
* Tidal capture efficiency factor (l = 2).
5
* ----------------------------------------
6
*
7
* Fitting functions for the second spherical harmonic index based on
8
* Ray et. al. (A & A 184, 164) and Lee & Ostriker (Ap. J. 310, 176).
9
* Developed at IOA by S. Portegies Zwart & T. Meinen (A & A 280, 174).
10
*
11
REAL*8
zeta,coeff(6),
efac2
12
*
13
*
14
* Select coefficients for a given polytropic index (ICASE = 1, 2, 3).
15
IF
(icase.EQ.1)
THEN
16
* Polytropic index 1.5.
17
coeff(1) =-0.397
18
coeff(2) = 1.678
19
coeff(3) = 1.277
20
coeff(4) =-12.42
21
coeff(5) = 9.446
22
coeff(6) =-5.550
23
ELSE
24
IF
(icase.EQ.2)
THEN
25
* Polytropic index 2.
26
coeff(1) =-0.517
27
coeff(2) =-0.906
28
coeff(3) = 23.88
29
coeff(4) =-93.49
30
coeff(5) = 112.3
31
coeff(6) =-44.15
32
ELSE
33
* Polytropic index 3.
34
coeff(1) =-1.124
35
coeff(2) = 0.877
36
coeff(3) =-13.37
37
coeff(4) = 21.55
38
coeff(5) =-16.48
39
coeff(6) = 4.124
40
END IF
41
END IF
42
*
43
* Obtain the tidal energy factor from fifth-order polynomial fit.
44
y = log10(zeta)
45
efac2
= ((((coeff(6)*y + coeff(5))*y + coeff(4))*y + coeff(3))*y
46
& + coeff(2))*y + coeff(1)
47
efac2
= 10.0**
efac2
48
*
49
RETURN
50
*
51
END
Nbody6
Ncode
efac2.f
Generated on Sat Nov 30 2013 19:04:36 for Nbody6 by
1.8.2