FORTRAN Generation
(/./ftp/cats/J/ApJS/251/26)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/251/26 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/ApJS/251/26      CO outflow cand. toward W3/4/5 complex. II.      (Li+, 2020)
*================================================================================
*CO outflow candidates toward the W3/4/5 complex.
*II. Feedback from candidate outflows.
*    Li Y., Xu Ye, Sun Y., Yang Ji
*   <Astrophys. J. Suppl. Ser., 251, 26 (2020)>
*   =2020ApJS..251...26L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Correlation among clouds and outflow candidates

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

C  Position composed of: GLON GLAT
      character*3   Reg         ! sub-region (1)
      character*6   Seq         ! Cloud Index
      real*8        GLON        ! (deg) [129.9/140]? Degree of Galactic Longitude
      real*8        GLAT        ! (deg) [-4.5/4.9]? Degree of Galactic Latitude
      real*4        Vlow        ! (km/s) [-66/0]? Lower velocity of integral interval
      real*4        Vhigh       ! (km/s) [-58/8]? Upper velocity of integral interval
      real*4        Mflow       ! (Msun) [0.03/186.3] Mass of outflow candidate (2)
      real*8        Pflow       ! (Msun.km/s) [0.05/1470.5] Momentum of outflow candidate (2)
      real*4        Eflow       ! (10-7J) [9.2e+41/1.4e+47] Kinetic energy of outflow
*                                 candidate in units of erg (2)
      real*4        Lflow       ! (10-7W) [2.9e+29/4.5e+34] Luminosity of outflow
*                                 candidate in units of erg/s (2)
      character*229 COI         ! Candidate outflow index(es)
*Note (1): Sub-regions as follows:
*   Per = Perseus arm (cloud indexes 1-136 and "Median")
*   Loc = Local arm (cloud indexes 137-260 and "Median")
*   Int = Interarm 1 (cloud indexes 261-265 and "Median")
*Note (2): The total outflow mass, momentum, kinetic energy and luminosity
*    corresponding to each cloud. They are from (Li+ 2019, J/ApJS/242/19),
*    but here we newly apply the corrections for inclination angle,
*    blending effect and optical depth. All values have been multiplied by
*    1.36 to take into consideration the mean molecular weight
*    (Brunt 2010A&A...513A..67B) instead of the weight of only molecular H2
*    given in (Li+ 2019, J/ApJS/242/19). The Median values are presented in
*    the end of these three sub-regions with index of "Median".

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

C  Declarations for 'table2.dat'	! Global physical properties of the clouds

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

      character*3   Reg_1       ! sub-region (G1)
      character*6   Seq_1       ! Cloud Index
      real*8        Mcloud      ! (Msun) [0.1/44942.5] Mass of cloud
      real*4        Rcloud      ! (pc) [0.04/17.12] Radius of cloud
      real*4        DVcloud     ! (km/s) [0.5/4] Line FWHM of cloud
      real*4        e_DVcloud   ! (km/s) [0.02/0.2] Uncertainty in DVcloud
      real*4        Tcloud      ! (K) [5.3/18.3] Excitation temperature of
*                                       cloud
      real*4        Taucloud    ! [0.1/0.7] Optical Depth of cloud
      real*4        Eturb       ! (10-7J) [6.6e+41/3.7e+48] Turbulent energy of
*                                       cloud in units of erg
      real*8        Pturb       ! (Msun.km/s) [0.1/129467.5] Turbulent momentum of
*                                       cloud
      real*4        tff         ! (Myr) [0.5/5.5] Free-fall timescale of cloud
      real*4        tdiss       ! (Myr) [0.2/7.1] Turbulent dissipation
*                                       timescale of cloud
      real*4        Lturb       ! (10-7W) [1.7e+28/1.4e+35] Turbulent dissipation
*                                       rate of cloud in units of erg/s
      real*4        vesc        ! (km/s) [0.1/4.8] Escape velocity of cloud
      real*4        Mesc        ! (Msun) [0.2/309.4] Escape mass of cloud
      real*4        Egrav       ! (10-7J) [1.3e+40/1e+49] Gravitational binding
*                                       energy of cloud in units of erg

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

C  Declarations for 'table3.dat'	! The effect of feedback of outflow candidates on
                              their parent cloud

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

      character*3   Reg_2       ! sub-region (G1)
      character*6   Seq_2       ! Cloud Index
      real*4        Ef_Et       ! [0.0008/47] Ratio of outflow to turbulent
*                                energies
      real*4        Pf_Pt       ! [0.00039/12] Ratio of outflow to turbulent
*                                momenta
      real*4        Lf_Lt       ! [0.003/350] Ratio of outflow to turbulent
*                                luminosities
      real*4        Ef_Eg       ! [0.0017/2400] Ratio of outflow to gravitational
*                                binding energies
      real*4        Me_Mc       ! [0.0005/84] Ratio of escape to total cloud
*                                masses
      real*4        Me_Mf       ! [1.6/34] Ratio of escape to outflow masses

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

C  Loading file 'table1.dat'	! Correlation among clouds and outflow candidates

C  Format for file interpretation

    1 format(
     +  A3,1X,A6,1X,F7.3,1X,F7.3,1X,F5.1,1X,F5.1,1X,F6.2,1X,F7.2,1X,
     +  E8.2,1X,E8.2,1X,A229)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,268
        read(1,'(A301)')ar__
        read(ar__,1)
     +  Reg,Seq,GLON,GLAT,Vlow,Vhigh,Mflow,Pflow,Eflow,Lflow,COI
        if(ar__(12:18) .EQ. '') GLON = rNULL__
        if(ar__(20:26) .EQ. '') GLAT = rNULL__
        if(ar__(28:32) .EQ. '') Vlow = rNULL__
        if(ar__(34:38) .EQ. '') Vhigh = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Reg,Seq,GLON,GLAT,Vlow,Vhigh,Mflow,Pflow,Eflow,Lflow,COI
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Global physical properties of the clouds

C  Format for file interpretation

    2 format(
     +  A3,1X,A6,1X,F7.1,1X,F5.2,1X,F3.1,1X,F4.2,1X,F4.1,1X,F3.1,1X,
     +  E7.1,1X,F8.1,1X,F3.1,1X,F3.1,1X,E7.1,1X,F3.1,1X,F5.1,1X,E7.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,277
        read(1,'(A93)')ar__1
        read(ar__1,2)
     +  Reg_1,Seq_1,Mcloud,Rcloud,DVcloud,e_DVcloud,Tcloud,Taucloud,
     +  Eturb,Pturb,tff,tdiss,Lturb,vesc,Mesc,Egrav
c    ..............Just test output...........
        write(6,2)
     +  Reg_1,Seq_1,Mcloud,Rcloud,DVcloud,e_DVcloud,Tcloud,Taucloud,
     +  Eturb,Pturb,tff,tdiss,Lturb,vesc,Mesc,Egrav
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! The effect of feedback of outflow candidates on
*                              their parent cloud

C  Format for file interpretation

    3 format(A3,1X,A6,1X,E7.1,1X,E7.1,1X,E7.1,1X,E7.1,1X,E7.1,1X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,277
        read(1,'(A55)')ar__2
        read(ar__2,3)Reg_2,Seq_2,Ef_Et,Pf_Pt,Lf_Lt,Ef_Eg,Me_Mc,Me_Mf
c    ..............Just test output...........
        write(6,3)Reg_2,Seq_2,Ef_Et,Pf_Pt,Lf_Lt,Ef_Eg,Me_Mc,Me_Mf
c    .......End.of.Just test output...........
      end do
      close(1)

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