-
Notifications
You must be signed in to change notification settings - Fork 2
/
const_pot.h
36 lines (27 loc) · 1014 Bytes
/
const_pot.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
c-----------------------------------------------------------------------
c This is 'const_pot.h'
c Constants related to galactic potential
REAL*8 Mb,Md,Mh,K
c general values
PARAMETER (R0 = 8.5d0, !solar radius in kpc
+ Usol = 10.00d0, ! solar velocity(Dehnen i Binney 1998)
+ Vsol = 5.25d0,
+ Wsol = 7.17d0)
c values of the potential
PARAMETER (ZMsol = 1.989d33,! Msol in gm
+ abu = 0.387d0,! en Kpc
+ ad = 5.3178d0,
+ ah = 12.0d0,
+ bd = 0.25d0,
+ Mb = 1.406d10,!in Msol
+ Md = 8.561d10,!in Msol
+ Mh = 1.071d11)!in Msol
c other variables:
PARAMETER(K=Gte*ZMsol/pc/1.d13,
+ abuabu=abu*abu,
+ bdbd=bd*bd,
+ rKMb=K*Mb,
+ q9=-3.d0*rKMb,
+ rKMd=K*Md,
+ q8=-3.d0*rKMd,
+ rKMh=K*Mh)