FORTRAN Generation
(/./ftp/cats/J/ApJ/831/155)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/831/155 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-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/831/155 VLA survey for faint compact radio sources in ONC (Sheehan+, 2016)
*================================================================================
*A VLA survey for faint compact radio sources in the Orion Nebula Cluster.
*    Sheehan P.D., Eisner J.A., Mann R.K., Williams J.P.
*   <Astrophys. J., 831, 155-155 (2016)>
*   =2016ApJ...831..155S    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Source detections, identifications, and fluxes

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      integer*4     Seq         ! Running sequence number
      character*8   v_RRS2008_  ! Proplyd name (from Ricci+, 2008, J/AJ/136/2136)
      integer*4     v_HC2000_   ! ? Hillenbrand & Carpenter 2000, J/ApJ/540/236
*                                  identifier
      character*2   GMR         ! Garay+ 1987ApJ...314..535G identifier
      integer*4     v_ZRK2004_  ! ? Zapata+ 2004AJ....127.2252Z identifier
      character*7   OName       ! Other Names
      integer*4     RAh         ! (h) [5] Hour of Right Ascension (J2000)
      integer*4     RAm         ! (min) [34/35] Minute of Right Ascension (J2000)
      real*4        RAs         ! (s) Second of Right Ascension (J2000)
      character*1   DE_         ! [-] Sign of the Declination (J2000)
      integer*4     DEd         ! (deg) [5] Degree of Declination (J2000)
      integer*4     DEm         ! (arcmin) [18/28] Arcminute of Declination (J2000)
      real*4        DEs         ! (arcsec) Arcsecond of Declination (J2000)
      real*4        F6          ! (mJy) [-0.4/36] The 6cm flux density
      real*4        e_F6        ! (mJy) [0.04/4] Uncertainty in F6
      real*4        F3_6        ! (mJy) [-1/37]? The 3.6cm flux density
      real*4        e_F3_6      ! (mJy) [0.03/4]? Uncertainty in F3.6
      real*4        F1_3a       ! (mJy) [-1.5/29]? Nov. 10, 2013 1.3cm flux density
      real*4        e_F1_3a     ! (mJy) [0.03/38]? Uncertainty in F1.3,1
      real*4        F1_3b       ! (mJy) [-1/41]? Mar. 3, 2014 1.3cm flux density
      real*4        e_F1_3b     ! (mJy) [0.03/25]? Uncertainty in F1.3,2
      real*4        F1_3c       ! (mJy) [-1.2/30]? Mar. 7, 2014 1.3cm flux density
      real*4        e_F1_3c     ! (mJy) [0.03/70]? Uncertainty in F1.3,3
      real*4        F1_3d       ! (mJy) [-1/23]? May 3, 2014 1.3cm flux density
      real*4        e_F1_3d     ! (mJy) [0.03/3]? Uncertainty in F1.3,4
      real*4        v_F1_3_     ! (mJy) [-0.9/19]? Mean 1.3cm flux density
      real*8        e__F1_3_    ! (mJy) [0.01/10335]? Uncertainty in <F1.3>

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

C  Declarations for 'table3.dat'	! Free-free emission model parameters and ALMA
                                 band fluxes for detected sources

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

      integer*4     Seq_1       ! [11/904] Index identifier
      character*1   l_nuturn    ! Limit flag on nuturn
      real*8        nuturn      ! (GHz) [5/300]? Frequency where the wind becomes
*                                 partially optically thick
      real*4        e_nuturn    ! (GHz) [0.9/60]? Uncertainty in {nu}turn [GHz]
      real*4        Fturn       ! (mJy) [0.05/24.5] Flux density at nuturn
      real*4        e_Fturn     ! (mJy) [0.01/1] Uncertainty in Fturn
      real*8        Fdust       ! (mJy) [0/274]? Dust flux density
      real*4        e_Fdust     ! (mJy) [0.06/18]? Uncertainty in Fdust
      real*4        S35         ! (mJy) [0.05/20.4] Flux density at 35GHz; Band 1
      real*4        S80         ! (mJy) [0.05/18.7] Flux density at 80GHz; Band 2
      real*4        S100        ! (mJy) [0.04/18.3] Flux density at 100GHz; Band 3
      real*4        S150        ! (mJy) [0.04/17.6] Flux density at 150GHz; Band 4
      real*4        S200        ! (mJy) [0.04/17.1] Flux density at 200GHz; Band 5
      real*4        S230        ! (mJy) [0.04/16.9] Flux density at 230GHz; Band 6
      real*4        S345        ! (mJy) [0.04/16.2] Flux density at 345GHz; Band 7
      real*4        S450        ! (mJy) [0.04/15.8] Flux density at 450GHz; Band 8
      real*4        S650        ! (mJy) [0.04/15.2] Flux density at 650GHz; Band 9
      real*4        S850        ! (mJy) [0.04/14.8] Flux density at 850GHz; Band 10

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

C  Declarations for 'table4.dat'	! Variability of ONC sources

      integer*4 nr__2
      parameter (nr__2=147)	! Number of records
      character*31 ar__2  	! Full-size record

      integer*4     Seq_2       ! [11/897] Index number
      character*1   v1_3Var     ! 1.3cm variable?
      integer*4     DelF1_3     ! (%) [19/905]? {Delta}F/F at 1.3cm
      integer*4     e_DelF1_3   ! (%) [3/574]? DelF1.3 uncertainty
      character*1   v3_6Var     ! 3.6cm variable?
      integer*4     DelF3_6     ! (%) [20/215]? {Delta}F/F at 3.6cm
      integer*4     e_DelF3_6   ! (%) [3/46]? DelF3.6 uncertainty
      character*1   v6Var       ! 6cm variable?
      integer*4     DelF6       ! (%) [47/102]? {Delta}F/F at 6cm
      integer*4     e_DelF6     ! (%) [3/11]? DelF6 uncertainty

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

C  Loading file 'table2.dat'	! Source detections, identifications, and fluxes

C  Format for file interpretation

    1 format(
     +  I3,1X,A8,1X,I3,1X,A2,1X,I2,1X,A7,1X,I1,1X,I2,1X,F5.2,1X,A1,I1,
     +  1X,I2,1X,F5.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F6.3,
     +  1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
     +  F9.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,909
        read(1,'(A151)')ar__
        read(ar__,1)
     +  Seq,v_RRS2008_,v_HC2000_,GMR,v_ZRK2004_,OName,RAh,RAm,RAs,DE_,
     +  DEd,DEm,DEs,F6,e_F6,F3_6,e_F3_6,F1_3a,e_F1_3a,F1_3b,e_F1_3b,
     +  F1_3c,e_F1_3c,F1_3d,e_F1_3d,v_F1_3_,e__F1_3_
        if(ar__(14:16) .EQ. '') v_HC2000_ = iNULL__
        if(ar__(21:22) .EQ. '') v_ZRK2004_ = iNULL__
        if(ar__(68:73) .EQ. '') F3_6 = rNULL__
        if(ar__(75:79) .EQ. '') e_F3_6 = rNULL__
        if(ar__(81:86) .EQ. '') F1_3a = rNULL__
        if(ar__(88:93) .EQ. '') e_F1_3a = rNULL__
        if(ar__(95:100) .EQ. '') F1_3b = rNULL__
        if(ar__(102:107) .EQ. '') e_F1_3b = rNULL__
        if(ar__(109:114) .EQ. '') F1_3c = rNULL__
        if(ar__(116:121) .EQ. '') e_F1_3c = rNULL__
        if(ar__(123:128) .EQ. '') F1_3d = rNULL__
        if(ar__(130:134) .EQ. '') e_F1_3d = rNULL__
        if(ar__(136:141) .EQ. '') v_F1_3_ = rNULL__
        if(ar__(143:151) .EQ. '') e__F1_3_ = rNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Seq,v_RRS2008_,v_HC2000_,GMR,v_ZRK2004_,OName,RAh,RAm,RAs,DE_,
     +  DEd,DEm,DEs,F6,e_F6,F3_6,e_F3_6,F1_3a,e_F1_3a,F1_3b,e_F1_3b,
     +  F1_3c,e_F1_3c,F1_3d,e_F1_3d,v_F1_3_,e__F1_3_
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Free-free emission model parameters and ALMA
*                                 band fluxes for detected sources

C  Format for file interpretation

    2 format(
     +  I3,1X,A1,F7.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F6.3,1X,F5.2,
     +  1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
     +  F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,175
        read(1,'(A107)')ar__1
        read(ar__1,2)
     +  Seq_1,l_nuturn,nuturn,e_nuturn,Fturn,e_Fturn,Fdust,e_Fdust,
     +  S35,S80,S100,S150,S200,S230,S345,S450,S650,S850
        if(ar__1(6:12) .EQ. '') nuturn = rNULL__
        if(ar__1(14:19) .EQ. '') e_nuturn = rNULL__
        if(ar__1(34:40) .EQ. '') Fdust = rNULL__
        if(ar__1(42:47) .EQ. '') e_Fdust = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Seq_1,l_nuturn,nuturn,e_nuturn,Fturn,e_Fturn,Fdust,e_Fdust,
     +  S35,S80,S100,S150,S200,S230,S345,S450,S650,S850
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Variability of ONC sources

C  Format for file interpretation

    3 format(I3,1X,A1,1X,I3,1X,I3,1X,A1,1X,I3,1X,I2,1X,A1,1X,I3,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,147
        read(1,'(A31)')ar__2
        read(ar__2,3)
     +  Seq_2,v1_3Var,DelF1_3,e_DelF1_3,v3_6Var,DelF3_6,e_DelF3_6,
     +  v6Var,DelF6,e_DelF6
        if(ar__2(7:9) .EQ. '') DelF1_3 = iNULL__
        if(ar__2(11:13) .EQ. '') e_DelF1_3 = iNULL__
        if(ar__2(17:19) .EQ. '') DelF3_6 = iNULL__
        if(ar__2(21:22) .EQ. '') e_DelF3_6 = iNULL__
        if(ar__2(26:28) .EQ. '') DelF6 = iNULL__
        if(ar__2(30:31) .EQ. '') e_DelF6 = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  Seq_2,v1_3Var,DelF1_3,e_DelF1_3,v3_6Var,DelF3_6,e_DelF3_6,
     +  v6Var,DelF6,e_DelF6
c    .......End.of.Just test output...........
      end do
      close(1)

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