FORTRAN Generation
(/./ftp/cats/J/ApJ/755/20)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/755/20 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-11
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/755/20     YSOs population of Sh2-294 with Spitzer/IRAC     (Samal+, 2012)
*================================================================================
*Star formation and young population of the H II complex Sh2-294.
*    Samal M.R., Pandey A.K., Ojha D.K., Chauhan N., Jose J., Pandey B.
*   <Astrophys. J., 755, 20 (2012)>
*   =2012ApJ...755...20S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Photometric data in IRAC bands
                            (97 different positions)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg       ! (deg) [109.0/109.2] Right Ascension (J2000)
      real*8        DEdeg       ! (deg) [-9.5/-9.3] Declination (J2000)
      real*4        v3_6mag     ! (mag) [9.6/15] Spitzer/IRAC 3.6 micron band magnitude
      real*4        e_3_6mag    ! (mag) Uncertainty in 3.6mag
      real*4        v4_5mag     ! (mag) [9.2/14.8] Spitzer/IRAC 4.5um band magnitude
      real*4        e_4_5mag    ! (mag) Uncertainty in 4.5mag
      real*4        v5_8mag     ! (mag) [7.8/14.1] Spitzer/IRAC 5.8um band magnitude
      real*4        e_5_8mag    ! (mag) Uncertainty in 5.8mag
      real*4        v8_0mag     ! (mag) [8.6/12.4]? Spitzer/IRAC 8.0um band magnitude
      real*4        e_8_0mag    ! (mag) [0.02/0.2]? Uncertainty in 8.0mag

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

C  Declarations for 'table2.dat'	! Inferred physical parameters from SED fits to YSOs

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

C  J2000 position composed of: RAdeg DEdeg
      character*2   ID          ! YSO identification (1)
      real*8        RAdeg_1     ! (deg) [109/109.2] Right ascension (J2000)
      real*8        DEdeg_1     ! (deg) [-9.5/-9.3] Declination (J2000)
      real*4        M           ! (Msun) [0.6/8.8] Stellar mass
      real*4        e_M         ! (Msun) [0.1/2] Mass uncertainty
      real*4        T           ! (kK) [3.4/16.1] Stellar temperature
      real*4        e_T         ! (kK) [0.1/5] Temperature uncertainty
      real*4        Age         ! (Myr) [0.04/8.1] Stellar age (t*)
      real*4        e_Age       ! (Myr) [0.04/3.3] Age uncertainty
      real*4        Mdisk       ! (Msun) [0/0.2] Disk mass
      real*4        e_Mdisk     ! (Msun) [0/0.3] Mdisk uncertainty
      real*4        dM_dt       ! (10-8Msun/yr) [0/390] Disk accretion rate (dM/dt_disk_)
      real*4        e_dM_dt     ! (10-8Msun/yr) [0/383] dM/dt uncertainty
      real*4        Av          ! (mag) [3.1/23.7] Foreground visual absorption
      real*4        e_Av        ! (mag) [0.1/6.7] Av uncertainty
      real*4        chi2        ! [0.8/35.9] Best fit {chi}^2^_min_
*Note (1): The IRAC-classified Class II, Class I/II, and Class I YSOs are
*  labeled with the numbers 1-25, T1-T5, and C1-C3, respectively (see Figure 3).

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

C  Loading file 'table1.dat'	! Photometric data in IRAC bands
*                            (97 different positions)

C  Format for file interpretation

    1 format(
     +  F9.5,2X,F8.5,2X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F6.3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,103
        read(1,'(A72)')ar__
        read(ar__,1)
     +  RAdeg,DEdeg,v3_6mag,e_3_6mag,v4_5mag,e_4_5mag,v5_8mag,
     +  e_5_8mag,v8_0mag,e_8_0mag
        if(ar__(61:66) .EQ. '') v8_0mag = rNULL__
        if(ar__(68:72) .EQ. '') e_8_0mag = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  RAdeg,DEdeg,v3_6mag,e_3_6mag,v4_5mag,e_4_5mag,v5_8mag,
     +  e_5_8mag,v8_0mag,e_8_0mag
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Inferred physical parameters from SED fits to YSOs

C  Format for file interpretation

    2 format(
     +  A2,1X,F10.6,1X,F9.6,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,F5.3,1X,F5.3,1X,F6.2,1X,F6.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,33
        read(1,'(A97)')ar__1
        read(ar__1,2)
     +  ID,RAdeg_1,DEdeg_1,M,e_M,T,e_T,Age,e_Age,Mdisk,e_Mdisk,dM_dt,
     +  e_dM_dt,Av,e_Av,chi2
c    ..............Just test output...........
        write(6,2)
     +  ID,RAdeg_1,DEdeg_1,M,e_M,T,e_T,Age,e_Age,Mdisk,e_Mdisk,dM_dt,
     +  e_dM_dt,Av,e_Av,chi2
c    .......End.of.Just test output...........
      end do
      close(1)

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