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
reflct.f
Go to the documentation of this file.
1
SUBROUTINE
reflct
(I,RI2)
2
*
3
*
4
* Boundary reflection.
5
* --------------------
6
*
7
include
'common6.h'
8
*
9
*
10
rdot = x(1,i)*xdot(1,i) + x(2,i)*xdot(2,i) + x(3,i)*xdot(3,i)
11
*
12
* Avoid multiple reflections for inward moving particles.
13
IF
(rdot.LT.0.0) go to 50
14
*
15
* Obtain radial velocity with respect to the inertial centre.
16
ri = sqrt(ri2)
17
rdot = rdot/ri
18
*
19
* Form radial velocity for mirror reflection.
20
rdot = -2.0*rdot
21
*
22
* Set reflected velocity at the boundary and initialize X0DOT.
23
DO
10 k = 1,3
24
xdot(k,i) = xdot(k,i) + rdot*x(k,i)/ri
25
x0dot(k,i) = xdot(k,i)
26
10
CONTINUE
27
*
28
* Predict current coordinates & velocities for all neighbours.
29
nnb = list(1,i)
30
CALL
xvpred
(i,nnb)
31
*
32
* Initialize force polynomials & time-steps.
33
CALL
fpoly1
(i,i,0)
34
CALL
fpoly2
(i,i,0)
35
*
36
ri2 = -1.0
37
* Return of negative argument indicates new polynomials & steps.
38
*
39
50
RETURN
40
*
41
END
Nbody6
Ncode
reflct.f
Generated on Sat Nov 30 2013 19:04:37 for Nbody6 by
1.8.2