Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/620/A87 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-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 (nr__) ! Object identifier according to CSIMon
* catalog (G1)
character*17 v2MASS (nr__) ! Object identifier according to 2MASS catalog
character*5 SpType (nr__) ! Object's spectral type (1)
real*4 Lbol (nr__) ! (Lsun) Stellar bolometric luminosity
real*4 Mass (nr__) ! (Msun) Stellar mass
real*4 Rad (nr__) ! (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 (nr__1) ! Object identifier according to CSIMon
* catalog (G1)
character*1 m_CSIMon (nr__1) ! [ab] Multiplicity index on CSIMon (1)
character*17 v2MASS_1 (nr__1) ! Object identifier according to 2MASS catalog
character*7 SpTlit (nr__1) ! Object's spectral type,
* from the literature (2)
character*5 SpTw (nr__1) ! Spectral type of the WTTS that was used as
* the photospheric template for this object (2)
real*4 Lbol_1 (nr__1) ! (Lsun) ? Stellar bolometric luminosity
real*4 Mass_1 (nr__1) ! (Msun) ? Stellar mass
real*4 Rad_1 (nr__1) ! (Rsun) ? Stellar radius
real*4 Macc (nr__1) ! (Msun/yr) ? Mass accretion rate (3)
real*4 Lacc (nr__1) ! (Lsun) ? Accretion luminosity (3)
real*4 Veil (nr__1) ! ? Veiling (ratio of excess flux to stellar
* flux) measured in our spectra
real*4 e_Veil (nr__1) ! ? Error in veiling
real*4 aIRAC (nr__1) ! ? alpha_IRAC_ index: measured as the slope of
* the star's SED between 3.6micron and 8micron
real*4 E_K_3_6 (nr__1) ! (mag) ? Infrared excess measured from the
* (K - 3.6micron) color index
real*4 Per (nr__1) ! (d) ? Stellar rotation period measured from the
* 2011 or 2008 CoRoT light curve
real*4 RV (nr__1) ! (km/s) ? Stellar radial velocity measured in our
* spectra
real*4 e_RV (nr__1) ! (km/s) ? Error in stellar radial velocity
real*4 vsini (nr__1) ! (km/s) ? vsini (stellar rotational velocity
* projected in our line-of-sight) measured in
* our spectra
real*4 e_vsini (nr__1) ! (km/s) ? Error in vsini
integer*4 inc (nr__1) ! (deg) ? System inclination derived from rotation
* properties
integer*4 E_inc (nr__1) ! (deg) ? Upper error bar in system inclination
integer*4 e_inc_1 (nr__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 (nr__2) ! Object identifier according to CSIMon
* catalog (G1)
real*4 EW_NLVC (nr__2) ! (0.1nm) ? Equivalent width of the narrow low-velocity
* component (NLVC)
real*4 vc_NLVC (nr__2) ! (km/s) ? Centroid velocity of the NLVC
real*4 e_vc_NLVC (nr__2) ! (km/s) ? Error on the centroid velocity of the NLVC
real*4 FW_NLVC (nr__2) ! (km/s) ? FWHM of the NLVC
real*4 e_FW_NLVC (nr__2) ! (km/s) ? Error on the FWHM of the NLVC
real*4 EW_BLVC (nr__2) ! (0.1nm) ? Equivalent width of the broad low-velocity
* component (BLVC)
real*4 vc_BLVC (nr__2) ! (km/s) ? Centroid velocity of the BLVC
real*4 e_vc_BLVC (nr__2) ! (km/s) ? Error on the centroid velocity of the BLVC
real*4 FW_BLVC (nr__2) ! (km/s) ? FWHM of the BLVC
real*4 e_FW_BLVC (nr__2) ! (km/s) ? Error on the FWHM of the BLVC
real*4 EW_BHVC (nr__2) ! (0.1nm) ? Equivalent width of the blueshifted
* high-velocity component (BHVC)
real*4 vc_BHVC (nr__2) ! (km/s) ? Centroid velocity of the BHVC
real*4 e_vc_BHVC (nr__2) ! (km/s) ? Error on the centroid velocity of the BHVC
real*4 FW_BHVC (nr__2) ! (km/s) ? FWHM of the BHVC
real*4 e_FW_BHVC (nr__2) ! (km/s) ? Error on the FWHM of the BHVC
real*4 EW_RHVC (nr__2) ! (0.1nm) ? Equivalent width of the redshifted
* high-velocity component (RHVC)
real*4 vc_RHVC (nr__2) ! (km/s) ? Centroid velocity of the RHVC
real*4 e_vc_RHVC (nr__2) ! (km/s) ? Error on the centroid velocity of the RHVC
real*4 FW_RHVC (nr__2) ! (km/s) ? FWHM of the RHVC
real*4 e_FW_RHVC (nr__2) ! (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 (nr__3) ! Object identifier according to CSIMon
* catalog (G1)
character*17 v2MASS_2 (nr__3) ! Object identifier according to 2MASS catalog
real*4 L_tot (nr__3) ! (Lsun) Luminosity of the total [OI]6300{AA} line
* profile
character*1 Upplim (nr__3) ! [*] indicates an upper limit at the
* detection threshold
real*4 L_NLVC (nr__3) ! (Lsun) ? Luminosity of the narrow low-velocity
* component of the [OI]6300{AA} line
real*4 L_BLVC (nr__3) ! (Lsun) ? Luminosity of the broad low-velocity
* component of the [OI]6300{AA} line
real*4 L_BHVC (nr__3) ! (Lsun) ? Luminosity of the blueshifted high-velocity
* component of the [OI]6300{AA} line
real*4 L_RHVC (nr__3) ! (Lsun) ? Luminosity of the redshifted high-velocity
* component of the [OI]6300{AA} line
integer*4 vjet_B (nr__3) ! (km/s) ? Jet velocity of the blue lobe, derived from
* the deprojected centroid velocity of the BHVC
integer*4 vjet_R (nr__3) ! (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(i__),v2MASS(i__),SpType(i__),Lbol(i__),Mass(i__),
+ Rad(i__)
c ..............Just test output...........
write(6,1)
+ CSIMon(i__),v2MASS(i__),SpType(i__),Lbol(i__),Mass(i__),
+ Rad(i__)
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(i__),m_CSIMon(i__),v2MASS_1(i__),SpTlit(i__),
+ SpTw(i__),Lbol_1(i__),Mass_1(i__),Rad_1(i__),Macc(i__),
+ Lacc(i__),Veil(i__),e_Veil(i__),aIRAC(i__),E_K_3_6(i__),
+ Per(i__),RV(i__),e_RV(i__),vsini(i__),e_vsini(i__),inc(i__),
+ E_inc(i__),e_inc_1(i__)
if(ar__1(48:52) .EQ. '') Lbol_1(i__) = rNULL__
if(ar__1(54:57) .EQ. '') Mass_1(i__) = rNULL__
if(ar__1(59:62) .EQ. '') Rad_1(i__) = rNULL__
if(ar__1(64:71) .EQ. '') Macc(i__) = rNULL__
if(ar__1(73:80) .EQ. '') Lacc(i__) = rNULL__
if(ar__1(82:85) .EQ. '') Veil(i__) = rNULL__
if(ar__1(87:90) .EQ. '') e_Veil(i__) = rNULL__
if(ar__1(92:96) .EQ. '') aIRAC(i__) = rNULL__
if(ar__1(98:102) .EQ. '') E_K_3_6(i__) = rNULL__
if(ar__1(104:108) .EQ. '') Per(i__) = rNULL__
if(ar__1(110:114) .EQ. '') RV(i__) = rNULL__
if(ar__1(116:119) .EQ. '') e_RV(i__) = rNULL__
if(ar__1(121:125) .EQ. '') vsini(i__) = rNULL__
if(ar__1(127:130) .EQ. '') e_vsini(i__) = rNULL__
if(ar__1(132:133) .EQ. '') inc(i__) = iNULL__
if(ar__1(135:136) .EQ. '') E_inc(i__) = iNULL__
if(ar__1(138:139) .EQ. '') e_inc_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ CSIMon_1(i__),m_CSIMon(i__),v2MASS_1(i__),SpTlit(i__),
+ SpTw(i__),Lbol_1(i__),Mass_1(i__),Rad_1(i__),Macc(i__),
+ Lacc(i__),Veil(i__),e_Veil(i__),aIRAC(i__),E_K_3_6(i__),
+ Per(i__),RV(i__),e_RV(i__),vsini(i__),e_vsini(i__),inc(i__),
+ E_inc(i__),e_inc_1(i__)
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(i__),EW_NLVC(i__),vc_NLVC(i__),e_vc_NLVC(i__),
+ FW_NLVC(i__),e_FW_NLVC(i__),EW_BLVC(i__),vc_BLVC(i__),
+ e_vc_BLVC(i__),FW_BLVC(i__),e_FW_BLVC(i__),EW_BHVC(i__),
+ vc_BHVC(i__),e_vc_BHVC(i__),FW_BHVC(i__),e_FW_BHVC(i__),
+ EW_RHVC(i__),vc_RHVC(i__),e_vc_RHVC(i__),FW_RHVC(i__),
+ e_FW_RHVC(i__)
if(ar__2(15:19) .EQ. '') EW_NLVC(i__) = rNULL__
if(ar__2(21:25) .EQ. '') vc_NLVC(i__) = rNULL__
if(ar__2(27:30) .EQ. '') e_vc_NLVC(i__) = rNULL__
if(ar__2(32:35) .EQ. '') FW_NLVC(i__) = rNULL__
if(ar__2(37:40) .EQ. '') e_FW_NLVC(i__) = rNULL__
if(ar__2(42:46) .EQ. '') EW_BLVC(i__) = rNULL__
if(ar__2(48:52) .EQ. '') vc_BLVC(i__) = rNULL__
if(ar__2(54:57) .EQ. '') e_vc_BLVC(i__) = rNULL__
if(ar__2(59:63) .EQ. '') FW_BLVC(i__) = rNULL__
if(ar__2(65:68) .EQ. '') e_FW_BLVC(i__) = rNULL__
if(ar__2(70:74) .EQ. '') EW_BHVC(i__) = rNULL__
if(ar__2(76:81) .EQ. '') vc_BHVC(i__) = rNULL__
if(ar__2(83:86) .EQ. '') e_vc_BHVC(i__) = rNULL__
if(ar__2(88:92) .EQ. '') FW_BHVC(i__) = rNULL__
if(ar__2(94:97) .EQ. '') e_FW_BHVC(i__) = rNULL__
if(ar__2(99:103) .EQ. '') EW_RHVC(i__) = rNULL__
if(ar__2(105:109) .EQ. '') vc_RHVC(i__) = rNULL__
if(ar__2(111:114) .EQ. '') e_vc_RHVC(i__) = rNULL__
if(ar__2(116:120) .EQ. '') FW_RHVC(i__) = rNULL__
if(ar__2(122:125) .EQ. '') e_FW_RHVC(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ CSIMon_2(i__),EW_NLVC(i__),vc_NLVC(i__),e_vc_NLVC(i__),
+ FW_NLVC(i__),e_FW_NLVC(i__),EW_BLVC(i__),vc_BLVC(i__),
+ e_vc_BLVC(i__),FW_BLVC(i__),e_FW_BLVC(i__),EW_BHVC(i__),
+ vc_BHVC(i__),e_vc_BHVC(i__),FW_BHVC(i__),e_FW_BHVC(i__),
+ EW_RHVC(i__),vc_RHVC(i__),e_vc_RHVC(i__),FW_RHVC(i__),
+ e_FW_RHVC(i__)
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(i__),v2MASS_2(i__),L_tot(i__),Upplim(i__),
+ L_NLVC(i__),L_BLVC(i__),L_BHVC(i__),L_RHVC(i__),vjet_B(i__),
+ vjet_R(i__)
if(ar__3(44:51) .EQ. '') L_NLVC(i__) = rNULL__
if(ar__3(53:60) .EQ. '') L_BLVC(i__) = rNULL__
if(ar__3(62:69) .EQ. '') L_BHVC(i__) = rNULL__
if(ar__3(71:78) .EQ. '') L_RHVC(i__) = rNULL__
if(ar__3(80:83) .EQ. '') vjet_B(i__) = iNULL__
if(ar__3(85:88) .EQ. '') vjet_R(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ CSIMon_3(i__),v2MASS_2(i__),L_tot(i__),Upplim(i__),
+ L_NLVC(i__),L_BLVC(i__),L_BHVC(i__),L_RHVC(i__),vjet_B(i__),
+ vjet_R(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end