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
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
|
subroutine synth
c******************************************************************************
c This program synthesizes a section of spectrum and compares it
c to an observation file.
c******************************************************************************
implicit real*8 (a-h,o-z)
include 'Atmos.com'
include 'Factor.com'
include 'Mol.com'
include 'Linex.com'
include 'Pstuff.com'
include 'Dummy.com'
c*****examine the parameter file
call params
c*****open the files for: standard output, raw spectrum depths, smoothed
c spectra, and (if desired) IRAF-style smoothed spectra
nf1out = 20
lscreen = 4
array = 'STANDARD OUTPUT'
nchars = 15
call infile ('output ',nf1out,'formatted ',0,nchars,
. f1out,lscreen)
nf2out = 21
lscreen = lscreen + 2
array = 'RAW SYNTHESIS OUTPUT'
nchars = 20
call infile ('output ',nf2out,'formatted ',0,nchars,
. f2out,lscreen)
if (plotopt .ne. 0) then
nf3out = 22
lscreen = lscreen + 2
array = 'SMOOTHED SYNTHESES OUTPUT'
nchars = 25
call infile ('output ',nf3out,'formatted ',0,nchars,
. f3out,lscreen)
if (f5out .ne. 'optional_output_file') then
nf5out = 26
lscreen = lscreen + 2
array = 'POSTSCRIPT PLOT OUTPUT'
nchars = 22
call infile ('output ',nf5out,'formatted ',0,nchars,
. f5out,lscreen)
endif
endif
if (iraf .ne. 0) then
nf4out = 23
lscreen = lscreen + 2
array = 'IRAF ("rtext") OUTPUT'
nchars = 24
call infile ('output ',nf4out,'formatted ',0,nchars,
. f4out,lscreen)
endif
c*****open and read the model atmosphere file
nfmodel = 30
lscreen = lscreen + 2
array = 'THE MODEL ATMOSPHERE'
nchars = 20
call infile ('input ',nfmodel,'formatted ',0,nchars,
. fmodel,lscreen)
call inmodel
c*****open the line list file and the strong line list file
nflines = 31
lscreen = lscreen + 2
array = 'THE LINE LIST'
nchars = 13
call infile ('input ',nflines,'formatted ',0,nchars,
. flines,lscreen)
if (dostrong .gt. 0) then
nfslines = 32
lscreen = lscreen + 2
array = 'THE STRONG LINE LIST'
nchars = 20
call infile ('input ',nfslines,'formatted ',0,nchars,
. fslines,lscreen)
endif
c*****do the syntheses
choice = '1'
do i=1,100
if (i .eq. 100) then
write (*,1002)
stop
endif
ncall = 1
call getsyns (lscreen,ncall)
c*****now either don't make a plot (plotopt = 0)
c plot the synthetic spectrum, (plotopt = 1)
c plot syntheses and observation (plotopt = 2)
c or just smooth the syntheses (plotopt = 3)
if (choice .eq. 'n') then
ncall = 2
else
ncall = 1
endif
if (plotopt .eq. 0) then
choice = 'q'
elseif (plotopt .eq. 1) then
call pltspec (lscreen,ncall)
elseif (plotopt .eq. 2) then
nfobs = 33
lscreen = lscreen + 2
array = 'THE OBSERVED SPECTRUM'
nchars = 21
if (specfileopt .eq. 1) then
call infile ('input ',nfobs,'unformatted',2880,nchars,
. fobs,lscreen)
else
call infile ('input ',nfobs,'formatted ',0,nchars,
. fobs,lscreen)
endif
call pltspec (lscreen,ncall)
elseif (plotopt .eq. 3) then
call smooth (-1,ncall)
choice = 'q'
else
write (*,1001)
stop
endif
if (choice .eq. 'q') then
call finish (0)
exit
endif
enddo
c*****format statements
1001 format ('for syntheses, parameter "plot" must be 0, 1, 2, or 3;',
. ' I QUIT!')
1002 format ('something wrong: max number (99) of synthesis ',
. 'cycles exceeded; I QUIT!')
end
|
