-
Notifications
You must be signed in to change notification settings - Fork 8
/
Copy pathFakeline.f
executable file
·91 lines (73 loc) · 2.46 KB
/
Fakeline.f
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
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
subroutine fakeline
c******************************************************************************
c This routine creates a representative atomic line with parameters
c that are arbitrarily set in this routine. Then this line is used to
c generate a curve-of-growth for use with the 'abfind' routine.
c******************************************************************************
implicit real*8 (a-h,o-z)
include 'Atmos.com'
include 'Linex.com'
include 'Quants.com'
include 'Dampdat.com'
c*****specify the input parameters of a representative "fake" line;
c*****here we use a real Fe I line because it has Barklem damping data
wave1(1) = 5006.126
atom1(1) = 26.0
e(1,1) = 2.833
e(1,2) = 5.308
gf(1) = 1.0e-3
iatom = atom1(1)
charge(1) = 1.0
amass(1) = xam(iatom)
chi(1,1) = xchi1(iatom)
chi(1,2) = xchi2(iatom)
chi(1,3) = xchi3(iatom)
gammabk(1) = -7.280
alphabk(1) = 0.238
gambark(1) = 10.**gammabk(1)
alpbark(1) = (1.-alphabk(1))/2.
if (dampingopt .eq. 0) then
damptype(1) = 'UNSLDc6'
elseif (dampingopt .eq. 1) then
damptype(1) = 'BKgamma'
elseif (dampingopt .eq. 2) then
damptype(1) = 'BLKWLc6'
else
damptype(1) = 'NEXTGEN'
endif
c*****do the partition function
call partfn (atom1(1),int(atom1(1)+0.001))
c*****now calculate the dependent quantities, such as doppler parameters,
c damping constant, line opacity at line center
lim1line = 1
lim2line = 1
nlines = 1
idump = nf2out
nf2out = 0
call nearly (3)
c*****setting some counting parameters and do a curve-of-growth
lim1 = 1
lim2 = 1
rwlow = -6.7
rwhigh = -3.7
rwstep = 0.15
call curve
nf2out = idump
c*****make an interpolated curve-of-growth at very small log gf steps
do i=2,ncurve-2
k = 30*(i-2)
l = -1
do m=k+1,k+30
l = l + 1
pp = 0.005/0.15*(l-1)
gftab(m) = gf1(2) + 0.005*(m-1)
rwtab(m) = w(i-1)*(-pp)*(pp-1.)*(pp-2.)/6. +
. w(i)*(pp*pp-1.)*(pp-2.)/2. +
. w(i+1)*(-pp)*(pp+1.)*(pp-2.)/2. +
. w(i+2)*pp*(pp*pp-1.)/6.
enddo
enddo
ntabtot = m - 1
c*****exit back to the abfind driver
return
end