4 August 2000 Programs for computing distribution functions This work is in collaboration with Fiorella Castelli in Trieste. The program DFSYNTHE is a modified version of the spectrum synthesis program SYNTHE. It can be run for arbitrary abundances (except that it has not yet been tested for pure H or pure He). Anyone who wants to use this program must go through the fortran to understand how it works. It is easy for users to modify DFSYNTHE for their own needs. This sample is set up for temperatures from 2000K to 200000K. If only temperatures for early type stars are needed, all of the molecules can be left out, for example. At one point in DFSYNTHE.FOR, for each temperature and pressure, there is an array that contains the whole 3507859 point line opacity spectrum. Instead of computing distribution functions, the user can keep that whole spectrum or process it any way she or he likes. Also, it would be possible to add the line and continuum absorption and scattering to produce total absorption and total scattering. At present the tabulation is in T and P. In order to do T and RHO it is necessary to compute the P that corresponds to the desired RHO. I will write a program to do that eventually. This sample run took 57 hours CPU/ 93 hours elapsed time on a 600Mhz Alpha Digital Personal Workstation running OpenVMS. The large df* files are saved here to illustrate the whole sequence of operations. They are normally deleted. Article on computing line opacity RAPIDOPACITY.PS;1 815/828 28-JUL-2000 23:08:29.18 Preprocess line data files diatomics edited down to 7771848 lines REPACKDI.COM;10 1/18 26-JUN-2000 17:21:46.19 REPACKDI.FOR;8 24/36 29-OCT-1997 00:01:47.21 REPACKDI.LOG;14 50/54 26-JUN-2000 17:21:54.82 H2O 48999850 lines REPACKH2O.COM;4 1/18 4-JAN-2000 17:05:04.88 REPACKH2O.FOR;11 5/18 22-JAN-2000 11:24:39.29 REPACKH2O.LOG;3 20/36 22-JAN-2000 11:24:57.52 high atomic ions 10057571 lines REPACKHI.COM;2 1/18 17-SEP-1997 15:03:44.73 REPACKHI.FOR;5 20/36 11-MAR-1998 16:36:17.01 REPACKHI.LOG;3 43/54 31-MAR-1998 14:35:37.68 low atomic ions 31108567 lines REPACKLOW.COM;1 1/18 17-SEP-1997 14:45:45.55 REPACKLOW.FOR;13 20/36 29-OCT-1997 12:23:43.42 REPACKLOW.LOG;5 44/54 28-OCT-1997 22:23:51.81 NLTE lines (in LTE) 38297 lines REPACKNLTE.COM;6 1/18 28-JUL-2000 22:25:24.84 REPACKNLTE.FOR;10 30/36 23-OCT-1997 17:35:54.60 REPACKNLTE.LOG;12 15613/15624 28-JUL-2000 22:25:35.15 NLTELINES.DAT;49 10250/10260 20-MAY-2000 21:33:30.36 TiO 36979284 lines REPACKTIO.COM;9 1/18 2-JAN-2000 19:09:48.53 REPACKTIO.FOR;11 5/18 4-JAN-2000 16:50:09.86 REPACKTIO.LOG;13 20/36 4-JAN-2000 15:59:53.96 Program for pretabulating number densities and continuum opacities XNFDF.FOR;52 1736/1746 28-JUL-2000 16:37:22.15 Sample run The temperatues and pressures are read in as artificial models Binary switches tell which depths to compute Abundances are the same as those in previous calculations at the San Diego Supercomputer Center XNFDFTEST.COM;46 281/288 25-JUL-2000 21:36:07.20 XNFDFTEST.LOG;2 49762/49770 28-JUL-2000 16:41:21.08 CONTINUA.DAT;6 9/18 3-FEB-1998 11:27:27.47 MOLECULES.DAT;9 13/18 8-APR-2000 21:04:37.73 Compute line opacity spectrum for 21 pressures at 56 temperatures and compute distribution functions from the spectrum DFP00.COM;63 3/18 30-JUL-2000 14:00:40.19 DFP00.LOG;46 3581/3582 30-JUL-2000 14:02:17.65 DFSYNTHE.FOR;412 292/306 28-JUL-2000 22:52:54.32 The output for temperature. Note elapsed times. DFT01.DAT;2 10107/10116 30-JUL-2000 14:02:19.97 DFT02.DAT;1 10107/10116 30-JUL-2000 23:15:26.24 DFT03.DAT;1 10107/10116 31-JUL-2000 06:27:26.60 DFT04.DAT;1 10107/10116 31-JUL-2000 12:56:30.62 DFT05.DAT;1 10107/10116 31-JUL-2000 18:49:41.98 DFT06.DAT;1 10107/10116 1-AUG-2000 00:11:54.71 DFT07.DAT;1 10107/10116 1-AUG-2000 05:08:00.76 DFT08.DAT;1 10107/10116 1-AUG-2000 09:38:32.28 DFT09.DAT;1 10107/10116 1-AUG-2000 13:45:11.76 DFT10.DAT;1 10107/10116 1-AUG-2000 17:27:08.08 DFT11.DAT;1 10107/10116 1-AUG-2000 20:46:30.41 DFT12.DAT;1 10107/10116 1-AUG-2000 23:43:50.89 DFT13.DAT;1 10107/10116 2-AUG-2000 02:17:10.31 DFT14.DAT;1 10107/10116 2-AUG-2000 04:25:38.28 DFT15.DAT;1 10107/10116 2-AUG-2000 06:10:42.57 DFT16.DAT;1 10107/10116 2-AUG-2000 07:37:43.84 DFT17.DAT;1 10107/10116 2-AUG-2000 08:51:14.97 DFT18.DAT;1 10107/10116 2-AUG-2000 09:54:56.76 DFT19.DAT;1 10107/10116 2-AUG-2000 10:51:43.22 DFT20.DAT;1 10107/10116 2-AUG-2000 11:43:12.28 DFT21.DAT;1 10107/10116 2-AUG-2000 12:22:43.36 DFT22.DAT;1 10107/10116 2-AUG-2000 13:02:50.17 DFT23.DAT;1 10107/10116 2-AUG-2000 13:45:25.54 DFT24.DAT;1 10107/10116 2-AUG-2000 14:32:14.95 DFT25.DAT;1 10107/10116 2-AUG-2000 15:22:19.32 DFT26.DAT;1 10107/10116 2-AUG-2000 16:14:49.70 DFT27.DAT;1 10107/10116 2-AUG-2000 17:11:49.95 DFT28.DAT;1 10107/10116 2-AUG-2000 18:10:12.64 DFT29.DAT;1 10107/10116 2-AUG-2000 19:10:02.87 DFT30.DAT;1 10107/10116 2-AUG-2000 20:09:44.43 DFT31.DAT;1 10107/10116 2-AUG-2000 21:08:31.13 DFT32.DAT;1 10107/10116 2-AUG-2000 22:04:45.90 DFT33.DAT;1 10107/10116 2-AUG-2000 22:57:40.20 DFT34.DAT;1 10107/10116 2-AUG-2000 23:48:02.93 DFT35.DAT;1 10107/10116 3-AUG-2000 00:37:43.30 DFT36.DAT;1 10107/10116 3-AUG-2000 01:28:39.85 DFT37.DAT;1 10107/10116 3-AUG-2000 02:20:15.53 DFT38.DAT;1 10107/10116 3-AUG-2000 03:08:32.14 DFT39.DAT;1 10107/10116 3-AUG-2000 03:53:06.87 DFT40.DAT;1 10107/10116 3-AUG-2000 04:33:29.24 DFT41.DAT;1 10107/10116 3-AUG-2000 05:04:40.17 DFT42.DAT;1 10107/10116 3-AUG-2000 05:29:37.41 DFT43.DAT;1 10107/10116 3-AUG-2000 05:52:53.08 DFT44.DAT;1 10107/10116 3-AUG-2000 06:14:02.20 DFT45.DAT;1 10107/10116 3-AUG-2000 06:34:49.92 DFT46.DAT;1 10107/10116 3-AUG-2000 06:55:32.01 DFT47.DAT;1 10107/10116 3-AUG-2000 07:16:09.09 DFT48.DAT;1 10107/10116 3-AUG-2000 07:36:42.61 DFT49.DAT;1 10107/10116 3-AUG-2000 07:57:01.32 DFT50.DAT;1 10107/10116 3-AUG-2000 08:17:09.15 DFT51.DAT;1 10107/10116 3-AUG-2000 08:37:06.87 DFT52.DAT;1 10107/10116 3-AUG-2000 08:56:58.76 DFT53.DAT;1 10107/10116 3-AUG-2000 09:16:37.88 DFT54.DAT;1 10107/10116 3-AUG-2000 09:36:04.49 DFT55.DAT;1 10107/10116 3-AUG-2000 09:55:10.29 DFT56.DAT;1 10107/10116 3-AUG-2000 10:14:00.81 Concatenates the above dft*.dat files into dfp00.dat Sorts dfp00.dat over various indices to put dfp00.srt in order of use DFSORT.COM;8 1/18 3-AUG-2000 15:46:27.66 DFSORT.LOG;14 14/18 3-AUG-2000 15:46:55.45 DFP00.DAT;2 565950/565956 3-AUG-2000 15:46:55.95 DFP00.SRT;2 565950/565956 3-AUG-2000 15:56:21.09 Divides dfp00.srt into 5 vturbs and into big and little wavelength intervals as in the previous calculations SEPARATEDF.COM;5 1/18 3-AUG-2000 22:39:20.66 SEPARATEDF.FOR;3 3/18 9-DEC-1997 14:41:58.15 SEPARATEDF.LOG;5 14/18 3-AUG-2000 22:39:47.97 NEWBDFP00BIG0.DAT;2 18081/18090 3-AUG-2000 22:39:51.03 NEWBDFP00BIG1.DAT;2 18081/18090 3-AUG-2000 22:42:28.66 NEWBDFP00BIG2.DAT;2 18081/18090 3-AUG-2000 22:45:19.50 NEWBDFP00BIG4.DAT;2 18081/18090 3-AUG-2000 22:49:06.86 NEWBDFP00BIG8.DAT;2 18081/18090 3-AUG-2000 22:53:12.71 NEWBDFP00LIT0.DAT;3 66812/66816 3-AUG-2000 22:40:20.29 NEWBDFP00LIT1.DAT;2 66812/66816 3-AUG-2000 22:43:04.45 NEWBDFP00LIT2.DAT;2 66812/66816 3-AUG-2000 22:45:56.86 NEWBDFP00LIT4.DAT;2 66812/66816 3-AUG-2000 22:49:59.05 NEWBDFP00LIT8.DAT;2 66812/66816 3-AUG-2000 22:54:00.19 I still have to update program KAPPA9 to generate Rosseland or Planck mean opacities using these distribution functions.