Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/153/174 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-20
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/AJ/153/174 SDSS-III/APOGEE. II. The multi-epoch ABE sample (Chojnowski+, 2017)
*================================================================================
*High-resolution H-band spectroscopy of Be stars with SDSS-III/APOGEE.
*II. Line profile and radial velocity variability.
* Chojnowski S.D., Wisniewski J.P., Whelan D.G., Labadie-Bartz J.,
* Borges Fernandes M., Lin C.-C., Majewski S.R., Stringfellow G.S.,
* Mennickent R.E., Roman-Lopes A., Tang B., Hearty F.R., Holtzman J.A.,
* Pepper J., Zasowski G.
* <Astron. J., 153, 174-174 (2017)>
* =2017AJ....153..174C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! *Additional new Be stars from APOGEE1
integer*4 nr__
parameter (nr__=3) ! Number of records
character*45 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 ABE ! Identifier of the APOGEE B-type emission-line
* (ABE) star
real*4 Hmag ! (mag) H-band magnitude
integer*4 Ns ! Number of spectra
character*8 Morph ! Morphology (1)
integer*4 RAh ! (h) Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the Declination
integer*4 DEd ! (deg) Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
*Note (1): "weak em." means double-peaked emission that is fully contained within
* broad photospheric absorption.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *APOGEE B-type emission-line (ABE) visit summary
integer*4 nr__1
parameter (nr__1=213) ! Number of records
character*101 ar__1 ! Full-size record
character*3 ABE_1 ! Identifier of the Apache Point Observatory
* Galactic Evolution Experiment (APOGEE)
* B-type emission-line (ABE) star (running
* number) (G1)
character*17 Name ! Star name
integer*4 Nobs ! [2/30] Number of APOGEE observations
integer*4 MJD0 ! (d) [55702/56812] Modified Julian Date (MJD) of
* first APOGEE observation (JD-2400000.5)
integer*4 dJD ! (d) [1/1862] Observational baseline ({Delta}JD)
integer*4 v_S_N_ ! [46/800] Mean Signal-to-Noise ratio
character*8 VType ! Variability classification (1)
integer*4 v_PSep_ ! (km/s) [0/629]? Mean Br11 peak separation (G2)
integer*4 PAmp ! (km/s) [0/128]? Scatter (maximum minus minimum)
* in Br11 peak separation
real*4 v_EW_ ! (0.1nm) [-17/6] Mean Br11 equivalent width W_Br11_
* (in {AA}) (G3)
real*4 e__EW_ ! (0.1nm) [0/2.93] Standard deviation ({sigma}) in <EW>
* (in {AA})
real*4 EWAmp ! (0.1nm) [0.01/8.37] Scatter (maximum minus minimum)
* in Br11 equivalent width (in {AA})
real*4 v_RV_ ! (km/s) [-93/76.6]? Mean radial velocity (G4)
real*4 e__RV_ ! (km/s) [0.2/14.4]? Mean standard deviation ({sigma})
* of line-by-line radial velocities
real*4 RVAmp ! (km/s) [0/112.9]? Scatter (maximum minus minimum)
* in radial velocity
real*4 o__RV_ ! [3/9]? Mean number of H-Br lines used for
* radial velocity determination (G5)
character*6 RVType ! Type of radial velocity measurement, as
* described in Section 3.4 (double, shell, or
* single) (G6)
*Note (1): Variability classification is defined as follows:
* RV = The scatter in radial velocity measurements exceeds twice the average
* standard deviation of individual H-Br lines measurements. A total of
* 37 stars are classified as RV-variable;
* W = Either the standard deviation of W_Br11_ measurements was greater
* than 0.65{AA}, or else that visual inspection clearly indicates
* variable emission strength despite a small standard deviation in
* multi-epoch WBr11 measurements. A total of 30 stars are classified as
* W_Br11_-variable;
* V/R = The ratio of intensities of V and R peaks is variable. The
* classification is assigned based on visual inspection of the spectra
* to confirm that the V/R ratio changed in multiple H-Br lines. A total
* of 16 stars are classified as V/R-variable;
* TD = 'Transient Disk', and indicates the APOGEE spectra include at least
* one mostly emission-less spectrum. The classification is assigned
* based on visual inspection of the spectra. A total of 18 stars are
* classified as having transient disks.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! *APOGEE B-type emission-line (ABE) star summary
integer*4 nr__2
parameter (nr__2=1101) ! Number of records
character*110 ar__2 ! Full-size record
character*3 ABE_2 ! Identifier of the Apache Point Observatory
* Galactic Evolution Experiment (APOGEE) B-type
* emission-line (ABE) star (running number) (G1)
character*17 Name_1 ! Star name
character*18 Visit ! Visit identifier (ap#-Plate-MJD-Fiber) (1)
character*21 Date ! UT date of mid-observation (ISO 8601 date time)
real*8 MJD ! (d) [55701.9/57701.6] Modified Julian Date (MJD) of
* mid-observation (JD-2400000.5)
integer*4 S_N ! [41/1427] APOGEE spectra Signal-to-Noise ratio
integer*4 PSep ! (km/s) [0/629]? Br11 peak separation (G2)
real*4 EW ! (0.1nm) [-18.4/6.4]? Br11 equivalent width W_Br11_
* (in {AA}) (G3)
real*4 RV ! (km/s) [-95.3/98.3]? Radial velocity derived from the
* H-Br lines (G4)
real*4 e_RV ! (km/s) [0.1/22.9]? Standard deviation ({sigma}) of
* line-by-line radial velocities
integer*4 o_RV ! [3/9]? Number of H-Br lines used for RV (G5)
character*6 RVType_1 ! Type of radial velocity measurement, as
* described in Section 3.4 (single, double, or
* shell) (G6)
*Note (1): The visit identifier, which begins with "ap1" for APOGEE1 spectra and
* "ap2" for APOGEE2 spectra, and which is followed by the plug-plate number,
* the observation MJD (JD-2400000), and the fiber number.
C=============================================================================
C Loading file 'table1.dat' ! *Additional new Be stars from APOGEE1
C Format for file interpretation
1 format(
+ A3,1X,F6.3,1X,I1,1X,A8,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
+ F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,3
read(1,'(A45)')ar__
read(ar__,1)ABE,Hmag,Ns,Morph,RAh,RAm,RAs,DE_,DEd,DEm,DEs
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)ABE,Hmag,Ns,Morph,RAh,RAm,RAs,DE_,DEd,DEm,DEs
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! *APOGEE B-type emission-line (ABE) visit summary
C Format for file interpretation
2 format(
+ A3,1X,A17,1X,I2,1X,I5,1X,I4,1X,I3,1X,A8,1X,I3,1X,I3,1X,F6.2,
+ 1X,F4.2,1X,F4.2,1X,F5.1,1X,F4.1,1X,F5.1,1X,F3.1,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,213
read(1,'(A101)')ar__1
read(ar__1,2)
+ ABE_1,Name,Nobs,MJD0,dJD,v_S_N_,VType,v_PSep_,PAmp,v_EW_,
+ e__EW_,EWAmp,v_RV_,e__RV_,RVAmp,o__RV_,RVType
if(ar__1(50:52) .EQ. '') v_PSep_ = iNULL__
if(ar__1(54:56) .EQ. '') PAmp = iNULL__
if(ar__1(75:79) .EQ. '') v_RV_ = rNULL__
if(ar__1(81:84) .EQ. '') e__RV_ = rNULL__
if(ar__1(86:90) .EQ. '') RVAmp = rNULL__
if(ar__1(92:94) .EQ. '') o__RV_ = rNULL__
c ..............Just test output...........
write(6,2)
+ ABE_1,Name,Nobs,MJD0,dJD,v_S_N_,VType,v_PSep_,PAmp,v_EW_,
+ e__EW_,EWAmp,v_RV_,e__RV_,RVAmp,o__RV_,RVType
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! *APOGEE B-type emission-line (ABE) star summary
C Format for file interpretation
3 format(
+ A3,1X,A17,1X,A18,1X,A21,1X,F11.5,1X,I4,1X,I3,1X,F6.2,1X,F5.1,
+ 1X,F4.1,1X,I1,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,1101
read(1,'(A110)')ar__2
read(ar__2,3)
+ ABE_2,Name_1,Visit,Date,MJD,S_N,PSep,EW,RV,e_RV,o_RV,RVType_1
if(ar__2(81:83) .EQ. '') PSep = iNULL__
if(ar__2(85:90) .EQ. '') EW = rNULL__
if(ar__2(92:96) .EQ. '') RV = rNULL__
if(ar__2(98:101) .EQ. '') e_RV = rNULL__
if(ar__2(103:103) .EQ. '') o_RV = iNULL__
c ..............Just test output...........
write(6,3)
+ ABE_2,Name_1,Visit,Date,MJD,S_N,PSep,EW,RV,e_RV,o_RV,RVType_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end