This example shows how to compute a synthetic spectrum in the region 2500-15000 A
for a model with parameters Teff=5750 K, logg=4.0, vturb=2 km/s [M/H]=+0.5.
The spectrum computed at a resolving power of 500000 is then degraded at
a resolution of 1 A.

To run the SYNTHE procedure type:

source r2500-10500p.com

where r2500-10500p.com is:

ln -s molecules.dat fort.2
ln -s continua.dat fort.17
./xnfpelsyn.exe < ap05t5750g40k2odfnew.mod >xnfpelsyn.out  !! ap05t5750g40k2odfnew.mod--> input model
mv fort.10 xnfp05t5750g40k2odfnew.dat
rm fort.*
./synbeg.exe<synbeg.out
AIR       250.0     1050.0    500000.   0.     0         -10 .001         1   00
AIRorVAC  WLBEG     WLEND     RESOLU    TURBV  IFNLTE LINOUT CUTOFF        NREAD
EOF
ln -s  gf0300.100  fort.11
./rgfalllinesnew.exe>gf0300.out
rm fort.11
ln -s  gf0400.100  fort.11
./rgfalllinesnew.exe>gf0400.out
rm fort.11
ln -s  gf0500.100  fort.11
./rgfalllinesnew.exe>gf0500.out
rm fort.11
ln -s  gf0600.100  fort.11
./rgfalllinesnew.exe>gf0600.out
rm fort.11
ln -s  gf0800.100  fort.11
./rgfalllinesnew.exe>gf0800.out
rm fort.11
ln -s  gf1200.100  fort.11
./rgfalllinesnew.exe>gf1200.out
rm fort.11
ln -s  fclowlines.bin  fort.11  !!lowlines.bin does not work with Linux
./rpredict.exe>predictlow.out
rm fort.11
ln -s  fchighlines.bin  fort.11  !!highlines.bin does not work with Linux
./rpredict.exe>predicthigh.out
rm fort.11
ln -s   chmasseron.asc  fort.11
./rmolecasc.exe>chmasseron.out
rm fort.11
ln -s  mgh.asc  fort.11
./rmolecasc.exe>mgh.out
rm fort.11
ln -s  nh.asc  fort.11
./rmolecasc.exe>nh.out
rm fort.11
ln -s   ohupdate.asc  fort.11
./rmolecasc.exe>oh.out
rm fort.11
ln -s   sihax.asc  fort.11
./rmolecasc.exe>sihax.out
rm fort.11
ln -s  h2.asc  fort.11
./rmolecasc.exe>h2.out
rm fort.11
ln -s  h2xx.asc  fort.11
./rmolecasc.exe>h2xx.out
rm fort.11
ln -s  hdxx.asc  fort.11
./rmolecasc.exe>hdxx.out
rm fort.11
ln -s  c2ax.asc  fort.11
./rmolecasc.exe>c2ax.out
rm fort.11
ln -s  c2ba.asc  fort.11
./rmolecasc.exe>c2ba.out
rm fort.11
ln -s   c2dabrookek.asc  fort.11
./rmolecasc.exe>c2da.out
rm fort.11
ln -s   c2ea.asc  fort.11
./rmolecasc.exe>c2ea.out
rm fort.11
ln -s  cnaxbrookek.asc  fort.11
./rmolecasc.exe>cnax.out
rm fort.11
ln -s  cnbxbrookek.asc  fort.11
./rmolecasc.exe>cnbx.out
rm fort.11
ln -s  cnxx12brooke.asc  fort.11
./rmolecasc.exe>cnxx12.out
rm fort.11
ln -s  coax.asc  fort.11
./rmolecasc.exe>coax.out
rm fort.11
ln -s  coxx.asc fort.11
./rmolecasc.exe>coxx.out
rm fort.11
ln -s  sioax.asc  fort.11
./rmolecasc.exe>sioax.out
rm fort.11
ln -s  sioex.asc  fort.11
./rmolecasc.exe>sioex.out
rm fort.11
ln -s  sioxx.asc  fort.11
./rmolecasc.exe>sioxx.out
rm fort.11
ln -s  tioschwenke.bin  fort.11
ln -s eschwenke.bin fort.48
./rschwenk.exe>rschwenk.out
rm fort.11
rm fort.48
ln -s  h2ofastfix.bin  fort.11
./rh2ofast.exe>h2ofastfix.out
rm fort.11
ln -s xnfp05t5750g40k2odfnew.dat fort.10
ln -s he1tables.dat fort.18
./synthe.exe>synthe.out
ln -s molecules.dat fort.2
cat <fort.25
0.0       0.        1.        0.        0.        0.        0.        0.
0.
RHOXJ     R1        R101      PH1       PC1       PSI1      PRDDOP    PRDPOW
EOF
./spectrv.exe<ap05t5750g40k2odfnew.mod>spectrv.out
mv fort.7 i2500-10500.dat
ln -s i2500-10500.dat fort.1
./rotate.exe<rotate.out
    1
0.
EOF
mv ROT1 f2500-10500vr0p05t5750g40k2.dat
rm fort.*
#!

THE SYNTHETIC SPECTRUM f2500-10500vr0t5750g40k2.dat 
IS IN BINARY FORM Hnu versus lambda in nm. Its resolution is 
500000.  TO DEGRADE IT TO 1 A RESOLUTION AND CONVERT IT TO ASCII
as Flamda versus lambda in A:


ln -s f2500-10500vr0p05t5750g40k2.dat fort.1
./fluxaverage1a_nmtoa.exe<f2500-10500p05t5750g40k2_1a.out
 250.0    1050.0
 250.0    250.1     250.05
EOF
mv fort.2 f2500-10500p05t5750g40k2_1ap.txt
rm fort.*


In the above command file there are the codes:

xnfpelsyn.for--xnfpelsyn is linked with atlas7v.
 It computes atomic and molecular number densities and continous opacities



The output model from ATLAS9 
ap05t5750g40k2odfnew.dat MUST BE MODIFIED in ap05t57500g40k2odfnew.mod

1) ADD the following control cards at the top of the 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



synbeg.for --initializes the synthetic spectrum computations

rgfalllinesnew.for--reads the atomic line lists gf****.100

rpredict.for -- reads lines from predicted levels from the binary file "fclowlines.bin" 
 
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 file)

rh2ofast.for--reads the H2O line list h2ofastfix.bin in binary format

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 --it computes the emergent flux for different rotational velocities vsini

broaden.for --broadens the emergent flux for macroturbolence or instrumental effects.


fluxaverage1a_nmtoa.for convets from binary to ascii, from nm to A, and degrades from 500000 resolving power to 1A resolution.