FORTRAN Generation
(/./ftp/cats/J/AJ/129/363)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/129/363 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-Mar-19
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/129/363    Spectroscopy in Orion Nebula Cluster    (Sicilia-Aguilar+, 2005)
*================================================================================
*Accretion, kinematics, and rotation in the Orion Nebula Cluster: initial results
*from Hectochelle.
*    Sicilia-Aguilar A., Hartmann L.W., Szentgyorgyi A.H., Fabricant D.G.,
*    Furesz G., Roll J., Conroy M.A., Calvet N., Tokarz S., Hernandez J.
*   <Astron. J., 129, 363-381 (2005)>
*   =2005AJ....129..363S
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! Results from 2003 December

      integer*4 nr__
      parameter (nr__=185)	! Number of records
      character*75 ar__   	! Full-size record

      integer*4     v_H97b_    (nr__) ! ? Identification number (G1)
      character*2   m__H97b_   (nr__) ! Multiple components
      character*16  v2MASS     (nr__) ! The 2MASS identification (HHMMSSss+DDMMSSs)
      character*14  SpType     (nr__) ! Hillenbrand et al. (1998, Cat. <J/AJ/116/1816>)
*                                  spectral type
      character*2   Type       (nr__) ! [CW: ] Type T Tauri star classification (G2)
      real*4        EW_Li      (nr__) ! (0.1nm) ? Equivalent width of Li (6707 Angstroms) in
*                                  Angstroms
      character*1   u_EW_Li    (nr__) ! Uncertainty flag in EW(Li)
      integer*4     EW_Ha      (nr__) ! (0.1nm) ? Equivalent width of H{alpha} in Angstroms
      character*22  Com        (nr__) ! Additional comments (1)
*Note (1): Proplyds named using their HST ID number as given in
*    Hillenbrand et al. 1998, Cat. <J/AJ/116/1816>) , Bally et al.
*    (1998AJ....116..293B). Winds (denoted PC, or P-Cygni profiles) 
*    and infall (denoted IPC, or inverse P-Cygni profiles) detected by
*    asymmetries in H{alpha} profile. Broad or very broad He (written as
*    bHe or very bHe) refers to the detection of the broad component of the
*    He emission line at 5876{AA}. Low signal to noise spectra are marked
*    with the label lS/N. Colons have been added to uncertain values. Stars
*    with near-IR excess in 2MASS are labeled as nIR. An additional comment
*    has been included to distinguish stars with systematically broad lines
*    as possible fast rotators. The label "(2004)" indicate stars that were
*    observed again in March 2004, which have the additional information
*    corresponding to the 2004 run in Table2.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! Results from March 2004

      integer*4 nr__1
      parameter (nr__1=153)	! Number of records
      character*92 ar__1  	! Full-size record

      integer*4     v_H97b__1  (nr__1) ! Identification number (G1)
      character*2   m__H97b__1 (nr__1) ! Multiple components
      character*16  v2MASS_1   (nr__1) ! The 2MASS identification (HHMMSSss+DDMMSSs)
      real*4        cz         (nr__1) ! (km/s) Heliocentric radial velocity
      real*4        e_cz       (nr__1) ! (km/s) Uncertainty in cz
      real*4        S_N        (nr__1) ! ? Signal-to-noise value of cross-correlation
*                                     peak (see text)
      real*4        vsini      (nr__1) ! (km/s) ? Rotational velocity
      real*4        e_vsini    (nr__1) ! (km/s) ? Uncertainty in vsini
      character*5   Type_1     (nr__1) ! T Tauri star classification (G2)
      real*4        EW_Ha_1    (nr__1) ! (0.1nm) ? Equivalent width of H{alpha} in Angstroms
      character*1   u_EW_Ha    (nr__1) ! Uncertainty flag in EW(Ha)
      character*12  SpType_1   (nr__1) ! Spectral type (1)
      character*14  Com_1      (nr__1) ! Additional comments (2)
*Note (1): From Hillenbrand (1997, Cat. <J/AJ/113/1733>) and subsequent updates
*          to electronically-available table except where noted:
*    P   = Prosser & Stauffer, unpublished;
*    Sam = A.E. Samuel, 1993 unpublished PhD thesis;
*    D   = Duncan (1993ApJ...406..172D); Strand (1958AJ.....63..152S);
*    Her = Herbig, quoted in Walker (1969ApJ...155..447W);
*    LR  = Luhman et al. 2000, Cat. <J/ApJ/540/1016>;
*    e   = Ca II infrared triplet emission;
*    :   = uncertain values.
*Note (2): Additional comments, defined as follows:
*    PC   = P Cygni profile;
*    IPC  = inverse P Cygni profile (see text);
*    SB2s = new double-lined binary;
*    nIR  = star with near-IR excess in 2MASS;
*    NM   = non-member from proper motions.

C=============================================================================

C  Loading file 'table1.dat'	! Results from 2003 December

C  Format for file interpretation

    1 format(I4,A2,1X,A16,1X,A14,1X,A2,1X,F4.2,A1,1X,I4,1X,A22)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,185
        read(1,'(A75)')ar__
        read(ar__,1)
     +  v_H97b_(i__),m__H97b_(i__),v2MASS(i__),SpType(i__),Type(i__),
     +  EW_Li(i__),u_EW_Li(i__),EW_Ha(i__),Com(i__)
        if(ar__(1:4) .EQ. '') v_H97b_(i__) = iNULL__
        if(ar__(43:46) .EQ. '') EW_Li(i__) = rNULL__
        if(ar__(49:52) .EQ. '') EW_Ha(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  v_H97b_(i__),m__H97b_(i__),v2MASS(i__),SpType(i__),Type(i__),
     +  EW_Li(i__),u_EW_Li(i__),EW_Ha(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2.dat'	! Results from March 2004

C  Format for file interpretation

    2 format(
     +  I4,A2,1X,A16,1X,F5.1,1X,F4.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A5,1X,
     +  F6.1,A1,1X,A12,1X,A14)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,153
        read(1,'(A92)')ar__1
        read(ar__1,2)
     +  v_H97b__1(i__),m__H97b__1(i__),v2MASS_1(i__),cz(i__),
     +  e_cz(i__),S_N(i__),vsini(i__),e_vsini(i__),Type_1(i__),
     +  EW_Ha_1(i__),u_EW_Ha(i__),SpType_1(i__),Com_1(i__)
        if(ar__1(36:39) .EQ. '') S_N(i__) = rNULL__
        if(ar__1(41:45) .EQ. '') vsini(i__) = rNULL__
        if(ar__1(47:50) .EQ. '') e_vsini(i__) = rNULL__
        if(ar__1(58:63) .EQ. '') EW_Ha_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  v_H97b__1(i__),m__H97b__1(i__),v2MASS_1(i__),cz(i__),
     +  e_cz(i__),S_N(i__),vsini(i__),e_vsini(i__),Type_1(i__),
     +  EW_Ha_1(i__),u_EW_Ha(i__),SpType_1(i__),Com_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end