FORTRAN Generation
(/./ftp/cats/VII/283)

Conversion of standardized ReadMe file for file /./ftp/cats/VII/283 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-Oct-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. VII/283            The Million Quasars (Milliquas) catalog (6.3)  (Flesch, 2019)
*================================================================================
*The Million Quasars (MILLIQUAS) catalogue, version 6.3.
*   Flesch E.W.
*   <Pub. Astron. Soc. Australia 32, 10 (2015)>
*   =2015PASA...32...10F
*   =2019yCat.7283....0F
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! The Milliquas catalogue, V6.3 (16 June 2019)

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg      (nr__) ! (deg) Right ascension J2000 (degrees)(1)
      real*8        DEdeg      (nr__) ! (deg) Declination J2000 (degrees)(1)
      character*25  Name       (nr__) ! ID from the literature, or J2000 (2)
      character*4   Cl_Ass     (nr__) ! [ABKQqNGL 2RX] Classification of object,
*                                  and associations (3)
      real*4        Rmag       (nr__) ! (mag) ?=0 Red optical magnitude (4)
      real*4        Bmag       (nr__) ! (mag) ?=0 Blue optical magnitude (4)
      character*3   Comment    (nr__) ! Comment on optical object (5)
      character*1   R          (nr__) ! Red optical PSF class (6)
      character*1   B          (nr__) ! Blue optical PSF class (6)
      real*4        z          (nr__) ! ? Redshift from the literature or estimated (7)
      character*6   rName      (nr__) ! Citation for name (8)
      character*6   rz         (nr__) ! Citation for redshift (8)
      integer*4     Qpct       (nr__) ! (pct) ? Probability that this object is a QSO (9)
      character*22  XName      (nr__) ! X-ray ID, if any (10)
      character*22  RName_1    (nr__) ! Radio ID, if any (10)
      character*22  Lobe1      (nr__) ! Radio lobe ID or extra R/X ID, if any (10)
      character*22  Lobe2      (nr__) ! Radio lobe ID or extra X-ray ID, if any (10)
*Note (1): These are to 7 decimals which suits Gaia-DR2 astrometry.
*Note (2): Nameless radio/X-ray associated objects here display the J2000
*  position in HHMMSS.SS+DDMMSS.S for the convenience of the user.  If needing
*  a name for it, just preface the J2000 with "MQ", e.g., MQ J000001.5-251706.
*Note (3): Legend of class/type:
*  Q = QSO, type-I broad-line core-dominated, 595972 of these.
*  A = AGN, type-I Seyferts/host-dominated, 27032 of these.
*  B = BL Lac object, 1724 of these.
*  q = photometric quasar candidate from SDSS or WISE, 1266417 of these.
*       ~1K are SDSS/LAMOST pipeline QSOs w/o subclass, with R/X associations.
*  L = lensed quasar extra image, only 61 of these in this optical data.
*  K = NLQSO, type-II narrow-line core-dominated, 5424 of these.
*  N = NLAGN, type-II Seyferts/host-dominated, 30625 of these.  Includes an
*       unknown number of legacy NELGs/ELGs/LINERs for completeness.
*  G = galaxy with double radio lobes (so calculated), 3581 of these.
*  R = radio association displayed.
*  X = X-ray association displayed.
*  2 = double radio lobes displayed (declared by data-driven algorithm).
*Note (4): Optical data is from the APM (http://www.ast.cam.ac.uk/~mike/apmcat),
*  USNO-A & USNO-B 
*  (https://www.usno.navy.mil/USNO/astrometry/optical-IR-prod/usno-b1.0), 
*  and the SDSS (http://sdss.org).
*  APM/USNO-A magnitudes have been recalibrated from the original values as
*  documented in QORG (2004,A&A,427,387), so such USNO-A magnitudes are often
*  used in preference to USNO-B.  APM galaxies < mag 17 are usually shown
*  too bright due to PSF modelling.  Integer magnitudes (e.g., 22.00) are
*  estimates if both bands are integer or one band empty.  Note: many SDSS
*  magnitudes are "extinction-corrected" ~0.3 mag brighter than observed.
*Note (5): Legend as follows:
*  p = optical magnitudes are POSS-I O (violet 4050A) and E (red 6400A).
*       These are preferred because O is well-offset from E, and those plates
*       were always taken on the same night, thus the red-blue color is
*       correct even for variable objects.
*  j = optical magnitudes are SERC J (Bj 4850A) and R (red 6400A) from the
*       POSS-II or UKST surveys.  Red-blue color is less reliable because the
*       red & blue plates were taken in different epochs, i.e., years apart.
*  b = blue magnitude is Vega 4400A (Johnson), red is 6400A (Cousins).
*  g = blue magnitude is SDSS green 4900A, red is SDSS r 6200A.
*  u = blue magnitude is SDSS ultraviolet 3850A.
*  v = red magnitude is visual, i.e., white, 5500A midpoint.
*  i = red magnitude is infrared 7500A.
*  z = red magnitude is infrared z 8600A.
*  r = red magnitude is r 6200A.
*  (blank) = red alone is 6400A (Cousins), else they are r 6200A & g 4900A.
*  G = Gaia-DR2 astrometry shown, precessed to J2000 by CDS.  If 'G' is alone
*       then the magnitudes are Gaia RP & BP, or Gaia G if red band only.
*  + = variability nominally detected in both red/blue over multi-epoch data.
*  m = proper motion nominally detected, from USNO-B.
*  a = object is host-dominated with faint nuclear activity, such as an SDSS
*       pipeline galaxy with an AGN subclass or AGN-classed elsewhere, see its
*       citation.  Milliquas class is 'A' if BROADLINE, else 'N'. (see note 2)
*Note (6): The APM, USNO-B, and SDSS provide PSF class, albeit using different
*  criteria. These are shown here as follows:
*  - = point source / stellar PSF (APM notation: -1, here truncated)
*  1 = fuzzy / galaxy shape       (APM notation: 1 and some 2)
*  n = no PSF available, whether borderline or too faint to tell, etc.
*  x = not seen in this band (fainter than plate depth, or confused, etc.)
*Note (7):  Spectroscopic/grism redshifts are required for objects classified as
*  Q/A/K/N/L; photometric redshifts for objects classified as 'q' are rounded
*  here to 0.1 z, except if pipeline redshifts are available from SDSS/LAMOST.
*  The XDQSO catalog does not provide photometric redshifts, so those are
*  provided either by NBCKDE or by this catalog using the four-color method
*  detailed in Appendix 2 of my HMQ paper (2015PASA...32...10F,PASA,32,10).
*  AllWISE photometric redshifts are similarly calculated using the four colors
*  B-R, R-W1, W1-W2 and W2-W3, and are displayed where calculated as 50%+
*  (usually 80%+) likely to be true within 0.5z of the displayed redshift value.
*Note (8): Legend (with counts of name and redshift) and references:
*  (for numerical reference code, see refs.dat file):
*  2dF (327,233)         : 2dF galaxy survey, Colless M. et al.,
*                           2001MNRAS.328.1039C, Cat. VII/250
*  2MRS (10,34)          : 2MASS, Huchra J.P. et al.,
*                           http://tdc-www.harvard.edu/2mrs
*  2QZ (27520,24161)     : Croom S.M. et al., 2004MNRAS.349.1397C, Cat. VII/241
*  2SLAQ (10365,8709)    : Croom S.M. et al., 2009, Cat. J/MNRAS/392/19
*  3FGL (16,16)          : Fermi cleanups, Paiano S. et al., 2017ApJ...851..135P
*  3FGL2 (22,22)         : Fermi cleanups II, Paiano S. et al.,
*                           2019ApJ...871..162P
*  3XLSS (25,25)         : The XXL Survey, Pierre M. et al., 2016A&A...592A...1P,
*                           Cat, IX/49
*  4LAC (366,199)        : Fermi AGN v4 + LL, Fermi-LAT collab.,
*                           2019,arXiv:1905.10771
*  6dF (300,277)         : 6dF galaxy survey, Jones D.H. et al.,
*                           2009MNRAS.399..683J, Cat. VII/259
*  AAOz (1491,1498)      : AAOmega XXL-South: Lidman C. et al.,
*                           2016PASA...33....1L
*  AGES (2046,2046)      : AGES survey, Kochanek C.S. et al.,
*                            2012, Cat. J/ApJS/200/8
*  AGNELL (2,3)          : DES lenses, Agnello A. et al., 2015MNRAS.454.1260A
*  AGNELA (4,4)          : SDSS J1433+6007 4-lens, Agnello A. et al.,
*                           2018MNRAS.474.3391A
*  AGNEL2 (13,13)        : VST-Gaia QSO pairs, Agnello A. et al.,
*                           2018MNRAS.475.2086A
*  AKARI (1,1)           : overlooked luminous quasar, Aoki K. et al.,
*                           2011PASJ...63S.457A
*  ALMA (4,4)            : ALMA hi-z, Roberto Decarli R. et al.,
*                           2018ApJ...854...97D
*  ANGUIT (1,1)          : COSMOS lens, Anguita T. et al., 2009A&A...507...35A
*  ATel (5,5)            : Astronomers Telegraph posts,
*                           http://www.astronomerstelegram.org
*  ATLAS (229,269)       : Mao M.Y. et al., 2012, Cat. J/MNRAS/426/3334
*  BAHM (24,24)          : dust-reddened QSOs, Banerji M. et al.,
*                           2015MNRAS.447.3368B
*  BASS (16,114)         : Swift-BAT AGN, Koss M. et al., 2017, Cat. J/ApJ/850/74
*  BERGHE (1,1)          : Pan-STARRS lens, Berghea C.T. et al.,
*                           2017ApJ...844...90B
*  BGGFC (4,4)           : COSMOS hi-z, Boutsia K. et al., 2018ApJ...869...20B
*  BQLS (17,17)          : BOSS QSO lenses & pairs, More A. et al.,
*                           2016MNRAS.456.1595M
*  BZCAT (5,4)           : Blazars catalog, Massaro E. et al.,
*                           http://www.asdc.asi.it/bzcat
*  C-COSM (180,180)      : Chandra COSMOS IDs, Marchesi S. et al.,
*                           2016ApJ...817...34M
*  ChaMP (191,187)       : Trichas M. et al., 2012ApJS..200...17T
*  DABAST (1,1)          : Diaz-Santos T. et al., 2018Sci...362.1034D
*  Dart (26,25)          : Heavily Obscured QSOs, Hviding R. et al.,
*                           2018MNRAS.474.1955H
*  DDC2 (14,14)          : variable AGN, De Cicco D. et al., 2019,
*                           arXiv:1905.10374
*  DEEP (143,139)        : DEEP2, Newman J. et al., 2013ApJS..208....5N;
*                           deep.ps.uci.edu/DR4
*  DES (1,1)             : Dark Energy hi-z, Reed S.L. et al.,
*                           2015MNRAS.454.3952R
*  DESQQ (26,26)         : STRIDES lenses, Anguita T. et al.,
*                           2018MNRAS.480.5017A
*  DESQQ2 (12,12)        : STRIDES lenses, Treu T. et al., 2018MNRAS.481.1041T
*  DPeake (683,656)      : Double-peaked NELGs, Ge J.-Q. et al.,
*                           2012, Cat. J/ApJS/201/31
*  DR12 (18,23)          : Alam S. et al., 2015ApJS..219...12A,
*                           http://sdss.org/dr12
*  DR12Q (129,153)       : SDSS-DR12Q, Paris I. et al., 2017A&A...597A..79P,
*                           Cat. VII/279
*  DR14 (42582,45125)    : Abolfathi B. et al., 2018ApJS..235...42A pipeline,
*                           data at
*                           https://data.sdss.org/sas/dr14/sdss/spectro/redux
*  DR14Q (513025,525752) : SDSS-DR14Q, Paris I. et al., 2018A&A...613A..51P
*                           data at http://data.sdss.org/sas/dr14/eboss/qso/DR14Q
*  DR7 (8,145)           : SDSS DR7, Abazajian K.N. et al., 2009ApJS..182..543A,
*                           files at
*                    http://classic.sdss.org/dr7/products/spectra/getspectra.html
*  DR7Q (2089,329)       : SDSS Quasar DR7, Schneider D. et al.,
*                           2010AJ....139.2360S, Cat. VII/26 data
*                http://classic.sdss.org/dr7/products/value_added/qsocat_dr7.html
*  DUHIZ (2,2)           : DECaLS-UKIRT hi-z, Wang F. et al., 2017ApJ...839...27W
*  Dusty (11,11)         : Dusty Starbursts, Rodighiero G. et al.,
*                           2019ApJ...877...38H
*  DUz6 (18,18)          : DESI & UKIRT hi-z, Wang F. et al.,
*                           2018,arXiv:1810.11926
*  eHAQ (82,81)          : Extended High AV, Krogager J.-K., 2016ApJ...832...49K
*  ELQ-PS (216,215)      : ELQS on PS1, Schindler J.-T. et al., 2019,
*                           arXiv:1905.04069
*  ELQS-N (38,38)        : ELQS in NGC, Schindler J.-T. et al.,
*                           2018ApJ...863..144S
*  ELQS-S (126,126)      : ELQS in SGC, Schindler J.-T. et al.,
*                           2019ApJ...871..258S
*  FISCBA (1,1)          : HST lens, Fischer/Schade/Barrientos
*                           1998ApJ...503L.127F
*  GAIA2 (263,263)       : Gaia DR2, Gaia Collaboration et al.,
*                           2018A&A...616A...1G, Cat. I/345 data as presented by
*                          SIMBAD, http://simbad.u-strasbg.fr/simbad
*  GEIER (1,1)           : Geier S.J. et al., 2019,arXiv:1904.01686
*  GLDD (1,1)            : Lensed QSO data-driven, Ostrovski F. et al.,
*                           2017MNRAS.465.4325O
*  GLIKMA (29,29)        : red WISE QSOs, Glikman E. et al., 2018ApJ...861...37G
*  GLRED (1,1)           : lensed red QSO, Glikman E. et al., 2018,
*                           arXiv:1807.05434
*  GUTI (1,1)            : not a ULX, Gutierrez C.M., 2013A&A...549A..81G
*  H-DOGs (16,16)        : Herschel DOGs, Riguccini L.A. et al.,
*                           2019A&A...625A...9D
*  HAQ (2,2)             : High AV serendipitous, Heintz K.E. et al.,
*                           2016AJ....152...13H
*  HAQC (1,1)            : High AV in COSMOS, Heintz K.E. et al.,
*                           2016A&A...595A..13H
*  HEINTZ (1,1)          : dusty absorbed QSO, Heintz K.E. et al.,
*                           2018A&A...615A..43H
*  HIZ7.5 (1,1)          : QSO z=7.5, Banados E. et al., 2018Natur.553..473B
*  HSTvar (42,42)        : variable AGN, Pouliasis E. et al., 2019,
*                           arXiv:1905.11995
*  IGMCP (10,10)         : IGM close pairs, Rorai A. et al., 2017Sci...356..418R
*  IKEDA (1,1)           : Ikeda H. et al., 2017ApJ...846...57I
*  IMS (1,1)             : IR medium-deep hi-z, Kim Y. et al.,
*                           2015ApJ...813...35H
*  IMS2 (10,11)          : IR medium-deep hi-z, Kim Y. et al.,
*                           2019ApJ...870...86K
*  KOVACS (1,1)          : bright QSO behind Milky Way, Kovacs T. et al.,
*                           2019, RNAAS, 3, 3
*  LAMDR4 (190,2249)     : LAMOST-DR4, pipeline, http://dr4.lamost.org
*  LAMQ1 (687,628)       : LAMOST QUASAR DR1, Ai Y.L. et al.,
*                           2016, Cat. J/AJ/151/24
*  LAMQ3 (6773,6663)     : LAMOST QUASAR DR3/DR2, Dong X.Y. et al.,
*                           2018, Cat. J/AJ/155/189
*  LAMQ5 (7998,7997)     : LAMOST QUASAR DR5/DR4, Yao S. et al.,
*                           2019ApJS..240....6Y
*  LEMON (40,40)         : 24 Gaia lenses, Lemon C. et al.,
*                           2018MNRAS.479.5060L
*  LEMON2 (30,30)        : 22 Gaia lenses, Lemon/Auger/McMahon,
*                           2019MNRAS.483.4242L
*  LGGS (11,11)          : M31/M33 area, Massey/Neugent/Levesque, 2019,
*                           arXiv:1904.07898
*  LIDMAN (1,1)          : SN Host Galaxy redshifts, Lidman C. et al.,
*                           2013PASA...30....1L
*  LIN (1,1)             : DES lens, Lin H. et al., 2017ApJ...838...15L
*  LIRAS (169,154)       : LoCuSS IR AGNs, Xu, L. et al.,
*                           2015, Cat. J/ApJS/219/18
*  LOZAGN (32,11009)     : Low-redshift AGN, Liu H.-Y. et al., 2019,
*                           arXiv:1906.05597
*  LSSA (2,2)            : 2 lenses, Lucey/Schechter/Smith/Anguita,
*                          2018MNRAS.476..927L
*  LUMIz5 (66,66)        : Luminous hi-z, Yang J. et al., 2018,
*                           arXiv:1810.11927
*  M31UV (1,1)           : UV flare QSO on M31, Meusinger H. et al.,
*                           2010A&A...512A...1M
*  MALS-N (70,70)        : MEERKAT QSOs, Krogager J.-K. et al.,
*                           2018ApJS..235...10K
*  MFJC (52,51)          : McGreer I.D., Fan X., Jiang L. & Cai Z.,
*                           2018, Cat. J/AJ/155/131
*  MORX (49324,0)        : Million Radio/X-ray Associations, Flesch E.,
*                           2016PASA...33...52F, Cat. V/148
*  MQ (6639,622652)      : MILLIQUAS, original data in this catalog,
*                           Flesch E., 2019
*  NBCKDE (40821,40930)  : Richards G.T. et al., 2009, Cat. J/ApJS/180/67
*  NBCKv3 (529572,563852): NBCKDE v3, Richards G.T. et al.,
*                           2015, Cat. J/ApJS/219/39
*  NED (3,3)             : NASA/IPAC Extragalactic Database,
*                           https://ned.ipac.caltech.edu
*  OGLE2 (2,2)           : OGLE quasars, Kozlowski S. et al., 2018,
*                           arXiv:1810.08622
*  OSTROV (1,1)          : lensed, Ostrovski et al., 2019, in preparation
*  OVRLAP (5,5)          : SDSS overlap hi-z QSOs, Jiang L. et al.,
*                           2015AJ....149..188J
*  OzDES (664,638)       : Dark Energy SN QSOs, Tie S.S. et al.,
*                           2017AJ....153..107T
*  P352-1 (1,1)          : P352-15, Banados E. et al., 2018ApJ...861L..14B
*  PETERS (12896,12914)  : photo special, Peters C.M. et al.,
*                           2015ApJ...811...95P
*  PGC (13623,19)        : Principal Galaxy Catalog, Paturel G. et al.,
*                           2003A&A...412...45P, Cat. VII/237
*  PHILLI (1,0)          : MERLIN lens, Phillips P.M. et al.,
*                           2000MNRAS.319L...7P
*  PS1 (63,63)           : PAN-STARRS1 hi-z, Banados E. et al.,
*                           2016,  Cat. J/ApJS/227/11
*  PS1hiz (1,1)          : Tang, Ji-Jia et al., 2017MNRAS.466.4568T
*  PS1MAZ (6,6)          : Mazzucchelli C. et al., 2017ApJ...849...91M
*  PSO (3,3)             : PAN-STARRS z-dropouts, Venemans B.P. et al.,
*                           2015ApJ...801L..11V
*  QLSV (20,19)          : QUEST-La Silla, Sanchez-Saez P. et al.,
*                           2019ApJS..242...10S
*  QPQ10 (70,70)         : Quasar pair DB, Findlay J.R. et al.,
*                           2018ApJS..236...44F
*  RBS (3,3)             : Laporte N. et al., 2017ApJ...851...40L
*  Redden (25,25)        : Reddened QSOs, Temple M.J. et al., 2019MNRAS.487.2594R
*  REQ4 (6,6)            : Reionization-Era quasars, Yang J. et al.,
*                           2019AJ....157..236Y
*  RLQ (4,3)             : Tuccillo D./Gonzalez-Serrano J.I./Benn C.R.,
*                           2015, Cat. J/MNRAS/449/2818
*  RSG (1,1)             : Dorn-Wallenstein T.Z. & Levesque E.,
*                           2017, IAUS, 329, 376
*  S82X (72,77)          : Stripe 82 AGN, LaMassa S.M. et al.,
*                           2019ApJ...876...50L
*  S82XRQ (8,8)          : Red Quasars, LaMassa S.M. et al., 2017ApJ...847..100L
*  SCULPT (2,3)          : Sculptor X-ray, Arnason R. M. et al.,
*                           2019MNRAS.485.2259A
*  SDLENS (3,3)          : SDSS Lenses, Williams P.R. et al.,
*                           2018MNRAS.477L..70W
*  SDSSHI (6,6)          : SDSS hi-z, Jiang L. et al., 2016ApJ...833..222J
*  SFM201 (1,1)          : Schulze S. et al., 2012A&A...546A..20S
*  SHELLQ (33,33)        : Subaru hi-z, Matsuoka Y. et al., 2018PASJ...70S..35M
*  SHELQS (30,30)        : SHELLQS hi-z, Matsuoka Y. et al., 2018ApJS..237....5M
*  SHELz7 (1,1)          : SHELLQS z=7, Matsuoka Y. et al., 2019ApJ...872L...2M
*  SMSSQ (2,2)           : SkyMapper hi-z, Zefeng Li Z. et al., 2018,
*                           arXiv:1805.03429
*  SPIN18 (1,1)          : KiDS-SQuaD lens, Spiniello C. et al.,
*                           2018MNRAS.480.1163S
*  SPIN19 (2,2)          : 2 lensed quasars, Spiniello C. et al.,
*                           2019MNRAS.485.5086S
*  SQLS (60,52)          : SDSS DR7 QSO Lens Search, Inada N. et al.,
*                           Cat. J/AJ/143/119
*  SQUAD (13,13)         : UVES DB DR1, Murphy M.T. et al., 2019MNRAS.482.3458M
*  SSLENS (3,3)          : South sky lenses, Spiniello C. et al.,
*                           2019MNRAS.483.3888S
*  SUV (21,21)           : SDSS-ULAS/VHS QSOs, Yang J. et al.,
*                           2017AJ....153..184Y
*  SXDF (39,39)          : Subaru-XMMDF redshifts, Simpson C. et al.,
*                           2012, Cat. J/MNRAS/421/3060
*  SXDS (308,307)        : Subaru-XMMDF spectra, Akiyama M. et al.,
*                           2015PASJ...67...82A
*  ULTRA (1,1)           : Ultraluminous hi-z, Wu, X.-B. et al.,
*                            2015Natur.518..512W
*  UVQS (435,503)        : UV QSOs, Monroe T.R. et al., 2016, Cat. J/AJ/152/25
*  VAHIZ (2,2)           : VST ATLAS hi-z, Carnall A.C. et al.,
*                           2015MNRAS.451L..16C
*  VAHIZ2 (1,1)          : bright z>6 QSOs, Chehade B. et al.,
*                           2018MNRAS.478.1649C
*  VAHIZ3 (1,1)          : VST-ATLAS lens, Schechter P.L. et al.,
*                           2018, RNAAS, 2b, 21
*  VAQL (11,11)          : VST-ATLAS quasar systems,Schechter P.L. et al.,
*                           2017AJ....153..219S
*  VDES (8,8)            : VISTA Dark Energy QSOs, Reed S.L. et al.,
*                           2017MNRAS.468.4702R
*  VDES2 (2,2)           : more VHS-DES quasars, Reed S.L. et al.,
*                           2019, arXiv:1901.07456
*  VIKING (4,4)          : VIKING IR, Venemans, G.A. et al., 2015MNRAS.453.2259V
*  VIPERS (241,284)      : VIPERS PDR-2, Scodeggio M. et al., 2018A&A...609A..84S
*  VMC (34,34)           : Magellanic IR QSOs, Ivanov V.D. et al.,
*                           2016A&A...588A..93I
*  WARSAW (3,3)          : OGLE lens, Kostrzewa-Rutkowska Z. et al.,
*                           2018MNRAS.476..663K
*  WERTZ (1,1)           : Gaia GraL lens, Wertz O. et al.,
*                           2018, arXiv:1810.02624
*  WGD (2,2)             : DES/Gaia lenses; Agnello A. et al.,
*                           2018MNRAS.479.4345A
*  WISEA (447228,0)      : AllWISE QSO candidates, Secrest N. et al.,
*                           2015, Cat. J/ApJS/221/12
*  WISEHI (72,70)        : Hi-z QSOs from WISE, Wang F. et al.,
*                           2016, Cat. J/ApJ/819/24
*  WOLF1 (1,1)           : most ultraluminous QSO, Wolf C. et al.,
*                           2018PASA...35...24W
*  XDQSO (235761,0)      : SDSS-XDQSO, Bovy J. et al., 2011ApJ...729..141B
*  XLSS (306,118)        : Stalin C.S. et al., 2010, Cat. J/MNRAS/401/294
*  XMM2 (12,12)          : 2XMM-Newton cross-search, Combi J.A. et al.,
*                           2011Ap&SS.331...53C, Cat. V/138
*  XMMSMC (6,6)          : SMC quasars, Maitral C. et al., 2019A&A...622A..29M
*  XMSS (182,148)        : Barcons X. et al., 2007, Cat. J/A+A/476/1191
*  XWAS (490,449)        : Esquej P. et al., 2013, Cat. J/A+A/557/A123
*  YQLF (25,25)          : deep CFHT QSOs, Yang J. et al.,
*                           2018, Cat. J/AJ/155/110
*  z6.51 (1,1)           : lensed quasar z=6.51, Fan X. et al.,
*                           2019ApJ...870L..11F
*  z7.02 (1,1)           : lensed quasar z=7.02, Wang F. et al.,
*                           2018ApJ...869L...9W
*
*  4-digit numeric citations are indexed in the HMQ (2015,PASA,32,10)
*  references list.  The citation for the classification (e.g., that the
*  object is a quasar) can be from either the name or redshift citation.
*Note (9):  For a QSO candidate (type starting with q/R/X/2), this shows the
*  percent chance that it is a QSO, based on photometric and/or radio/X-ray
*  association analysis. Included candidates are those of 80%+ likelihood
*  of being true quasars. For a known QSO (type=Q/A), Bl Lac type (type=B),
*  type-II (type=K/N) or active galaxy (type=G), this shows the percent
*  chance that the shown radio/X-ray detection(s) is truly associated to it.
*
*  Candidates (objects without spectroscopic confirmation) originate from
*  three types of sources:
*
*  (1) Photometric quasars from the SDSS-based NBCKDE/NBCKDE-v3/XDQSO/Peters
*       quasar candidate catalogs.  Those catalogs give pQSOs (calculated QSO
*       probabilities) for their objects which however are not reported here;
*       instead, I have calibrated those pQSOs against SDSS-DR12Q classified
*       objects to produce these QSO likelihoods -- see appendix 1 of my HMQ
*       paper (Flesch 2015PASA...32...10F) for details of the calibration method.
*
*  (2) WISE-sourced candidates (from Secrest et al. 2015ApJS..221...12S,
*       Cat. J/ApJS/221/12) are presented by their authors with a bulk pQSO as
*       an average figure, with no individual pQSOs given. These I have
*       generated by matching WISE candidates to optical objects (within a two
*       arcsecond radius) and then binning them into four-color (colors B-R,
*       R-W1, W1-W2, and W2-W3) subsets which are then calibrated against
*       SDSS-DR12Q classified objects to yield the pQSO for each four-color
*       subset.
*
*  (3) 56182 radio/X-ray associated optical objects are presented only here in
*       Milliquas, with pQSO calculated as described in the QORG paper (Flesch
*       & Hardcastle, 2004A&A...427..387F, Cat. J/A+A/427/387), with no redshift
*       displayed. However, an additional 349 such objects show pipeline
*       redshifts from SDSS-DR14.
*
*  71543 photometric quasars are also radio/X-ray associated, and the
*  displayed probability figure combines the calibrated photometric QSO
*  likelihood P1 and the radio/X-ray derived QSO likelihood P2 as
*  P = 1/(1+((1-P1)*(1-P2))/(P1*P2)).  Over all objects, using the probability
*  as expected yield, the 1945385 type-I Milliquas objects will yield 1890994
*  actual type-I quasars/AGN.
*Note (10): Four columns of Radio/X-ray detections are presented:
*   * 1st column: best X-ray detection (i.e. highest probability association).
*   * 2nd column: best core Radio detection.
*   * 3rd column: a radio lobe if the description (see note 2) shows a "2",
*                 otherwise this is an additional radio or X-ray detection.
*   * 4th column: a radio lobe if the description (see note 2) shows a "2",
*                 otherwise this is an additional X-ray detection.
*
*   Legend of Radio/X-ray detection prefixes and catalog home pages:
*
*   FIRST: VLA FIRST survey, 13Jun05 version, http://sundog.stsci.edu
*   NVSS: NRAO VLA sky survey, http://www.cv.nrao.edu/nvss
*   SUMSS: Sydney U. Molonglo, http://www.physics.usyd.edu.au/sifa/Main/SUMSS
*   MGPS: Molonglo galactic plane survey, same attribution as SUMSS
*   ROSAT catalogs home page: www.mpe.mpg.de/xray/wave/rosat/catalogue , for:
*    - 1RXH: ROSAT HRI (high resolution imager)
*    - 2RXP/2RXF: ROSAT PSPC (position sensitive proportional counter)
*    - 1RXS: ROSAT RASS (all-sky survey, both bright & faint)
*   2RXS: 2nd RASS source catalog, Boller Th. et al., 2016,A&A,588,103
*   1WGA: White, Giommi & Angelini, wgacat.gsfc.nasa.gov/wgacat/wgacat.html
*   CXO: Chandra Source Catalog v1.1, http://cxc.cfa.harvard.edu/csc
*   CXOG: Chandra ACIS source catalog, Wang S. et al., 2016,ApJS,224,40
*   CXOX: XAssist Chandra source list, http://xassist.pha.jhu.edu/
*   2XMM/2XMMi: XMM-Newton DR3, http://xmm.esac.esa.int/xsa/versions.shtml
*   3XMM: XMM-Newton DR8, http://www.cosmos.esa.int/web/xmm-newton/xsa
*   XMMSL: XMM-Newton Slew Survey Release 2.0, same attribution as 3XMM
*   XMMX: XAssist XMM-Newton source list, http://xassist.pha.jhu.edu/
*   1SXPS: Swift X-ray Point Source catalog, http://www.swift.ac.uk/1SXPS
*
*   Optical field solutions are calculated from the raw source positions of all
*   these catalogs as described in my MORX paper 2016PASA...33...52F.

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

C  Declarations for 'refs.dat'	! References (for numerical references)

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

      integer*4     Ref        (nr__1) ! Reference number
      integer*4     N1         (nr__1) ! ? Number of this reference for Name
      integer*4     N2         (nr__1) ! ? Number of this reference for redshift
      character*19  Bibcode    (nr__1) ! BibCode
      character*31  Aut        (nr__1) ! Author's name
      character*90  Com        (nr__1) ! Comments

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

C  Loading file 'catalog.dat'	! The Milliquas catalogue, V6.3 (16 June 2019)

C  Format for file interpretation

    1 format(
     +  F11.7,1X,F11.7,2X,A25,1X,A4,1X,F5.2,1X,F5.2,1X,A3,1X,A1,1X,A1,
     +  1X,F6.3,1X,A6,1X,A6,1X,I3,1X,A22,1X,A22,1X,A22,1X,A22)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,1986800
        read(1,'(A192)')ar__
        read(ar__,1)
     +  RAdeg(i__),DEdeg(i__),Name(i__),Cl_Ass(i__),Rmag(i__),
     +  Bmag(i__),Comment(i__),R(i__),B(i__),z(i__),rName(i__),
     +  rz(i__),Qpct(i__),XName(i__),RName_1(i__),Lobe1(i__),
     +  Lobe2(i__)
        if(ar__(77:82) .EQ. '') z(i__) = rNULL__
        if(ar__(98:100) .EQ. '') Qpct(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  RAdeg(i__),DEdeg(i__),Name(i__),Cl_Ass(i__),Rmag(i__),
     +  Bmag(i__),Comment(i__),R(i__),B(i__),z(i__),rName(i__),
     +  rz(i__),Qpct(i__),XName(i__),RName_1(i__),Lobe1(i__),
     +  Lobe2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References (for numerical references)

C  Format for file interpretation

    2 format(I4,1X,I6,1X,I6,2X,A19,1X,A31,1X,A90)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,2171
        read(1,'(A162)')ar__1
        read(ar__1,2)
     +  Ref(i__),N1(i__),N2(i__),Bibcode(i__),Aut(i__),Com(i__)
        if(ar__1(6:11) .EQ. '') N1(i__) = iNULL__
        if(ar__1(13:18) .EQ. '') N2(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  Ref(i__),N1(i__),N2(i__),Bibcode(i__),Aut(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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