FORTRAN Generation
(/./ftp/cats/J/ApJ/772/45)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/772/45 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-07
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/772/45   Hi-GAL obs.: star formation in the third quadrant   (Elia+, 2013)
*================================================================================
*The first Hi-GAL observations of the outer galaxy: a look at star formation in
*the third galactic quadrant in the longitude range 216.5{deg}<~l<~225.5{deg}.
*    Elia D., Molinari S., Fukui Y., Schisano E., Olmi L., Veneziani M.,
*    Hayakawa T., Pestalozzi M., Schneider N., Benedettini M., di Giorgio A.M.,
*    Ikhenaode D., Mizuno A., Onishi T., Pezzuto S., Piazzo L., Polychroni D.,
*    Rygl K.L.J., Yamamoto H., Maruccia Y.
*   <Astrophys. J., 772, 45 (2013)>
*   =2013ApJ...772...45E
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! List of the clumps detected in the CO(1-0) line

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

C  Position composed of: GLON GLAT
      integer*4     Seq         ! [1/321] Clump running sequence number
      real*8        GLON        ! (deg) Galactic longitude
      real*8        GLAT        ! (deg) Galactic latitude
      real*4        VLSR        ! (km/s) [0/65] Local Standard of Rest velocity
      real*4        sigma       ! (km/s) [0.4/3] Velocity dispersion (1)
      real*4        Dist        ! (kpc) [0.2/7.5] Heliocentric distance (2)
      real*4        Rad         ! (pc) [0.3/20] Deconvolved radius (3)
      integer*4     Mass        ! (Msun) [4/47071] Clump mass (4)
      integer*4     Mvir        ! (Msun) [283/72960] Virial mass (5)
      real*4        Lum         ! (K.km/s.pc2) [0.6/6139] The CO(1-0) line luminosity (6)
*Note (1): Velocity dispersion is derived as the second moment of the intensity
*          distribution along the velocity dimension.
*Note (2): To derive the heliocentric distances, the Galactic rotation curve of
*          Brunthaler et al. (2011AN....332..461B, R_0_=8.3kpc,
*          {theta}_0_=239km/s) has been used instead of the default one of
*          Reid et al. (2009ApJ...700..137R) implemented in CLOUDPROPS.
*Note (3): The deconvolved radius r is estimated as the geometric mean of the
*          deconvolved major and minor axes of the clump.
*Note (4): The mass of the clump is obtained as in Section 3.1, i.e.,
*          integrating the column density over the solid angle of the clump.
*          See also section 4.2 for further details.
*Note (5): The virial mass is calculated as:
*          M_vir_[M_{sun}_]=1040{sigma}_v_^2^[km/s]r[pc], under the assumption
*          that each clump is spherical and virialized with a density profile of
*          the form {rho}{propto}r^-1^ (Solomon et al. 1987ApJ...319..730S).
*          Density power laws with exponents 0 and -2 would require, instead, a
*          constant of 1150 and 690M_{sun}_(km/s)^-1^/pc, respectively.
*Note (6): The conversion factor to solar luminosities for this transition is
*          4.9x10^-5^L_{sun}_[K.km/s.pc^2^]^-1^
*          (Solomon et al. 1992ApJ...398L..29S).

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

C  Declarations for 'table3.dat'	! List of the analyzed Hi-GAL sources

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

C  Position composed of: GLON GLAT
      integer*4     Seq_1       ! [1/943] Source running sequence number
      real*8        GLON_1      ! (deg) Galactic longitude
      real*8        GLAT_1      ! (deg) Galactic latitude
      real*4        Dist_1      ! (kpc) [0.4/7.5] Distance
      real*4        F70         ! (Jy) ? Herschel/PACS 70 micron band flux density
      real*4        e_F70       ! (Jy) ? Uncertainty in F70
      real*4        F160        ! (Jy) ? Herschel/PACS 160 micron band flux density
      real*4        e_F160      ! (Jy) ? Uncertainty in F160
      real*4        F250        ! (Jy) [0.1/223] Herschel/SPIRE 250um band flux density
      real*4        e_F250      ! (Jy) [0.02/38] Uncertainty in F250
      real*4        F350        ! (Jy) Herschel/SPIRE 350 micron band flux density
      real*4        e_F350      ! (Jy) Uncertainty in F350
      real*4        F500        ! (Jy) ? Herschel/SPIRE 500 micron band flux density
      real*4        e_F500      ! (Jy) ? Uncertainty in F500
      real*4        Rad_1       ! (arcsec) [1.2/79] Angular radius
      real*4        Mass_1      ! (Msun) [0/729] Mass
      real*4        Temp        ! (K) [7.1/44] Temperature
      real*4        F22         ! (Jy) [0.01/475]? WISE 22 micron band flux density
      real*4        Lbol        ! (Lsun) [0/9680] Bolometric luminosity
      integer*4     Flag        ! [0/2]? Source flag (1)
*Note (1): Flag as follows:
*    0 = unbound;
*    1 = pre-stellar;
*    2 = proto-stellar.

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

C  Loading file 'table2.dat'	! List of the clumps detected in the CO(1-0) line

C  Format for file interpretation

    1 format(
     +  I3,1X,F9.5,1X,F8.5,1X,F4.1,1X,F3.1,1X,F3.1,1X,F4.1,1X,I5,1X,
     +  I5,1X,F6.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,321
        read(1,'(A59)')ar__
        read(ar__,1)Seq,GLON,GLAT,VLSR,sigma,Dist,Rad,Mass,Mvir,Lum
c    ..............Just test output...........
        write(6,1)Seq,GLON,GLAT,VLSR,sigma,Dist,Rad,Mass,Mvir,Lum
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! List of the analyzed Hi-GAL sources

C  Format for file interpretation

    2 format(
     +  I3,1X,F9.5,1X,F8.5,1X,F3.1,1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,
     +  F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.1,1X,F5.1,
     +  1X,F4.1,1X,F6.2,1X,F6.1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,943
        read(1,'(A119)')ar__1
        read(ar__1,2)
     +  Seq_1,GLON_1,GLAT_1,Dist_1,F70,e_F70,F160,e_F160,F250,e_F250,
     +  F350,e_F350,F500,e_F500,Rad_1,Mass_1,Temp,F22,Lbol,Flag
        if(ar__1(28:33) .EQ. '') F70 = rNULL__
        if(ar__1(35:38) .EQ. '') e_F70 = rNULL__
        if(ar__1(40:45) .EQ. '') F160 = rNULL__
        if(ar__1(47:50) .EQ. '') e_F160 = rNULL__
        if(ar__1(77:81) .EQ. '') F500 = rNULL__
        if(ar__1(83:87) .EQ. '') e_F500 = rNULL__
        if(ar__1(105:110) .EQ. '') F22 = rNULL__
        if(ar__1(119:119) .EQ. '') Flag = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  Seq_1,GLON_1,GLAT_1,Dist_1,F70,e_F70,F160,e_F160,F250,e_F250,
     +  F350,e_F350,F500,e_F500,Rad_1,Mass_1,Temp,F22,Lbol,Flag
c    .......End.of.Just test output...........
      end do
      close(1)

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