This example shows how to compute a synthetic spectrum in the region 7000-7210 A
for a model with parameters Teff=4970 K, logg=4.6, vturb=1 km/s [M/H]=+0.4
The last release of the Kurucz atomic line data is availabe at:
http://kurucz.harvard.edu/linelists/gfnew/
The example does not use the last release !!!!
To run the SYNTHE procedure type:
source r7000-7210.com
where r7000-7210.com is:
ln -s molecules.dat fort.2
ln -s continua.dat fort.17
./xnfpelsyn.exe< ap00t4950g46k1at12.mod xnfpelsyn.out
mv fort.10 xnft4950g46k1at12.dat
rm fort.*
./synbeg.exe<<"EOF">synbeg.out
AIR 700.0 721.0 500000. 0. 0 10 .001 0 00
AIRorVAC WLBEG WLEND RESOLU TURBV IFNLTE LINOUT CUTOFF NREAD
EOF
ln -s gf0800.100 fort.11
./rgfalllinewnew.exe>gf0800.out
rm fort.11
ln -s gf1200.100 fort.11
./rgfalllinesnew.exe>gf1200.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/ch/ chmasseron.asc fort.11
./rmolecasc.exe>chmasseron.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/ mgh.asc fort.11
./rmolecasc.exe>mgh.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/ nh.asc fort.11
./rmolecasc.exe>nh.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/oh/ ohupdate.asc fort.11
./rmolecasc.exe>oh.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/sih/ sihax.asc fort.11
./rmolecasc.exe>sihax.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/ h2.asc fort.11
./rmolecasc.exe>h2.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/h2/ h2xx.asc fort.11
./rmolecasc.exe>h2xx.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/h2/ hdxx.asc fort.11
./rmolecasc.exe>hdxx.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/ c2ax.asc fort.11
./rmolecasc.exe>c2ax.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/ c2ba.asc fort.11
./rmolecasc.exe>c2ba.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/c2/ c2dabrookek.asc fort.11
./rmolecasc.exe>c2da.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/ c2ea.asc fort.11
./rmolecasc.exe>c2ea.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/cn/ cnaxbrookek.asc fort.11
./rmolecasc.exe>cnax.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/cn/ cnbxbrookek.asc fort.11
./rmolecasc.exe>cnbx.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/cn/ cnxx12brooke.asc fort.11
./rmolecasc.exe>cnxx12.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/co/ coax.asc fort.11
./rmolecasc.exe>coax.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/coxx.asc fort.11
./rmolecasc.exe>coxx.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/sioax.asc fort.11
./rmolecasc.exe>sioax.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/sioex.asc fort.11
./rmolecasc.exe>sioex.out
rm fort.11
ln -s http://kurucz.harvard.edu/linelists/linesmol/sioxx.asc fort.11
./rmolecasc.exe>sioxx.out
rm fort.11
ln -s http://kurucz.harvard.edu/molecules/tio/tioschwenke.bin fort.11
ln -s eschwenke.bin fort.48
./rschwenk.exe>rschwenk.out
rm fort.11
rm fort.48
ln -s http://kurucz.harvard.edu/molecules/h2o/h2ofastfix.bin fort.11
./rh2ofast.exe>h2ofastfix.out
rm fort.11
ln xnft4950g46k1at12.dat fort.10
ln -s he1tables.dat fort.18
./synthe.exe>synthe.out
ln -s molecules.dat fort.2
cat <<"EOF" >fort.25
0.0 0. 1. 0. 0. 0. 0. 0.
0.
RHOXJ R1 R101 PH1 PC1 PSI1 PRDDOP PRDPOW
EOF
./spectrv.exe<"ap04t4970g46k1at12.mod">spectrv.out
mv fort.7 i7000-7210.dat
ln -s i7000-7210.dat fort.1
./rotate.exe<<"EOF">rotate.out
1
2.
EOF
mv ROT1 f7000-7210vr2.dat
ln -s f7000-7210vr2.dat fort.21
./broaden.exe<<"EOF">broaden.out
GAUSSIAN 48000. RESOLUTION
EOF
mv fort.22 f7000-7210vr2br48000ap04t4970g46k1at12.bin
rm fort.*
THE SYNTHETIC SPECTRUM f7000-7210vr2br48000ap04t4970g46k1at12.bin
IS IN BINARY FORM Hnu versus lambda in nm. TO CONVERT IT TO ASCII
as Flamda versus lambda in A:
ln -s f7000-7210vr2br48000ap04t4970g46k1at12.dat fort.1
./converfsynnmtoa.exe
mv fort.2 f7000-7210vr2br48000ap04t4970g46k1at12.asc
rm fort.*
In the above command file there are the codes:
xnfpelsyn which is linked with atlas7v.
It computes atomic and molecular number densities and continous opacities
Input data are:
molecules.dat
continua.dat
The output model from ATLAS12 ap04t4970g46k1at12.dat MUST BE MODIFIED in ap04t4970g46k1at12.mod
1) ADD the following control cards at the top of the model (this holds also for an ATLAS9 input model):
SURFACE INTENSI 17 1.,.9,.8,.7,.6,.5,.4,.3,.25,.2,.15,.125,.1,.075,.05,.025,.01
ITERATIONS 1 PRINT 2 PUNCH 2
CORRECTION OFF
PRESSURE OFF
READ MOLECULES
MOLECULES ON
2) DROP the cards of the ATLAS12 model:
ABUNDANCE TABLE
1H XXX 2He XXX
3Li XXX XX 4Be XXX XX etc..
up to
98Cf XXX CC 99Es XXX XX
This DROPPING is no more needed in the last SYNTHE version (Kurucz, 2015)
synbeg.for--initializes the synthetic spectrum computations
rgfalllinesnew.for--reads the atomic line lists gf****.100 which are available at:
http://kurucz.harvard.edu/LINELISTS/GF100/
More recent data are availabe at:
http://kurucz.harvard.edu/linelists/gfnew/
Use http://kurucz.harvard.edu/linelists/gfnew/splitgf100.for to split the big file gfall21oct16.dat in the gfxxxx.100 files
Modified line lists are available at http://wwwuser.oats.inaf.it/castelli/linelists
rmolecasc.for--reads the molecular line lists in ASCII format.
rschwenk.for--reads the TiO line list tioschwenke.bin (binary file)
and the TiO energy levels eschwenke.bin (binary)
rh2ofast.for--reads the H2O line list h2ofastfix.bin (binary)
synthe.for--computes the line opacity
he1tables.dat--Stark broadening for HeI 4026, 4387, 4471, 4922 A
spectrv.for--spectrv is linked with atlas7v.
It computes in the wavelength range WLBEG--WLEND fixed by synbeg the emergent intensities
at 17 angles or the emergent flux at the resolution RESOLU fixed by synbeg
rotate.for-- computes the emergent flux for different rotational
velocities vsini
broaden.for--broadens the emergent flux for macroturbolence
or instrumental effects.
converfsynnmtoa.for converts the computed flux from binary to ascii