FORTRAN Generation
(/./ftp/cats/J/ApJ/673/331)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/673/331 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 2024-Mar-28
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/673/331        X-ray observations of CG 12               (Getman+, 2008)
*================================================================================
*The stellar population and origin of the mysterious high-latitude star-forming
*cloud CG 12.
*    Getman K.V., Feigelson E.D., Lawson W.A., Broos P.S., Garmire G.P.
*   <Astrophys. J., 673, 331-353 (2008)>
*   =2008ApJ...673..331G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Basic source properties

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

C  J2000 position composed of: RAdeg DEdeg
      character*3   Fld         ! Chandra ACIS-I fields surveyed in the CG 12
*                                    region (G1)
      integer*4     Src         ! X-ray source number within the field,
*                                    [GFL2008] "Field" "Source" in Simbad
      character*18  CXOU        ! CXOUJ designation (HHMMSS.dd+DDMMSS.s)
      real*8        RAdeg       ! (deg) Right ascension in decimal degrees (J2000)
      real*8        DEdeg       ! (deg) Declination in decimal degrees (J2000)
      real*4        e_Pos       ! (arcsec) Positional uncertainty
      real*4        Theta       ! (arcmin) Off-axis angle
      real*4        Fnet        ! (ct) Estimated net counts from extracted area in
*                                    total energy band 0.5-8.0keV
      real*4        e_Fnet      ! (ct) rms uncertainty on ctNetF
      real*4        Fbg         ! (ct) Estimated background counts from extracted
*                                    area in total energy band
      real*4        Hnet        ! (ct) Estimated net counts from extracted area in
*                                    hard energy band 2-8keV
      real*4        Frac        ! Fraction of the PSF at the fiducial energy of
*                                    1.497keV enclosed within the extracted
*                                    area (2)
      real*4        Signi       ! Source significance in {sigma}
      character*2   Anom        ! [ge ] Source anomalies (3)
      character*1   Var         ! [abc] Variability characterization based on
*                                         K-S statistics (4)
      real*4        Texp        ! (ks) Source effective exposure
      real*4        Emed        ! (keV) Background-corrected median photon energy in
*                                    the total 0.5-8keV energy band
      real*4        logPB       ! ([-]) Logarithm of K-S probability that background
*                                    and source spectra are drawn from the same
*                                    distribution (5)
*Note (2): Note that a reduced PSF fraction (significantly below 90%) may
*     indicate that the source is in a crowded region or on the edge of the
*     field.
*Note (3): Source anomalies flags as follows:
*      g = fractional time that source was on a detector (FRACEXPO from mkarf)
*          is <0.9
*      e = source on field edge
*      p = source piled up
*      s = source on readout streak
*Note (4): Source variability flags as follows:
*      a = no evidence for variability
*      b = possibly variable
*      c = definitely variable.
*  blank = No test is performed for sources with fewer than 4 total full-band
*          counts. No value is reported for sources in chip gaps
*          or on field edges.
*Note (5): The value of logPB has been truncated at -4.0 for sources with
*     significantly different background and source spectra.

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

C  Declarations for 'table3.dat'	! X-ray spectroscopy for photometrically selected
                              sources: thermal plasma fits

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

      character*3   Fld_1       ! Chandra ACIS-I fields surveyed in the CG 12
*                                   region (G1)
      integer*4     Src_1       ! X-ray source number within the field
      character*18  CXOU_1      ! CXOUJ designation (HHMMSS.dd+DDMMSS.s)
      real*4        Fnet_1      ! (ct) Estimated net counts from extracted area in
*                                   total energy band 0.5-8.0keV
      real*4        Emed_1      ! (keV) Background-corrected median photon energy in
*                                   the total 0.5-8keV energy band
      real*4        logNH       ! ([cm-2]) Best-fit value or the column density (1)
      real*4        E_logNH     ! ([cm-2]) ? Error on logNH (upper limit)
      real*4        e_logNH_1   ! ([cm-2]) ? Error on logNH (lower limit)
      real*4        kT          ! (keV) Best-fit value for plasma temperature (1)
      real*4        E_kT        ! (keV) ? Error on kT (upper limit)
      real*4        e_kT_1      ! (keV) ? Error on kT (lower limit)
      real*4        logLh       ! ([10-7W]) Hard band (2-8keV) X-ray luminosity (2)
      real*4        e_logLh     ! ([10-7W]) rms uncertainty on logLh
      real*4        logLhc      ! ([10-7W]) Hard band (2-8keV) absorption-corrected
*                                   X-ray luminosity (2)
      real*4        e_logLhc    ! ([10-7W]) rms uncertainty on logLhc
      real*4        logLt       ! ([10-7W]) Total band (-0.5-8keV) X-ray luminosity (2)
      real*4        e_logLt     ! ([10-7W]) rms uncertainty on logLt
      real*4        logLtc      ! ([10-7W]) Total band (-0.5-8keV) absorption-corrected
*                                   X-ray luminosity (2)
      real*4        e_logLtc    ! ([10-7W]) rms uncertainty on logLtc
*Note (1): All fits used the "wabs(apec)" model in XSPEC and assumed
*     0.3Z_{sun}_ abundances. Uncertainties represent 90% confidence
*     intervals. More significant digits are used for uncertainties <0.1 in
*     order to avoid large rounding errors; for consistency, the same number
*     of significant digits is used for both lower and upper uncertainties.
*     Uncertainties are missing when XSPEC was unable to compute them or
*     when their values were so large that the parameter is effectively
*     unconstrained. Fitted values with logNH<20.0cm^-2^ are truncated at
*     20.0 because ACIS-I spectra are insensitive to differences in very 
*     low column densities. Fitted values with kT>15keV are truncated at
*     15keV because the data often cannot discriminate between very high
*     temperature values.
*Note (2): X-ray luminosities are calculated assuming a distance of 550pc

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

C  Declarations for 'table4.dat'	! Optical and NIR counterparts

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

      character*3   Fld_2       ! Chandra ACIS-I fields surveyed in the CG 12
*                                    region (G1)
      integer*4     Src_2       ! X-ray source number within the field
      character*18  CXOU_2      ! CXOUJ designation (HHMMSS.dd+DDMMSS.s)
      character*12  NOMAD       ! NOMAD identifier (NNNN-NNNNNNN)
      real*4        Sep1        ! (arcsec) ? Chandra-NOMAD positional offset
      character*5   v_MMB2004b_ ! Source identifier from Maheswar et al.
*                                    (2004, Cat. <J/MNRAS/355/1272>)
      real*4        Vmag        ! (mag) ? V magnitude (1)
      real*4        e_Vmag      ! (mag) ? rms uncertainty on Vmag
      real*4        Rmag        ! (mag) ? R magnitude (1)
      real*4        e_Rmag      ! (mag) ? rms uncertainty on Rmag
      real*4        pmRA        ! (mas/yr) ? Best proper motion along RA
*                                      (pmRA*cosDE) (2)
      real*4        e_pmRA      ! (mas/yr) ? rms uncertainty on pmRA
      real*4        pmDE        ! (mas/yr) ? Best proper motion along DE (2)
      real*4        e_pmDE      ! (mas/yr) ? rms uncertainty on pmDE
      character*16  v2MASS      ! 2MASS designation (HHMMSSss+DDMMSSs)
      real*4        Sep2        ! (arcsec) ? Chandra-2MASS positional offset
      real*4        Jmag        ! (mag) ? 2MASS J magnitude
      real*4        e_Jmag      ! (mag) ? rms uncertainty on Jmag
      real*4        Hmag        ! (mag) ? 2MASS H magnitude
      real*4        e_Hmag      ! (mag) ? rms uncertainty on Hmag
      real*4        Ksmag       ! (mag) ? 2MASS Ks magnitude
      real*4        e_Ksmag     ! (mag) ? rms uncertainty on Ksmag
      character*6   Flags       ! 2MASS photometry quality and
*                                    confusion-contamination flag
*Note (1): When formal statistical errors are given, V and R magnitudes are
*     from Maheswar et al. (2004, Cat. <J/MNRAS/355/1272>). V magnitudes for
*     three B-type stars, I 57, I 58, and I 74, are from Williams et al.
*     (1977MNRAS.180..709W). Otherwise, V and R magnitudes are from the
*     NOMAD catalog, in which case all V magnitudes except for the Chandra
*     source I 81 are from the YB6 (=USNO unpublished) catalog, and all R
*     magnitudes except for the Chandra source I 57 are from the USNO-B1.0
*     catalog (R magnitudes from USNO-B1.0 are second red magnitudes R_2_,
*     except for the Chandra sources I 70 and II 9, for which it is a first
*     red magnitude R_1_). For the Chandra source I 81 V magnitude is from
*     the Tycho-2 catalog, and for the Chandra source I 57 R magnitude is
*     from the UCAC2 catalog.
*Note (2): "Best" proper motions and their mean errors from catalogs
*     compiled in NOMAD (Cat. I/297)

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

C  Loading file 'table2.dat'	! Basic source properties

C  Format for file interpretation

    1 format(
     +  A3,1X,I2,1X,A18,1X,F10.6,1X,F10.6,1X,F4.2,1X,F5.2,1X,F5.1,1X,
     +  F4.1,1X,F4.1,1X,F5.1,1X,F4.2,1X,F4.1,1X,A2,1X,A1,1X,F4.1,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,128
        read(1,'(A111)')ar__
        read(ar__,1)
     +  Fld,Src,CXOU,RAdeg,DEdeg,e_Pos,Theta,Fnet,e_Fnet,Fbg,Hnet,
     +  Frac,Signi,Anom,Var,Texp,Emed,logPB
c    ..............Just test output...........
        write(6,1)
     +  Fld,Src,CXOU,RAdeg,DEdeg,e_Pos,Theta,Fnet,e_Fnet,Fbg,Hnet,
     +  Frac,Signi,Anom,Var,Texp,Emed,logPB
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! X-ray spectroscopy for photometrically selected
*                              sources: thermal plasma fits

C  Format for file interpretation

    2 format(
     +  A3,1X,I2,1X,A18,1X,F5.1,1X,F4.2,1X,F4.1,1X,F3.1,1X,F3.1,1X,
     +  F5.2,F3.1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,45
        read(1,'(A107)')ar__1
        read(ar__1,2)
     +  Fld_1,Src_1,CXOU_1,Fnet_1,Emed_1,logNH,E_logNH,e_logNH_1,kT,
     +  E_kT,e_kT_1,logLh,e_logLh,logLhc,e_logLhc,logLt,e_logLt,
     +  logLtc,e_logLtc
        if(ar__1(43:45) .EQ. '') E_logNH = rNULL__
        if(ar__1(47:49) .EQ. '') e_logNH_1 = rNULL__
        if(ar__1(56:58) .EQ. '') E_kT = rNULL__
        if(ar__1(60:63) .EQ. '') e_kT_1 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Fld_1,Src_1,CXOU_1,Fnet_1,Emed_1,logNH,E_logNH,e_logNH_1,kT,
     +  E_kT,e_kT_1,logLh,e_logLh,logLhc,e_logLhc,logLt,e_logLt,
     +  logLtc,e_logLtc
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Optical and NIR counterparts

C  Format for file interpretation

    3 format(
     +  A3,1X,I2,1X,A18,1X,A12,1X,F3.1,1X,A5,1X,F5.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,F5.1,1X,F4.1,1X,F5.1,1X,F4.1,1X,A16,1X,F3.1,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,128
        read(1,'(A153)')ar__2
        read(ar__2,3)
     +  Fld_2,Src_2,CXOU_2,NOMAD,Sep1,v_MMB2004b_,Vmag,e_Vmag,Rmag,
     +  e_Rmag,pmRA,e_pmRA,pmDE,e_pmDE,v2MASS,Sep2,Jmag,e_Jmag,Hmag,
     +  e_Hmag,Ksmag,e_Ksmag,Flags
        if(ar__2(40:42) .EQ. '') Sep1 = rNULL__
        if(ar__2(50:54) .EQ. '') Vmag = rNULL__
        if(ar__2(56:59) .EQ. '') e_Vmag = rNULL__
        if(ar__2(61:65) .EQ. '') Rmag = rNULL__
        if(ar__2(67:70) .EQ. '') e_Rmag = rNULL__
        if(ar__2(72:76) .EQ. '') pmRA = rNULL__
        if(ar__2(78:81) .EQ. '') e_pmRA = rNULL__
        if(ar__2(83:87) .EQ. '') pmDE = rNULL__
        if(ar__2(89:92) .EQ. '') e_pmDE = rNULL__
        if(ar__2(111:113) .EQ. '') Sep2 = rNULL__
        if(ar__2(115:119) .EQ. '') Jmag = rNULL__
        if(ar__2(121:124) .EQ. '') e_Jmag = rNULL__
        if(ar__2(126:130) .EQ. '') Hmag = rNULL__
        if(ar__2(132:135) .EQ. '') e_Hmag = rNULL__
        if(ar__2(137:141) .EQ. '') Ksmag = rNULL__
        if(ar__2(143:146) .EQ. '') e_Ksmag = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Fld_2,Src_2,CXOU_2,NOMAD,Sep1,v_MMB2004b_,Vmag,e_Vmag,Rmag,
     +  e_Rmag,pmRA,e_pmRA,pmDE,e_pmDE,v2MASS,Sep2,Jmag,e_Jmag,Hmag,
     +  e_Hmag,Ksmag,e_Ksmag,Flags
c    .......End.of.Just test output...........
      end do
      close(1)

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