FORTRAN Generation
(/./ftp/cats/J/ApJ/833/229)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/833/229 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-Aug-16
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/833/229  Star forming cloud-giant molecular cloud complexes  (Lee+, 2016)
*================================================================================
*Observational evidence of dynamic star formation rate in Milky Way giant
*molecular clouds.
*    Lee E.J., Miville-Deschenes M.-A., Murray N.W.
*   <Astrophys. J., 833, 229-229 (2016)>
*   =2016ApJ...833..229L    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! *MML16-SFC match (sorted by SFC luminosities)

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

C  Position composed of: GLONg GLATg
      integer*4     GMCno      (nr__) ! [13/19689] GMC record number (1)
      real*8        GLONg      (nr__) ! (deg) Galactic longitude of GMC
      real*4        Sigl       (nr__) ! (deg) [0.08/0.7] RMS dispersion in galactic longitude
*                                 of GMC ({sigma}_l_)
      real*8        GLATg      (nr__) ! (deg) Galactic latitude of GMC
      real*4        Sigb       (nr__) ! (deg) [0.08/0.6] RMS dispersion in galactic latitude
*                                 of GMC ({sigma}_b_)
      real*8        Vlsrg      (nr__) ! (km/s) [-122/113] LSR velocity centroid of GMC
      real*4        SigV       (nr__) ! (km/s) [0.8/16] RMS velocity dispersion of GMC
*                                 ({sigma}_V_)
      real*4        Distg      (nr__) ! (kpc) [0.4/16] Kinematic distance to GMC
      integer*4     SFCno      (nr__) ! [0/277] Matched SFC record number (1)
      real*8        GLONs      (nr__) ! (deg) Galactic longitude of SFC
      real*8        GLATs      (nr__) ! (deg) Galactic latitude of SFC
      real*4        Rad        (nr__) ! (deg) [0.01/2.2] Geometric mean radius of SFC
      real*8        Vlsrs      (nr__) ! (km/s) [-122/112] LSR mean velocity of SFC
      real*4        ExpV       (nr__) ! (km/s) [2.7/49] Expansion velocity of SFC
      real*4        Dists      (nr__) ! (kpc) [0.4/16] Kinematic distance to SFC
*Note (1): Coordinates and velocities of star-forming complexes (SFCs;
*          Lee+, 2012, J/ApJ/752/146) and MML16
*          (Miville-Deschenes+ 2017, J/ApJ/834/57) giant molecular clouds (GMCs)
*          matched to each other.

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

C  Declarations for 'table2.dat'	! SFC-GMC complexes (sorted by luminosities)

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

C  Position composed of: GLON GLAT
      integer*4     SFCno_1    (nr__1) ! [0/277] SFC record number (1)
      real*4        GLON       (nr__1) ! (deg) Mass-weighted mean galactic longitude
      real*4        GLAT       (nr__1) ! (deg) Mass-weighted mean galactic latitude
      real*8        Vmean      (nr__1) ! (km/s) [-122/110] Mass-weighted mean velocity
      real*4        SigV_1     (nr__1) ! (km/s) [1.2/22.2] Mass-weighted velocity dispersion
      real*4        Rang       (nr__1) ! (deg) [0.1/1] Geometric mean angular radius (R_ang_) (2)
      real*4        MaxAxis    (nr__1) ! (deg) [0.2/1.3] Maximum eigenvalue of
*                               moment of inertia matrix (R_max_)
      real*4        MinAxis    (nr__1) ! (deg) [0.1/1] Minimum eigenvalue of
*                               moment of inertia matrix (R_min_)
      real*4        Rgal       (nr__1) ! (kpc) [0.02/18] Mass-weighted mean galactocentric
*                               radius (3)
      real*4        Dist       (nr__1) ! (kpc) [0.4/16] Mass-weighted mean kinematic distance (3)
      real*4        Wco        (nr__1) ! (K.km/s) [23.9/9640] Total CO flux
      real*4        Mgas       (nr__1) ! (Msun) [1580/] Total gas mass
*Note (1): The order of SFCno appear shuffled compared to Table 1 because we
*          adopt the mass-weighted distances of matched GMCs for each SFC here.
*Note (2): The geometric mean angular radius is defined as
*           (MaxAxis*MinAxis^2^)^{1/3}^
*Note (3): We detail the calculations of the mass-weighted parameters in
*          Section 2.4.

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

C  Declarations for 'table3.dat'	! Star formation properties of SFC-GMC complexes
                                 and unmatched clouds more massive than
                                 10^4.75^Msun (sorted by ionizing luminosities)

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

      integer*4     SFCno_2    (nr__2) ! [0/277] SFC record number (SFC-GMC complex)
      real*4        SigV_2     (nr__2) ! (km/s) [1/22.2] Velocity dispersion
      real*4        Dist_1     (nr__2) ! (kpc) [0.4/16] Kinematic distance
      real*4        Rad_1      (nr__2) ! (pc) [2/166] Physical radius (1)
      real*4        Q          (nr__2) ! (s-1) Ionizing luminosity
      real*4        Mgas_1     (nr__2) ! (Msun) [1580/] Gas mass
      real*4        Sigma      (nr__2) ! (Msun/pc2) [4/696] Gas column density ({Sigma}_g_) (2)
      real*4        SFEbr      (nr__2) ! [/0.4] Star formation efficiency
*                                 ({epsilon}_br_) (3)
      real*4        SFRbr      (nr__2) ! [0.0001/1.3] Star formation rate per
*                                 free-fall time ({epsilon}_ff,br_) (3)
      real*4        tff        (nr__2) ! (10+6yr) [1.4/40] Free-fall time ((3{pi}/32G{rho})^0.5^)
      real*4        aVir       (nr__2) ! [0.2/26] Virial parameter ({alpha}_vir_) (4)
*Note (1): The physical radius is defined as Dist*tan(R_ang_).
*          The angular radius is defined in Section 2.4.
*Note (2): The gas surface density ({Sigma}_g_) is defined as
*          Mgas/(Dist({pi}*MaxAxis*MinAxis))^2^
*          where MaxAxis and MinAxis are the maximum and minimum eigenvalues of
*          the moment of inertia matrix for each cloud (see Section 2.4).
*Note (3): Star formation efficiency {epsilon} is defined as the ratio of the
*          mass in protostars to the total mass in a given star-forming region is
*          defined in Equation (1) as {epsilon}=M*/(Mgas+M*) where the
*          star-forming region may be a GMC or a smaller sub-region of a GMC.
*          The star formation rate per freefall time epsilonff is defined in
*          Equation (2): {epsilon}_ff_={epsilon}*(tff/{tau}*) where {tau}* is
*          the lifetime of the (proto-)stellar object in question, and tff is
*          the freefall time of the star-forming region (a GMC or its
*          sub-region), which is assumed to have a mean density {rho}.
*          The subscript "br" stands for bremsstrahlung (see Figure 3).
*Note (4): The virial parameter is defined as 5*SigV^2^*Rad/G*Mgas where
*          G is the gravitational constant.

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

C  Loading file 'table1.dat'	! *MML16-SFC match (sorted by SFC luminosities)

C  Format for file interpretation

    1 format(
     +  I5,1X,F8.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.2,1X,F5.2,1X,F5.2,1X,
     +  I3,1X,F8.4,1X,F7.4,1X,F4.2,1X,F7.2,1X,F5.2,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,389
        read(1,'(A103)')ar__
        read(ar__,1)
     +  GMCno(i__),GLONg(i__),Sigl(i__),GLATg(i__),Sigb(i__),
     +  Vlsrg(i__),SigV(i__),Distg(i__),SFCno(i__),GLONs(i__),
     +  GLATs(i__),Rad(i__),Vlsrs(i__),ExpV(i__),Dists(i__)
c    ..............Just test output...........
        write(6,1)
     +  GMCno(i__),GLONg(i__),Sigl(i__),GLATg(i__),Sigb(i__),
     +  Vlsrg(i__),SigV(i__),Distg(i__),SFCno(i__),GLONs(i__),
     +  GLATs(i__),Rad(i__),Vlsrs(i__),ExpV(i__),Dists(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! SFC-GMC complexes (sorted by luminosities)

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,191
        read(1,'(A75)')ar__1
        read(ar__1,2)
     +  SFCno_1(i__),GLON(i__),GLAT(i__),Vmean(i__),SigV_1(i__),
     +  Rang(i__),MaxAxis(i__),MinAxis(i__),Rgal(i__),Dist(i__),
     +  Wco(i__),Mgas(i__)
c    ..............Just test output...........
        write(6,2)
     +  SFCno_1(i__),GLON(i__),GLAT(i__),Vmean(i__),SigV_1(i__),
     +  Rang(i__),MaxAxis(i__),MinAxis(i__),Rgal(i__),Dist(i__),
     +  Wco(i__),Mgas(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Star formation properties of SFC-GMC complexes
*                                 and unmatched clouds more massive than
*                                 10^4.75^Msun (sorted by ionizing luminosities)

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,191
        read(1,'(A82)')ar__2
        read(ar__2,3)
     +  SFCno_2(i__),SigV_2(i__),Dist_1(i__),Rad_1(i__),Q(i__),
     +  Mgas_1(i__),Sigma(i__),SFEbr(i__),SFRbr(i__),tff(i__),
     +  aVir(i__)
c    ..............Just test output...........
        write(6,3)
     +  SFCno_2(i__),SigV_2(i__),Dist_1(i__),Rad_1(i__),Q(i__),
     +  Mgas_1(i__),Sigma(i__),SFEbr(i__),SFRbr(i__),tff(i__),
     +  aVir(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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