FORTRAN Generation
(/./home/cats/J/A_A/618/A144)

Conversion of standardized ReadMe file for file /./home/cats/J/A_A/618/A144 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-May-19
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/A+A/618/A144     QSO candidates catalog with APOP & ALLWISE (QCC) (Guo+, 2018)
*================================================================================
*Identifying quasars with astrometric and photometric methods from
*APOP and ALLWISE.
*    Guo S., Qi Z., Liao S., Cao Z., Lattanzi M.G., Bucciarelli B., Tang Z.,
*    Yan Q.-Z.
*    <Astron. Astrophys. 618, A144 (2018)>
*    =2018A&A...618A.144G        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'qcc.dat'	! Quasar candidates

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

C  J2000.0 position composed of: RAdeg DEdeg
      integer*4     ID1        (nr__) ! ID1 of the unique union id (ID1,ID2) (1)
      integer*4     ID2        (nr__) ! ID2 of the unique union id (ID1,ID2) (1)
      real*8        RAdeg      (nr__) ! (deg) Right ascension (J2000.0) (2)
      real*8        DEdeg      (nr__) ! (deg) Declination (J2000.0) (2)
      integer*4     e_RAdeg    (nr__) ! (mas) Error of right ascension
      integer*4     e_DEdeg    (nr__) ! (mas) Error of declination
      real*4        pmRA       (nr__) ! (mas/yr) Proper motion in RA direction
      real*4        pmDE       (nr__) ! (mas/yr) Proper motion in DEcdirection
      real*4        e_pmRA     (nr__) ! (mas/yr) Error of proper motion in ra direction
      real*4        e_pmDE     (nr__) ! (mas/yr) Error of proper motion in dec direction
      real*8        Rmag       (nr__) ! (mag) ?=-32.768 R photographic magnitude (3)
      real*8        Bmag       (nr__) ! (mag) ?=-32.768 B photographic magnitude (3)
      real*8        Vmag       (nr__) ! (mag) ?=-32.768 V photographic magnitude (3)
      real*8        Imag       (nr__) ! (mag) ?=-32.768 I photographic magnitude (3)
      real*8        Jmag       (nr__) ! (mag) ?=-32.768 J photographic magnitude (3)
      real*8        Hmag       (nr__) ! (mag) ?=-32.768 H photographic magnitude (3)
      real*8        Kmag       (nr__) ! (mag) ?=-32.768 K photographic magnitude (3)
      real*8        W1mag      (nr__) ! (mag) ?=-32.768 W1 photographic magnitude (3)
      real*8        W2mag      (nr__) ! (mag) ?=-32.768 W2 photographic magnitude (3)
      real*8        W3mag      (nr__) ! (mag) ?=-32.768 W3 photographic magnitude (3)
      real*8        W4mag      (nr__) ! (mag) ?=-32.768 W4 photographic magnitude (3)
      character*1   QSOflag    (nr__) ! [sgo] QSO in SDSS (s) or in LAMOST (g) (4)
      real*8        z          (nr__) ! ?=-9.9 Redshift from SDSS and LAMOST,
      character*1   zflag      (nr__) ! [sgo] Original of redshift:
*                                    SDSS (s) or LAMOST (g) (4)
      character*1   Morflag    (nr__) ! [pe] The classification codes are
*                                     'p' = point-like or 'e' = extend source (5)
*Note (1): The combination of ID1 and ID2 is the source unique identifier for
*  object in the QCC.
*  Here the ID1 is MasterID and the ID2 is Running_Number in HEALPix
*  (the Hierarchical, Equal Area, and iso-Latitude Pixelisation of the sphere)
*  system.
*Note (2): The positions and proper motions are extracted directly from the
*  catalog of absolute proper motions outside The galactic plane (APOP),
*  (APOP; Qi et al. 2012, J/AJ/150/A137).
*Note (3): The photographic magnitudes (R, B, V, I, J, H, K) are taken from APOP.
* The W1, W2, W3, W4 magnitudes are extracted from ALLWISE (AllWISE Data
*  Release, Cutri R.M., et al., IPAC/Caltech (2013, Cat. II/328).
*  Default values for other cases(null) = -32.726.
*Note (4): Codes as follows:
*  s = SDSS
*  g = LAMOST
*  o = default values for other cases
*Note (5): The Morflag in APOP is extracted from GSC2.3, as follows:
*  p = stellar objects in GSC2.3
*  e = nonstellar objects in GSC2.3
*
*  There are many faint point-like sources in nonstellar classification.
*  This flag should be used carefully.

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

C  Loading file 'qcc.dat'	! Quasar candidates

C  Format for file interpretation

    1 format(
     +  I6,1X,I6,1X,F14.10,1X,F14.10,1X,I3,1X,I3,1X,F6.2,1X,F6.2,2X,
     +  F5.2,2X,F5.2,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,
     +  1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,A1,2X,F7.4,1X,A1,
     +  1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'qcc.dat')
      write(6,*) '....Loading file: qcc.dat'
      do i__=1,662753
        read(1,'(A182)')ar__
        read(ar__,1)
     +  ID1(i__),ID2(i__),RAdeg(i__),DEdeg(i__),e_RAdeg(i__),
     +  e_DEdeg(i__),pmRA(i__),pmDE(i__),e_pmRA(i__),e_pmDE(i__),
     +  Rmag(i__),Bmag(i__),Vmag(i__),Imag(i__),Jmag(i__),Hmag(i__),
     +  Kmag(i__),W1mag(i__),W2mag(i__),W3mag(i__),W4mag(i__),
     +  QSOflag(i__),z(i__),zflag(i__),Morflag(i__)
c    ..............Just test output...........
        write(6,1)
     +  ID1(i__),ID2(i__),RAdeg(i__),DEdeg(i__),e_RAdeg(i__),
     +  e_DEdeg(i__),pmRA(i__),pmDE(i__),e_pmRA(i__),e_pmDE(i__),
     +  Rmag(i__),Bmag(i__),Vmag(i__),Imag(i__),Jmag(i__),Hmag(i__),
     +  Kmag(i__),W1mag(i__),W2mag(i__),W3mag(i__),W4mag(i__),
     +  QSOflag(i__),z(i__),zflag(i__),Morflag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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