Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/620/A87 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 2026-Oct-12
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/A+A/620/A87 Ejection processes in NGC 2264 (McGinnis+, 2018)
*================================================================================
*Ejection processes in the young open cluster NGC 2264.
*A study of the [OI]{lambda}6300{AA} emission line.
* McGinnis P., Dougados C., Alencar S.H.P., Bouvier J., Cabrit S.
* <Astron. Astrophys. 620, A87 (2018)>
* =2018A&A...620A..87M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'wtts.dat' ! Stellar parameters of weak-lined T Tauri stars
in NGC 2264
integer*4 nr__
parameter (nr__=34) ! Number of records
character*52 ar__ ! Full-size record
character*13 CSIMon ! Object identifier according to CSIMon
* catalog (G1)
character*17 v2MASS ! Object identifier according to 2MASS catalog
character*5 SpType ! Object's spectral type (1)
real*4 Lbol ! (Lsun) Stellar bolometric luminosity
real*4 Mass ! (Msun) Stellar mass
real*4 Rad ! (Rsun) Stellar radius
*Note (1): Spectral types were taken from Dahm & Simon (2005AJ....129..829D,
* Cat. J/AJ/129/829), Rebull et al. (2002AJ....123.1528R, Cat. J/AJ/123/1528),
* Walker (1956ApJS....2..365W) and Venuti et al. (2014A&A...570A..82V,
* Cat. J/A+A/570/A82), in that order of preference.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ctts.dat' ! Stellar parameters of classical T Tauri stars
in NGC 2264
integer*4 nr__1
parameter (nr__1=186) ! Number of records
character*139 ar__1 ! Full-size record
character*13 CSIMon_1 ! Object identifier according to CSIMon
* catalog (G1)
character*1 m_CSIMon ! [ab] Multiplicity index on CSIMon (1)
character*17 v2MASS_1 ! Object identifier according to 2MASS catalog
character*7 SpTlit ! Object's spectral type,
* from the literature (2)
character*5 SpTw ! Spectral type of the WTTS that was used as
* the photospheric template for this object (2)
real*4 Lbol_1 ! (Lsun) ? Stellar bolometric luminosity
real*4 Mass_1 ! (Msun) ? Stellar mass
real*4 Rad_1 ! (Rsun) ? Stellar radius
real*4 Macc ! (Msun/yr) ? Mass accretion rate (3)
real*4 Lacc ! (Lsun) ? Accretion luminosity (3)
real*4 Veil ! ? Veiling (ratio of excess flux to stellar
* flux) measured in our spectra
real*4 e_Veil ! ? Error in veiling
real*4 aIRAC ! ? alpha_IRAC_ index: measured as the slope of
* the star's SED between 3.6micron and 8micron
real*4 E_K_3_6 ! (mag) ? Infrared excess measured from the
* (K - 3.6micron) color index
real*4 Per ! (d) ? Stellar rotation period measured from the
* 2011 or 2008 CoRoT light curve
real*4 RV ! (km/s) ? Stellar radial velocity measured in our
* spectra
real*4 e_RV ! (km/s) ? Error in stellar radial velocity
real*4 vsini ! (km/s) ? vsini (stellar rotational velocity
* projected in our line-of-sight) measured in
* our spectra
real*4 e_vsini ! (km/s) ? Error in vsini
integer*4 inc ! (deg) ? System inclination derived from rotation
* properties
integer*4 E_inc ! (deg) ? Upper error bar in system inclination
integer*4 e_inc_1 ! (deg) ? Lower error bar in system inclination
*Note (1): The characters 'a' and 'b' following the identifier indicate that
* this object is a binary.
*Note (2): Spectral types were taken from Dahm & Simon (2005AJ....129..829D,
* Cat. J/AJ/129/829), Rebull et al. (2002AJ....123.1528R, Cat. J/AJ/123/1528),
* Walker (1956ApJS....2..365W) and Venuti et al. (2014A&A...570A..82V,
* Cat. J/A+A/570/A82), in that order of preference.
*Note (3): Mass accretion rates and accretion luminosities were taken
* preferentially from Venuti et al. (2014A&A...570A..82V, Cat. J/AJ/123/1528),
* who derived them from (u-r) color excesses, or else from the luminosity of
* the Halpha line.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'oi_comp.dat' ! Equivalent widths, centroid velocities and FWHM
of the components of the [OI]6300{AA} line
integer*4 nr__2
parameter (nr__2=110) ! Number of records
character*125 ar__2 ! Full-size record
character*13 CSIMon_2 ! Object identifier according to CSIMon
* catalog (G1)
real*4 EW_NLVC ! (0.1nm) ? Equivalent width of the narrow low-velocity
* component (NLVC)
real*4 vc_NLVC ! (km/s) ? Centroid velocity of the NLVC
real*4 e_vc_NLVC ! (km/s) ? Error on the centroid velocity of the NLVC
real*4 FW_NLVC ! (km/s) ? FWHM of the NLVC
real*4 e_FW_NLVC ! (km/s) ? Error on the FWHM of the NLVC
real*4 EW_BLVC ! (0.1nm) ? Equivalent width of the broad low-velocity
* component (BLVC)
real*4 vc_BLVC ! (km/s) ? Centroid velocity of the BLVC
real*4 e_vc_BLVC ! (km/s) ? Error on the centroid velocity of the BLVC
real*4 FW_BLVC ! (km/s) ? FWHM of the BLVC
real*4 e_FW_BLVC ! (km/s) ? Error on the FWHM of the BLVC
real*4 EW_BHVC ! (0.1nm) ? Equivalent width of the blueshifted
* high-velocity component (BHVC)
real*4 vc_BHVC ! (km/s) ? Centroid velocity of the BHVC
real*4 e_vc_BHVC ! (km/s) ? Error on the centroid velocity of the BHVC
real*4 FW_BHVC ! (km/s) ? FWHM of the BHVC
real*4 e_FW_BHVC ! (km/s) ? Error on the FWHM of the BHVC
real*4 EW_RHVC ! (0.1nm) ? Equivalent width of the redshifted
* high-velocity component (RHVC)
real*4 vc_RHVC ! (km/s) ? Centroid velocity of the RHVC
real*4 e_vc_RHVC ! (km/s) ? Error on the centroid velocity of the RHVC
real*4 FW_RHVC ! (km/s) ? FWHM of the RHVC
real*4 e_FW_RHVC ! (km/s) ? Error on the FWHM of the RHVC
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'oi_lum.dat' ! Luminosities of the components of the
[OI]6300{AA} line and jet velocities
integer*4 nr__3
parameter (nr__3=184) ! Number of records
character*88 ar__3 ! Full-size record
character*13 CSIMon_3 ! Object identifier according to CSIMon
* catalog (G1)
character*17 v2MASS_2 ! Object identifier according to 2MASS catalog
real*4 L_tot ! (Lsun) Luminosity of the total [OI]6300{AA} line
* profile
character*1 Upplim ! [*] indicates an upper limit at the
* detection threshold
real*4 L_NLVC ! (Lsun) ? Luminosity of the narrow low-velocity
* component of the [OI]6300{AA} line
real*4 L_BLVC ! (Lsun) ? Luminosity of the broad low-velocity
* component of the [OI]6300{AA} line
real*4 L_BHVC ! (Lsun) ? Luminosity of the blueshifted high-velocity
* component of the [OI]6300{AA} line
real*4 L_RHVC ! (Lsun) ? Luminosity of the redshifted high-velocity
* component of the [OI]6300{AA} line
integer*4 vjet_B ! (km/s) ? Jet velocity of the blue lobe, derived from
* the deprojected centroid velocity of the BHVC
integer*4 vjet_R ! (km/s) ? Jet velocity of the red lobe, derived from
* the deprojected centroid velocity of the RHVC
C=============================================================================
C Loading file 'wtts.dat' ! Stellar parameters of weak-lined T Tauri stars
* in NGC 2264
C Format for file interpretation
1 format(A13,1X,A17,1X,A5,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'wtts.dat')
write(6,*) '....Loading file: wtts.dat'
do i__=1,34
read(1,'(A52)')ar__
read(ar__,1)CSIMon,v2MASS,SpType,Lbol,Mass,Rad
c ..............Just test output...........
write(6,1)CSIMon,v2MASS,SpType,Lbol,Mass,Rad
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ctts.dat' ! Stellar parameters of classical T Tauri stars
* in NGC 2264
C Format for file interpretation
2 format(
+ A13,A1,1X,A17,1X,A7,1X,A5,1X,F5.2,1X,F4.2,1X,F4.2,1X,E8.2,1X,
+ E8.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.1,1X,F4.1,
+ 1X,F5.1,1X,F4.1,1X,I2,1X,I2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'ctts.dat')
write(6,*) '....Loading file: ctts.dat'
do i__=1,186
read(1,'(A139)')ar__1
read(ar__1,2)
+ CSIMon_1,m_CSIMon,v2MASS_1,SpTlit,SpTw,Lbol_1,Mass_1,Rad_1,
+ Macc,Lacc,Veil,e_Veil,aIRAC,E_K_3_6,Per,RV,e_RV,vsini,e_vsini,
+ inc,E_inc,e_inc_1
if(ar__1(48:52) .EQ. '') Lbol_1 = rNULL__
if(ar__1(54:57) .EQ. '') Mass_1 = rNULL__
if(ar__1(59:62) .EQ. '') Rad_1 = rNULL__
if(ar__1(64:71) .EQ. '') Macc = rNULL__
if(ar__1(73:80) .EQ. '') Lacc = rNULL__
if(ar__1(82:85) .EQ. '') Veil = rNULL__
if(ar__1(87:90) .EQ. '') e_Veil = rNULL__
if(ar__1(92:96) .EQ. '') aIRAC = rNULL__
if(ar__1(98:102) .EQ. '') E_K_3_6 = rNULL__
if(ar__1(104:108) .EQ. '') Per = rNULL__
if(ar__1(110:114) .EQ. '') RV = rNULL__
if(ar__1(116:119) .EQ. '') e_RV = rNULL__
if(ar__1(121:125) .EQ. '') vsini = rNULL__
if(ar__1(127:130) .EQ. '') e_vsini = rNULL__
if(ar__1(132:133) .EQ. '') inc = iNULL__
if(ar__1(135:136) .EQ. '') E_inc = iNULL__
if(ar__1(138:139) .EQ. '') e_inc_1 = iNULL__
c ..............Just test output...........
write(6,2)
+ CSIMon_1,m_CSIMon,v2MASS_1,SpTlit,SpTw,Lbol_1,Mass_1,Rad_1,
+ Macc,Lacc,Veil,e_Veil,aIRAC,E_K_3_6,Per,RV,e_RV,vsini,e_vsini,
+ inc,E_inc,e_inc_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'oi_comp.dat' ! Equivalent widths, centroid velocities and FWHM
* of the components of the [OI]6300{AA} line
C Format for file interpretation
3 format(
+ A13,1X,F5.3,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F5.3,1X,F5.1,
+ 1X,F4.1,1X,F5.1,1X,F4.1,1X,F5.3,1X,F6.1,1X,F4.1,1X,F5.1,1X,
+ F4.1,1X,F5.3,1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'oi_comp.dat')
write(6,*) '....Loading file: oi_comp.dat'
do i__=1,110
read(1,'(A125)')ar__2
read(ar__2,3)
+ CSIMon_2,EW_NLVC,vc_NLVC,e_vc_NLVC,FW_NLVC,e_FW_NLVC,EW_BLVC,
+ vc_BLVC,e_vc_BLVC,FW_BLVC,e_FW_BLVC,EW_BHVC,vc_BHVC,e_vc_BHVC,
+ FW_BHVC,e_FW_BHVC,EW_RHVC,vc_RHVC,e_vc_RHVC,FW_RHVC,e_FW_RHVC
if(ar__2(15:19) .EQ. '') EW_NLVC = rNULL__
if(ar__2(21:25) .EQ. '') vc_NLVC = rNULL__
if(ar__2(27:30) .EQ. '') e_vc_NLVC = rNULL__
if(ar__2(32:35) .EQ. '') FW_NLVC = rNULL__
if(ar__2(37:40) .EQ. '') e_FW_NLVC = rNULL__
if(ar__2(42:46) .EQ. '') EW_BLVC = rNULL__
if(ar__2(48:52) .EQ. '') vc_BLVC = rNULL__
if(ar__2(54:57) .EQ. '') e_vc_BLVC = rNULL__
if(ar__2(59:63) .EQ. '') FW_BLVC = rNULL__
if(ar__2(65:68) .EQ. '') e_FW_BLVC = rNULL__
if(ar__2(70:74) .EQ. '') EW_BHVC = rNULL__
if(ar__2(76:81) .EQ. '') vc_BHVC = rNULL__
if(ar__2(83:86) .EQ. '') e_vc_BHVC = rNULL__
if(ar__2(88:92) .EQ. '') FW_BHVC = rNULL__
if(ar__2(94:97) .EQ. '') e_FW_BHVC = rNULL__
if(ar__2(99:103) .EQ. '') EW_RHVC = rNULL__
if(ar__2(105:109) .EQ. '') vc_RHVC = rNULL__
if(ar__2(111:114) .EQ. '') e_vc_RHVC = rNULL__
if(ar__2(116:120) .EQ. '') FW_RHVC = rNULL__
if(ar__2(122:125) .EQ. '') e_FW_RHVC = rNULL__
c ..............Just test output...........
write(6,3)
+ CSIMon_2,EW_NLVC,vc_NLVC,e_vc_NLVC,FW_NLVC,e_FW_NLVC,EW_BLVC,
+ vc_BLVC,e_vc_BLVC,FW_BLVC,e_FW_BLVC,EW_BHVC,vc_BHVC,e_vc_BHVC,
+ FW_BHVC,e_FW_BHVC,EW_RHVC,vc_RHVC,e_vc_RHVC,FW_RHVC,e_FW_RHVC
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'oi_lum.dat' ! Luminosities of the components of the
* [OI]6300{AA} line and jet velocities
C Format for file interpretation
4 format(
+ A13,1X,A17,1X,E8.2,1X,A1,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,
+ I4,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'oi_lum.dat')
write(6,*) '....Loading file: oi_lum.dat'
do i__=1,184
read(1,'(A88)')ar__3
read(ar__3,4)
+ CSIMon_3,v2MASS_2,L_tot,Upplim,L_NLVC,L_BLVC,L_BHVC,L_RHVC,
+ vjet_B,vjet_R
if(ar__3(44:51) .EQ. '') L_NLVC = rNULL__
if(ar__3(53:60) .EQ. '') L_BLVC = rNULL__
if(ar__3(62:69) .EQ. '') L_BHVC = rNULL__
if(ar__3(71:78) .EQ. '') L_RHVC = rNULL__
if(ar__3(80:83) .EQ. '') vjet_B = iNULL__
if(ar__3(85:88) .EQ. '') vjet_R = iNULL__
c ..............Just test output...........
write(6,4)
+ CSIMon_3,v2MASS_2,L_tot,Upplim,L_NLVC,L_BLVC,L_BHVC,L_RHVC,
+ vjet_B,vjet_R
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end