Conversion of standardized ReadMe file for
file /./home/cats/J/ApJ/752/59 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 2026-May-15
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/ApJ/752/59 Low-mass population in L1641 (Hsu+, 2012)
*================================================================================
*The low-mass stellar population in L1641: evidence for environmental
*dependence of the stellar initial mass function.
* Hsu W.-H., Hartmann L., Allen L., Hernandez J., Megeath S.T., Mosby G.,
* Tobin J.J., Espaillat C.
* <Astrophys. J., 752, 59 (2012)>
* =2012ApJ...752...59H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Confirmed members in L1641
integer*4 nr__
parameter (nr__=865) ! Number of records
character*104 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq (nr__) ! Running sequence number
real*8 RAdeg (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
real*8 Vmag (nr__) ! (mag) [12/23]?=-100 V band magnitude
real*8 Imag (nr__) ! (mag) [11/20]?=-100 I band magnitude
real*8 Jmag (nr__) ! (mag) [9/18]?=-100 2MASS J band magnitude
real*8 Hmag (nr__) ! (mag) [8/17]?=-100 2MASS H band magnitude
real*8 Kmag (nr__) ! (mag) [7/17]?=-100 2MASS Ks band magnitude
character*4 SpT (nr__) ! MK spectral type
real*4 e_SpT (nr__) ! ? Uncertainty in SpT
character*1 r_SpT (nr__) ! Origin of SpT (1)
character*1 IR (nr__) ! [Y/N] IR excess? (2)
real*8 EWHa (nr__) ! (0.1nm) H{alpha} equivalent width
real*4 EWLi (nr__) ! (0.1nm) ? Lithium equivalent width
character*11 Date (nr__) ! ("Y/M/D") ? Date of the observation (3)
*Note (1): Flag as follows:
* c = Data from Hectospec spectra.
* d = Data from IMACS spectra.
* e = Data from Fang et al. (2009, Cat. J/A+A/504/461). No spectral-typing
* error estimate.
*Note (2): Criteria for IR excess are defined in Gutermuth et al. (2009,
* Cat. J/ApJS/184/18) and Megeath et al. (2012AJ....144..192M).
*Note (3): Local date at the beginning of the night.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Probable members in L1641
integer*4 nr__1
parameter (nr__1=102) ! Number of records
character*104 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq_1 (nr__1) ! Running sequence number
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000)
real*8 Vmag_1 (nr__1) ! (mag) [12/23]?=-100 V band magnitude
real*8 Imag_1 (nr__1) ! (mag) [11/20]?=-100 I band magnitude
real*8 Jmag_1 (nr__1) ! (mag) [9/18]?=-100 2MASS J band magnitude
real*8 Hmag_1 (nr__1) ! (mag) [8/17]?=-100 2MASS H band magnitude
real*8 Kmag_1 (nr__1) ! (mag) [7/17]?=-100 2MASS Ks band magnitude
character*4 SpT_1 (nr__1) ! MK spectral type
real*4 e_SpT_1 (nr__1) ! ? Uncertainty in SpT
character*1 r_SpT_1 (nr__1) ! Origin of SpT (1)
character*1 IR_1 (nr__1) ! [Y/N] IR excess? (2)
real*8 EWHa_1 (nr__1) ! (0.1nm) H{alpha} equivalent width
real*4 EWLi_1 (nr__1) ! (0.1nm) ? Lithium equivalent width
character*11 Date_1 (nr__1) ! ("Y/M/D") ? Date of the observation (3)
*Note (1): Flag as follows:
* c = Data from Hectospec spectra.
* d = Data from IMACS spectra.
* e = Data from Fang et al. (2009, Cat. J/A+A/504/461). No spectral-typing
* error estimate.
*Note (2): Criteria for IR excess are defined in Gutermuth et al. (2009,
* Cat. J/ApJS/184/18) and Megeath et al. (2012AJ....144..192M).
*Note (3): Local date at the beginning of the night.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! YSOs with strong emission lines
integer*4 nr__2
parameter (nr__2=161) ! Number of records
character*108 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq_2 (nr__2) ! [6/861] Running sequence number
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination in decimal degrees (J2000)
real*4 W5577 (nr__2) ! (0.1nm) ? Equivalent width of [OI] 5577{AA}
real*4 W6300 (nr__2) ! (0.1nm) ? Equivalent width of [OI] 6300{AA}
real*4 W6363 (nr__2) ! (0.1nm) ? Equivalent width of [OI] 6363{AA}
real*4 W6538 (nr__2) ! (0.1nm) ? Equivalent width of [NII] 6548{AA}
real*4 W6583 (nr__2) ! (0.1nm) ? Equivalent width of [NII] 6583{AA}
real*4 W6678 (nr__2) ! (0.1nm) ? Equivalent width of HeI 6678{AA}
real*4 W6716 (nr__2) ! (0.1nm) ? Equivalent width of [SII] 6716{AA}
real*4 W6731 (nr__2) ! (0.1nm) ? Equivalent width of [SII] 6731{AA}
real*4 W7773 (nr__2) ! (0.1nm) ? Equivalent width of [OI] 7773{AA}
real*4 W8446 (nr__2) ! (0.1nm) ? Equivalent width of [OI] 8446{AA}
real*4 W8498 (nr__2) ! (0.1nm) ? Equivalent width of [CaII] 8498{AA}
real*4 W8542 (nr__2) ! (0.1nm) ? Equivalent width of [CaII] 8542{AA}
real*4 W8662 (nr__2) ! (0.1nm) ? Equivalent width of [CaII] 8662{AA}
C=============================================================================
C Loading file 'table2.dat' ! Confirmed members in L1641
C Format for file interpretation
1 format(
+ I3,1X,F8.5,1X,F8.5,1X,F8.3,1X,F8.3,1X,F8.3,1X,F8.3,1X,F8.3,1X,
+ A4,1X,F3.1,1X,A1,1X,A1,1X,F8.3,1X,F3.1,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,865
read(1,'(A104)')ar__
read(ar__,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),Vmag(i__),Imag(i__),Jmag(i__),
+ Hmag(i__),Kmag(i__),SpT(i__),e_SpT(i__),r_SpT(i__),IR(i__),
+ EWHa(i__),EWLi(i__),Date(i__)
if(ar__(73:75) .EQ. '') e_SpT(i__) = rNULL__
if(ar__(90:92) .EQ. '') EWLi(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),RAdeg(i__),DEdeg(i__),Vmag(i__),Imag(i__),Jmag(i__),
+ Hmag(i__),Kmag(i__),SpT(i__),e_SpT(i__),r_SpT(i__),IR(i__),
+ EWHa(i__),EWLi(i__),Date(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Probable members in L1641
C Format for file interpretation
2 format(
+ I3,1X,F8.5,1X,F8.5,1X,F8.3,1X,F8.3,1X,F8.3,1X,F8.3,1X,F8.3,1X,
+ A4,1X,F3.1,1X,A1,1X,A1,1X,F8.3,1X,F3.1,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,102
read(1,'(A104)')ar__1
read(ar__1,2)
+ Seq_1(i__),RAdeg_1(i__),DEdeg_1(i__),Vmag_1(i__),Imag_1(i__),
+ Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),SpT_1(i__),e_SpT_1(i__),
+ r_SpT_1(i__),IR_1(i__),EWHa_1(i__),EWLi_1(i__),Date_1(i__)
if(ar__1(73:75) .EQ. '') e_SpT_1(i__) = rNULL__
if(ar__1(90:92) .EQ. '') EWLi_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),RAdeg_1(i__),DEdeg_1(i__),Vmag_1(i__),Imag_1(i__),
+ Jmag_1(i__),Hmag_1(i__),Kmag_1(i__),SpT_1(i__),e_SpT_1(i__),
+ r_SpT_1(i__),IR_1(i__),EWHa_1(i__),EWLi_1(i__),Date_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! YSOs with strong emission lines
C Format for file interpretation
3 format(
+ I3,1X,F8.5,1X,F8.5,1X,F6.2,1X,F6.2,1X,F6.2,1X,F4.1,1X,F6.2,1X,
+ F5.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,161
read(1,'(A108)')ar__2
read(ar__2,3)
+ Seq_2(i__),RAdeg_2(i__),DEdeg_2(i__),W5577(i__),W6300(i__),
+ W6363(i__),W6538(i__),W6583(i__),W6678(i__),W6716(i__),
+ W6731(i__),W7773(i__),W8446(i__),W8498(i__),W8542(i__),
+ W8662(i__)
if(ar__2(23:28) .EQ. '') W5577(i__) = rNULL__
if(ar__2(30:35) .EQ. '') W6300(i__) = rNULL__
if(ar__2(37:42) .EQ. '') W6363(i__) = rNULL__
if(ar__2(44:47) .EQ. '') W6538(i__) = rNULL__
if(ar__2(49:54) .EQ. '') W6583(i__) = rNULL__
if(ar__2(56:60) .EQ. '') W6678(i__) = rNULL__
if(ar__2(62:67) .EQ. '') W6716(i__) = rNULL__
if(ar__2(69:74) .EQ. '') W6731(i__) = rNULL__
if(ar__2(76:80) .EQ. '') W7773(i__) = rNULL__
if(ar__2(82:87) .EQ. '') W8446(i__) = rNULL__
if(ar__2(89:94) .EQ. '') W8498(i__) = rNULL__
if(ar__2(96:101) .EQ. '') W8542(i__) = rNULL__
if(ar__2(103:108) .EQ. '') W8662(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Seq_2(i__),RAdeg_2(i__),DEdeg_2(i__),W5577(i__),W6300(i__),
+ W6363(i__),W6538(i__),W6583(i__),W6678(i__),W6716(i__),
+ W6731(i__),W7773(i__),W8446(i__),W8498(i__),W8542(i__),
+ W8662(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end