Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/230/16 into FORTRAN code for reading data files line by line.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-19
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/ApJS/230/16 05 through L3 empirical stellar spectra from SDSS (Kesseli+, 2017)
*================================================================================
*An empirical template library of stellar spectra for a wide range of spectral
*classes, luminosity classes, and metallicities using SDSS BOSS spectra.
* Kesseli A.Y., West A.A., Veyette M., Harrison B., Feldman D.,
* Bochanski J.J.
* <Astrophys. J. Suppl. Ser., 230, 16-16 (2017)>
* =2017ApJS..230...16K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! The template spectra
integer*4 nr__
parameter (nr__=323) ! Number of records
character*198 ar__ ! Full-size record
character*18 FileName ! Name of the template file in subdirectory fits
integer*4 Nsp ! [1/590] Number of coadded spectra
real*4 VO7912 ! Spectral index, flux in 7900.0-7980.0{AA}
* divided by 8100.0-8150.0{AA}
real*4 NaD ! Spectral index, flux in 5880.0-5905.0{AA}
* divided by 5910.0-5935.0{AA}
real*4 CaH2 ! Spectral index, flux in 6814.0-6845.0{AA}
* divided by 7042.0-7046.0{AA}
real*4 MgI ! Spectral index, flux in 5152.7-5192.7{AA}
* divided by 5100.0-5150.0{AA}
real*4 Gband ! Spectral index, flux in 4285.0-4315.0{AA}
* divided by 4260.0-4285.0{AA}
real*4 Hgam ! Spectral index, flux in 4332.5-4347.5{AA}
* divided by 4355.0-4370.0{AA}
real*4 FeH ! Spectral index, flux in 9880.0-10000.0{AA}
* divided by 9820.0-9860.0{AA}
real*4 Rb_B ! Spectral index, flux in 7922.6-7932.6 and
* 7962.6-7972.6{AA} divided by 7942.6-7952.6{AA}
real*4 FeI8689 ! Spectral index, flux in 8684.0-8694.0{AA}
* divided by 8664.0-8674.0{AA}
real*4 VO7434 ! Spectral index, flux in 7430.0-7470.0{AA}
* divided by 7550.0-7570.0{AA}
real*4 Hbeta ! Spectral index, flux in 4847.0-4877.0{AA}
* divided by 4817.0-4847.0{AA}
real*4 FeI4383 ! Spectral index, flux in 4378.6-4388.6{AA}
* divided by 4355.0-4370.0{AA}
real*4 TiO5 ! Spectral index, flux in 7126.0-7135.0{AA}
* divided by 7042.0-7046.0{AA}
real*4 CrH_A ! Spectral index, flux in 8580.0-8600.0{AA}
* divided by 8621.0-8641.0{AA}
real*4 Halpha ! Spectral index, flux in 6548.0-6578.0{AA}
* divided by 6583.0-6613.0{AA}
real*4 VO_B ! Spectral index, flux in 7860.0-7880.0 and
* 8080.0-8100.0{AA} divided by 7960.0-8000.0{AA}
real*4 CaII8662 ! Spectral index, flux in 8650.0-8675.0{AA}
* divided by 8625.0-8650.0{AA}
real*4 CaH3 ! Spectral index, flux in 6960.0-6990.0{AA}
* divided by 7042.0-7046.0{AA}
real*4 CaI6162 ! Spectral index, flux in 6150.0-6175.0{AA}
* divided by 6120.0-6145.0{AA}
real*4 TiO8440 ! Spectral index, flux in 8440.0-8470.0{AA}
* divided by 8400.0-8420.0{AA}
real*4 CaII8498 ! Spectral index, flux in 8483.0-8513.0{AA}
* divided by 8513.0-8543.0{AA}
real*4 Hdel ! Spectral index, flux in 4086.7-4116.7{AA}
* divided by 4136.7-4176.7{AA}
real*4 VO7445 ! Spectral index, flux in 7350.0-7400.0 and
* 7510.0-7560.0{AA} divided by 7420.0-7470.0{AA}
real*4 NaI ! Spectral index, flux in 8177.0-8201.0{AA}
* divided by 8151.0-8175.0{AA}
real*4 CaK ! Spectral index, flux in 3923.7-3943.7{AA}
* divided by 3943.7-3953.7{AA}
real*4 CaI4227 ! Spectral index, flux in 4216.7-4236.7{AA}
* divided by 4236.7-4256.7{AA}
real*4 FeI4404 ! Spectral index, flux in 4399.8-4409.8{AA}
* divided by 4414.8-4424.8{AA}
real*4 Cs_A ! Spectral index, flux in 8496.1-8506.1 and
* 8536.1-8546.1, 0.5{AA} divided
* by 8516.1-8526.1{AA}
real*4 TiO8 ! Spectral index, flux in 8400.0-8415.0{AA}
* divided by 8455.0-8470.0{AA}
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Photometry of templates
integer*4 nr__1
parameter (nr__1=329) ! Number of records
character*99 ar__1 ! Full-size record
character*2 SpT ! MK spectral type (O5 through L6)
real*4 v_Fe_H_ ! [-2/1]? Metallicity bin (1)
character*3 Lum ! Luminosity class (III or V)
real*4 u_g ! (mag) [-0.5/3]? template u-g color index
real*8 u_grms ! (mag) [0/0.3]?=-999 RMS scatter in u-g template
* color index
real*8 e_u_g ! (mag) [0/0.2]?=-999 Uncertainty in u-g
real*4 g_r ! (mag) [-0.6/1.7]? template g-r color index
real*8 g_rrms ! (mag) [0/0.3]?=-999 RMS scatter in g-r template
* color index
real*8 e_g_r ! (mag) [0/0.05]?=-999 Uncertainty in g-r
real*4 r_i ! (mag) [-0.4/2.7]? template r-i color index
real*8 r_irms ! (mag) [0/0.2]?=-999 RMS scatter in r-i template
* color index
real*8 e_r_i ! (mag) [0/0.06]?=-999 Uncertainty in r-i
real*4 i_z ! (mag) [-0.4/2.3]? template i-z color index
real*8 i_zrms ! (mag) [0/0.08]?=-999 RMS scatter in i-z template
* color index
real*8 e_i_z ! (mag) [0/0.09]?=-999 Uncertainty in i-z
*Note (1): All bins are 0.5dex with the number reported denoting the highest
* value in the bin (e.g., 0.0 means the metallicity is between -0.5 and
* 0.0dex compared to solar metallicity).
C=============================================================================
C Loading file 'table2.dat' ! The template spectra
C Format for file interpretation
1 format(
+ A18,1X,I3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,323
read(1,'(A198)')ar__
read(ar__,1)
+ FileName,Nsp,VO7912,NaD,CaH2,MgI,Gband,Hgam,FeH,Rb_B,FeI8689,
+ VO7434,Hbeta,FeI4383,TiO5,CrH_A,Halpha,VO_B,CaII8662,CaH3,
+ CaI6162,TiO8440,CaII8498,Hdel,VO7445,NaI,CaK,CaI4227,FeI4404,
+ Cs_A,TiO8
c ..............Just test output...........
write(6,1)
+ FileName,Nsp,VO7912,NaD,CaH2,MgI,Gband,Hgam,FeH,Rb_B,FeI8689,
+ VO7434,Hbeta,FeI4383,TiO5,CrH_A,Halpha,VO_B,CaII8662,CaH3,
+ CaI6162,TiO8440,CaII8498,Hdel,VO7445,NaI,CaK,CaI4227,FeI4404,
+ Cs_A,TiO8
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Photometry of templates
C Format for file interpretation
2 format(
+ A2,1X,F4.1,1X,A3,1X,F5.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,
+ F7.2,1X,F5.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,F7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,329
read(1,'(A99)')ar__1
read(ar__1,2)
+ SpT,v_Fe_H_,Lum,u_g,u_grms,e_u_g,g_r,g_rrms,e_g_r,r_i,r_irms,
+ e_r_i,i_z,i_zrms,e_i_z
if(ar__1(4:7) .EQ. '') v_Fe_H_ = rNULL__
if(ar__1(13:17) .EQ. '') u_g = rNULL__
if(ar__1(35:39) .EQ. '') g_r = rNULL__
if(ar__1(57:61) .EQ. '') r_i = rNULL__
if(ar__1(79:83) .EQ. '') i_z = rNULL__
c ..............Just test output...........
write(6,2)
+ SpT,v_Fe_H_,Lum,u_g,u_grms,e_u_g,g_r,g_rrms,e_g_r,r_i,r_irms,
+ e_r_i,i_z,i_zrms,e_i_z
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end