Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/121/2124 into FORTRAN code for loading all data files into arrays.
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-16
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/121/2124 {lambda} Ori pre-main-sequence stars photometry (Dolan+, 2001)
*================================================================================
*The spatial distribution of the {lambda} Orionis pre-main-sequence population.
* Dolan C.J., Mathieu R.D.
* <Astron. J. 121, 2124 (2001)>
* =2001AJ....121.2124D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Improved data for Paper I (Dolan & Mathieu,
1999, Cat. <J/AJ/118/2409>) stars
integer*4 nr__
parameter (nr__=72) ! Number of records
character*72 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 v_DM99_ (nr__) ! [1/72] The {lambda} Ori running number
character*16 JName (nr__) ! The {lambda} Ori designation
integer*4 RAh (nr__) ! (h) Right Ascension from USNO-A 2.0 catalog (J2000)
integer*4 RAm (nr__) ! (min) Right Ascension from USNO-A 2.0 catalog (J2000)
real*4 RAs (nr__) ! (s) Right Ascension from USNO-A 2.0 catalog (J2000)
character*1 DE_ (nr__) ! Declination sign (J2000)
integer*4 DEd (nr__) ! (deg) Declination from USNO-A 2.0 catalog (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination from USNO-A 2.0 catalog (J2000)
real*4 DEs (nr__) ! (arcsec) Declination from USNO-A 2.0 catalog (J2000)
real*4 Vmag (nr__) ! (mag) The V band magnitude
real*4 V_R (nr__) ! (mag) The (V-R) color
real*4 R_I (nr__) ! (mag) The (R-I) color
real*4 LogAge (nr__) ! ([yr]) Log of the Age (1)
real*4 Mass (nr__) ! (solMass) Mass (1)
character*1 f__DM99_ (nr__) ! [b] Indicates star does not exist in the
* USNO-A2.0 database. The coordinates (A1.0) and
* magnitudes listed are from Paper I,
* Dolan & Mathieu, 1999AJ....118.2409D.
*Note (1): Interpolated from the Palla & Stahler (1999ApJ...525..772P) stellar
* evolution models.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Pre-main-sequence stellar data
integer*4 nr__1
parameter (nr__1=204) ! Number of records
character*106 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 v_DM99__1 (nr__1) ! [73/266] The {lambda} Ori running number
character*16 JName_1 (nr__1) ! The {lambda} Ori designation
integer*4 RAh_1 (nr__1) ! (h) Right Ascension from USNO-A 2.0 catalog (J2000)
integer*4 RAm_1 (nr__1) ! (min) Right Ascension from USNO-A 2.0 catalog (J2000)
real*4 RAs_1 (nr__1) ! (s) Right Ascension from USNO-A 2.0 catalog (J2000)
character*1 DE__1 (nr__1) ! Declination sign (J2000)
integer*4 DEd_1 (nr__1) ! (deg) Declination from USNO-A 2.0 catalog (J2000)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination from USNO-A 2.0 catalog (J2000)
real*4 DEs_1 (nr__1) ! (arcsec) Declination from USNO-A 2.0 catalog (J2000)
real*4 Vmag_1 (nr__1) ! (mag) The V band magnitude
character*1 f_Vmag (nr__1) ! [c] Indicates photometry from the
* lower-quality 1998 data
real*4 V_R_1 (nr__1) ! (mag) The (V-R) color
character*1 f_V_R (nr__1) ! [c] Indicates photometry from the
* lower-quality 1998 data
real*4 R_I_1 (nr__1) ! (mag) The (R-I) color
character*1 f_R_I (nr__1) ! [c] Indicates photometry from the
* lower-quality 1998 data
real*4 EqWidLi (nr__1) ! (0.1nm) ? Equivalent width of Li (1)
character*1 f_EqWidLi (nr__1) ! [gh] Flag on EqWidLi (2)
real*4 EqWidHa (nr__1) ! (0.1nm) ? Equivalent width of H{alpha} (3)
character*1 f_EqWidHa (nr__1) ! [d] Indicates star shows a strong
* P Cygni profile
real*4 RVel (nr__1) ! (km/s) Stellar radial velocity
character*1 f_RVel (nr__1) ! [fij] Flag on RVel (4)
real*4 LogAge_1 (nr__1) ! ([yr]) Log of the Age (5)
real*4 Mass_1 (nr__1) ! (solMass) Mass (5)
character*1 f_Mass (nr__1) ! [e] Indicates star is outside the age/mass grid
*Note (1): Positive values indicate absorption.
*Note (2):
* g = These stars have one EqWidLi above 0.2{AA} and one below.
* h = EqWidLi measured from the combined light of both components.
*Note (3): Positive values indicate emission.
*Note (4):
* f = Double-lined spectrum.
* i = Two observations of a blended double-line spectrum. The first pair are
* JD2451235.67 while the second pair are two days later (JD2451237.68).
* j = The first pair are from a blended double-line spectrum (JD2451564.82).
* A second spectrum about three weeks later (JD2451586.65) was not
* resolvable as two spectra.
*Note (5): Interpolated from the Palla & Stahler (1999ApJ...525..772P)
* stellar evolution models.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Massive star photometry
integer*4 nr__2
parameter (nr__2=20) ! Number of records
character*88 ar__2 ! Full-size record
integer*4 HD (nr__2) ! The HD (Cat. <III/135>) number
integer*4 PPM (nr__2) ! ? The Positions & Proper Motions Catalog
* (Cat. <I/146>) number
character*2 SpType (nr__2) ! Spectral type (1)
character*1 f_SpType (nr__2) ! [b] Flag on SpType (2)
real*4 Vmag_2 (nr__2) ! (mag) The V band magnitude
real*4 e_Vmag (nr__2) ! (mag) Uncertainty in Vmag
real*4 b_y (nr__2) ! (mag) The (b-y) Stroemgren color
real*4 e_b_y (nr__2) ! (mag) Uncertainty in (b-y)
real*4 m1 (nr__2) ! (mag) The m_1_ band magnitude
real*4 e_m1 (nr__2) ! (mag) Uncertainty in m1
real*4 c1 (nr__2) ! (mag) The c_1_ band magnitude
real*4 e_c1 (nr__2) ! (mag) Uncertainty in c1
integer*4 o_b_y (nr__2) ! Number of Stroemgren uvby observations
real*4 beta (nr__2) ! (mag) The {beta} band magnitude
real*4 e_beta (nr__2) ! (mag) Uncertainty in beta
integer*4 o_beta (nr__2) ! Number of beta observations
integer*4 r_HD (nr__2) ! Reference on HD
*Note (1): Spectral types are reported from the PPM catalog, except for
* HD 245185 which is from the HD catalog.
*Note (2): Terranegra et al. (1994, Cat. <J/A+AS/104/557>) find a photometric
* spectral type of B9.5III-V for this Herbig Ae-Be star.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Massive star derived data
integer*4 nr__3
parameter (nr__3=20) ! Number of records
character*35 ar__3 ! Full-size record
integer*4 HD_1 (nr__3) ! The HD (Cat. <III/135>) number
real*4 E_B_V (nr__3) ! (mag) The E(B-V) reddening
real*4 Teff (nr__3) ! ([K]) Log of the effective temperature
real*4 Lum (nr__3) ! ([solLum]) Log of the luminosity in solar units
real*4 Mass_2 (nr__3) ! ([solMass]) Log of the mass in solar units
integer*4 Dist (nr__3) ! (pc) ? Distance
integer*4 e_Dist (nr__3) ! (%) ? The parallax fractional error
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Follow-up data for DIL (Duerr, Imhoff & Lada,
1982ApJ...261..135D) stars
integer*4 nr__4
parameter (nr__4=99) ! Number of records
character*48 ar__4 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 (nr__4) ! (deg) Right Ascension J2000
real*8 DEdeg_2 (nr__4) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 DIL (nr__4) ! The Duerr, Imhoff & Lada (1982ApJ...261..135D)
* survey number
character*13 USNO_A2_0 (nr__4) ! The USNO A2.0 identification
character*1 f_USNO_A2_0(nr__4) ! [abcde] Flag on USNO-A2.0 (1)
integer*4 RAh_2 (nr__4) ! (h) Right Ascension (J2000)
integer*4 RAm_2 (nr__4) ! (min) Right Ascension (J2000)
real*4 RAs_2 (nr__4) ! (s) Right Ascension (J2000)
character*1 DE__2 (nr__4) ! Declination sign (J2000)
integer*4 DEd_2 (nr__4) ! (deg) Declination (J2000)
integer*4 DEm_2 (nr__4) ! (arcmin) Declination (J2000)
real*4 DEs_2 (nr__4) ! (arcsec) Declination (J2000)
integer*4 v_DM99__2 (nr__4) ! ? This paper's identification number
*Note (1):
* a = DIL 27 appears to be part of a very close pair. We have indicated the
* likely match, but the other is at 05:31:24.660 12:12:09.71.
* b = DIL 42 is not marked on the DIL finding charts. We found the star in
* Haro, Iriarte & Chavira (1953BOTT....1h...3H), but there are too few
* stars visible in their finding chart to be absolutely sure this is the
* right one.
* c = DIL 49 and 50 have reversed names in the DIL catalog. DIL 49 should
* be ``H 6-70'' while DIL 50 should be ``H 6-69''. Column 8 of the DIL
* catalog should be reversed for these stars as well.
* d = DIL 70 is unclearly marked on both the DIL and Manova
* (1959SvA.....3..188M) finding charts.
* We have indicated the likely match, but another possibility is at
* 05:41:29.597 09:07:11.21.
* e = DIL 80 is unclearly marked on both the DIL and Manova finding charts.
* We have indicated the likely match, but another possibility is at
* 05:43:57.008 09:16:25.65.
C=============================================================================
C Loading file 'table1.dat' ! Improved data for Paper I (Dolan & Mathieu,
* 1999, Cat. <J/AJ/118/2409>) stars
C Format for file interpretation
1 format(
+ I2,1X,A16,1X,I2,1X,I2,1X,F6.3,A1,I2,1X,I2,1X,F5.2,1X,F5.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,72
read(1,'(A72)')ar__
read(ar__,1)
+ v_DM99_(i__),JName(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vmag(i__),V_R(i__),R_I(i__),
+ LogAge(i__),Mass(i__),f__DM99_(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ v_DM99_(i__),JName(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vmag(i__),V_R(i__),R_I(i__),
+ LogAge(i__),Mass(i__),f__DM99_(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Pre-main-sequence stellar data
C Format for file interpretation
2 format(
+ I3,1X,A16,1X,I2,1X,I2,1X,F6.3,A1,I2,1X,I2,1X,F5.2,1X,F5.2,1X,
+ A1,1X,F4.2,1X,A1,1X,F4.2,1X,A1,1X,F4.2,1X,A1,1X,F6.2,1X,A1,1X,
+ F6.2,1X,A1,1X,F4.2,1X,F4.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,204
read(1,'(A106)')ar__1
read(ar__1,2)
+ v_DM99__1(i__),JName_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Vmag_1(i__),
+ f_Vmag(i__),V_R_1(i__),f_V_R(i__),R_I_1(i__),f_R_I(i__),
+ EqWidLi(i__),f_EqWidLi(i__),EqWidHa(i__),f_EqWidHa(i__),
+ RVel(i__),f_RVel(i__),LogAge_1(i__),Mass_1(i__),f_Mass(i__)
if(ar__1(69:72) .EQ. '') EqWidLi(i__) = rNULL__
if(ar__1(76:81) .EQ. '') EqWidHa(i__) = rNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ v_DM99__1(i__),JName_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Vmag_1(i__),
+ f_Vmag(i__),V_R_1(i__),f_V_R(i__),R_I_1(i__),f_R_I(i__),
+ EqWidLi(i__),f_EqWidLi(i__),EqWidHa(i__),f_EqWidHa(i__),
+ RVel(i__),f_RVel(i__),LogAge_1(i__),Mass_1(i__),f_Mass(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Massive star photometry
C Format for file interpretation
3 format(
+ I6,1X,I6,1X,A2,1X,A1,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,I1,1X,F5.3,1X,F5.3,1X,I1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,20
read(1,'(A88)')ar__2
read(ar__2,3)
+ HD(i__),PPM(i__),SpType(i__),f_SpType(i__),Vmag_2(i__),
+ e_Vmag(i__),b_y(i__),e_b_y(i__),m1(i__),e_m1(i__),c1(i__),
+ e_c1(i__),o_b_y(i__),beta(i__),e_beta(i__),o_beta(i__),
+ r_HD(i__)
if(ar__2(8:13) .EQ. '') PPM(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ HD(i__),PPM(i__),SpType(i__),f_SpType(i__),Vmag_2(i__),
+ e_Vmag(i__),b_y(i__),e_b_y(i__),m1(i__),e_m1(i__),c1(i__),
+ e_c1(i__),o_b_y(i__),beta(i__),e_beta(i__),o_beta(i__),
+ r_HD(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Massive star derived data
C Format for file interpretation
4 format(I6,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.1,1X,I3,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,20
read(1,'(A35)')ar__3
read(ar__3,4)
+ HD_1(i__),E_B_V(i__),Teff(i__),Lum(i__),Mass_2(i__),Dist(i__),
+ e_Dist(i__)
if(ar__3(28:30) .EQ. '') Dist(i__) = iNULL__
if(ar__3(32:35) .EQ. '') e_Dist(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ HD_1(i__),E_B_V(i__),Teff(i__),Lum(i__),Mass_2(i__),Dist(i__),
+ e_Dist(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Follow-up data for DIL (Duerr, Imhoff & Lada,
* 1982ApJ...261..135D) stars
C Format for file interpretation
5 format(
+ I2,1X,A13,1X,A1,1X,I2,1X,I2,1X,F6.3,A1,I2,1X,I2,1X,F5.2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,99
read(1,'(A48)')ar__4
read(ar__4,5)
+ DIL(i__),USNO_A2_0(i__),f_USNO_A2_0(i__),RAh_2(i__),
+ RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),
+ DEs_2(i__),v_DM99__2(i__)
if(ar__4(45:47) .EQ. '') v_DM99__2(i__) = iNULL__
RAdeg_2(i__) = rNULL__
DEdeg_2(i__) = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c ..............Just test output...........
write(6,5)
+ DIL(i__),USNO_A2_0(i__),f_USNO_A2_0(i__),RAh_2(i__),
+ RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),
+ DEs_2(i__),v_DM99__2(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end