|The programs described on this page are IDL utilities to provide a useful interface to the library of TLUSTY runs performed on the APS nonLTE PC cluster. Before using these programs you should understand those model grids, and their limitations.|
|Conditions of use|
|The programs and model atmosphere results that compose this suite of tools are Copyright Robert Ryans, and the Queen's University of Belfast. Any use of these tools must be by agreement with the author, and subject to suitable recognition in any publication which makes use of them.|
To run the programs you need to set some environment variables so that IDL knows
where to find the data and program source codes. Add the lines below to your
setenv TLABUND_DATADIR /scratch/rsir/nlte/tlabund
tlgui is an easy to use interactive GUI tool for accessing the grids of
nonLTE line equivalent widths, much in the spirit of the old Xabund tool. It is
suitable for interactive analysis of a small number of lines, exploring the
datasets, and some debugging of the datasets themselves.
To start the program, type tlgui in an IDL terminal window. A GUI control panel will appear as shown below:
The layout of the panel should be intuitive to people familar with model atmosphere analysis work. In summary:
When you start tlgui, and run a calculation, you will probably see output like this in the terminal window:
IDL> tlgui TLSetup complete. T=20000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew=150.0 grid=NLGAL Linear abund: 6.96 Poly. abund: 6.96
The output summarises the model that has been run, giving model atmosphere paramters, the line selected, equivalent width, and grid. The second line gives two abundance measurements, one linear and one polynomial. The codes use two types of interpolation (or extrapolation) as a sanity check - if these do not agree, there is a serious problem with your calculation, and you'll probably have seen one of the error messages below.
tlcli is a command line interactive interface to the model codes, similar
in spirit to the LTE spectrum code.
Start the program by typing tlcli - which will produce output like
IDL> tlcli TLSetup complete. Initialising line lists. Using terminal input. Type exit to quit. Enter element, wavelength, ew. Example: si 4128 55 >>>From here you can then start to set the model atmosphere parameters, select the metallicity grid and type, and - of course - enter lines to analyse. Commands in tlcli are:
>>> si 4128 25 You entered 4128.00, setting line wavelength to 4128.05 T=25000 log g=3.50 vt=10 Atom: Si II 4128.05 ew= 25.0 grid=NLGAL Linear abund: 7.90 Poly. abund: 7.95 >>> si 4200 25 *** No line found within 1.0 angstroms of specified wavelength 4200.00 The nearest is 4212.00
tlbatch is a batch mode program, which runs a list of lines through a
specified model atmosphere. Output can be logged to a file for later
The format of the input file is quite simple, just one line per line of the
file, specified the same way as in tlcli, i.e.
si 4128 55 si 4131 67 mg 4481 250This is saved to a simple plain text file.
The syntax of the tlbatch command is:
The name of the model grid is selected from nlgal, nllmc, nlsmc, nlmcb, ltgal,
ltlmc, ltsmc, ltmcb, as in tlcli. input_file is the name of the
file the linelist is saved to. An example:
IDL> tlbatch, 20000, 4.00, 5, 'nlgal', 'test.lis' TLSetup complete. Initialising line lists. Using input from file: test.lis You entered 4128.00, setting line wavelength to 4128.05 T=20000 log g=4.00 vt= 5 Atom: Si II 4128.05 ew= 55.0 grid=NLGAL Linear abund: 7.47 Poly. abund: 7.47 You entered 4131.00, setting line wavelength to 4130.89 T=20000 log g=4.00 vt= 5 Atom: Si II 4130.89 ew= 67.0 grid=NLGAL Linear abund: 7.44 Poly. abund: 7.44 You entered 4481.00, setting line wavelength to 4481.00 T=20000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew=250.0 grid=NLGAL Linear abund: 7.57 Poly. abund: 7.57Note the quoting used for the two final strings on the command line.
Extra command options for tlbatch are:
IDL> tlbatch, 20000, 4.00, 5, 'nlgal', 'test.lis', /extrap, log='log.txt'
|To be added.|
The error messages you're most likely to see involve the model you are trying to
run lying outside the bounds of the grids that have been calculated. There are
two cases - outside the limits of the equivalent widths, or outside the T/g
*** EW is more than maximum in model. EW=400.0 maxEW=334.9 *** To calculate this value, run with the extrapolate flag set! T=20000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew=400.0 grid=NLGAL *** No abundance returned *** EW is less than minimum in model. EW= 60.0 minEW= 76.7 *** To calculate this value, run with the extrapolate flag set! T=20000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew= 60.0 grid=NLGAL *** No abundance returnedIf you absolutely need the result, you can click the Extrapolate button, which will force the calculation through. In these cases it is vital to check for consistency between the linear and polynomial abundances. You should also check how far away from the grid values you are by comparing your EW to the grid limits displayed by the code when the extrapolate button is pressed. Clearly, the further you go from the models, the riskier.
*** EW is more than maximum in model. Extrapolating. EW=400.0 maxEW=334.9 T=20000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew=400.0 grid=NLGAL Linear abund: 8.42 Poly. abund: 8.40 *** EW is less than minimum in model. Extrapolating. EW= 60.0 minEW= 76.7 T=20000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew= 60.0 grid=NLGAL Linear abund: 6.33 Poly. abund: 6.32
*** Temperature outside bounds of grid. T, Tmin, Tmax = 36000 13500 35000 *** Bailing out at lower VT - grid bounds exceeded or file not found. T=36000 log g=4.00 vt= 5 Atom: Mg II 4481.00 ew=100.0 grid=NLGAL *** No abundance returned *** Gravity outside grid limits. g, gmin, gmax = 2.00 2.75 4.50 *** Unable to complete due to log g exceeding grid bounds, or file error *** Bailing out at lower VT - grid bounds exceeded or file not found. T=25000 log g=2.00 vt= 5 Atom: Mg II 4481.00 ew=100.0 grid=NLGAL *** No abundance returned
File not found
Features in emission
*** Warning - some features in emission! *** EW is more than maximum in model. EW= 50.0 maxEW= 0.0 *** To calculate this value, run with the extrapolate flag set! T=32500 log g=3.25 vt=10 Atom: Si IV 4631.00 ew= 50.0 grid=NLGAL *** No abundance returnedIn the example above no abundance is returned, but this is not always the case. Each model grid point is calculated at 5 abundances, and those 'in emission' are automatically set to zero. If sufficient other points exist to allow the calculation to be performed, a result may be returned. However, it would be prudent in such a case to use the verbose mode to make sure it is a correct return. Clearly, you really do need to know what you're doing in these cases...
Multiple error messages