/
cos_quasarcomposite.pro
59 lines (57 loc) · 1.41 KB
/
cos_quasarcomposite.pro
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
; docformat = 'rst'
;
;+
;
;
; :Categories:
; IFSFIT
;
; :Returns:
; The best fit continuum spectrum (over all wavelengths, not just
; those fit).
;
; :Params:
; x: in, required, type=dblarr(N)
; Wavelength array
; p: in, required, type=dblarr(N)
; Parameter array
; ymod: out, required, type=dblarr(N)
; Model flux
;
; :Keywords:
;
; :Author:
; David S. N. Rupke::
; Rhodes College
; Department of Physics
; 2000 N. Parkway
; Memphis, TN 38104
; drupke@gmail.com
;
; :History:
; Change History::
; 2018feb16, DSNR, created
;
; :Copyright:
; Copyright (C) 2018 David S. N. Rupke
;
; This program is free software: you can redistribute it and/or
; modify it under the terms of the GNU General Public License as
; published by the Free Software Foundation, either version 3 of
; the License or any later version.
;
; This program is distributed in the hope that it will be useful,
; but WITHOUT ANY WARRANTY; without even the implied warranty of
; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
; General Public License for more details.
;
; You should have received a copy of the GNU General Public License
; along with this program. If not, see
; http://www.gnu.org/licenses/.
;
;-
pro cos_quasarcomposite,x,p,ymod,twave=twave,tflux=tflux
ymod = tflux
ymod *= p[0]*twave^p[1]
ymod += (p[2] + p[3]*twave)
end