Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/147/305 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.5, on 2013-May-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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJS/147/305 IUE atlas of pre-main-sequence stars (Valenti+, 2003)
*================================================================================
*An IUE atlas of pre-main-sequence stars.
*III. Co-added final archive spectra from the long-wavelength cameras.
* Valenti J.A., Fallon A.A., Johns-Krull C.M.
* <Astrophys. J. Suppl. Ser. 147, 305 (2003)>
* =2003ApJS..147..305V
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Pre-main-sequence stars observed by IUE with
LW cameras
integer*4 nr__
parameter (nr__=238) ! Number of records
character*73 ar__ ! Full-size record
integer*4 Seq (nr__) ! Internal identification number
character*1 u_Seq (nr__) ! [)] ) if none of the LW spectra are useful
integer*4 VJL (nr__) ! ? Identification number from Paper I,
* Valenti et al., 2000ApJS..129..399V
character*1 u_VJL (nr__) ! [)] ) if none of the SW spectra are useful
integer*4 HBC (nr__) ! ? Number of Herbig & Bell, 1988LicOB1111....1H
character*15 Name (nr__) ! Name of source
character*16 OName (nr__) ! Other name
character*13 SpType (nr__) ! Spectral type (1)
character*6 Cat (nr__) ! PMS category
character*1 n_Cat (nr__) ! [h] Note on Cat (2)
integer*4 N (nr__) ! Number of useful IUE LW spectra
character*1 sp (nr__) ! [-] when there is no spectrum associated.
*Note (1): Spectral types preferentially from Herbig & Bell (1988LicOB1111....1H)
* for TTS and The et al. (1994, Cat. <J/A+AS/104/315> for HAEBE stars.
*Note (2): h: We analyze as a CTTS, despite a HAEBE classification by
* Gregorio-Hetem et al. (1992AJ....103..549G).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Additional classifications from the literature
integer*4 nr__1
parameter (nr__1=22) ! Number of records
character*76 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! Internal identification number
character*13 Name_1 (nr__1) ! Name of source
character*6 Cat_1 (nr__1) ! PMS category
character*1 n_Cat_1 (nr__1) ! [b] Note on Cat (1)
integer*4 r_Cat (nr__1) ! References, in refs.dat file
character*47 Reasons (nr__1) ! Reasons of the classification
*Note (1): b: Status is uncertain as Hillenbrand et al. (1998,
* Cat. <J/AJ/116/1816> detect weak emission in the Ca II 8542{AA} line,
* possibly indicating accretion.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__2
parameter (nr__2=17) ! Number of records
character*78 ar__2 ! Full-size record
integer*4 Ref (nr__2) ! Reference number
character*19 BibCode (nr__2) ! BibCode
character*21 Aut (nr__2) ! Author's name
character*31 Com (nr__2) ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Measured Flux Variations
integer*4 nr__3
parameter (nr__3=126) ! Number of records
character*44 ar__3 ! Full-size record
integer*4 Seq_2 (nr__3) ! Identification number
character*12 Name_2 (nr__3) ! Name
character*5 Cat_2 (nr__3) ! PMS category
integer*4 N_1 (nr__3) ! Number of good IUE spectra used in analysis.
real*8 chi2r (nr__3) ! {chi}^2^_r_ value
character*1 l_Amp (nr__3) ! Limit flag on Amp
real*4 Amp (nr__3) ! (%) Relative amplitude of fluctuation
real*4 e_Amp (nr__3) ! (%) rms uncertainty on Amp
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Measured Fluxes for TTS
integer*4 nr__4
parameter (nr__4=107) ! Number of records
character*139 ar__4 ! Full-size record
integer*4 Seq_3 (nr__4) ! Identification number
character*1 l_ct2257 (nr__4) ! Upper limit flag on ct2257 (1)
real*8 ct2257 (nr__4) ! ([10mW/m2/nm]) Log of the continuum flux
* near 2257{AA} (2)
real*4 e_ct2257 (nr__4) ! ([10mW/m2/nm]) ? Error in ct2257 (2)
real*8 ct3060 (nr__4) ! ([10mW/m2/nm]) Log of the continuum flux
* near 3060{AA} (2)
real*4 e_ct3060 (nr__4) ! ([10mW/m2/nm]) Error in ct3060 (2)
character*1 l_CII2330 (nr__4) ! Upper limit flag on CII2330 (1)
real*4 CII2330 (nr__4) ! (10-16W/m2) Integrated C II (2330) line flux (3)
real*4 e_CII2330 (nr__4) ! (10-16W/m2) ? Error in CII2330 (3)
character*1 l_FeII2392 (nr__4) ! Upper limit flag on FeII2392 (1)
real*8 FeII2392 (nr__4) ! (10-16W/m2) Integrated Fe II (2392) line flux (3)
real*4 e_FeII2392 (nr__4) ! (10-16W/m2) ? Error in FeII2392 (3)
character*1 l_FeII2508 (nr__4) ! Upper limit flag on FeII2508 (1)
real*4 FeII2508 (nr__4) ! (10-16W/m2) Integrated Fe II (2508) line flux (3)
real*4 e_FeII2508 (nr__4) ! (10-16W/m2) ? Error in FeII2508 (3)
character*1 l_FeII2600 (nr__4) ! Upper limit flag on FeII2600 (1)
real*8 FeII2600 (nr__4) ! (10-16W/m2) Integrated Fe II (2600) line flux (3)
real*4 e_FeII2600 (nr__4) ! (10-16W/m2) ? Error in FeII2600 (3)
character*1 l_FeII2742 (nr__4) ! Upper limit flag on FeII2742 (1)
real*4 FeII2742 (nr__4) ! (10-16W/m2) Integrated Fe II (2742) line flux (3)
real*4 e_FeII2742 (nr__4) ! (10-16W/m2) ? Error in FeII2742 (3)
character*1 l_MgII2798 (nr__4) ! Upper limit flag on MgII2798 (1)
real*8 MgII2798 (nr__4) ! (10-16W/m2) Integrated Mg II (2798) line flux (3)
real*4 e_MgII2798 (nr__4) ! (10-16W/m2) ? Error in MgII2798 (3)
character*1 l_FeII2995 (nr__4) ! Upper limit flag on FeII2995 (1)
real*8 FeII2995 (nr__4) ! (10-16W/m2) Integrated Fe II (2995) line flux (3)
real*4 e_FeII2995 (nr__4) ! (10-16W/m2) ? Error in FeII2995 (3)
*Note (1): Tabulated upper limits are twice the measured uncertainty.
*Note (2): In units of erg/s/cm^2^/{AA}.
*Note (3): In units of 10^-13^erg/s/cm^2^. Line fluxes in several columns
* include significant contributions from more than one species.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Main-sequence template stars
integer*4 nr__5
parameter (nr__5=32) ! Number of records
character*28 ar__5 ! Full-size record
integer*4 HD (nr__5) ! HD number
character*5 SpType_1 (nr__5) ! MK Spectral type
integer*4 N_2 (nr__5) ! Number of IUE spectra
real*4 AV (nr__5) ! (mag) Visual extinction
* (from Neckel et al., 1980, Cat. <II/62>
real*4 v_3000_V (nr__5) ! (mag) Dereddened (3000-V) colour index
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Continuum Fluxes and Excess Line Emission
for HAEBE Stars
integer*4 nr__6
parameter (nr__6=97) ! Number of records
character*87 ar__6 ! Full-size record
integer*4 Seq_4 (nr__6) ! Identification number
character*1 l_ct2257_1 (nr__6) ! Upper limit flag on ct2257 (1)
real*8 ct2257_1 (nr__6) ! ([10mW/m2/nm]) Log of the continuum flux
* near 2257{AA} (2)
real*4 e_ct2257_1 (nr__6) ! ([10mW/m2/nm]) ? Error in ct2257 (2)
real*8 ct3060_1 (nr__6) ! ([10mW/m2/nm]) Log of the continuum flux
* near 3060{AA} (2)
real*4 e_ct3060_1 (nr__6) ! ([10mW/m2/nm]) Error in ct3060 (2)
character*1 l_FeII2392_1(nr__6) ! Upper limit flag on FeII2392 (1)
real*4 FeII2392_1 (nr__6) ! (10-16W/m2) ? Integrated Fe II (2392) line flux (3)
real*4 e_FeII2392_1(nr__6) ! (10-16W/m2) ? Error in FeII2392 (3)
character*1 l_FeII2600_1(nr__6) ! Upper limit flag on FeII2600 (1)
real*4 FeII2600_1 (nr__6) ! (10-16W/m2) ? Integrated Fe II (2600) line flux (3)
real*4 e_FeII2600_1(nr__6) ! (10-16W/m2) ? Error in FeII2600 (3)
character*1 u_FeII2600 (nr__6) ! [:] Uncertainty flag on FeII2600
character*1 l_FeII2742_1(nr__6) ! Upper limit flag on FeII2742 (1)
real*4 FeII2742_1 (nr__6) ! (10-16W/m2) ? Integrated Fe II (2742) line flux (3)
real*4 e_FeII2742_1(nr__6) ! (10-16W/m2) ? Error in FeII2742 (3)
character*1 u_FeII2742 (nr__6) ! [:] Uncertainty flag on FeII2742
character*1 l_MgII2798_1(nr__6) ! Upper limit flag on MgII2798 (1)
real*4 MgII2798_1 (nr__6) ! (10-16W/m2) ? Integrated Mg II (2798) line flux (3)
real*4 e_MgII2798_1(nr__6) ! (10-16W/m2) ? Error in MgII2798 (3)
*Note (1): Tabulated upper limits are twice the measured uncertainty.
*Note (2): In units of erg/s/cm^2^/Angstrom.
*Note (3): In units of 10^-13^ erg/s/cm^2^. Line fluxes in several columns
* include significant contributions from more than one species.
C=============================================================================
C Loading file 'table1.dat' ! Pre-main-sequence stars observed by IUE with
* LW cameras
C Format for file interpretation
1 format(
+ I3,A1,1X,I3,A1,1X,I3,1X,A15,1X,A16,1X,A13,1X,A6,A1,1X,I2,1X,A1)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,238
read(1,'(A73)')ar__
read(ar__,1)
+ Seq(i__),u_Seq(i__),VJL(i__),u_VJL(i__),HBC(i__),Name(i__),
+ OName(i__),SpType(i__),Cat(i__),n_Cat(i__),N(i__),sp(i__)
if(ar__(6:8) .EQ. '') VJL(i__) = iNULL__
if(ar__(11:13) .EQ. '') HBC(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Seq(i__),u_Seq(i__),VJL(i__),u_VJL(i__),HBC(i__),Name(i__),
+ OName(i__),SpType(i__),Cat(i__),n_Cat(i__),N(i__),sp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Additional classifications from the literature
C Format for file interpretation
2 format(I3,1X,A13,1X,A6,A1,1X,I2,1X,A47)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,22
read(1,'(A76)')ar__1
read(ar__1,2)
+ Seq_1(i__),Name_1(i__),Cat_1(i__),n_Cat_1(i__),r_Cat(i__),
+ Reasons(i__)
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),Name_1(i__),Cat_1(i__),n_Cat_1(i__),r_Cat(i__),
+ Reasons(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
3 format(I2,1X,A19,1X,A21,3X,A31)
C Effective file loading
open(unit=1,file='refs.dat', status='old')
write(6,*) '....Loading file: refs.dat'
do i__=1,17
read(1,'(A78)')ar__2
read(ar__2,3)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,3)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Measured Flux Variations
C Format for file interpretation
4 format(I3,1X,A12,1X,A5,1X,I2,1X,F7.1,1X,A1,F4.1,1X,F4.1)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,126
read(1,'(A44)')ar__3
read(ar__3,4)
+ Seq_2(i__),Name_2(i__),Cat_2(i__),N_1(i__),chi2r(i__),
+ l_Amp(i__),Amp(i__),e_Amp(i__)
c ..............Just test output...........
write(6,4)
+ Seq_2(i__),Name_2(i__),Cat_2(i__),N_1(i__),chi2r(i__),
+ l_Amp(i__),Amp(i__),e_Amp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Measured Fluxes for TTS
C Format for file interpretation
5 format(
+ I3,1X,A1,F8.4,1X,F6.4,1X,F8.4,1X,F6.4,1X,A1,F6.3,1X,F5.3,1X,
+ A1,F7.3,1X,F5.3,1X,A1,F6.3,1X,F5.3,1X,A1,F7.3,1X,F5.3,1X,A1,
+ F6.3,1X,F5.3,1X,A1,F7.3,1X,F5.3,1X,A1,F7.3,1X,F5.3)
C Effective file loading
open(unit=1,file='table6.dat', status='old')
write(6,*) '....Loading file: table6.dat'
do i__=1,107
read(1,'(A139)')ar__4
read(ar__4,5)
+ Seq_3(i__),l_ct2257(i__),ct2257(i__),e_ct2257(i__),
+ ct3060(i__),e_ct3060(i__),l_CII2330(i__),CII2330(i__),
+ e_CII2330(i__),l_FeII2392(i__),FeII2392(i__),e_FeII2392(i__),
+ l_FeII2508(i__),FeII2508(i__),e_FeII2508(i__),l_FeII2600(i__),
+ FeII2600(i__),e_FeII2600(i__),l_FeII2742(i__),FeII2742(i__),
+ e_FeII2742(i__),l_MgII2798(i__),MgII2798(i__),e_MgII2798(i__),
+ l_FeII2995(i__),FeII2995(i__),e_FeII2995(i__)
if(ar__4(15:20) .EQ. '') e_ct2257(i__) = rNULL__
if(ar__4(46:50) .EQ. '') e_CII2330(i__) = rNULL__
if(ar__4(61:65) .EQ. '') e_FeII2392(i__) = rNULL__
if(ar__4(75:79) .EQ. '') e_FeII2508(i__) = rNULL__
if(ar__4(90:94) .EQ. '') e_FeII2600(i__) = rNULL__
if(ar__4(104:108) .EQ. '') e_FeII2742(i__) = rNULL__
if(ar__4(119:123) .EQ. '') e_MgII2798(i__) = rNULL__
if(ar__4(134:138) .EQ. '') e_FeII2995(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Seq_3(i__),l_ct2257(i__),ct2257(i__),e_ct2257(i__),
+ ct3060(i__),e_ct3060(i__),l_CII2330(i__),CII2330(i__),
+ e_CII2330(i__),l_FeII2392(i__),FeII2392(i__),e_FeII2392(i__),
+ l_FeII2508(i__),FeII2508(i__),e_FeII2508(i__),l_FeII2600(i__),
+ FeII2600(i__),e_FeII2600(i__),l_FeII2742(i__),FeII2742(i__),
+ e_FeII2742(i__),l_MgII2798(i__),MgII2798(i__),e_MgII2798(i__),
+ l_FeII2995(i__),FeII2995(i__),e_FeII2995(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Main-sequence template stars
C Format for file interpretation
6 format(I6,2X,A5,I3,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='table7.dat', status='old')
write(6,*) '....Loading file: table7.dat'
do i__=1,32
read(1,'(A28)')ar__5
read(ar__5,6)
+ HD(i__),SpType_1(i__),N_2(i__),AV(i__),v_3000_V(i__)
c ..............Just test output...........
write(6,6)
+ HD(i__),SpType_1(i__),N_2(i__),AV(i__),v_3000_V(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Continuum Fluxes and Excess Line Emission
* for HAEBE Stars
C Format for file interpretation
7 format(
+ I3,1X,A1,F8.4,1X,F6.4,1X,F8.4,1X,F6.4,1X,A1,F5.2,1X,F4.2,1X,
+ A1,F5.2,1X,F4.2,A1,1X,A1,F5.2,1X,F4.2,A1,1X,A1,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,file='table9.dat', status='old')
write(6,*) '....Loading file: table9.dat'
do i__=1,97
read(1,'(A87)')ar__6
read(ar__6,7)
+ Seq_4(i__),l_ct2257_1(i__),ct2257_1(i__),e_ct2257_1(i__),
+ ct3060_1(i__),e_ct3060_1(i__),l_FeII2392_1(i__),
+ FeII2392_1(i__),e_FeII2392_1(i__),l_FeII2600_1(i__),
+ FeII2600_1(i__),e_FeII2600_1(i__),u_FeII2600(i__),
+ l_FeII2742_1(i__),FeII2742_1(i__),e_FeII2742_1(i__),
+ u_FeII2742(i__),l_MgII2798_1(i__),MgII2798_1(i__),
+ e_MgII2798_1(i__)
if(ar__6(15:20) .EQ. '') e_ct2257_1(i__) = rNULL__
if(ar__6(39:43) .EQ. '') FeII2392_1(i__) = rNULL__
if(ar__6(45:48) .EQ. '') e_FeII2392_1(i__) = rNULL__
if(ar__6(51:55) .EQ. '') FeII2600_1(i__) = rNULL__
if(ar__6(57:60) .EQ. '') e_FeII2600_1(i__) = rNULL__
if(ar__6(64:68) .EQ. '') FeII2742_1(i__) = rNULL__
if(ar__6(70:73) .EQ. '') e_FeII2742_1(i__) = rNULL__
if(ar__6(77:81) .EQ. '') MgII2798_1(i__) = rNULL__
if(ar__6(83:86) .EQ. '') e_MgII2798_1(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ Seq_4(i__),l_ct2257_1(i__),ct2257_1(i__),e_ct2257_1(i__),
+ ct3060_1(i__),e_ct3060_1(i__),l_FeII2392_1(i__),
+ FeII2392_1(i__),e_FeII2392_1(i__),l_FeII2600_1(i__),
+ FeII2600_1(i__),e_FeII2600_1(i__),u_FeII2600(i__),
+ l_FeII2742_1(i__),FeII2742_1(i__),e_FeII2742_1(i__),
+ u_FeII2742(i__),l_MgII2798_1(i__),MgII2798_1(i__),
+ e_MgII2798_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end