Conversion of standardized ReadMe file for
file /./ftp/cats/VII/280 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. VII/280 The Million Quasars (Milliquas) catalog (V5.2) (Flesch, 2017) *================================================================================ *The Million Quasars (Milliquas) catalogue, version 5.2. * Flesch E.W. * <Pub. Astron. Soc. Australia 32, 10 (2015)> * =2015PASA...32...10F * =2017yCat.7280....0F C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'catalog.dat' ! The Milliquas catalogue, V5.2 (5 August 2017) integer*4 nr__ parameter (nr__=1998464) ! Number of records character*192 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 RAdeg ! (deg) Right ascension (J2000) (1) real*8 DEdeg ! (deg) Declination (J2000) (1) character*25 Name ! ID from the literature, or J2000 (2) character*1 Cl ! [ABKQqN] classification of object (3) character*3 Assoc ! [2RX ] Summary of associations for object (4) real*4 Rmag ! (mag) Red optical magnitude (5) real*4 Bmag ! (mag) Blue optical magnitude (5) character*3 Comment ! Comment on optical object (6) character*1 R ! Red optical PSF class (7) character*1 B ! Blue optical PSF class (7) real*4 z ! ? Redshift from the literature or estimated (8) character*6 r_Name ! Citation for name (9) character*6 r_z ! Citation for redshift (9) integer*4 Qpct ! (pct) ? Probability that this object is a QSO (10) character*22 XName ! X-ray ID, if any (11) character*22 RName ! Radio ID, if any (11) character*22 Lobe1 ! Radio lobe ID or extra R/X ID, if any (11) character*22 Lobe2 ! Radio lobe ID or extra X-ray ID, if any (11) *Note (1): These are to 7 decimals which is too precise, but which accomodates a * miniscule round-up which prevents inadvertent truncation by the user when * converting to sexagesimal. *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): Classification of object as follows: * Q = QSO, type I broad-line unresolved, 577146 of these. * A = AGN, type I Seyferts/host-dominated, 30062 of these. * B = BL Lac object, 1615 of these. * q = photometric quasar candidate from SDSS or WISE, 1297111 of these. * N = type-II narrow-line Seyfert/NELG galaxy, 38110 of these. This class * includes an unknown number of LINERs, presented for completeness. * K = type-II narrow-line NLAGN, quasar-like high luminosity, 1355 of these. *Note (4): Summary of associations for object as follows: * R = Radio association displayed. * X = X-ray association displayed. * 2 = Double radio lobes displayed (declared by data-driven algorithm). *Note (5): 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) can be * estimates if both bands are integer or one band empty. Note: many SDSS * magnitudes are "extinction-corrected" ~0.3 mag brighter than observed. *Note (6): Comment on optical object legend as follows: * p = optical magnitudes are POSS-I O (violet 4050{AA}) and E (red 6400{AA}). * These are preferred because O is well-offset from E, and these plates * were always taken on the same night, thus the red-blue color is * correct even for variable objects. * j = blue magnitude is SERC J (Bj 4850A blue-green) from the POSS-II or * UKST surveys. Red-blue color is less reliable because the plates were * taken in different epochs, i.e. years apart. * g = blue magnitude is SDSS green 4900A, red is SDSS r 6200A. * u = blue magnitude is SDSS ultraviolet 3850A. * b = blue magnitude is Vega 4400A (Johnson). * v = red magnitude is visual., ie, white, 5500A midpoint. * i = red magnitude is infrared 7500A. * z = red magnitude is infrared z 8500A. * (not v/i/z) = standard red 6400A (Cousins). * + = variability nominally detected in both red/blue, over multi-epoch data. * m = proper motion nominally detected, from USNO-B. * a = object is an SDSS galaxy with AGN subclass. If also BROADLINE then * Milliquas class is 'A', otherwise 'N' (see note 3). * ? = identification uncertain (quasar may be located elsewhere, 4 of these). *Note (7): 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 color (fainter than plate depth, or confused, etc.) *Note (8): Photometric redshifts (for objects classified as 'q') are here rounded * to 0.1 z. The XDQSO catalog does not provide photometric redshifts, so * those are provided either by NBCKDE or by this catalog using the method * detailed in Appendix 2 of my HMQ paper (2015PASA...32...10F). 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%+ likely to * be true within 0.5z of the displayed redshift value -- usual likelihood is * more like 80%. *Note (9): Legend (with counts of name and redshift) and references: * 2dF (366,232) : 2dF galaxy survey, Colless M. et al., * 2001MNRAS.328.1039C, Cat. VII/250 * 2QZ (27515,24161) : Croom S.M. et al., 2004MNRAS.349.1397C, Cat. VII/241 * 2SLAQ (10388,8717) : Croom S.M. et al., 2009, Cat. J/MNRAS/392/19 * 3LAC (30,28) : FERMI AGN v3, Ackermann M. et al., * 2015, Cat. J/ApJ/810/14 * 6dF (287,227) : 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 (2047,2047) : AGES survey, Kochanek C.S. et al., * 2012, Cat. J/ApJS/200/8 * ATLAS (229,269) : Mao M.Y. et al., 2012, Cat. J/MNRAS/426/3334 * BASS (15,124) : Swift-BAT AGN, Koss M. et al., * 2017,arXiv:1707.08123, ApJ * BQLS (9,9) : BOSS QSO lenses & pairs, More A. et al., * 2016MNRAS.456.1595M * C-COSM (180,180) : Chandra COSMOS IDs, Marchesi S. et al., * 2016ApJ...817...34M * ChaMP (191,188) : Trichas M. et al., 2012ApJS..200...17T * DES (1,1) : Dark Energy hi-z, Reed S.L. et al., * 2015MNRAS.454.3952R * DPeake (7,3450) : Double-peaked NELGs, Ge J.-Q. et al., * 2012, Cat. J/ApJS/201/31 * DR12 (1628,5332) : Alam S. et al., 2015ApJS..219...12A, * http://sdss.org/dr12, Cat. V/147 * DR12Q (233,415) : SDSS-DR12Q, Paris I. et al., 2017A&A...597A..79P, * Cat. VII/279 * DR14 (32900,40157) : Abolfathi B. et al., 2017,arXiv:1707.09322, data at * https://data.sdss.org/sas/dr14/sdss/spectro/redux * DR14Q (512359,525951) : SDSS-DR14Q, Paris I. et al., 2017, in preparation, * data at http://data.sdss.org/sas/dr14/eboss/qso/DR14Q * DR7 (6378,4872) : SDSS DR7, Abazajian K.N. et al., 2009ApJS..182..543A, * files at * http://classic.sdss.org/dr7/products/spectra/getspectra.html * DR7Q (2128,340) : SDSS Quasar DR7, Schneider D. et al., * 2010AJ....139.2360S, Cat. VII/260 * DUHIZ (2,2) : DECaLS-UKIRT hi-z, Wang F. et al., 2017ApJ...839...27W * eHAQ (82,81) : Extended High AV, Krogager J.-K., 2016ApJ...832...49K * GLDD (1,1) : Lensed QSO data-driven, Ostrovski F. et al., * 2017MNRAS.465.4325O * 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 * HMQ (28638,22897) : Half Million Quasars catalog, Flesch E., * 2015PASA...32...10F, Cat. VII/273, references * therein to the original authors & papers. * IGMCP (10,10) : IGM close pairs, Rorai A. et al., 2017Sci...356..418R * IMS (1, 1) : IR medium-deep hi-z, Kim Y. et al., * 2015ApJ...813...35H * IKEDA (1,1) : Ikeda H. et al., 2017, arXiv:1708.00314, ApJ accepted * LAMOST (4314,3786) : LAMOST-DR3, http://dr3.lamost.org * LIDMAN (1,1) : SN Host Galaxy redshifts, Lidman C. et al., * 2013PASA...30....1L * LIRAS (168,158) : LoCuSS IR AGNs, Xu, L. et al., * 2015, Cat. J/ApJS/219/18 * MORX (49916,0) : Million Radio/X-ray Associations, Flesch E., * 2016PASA...33...52F, Cat. V/148 * MQ (3144,635633) : MILLIQUAS, original data in this catalog, * Flesch E.,2017 * MZZ (9,9) : Marano B., Zamorani G., Zitelli V., * 1988MNRAS.232..111M * NBCKDE (40679,40679) : Richards G.T. et al., 2009, Cat. J/ApJS/180/67 * NBCKv3 (576782,576782): NBCKDE v3, Richards G.T. et al., * 2015, Cat. J/ApJS/219/39 * 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 * PETERS (12952,12952) : photo special, Peters C.M. et al., * 2015ApJ...811...95P * PGC (13548,8) : Principal Galaxy Catalogue, Paturel G. et al., * 2003A&A...412...45P, Cat. VII/237 * PS1 (63,63) : PAN-STARRS1 hi-z, Banados E. et al., * Cat. J/ApJS/227/11 * PSO (3,3) : PAN-STARRS z-dropouts, Venemans B.P. et al., * 2015ApJ...801L..11V * 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,arXiv:1701.07888 * 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., * 2017,arXiv:1704.05854, PASJ * SUV (21,21) : SDSS-ULAS/VHS QSOs, Yang J. et al., 2017,AJ,153,184 * SXDF (39,39) : Subaru-XMMDF redshifts, Simpson C. et al., * 2012, Cat. J/MNRAS/421/3060 * SXDS (310,309) : 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 (433,517) : UV QSOs, Monroe T.R. et al., 2016,AJ,152,25 * VAHIZ (2,2) : VST ATLAS hi-z, Carnall A.C. et al., * 2015MNRAS.451L..16C * VDES (8,8) : VISTA Dark Energy Qsos, Reed S.L. et al., * 2017MNRAS.468.4702R * VIKING (4,4) : VIKING IR, Venemans, G.A. et al., 2015MNRAS.453.2259V * VIPERS (241,285) : VIPERS PDR-2, Scodeggio M. et al., * 2016, arXiv:1611.07048 * VMC (34,34) : Magellanic IR QSOs, Ivanov V.D. et al., * 2016A&A...588A..93I * WARREN (1,1) : Warren S.J. et al., 2017, in preparation * WISEA (430999,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 * XDQSO (235694,0) : SDSS-XDQSO, Bovy J. et al., 2011ApJ...729..141B * XLSS (307,120) : Stalin C.S. et al., 2010, Cat. J/MNRAS/401/294 * XMSS (182,150) : 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 * * The citation for the classification (e.g., that the object is a quasar) can * be from either the name or redshift citation. *Note (10): 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) * or type-II (type=K/N), this shows the percent chance that the shown * radio/X-ray detection(s) is truly associated to it. * * Candidates (objects without spectroscopic confirmation) consist of three * separate classes of objects: * * (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 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 dividing them into four-color (the colors * being B-R, R-W1, W1-W2 & W2-W3) subsets which are then calibrated * against SDSS-DR12Q classified objects to yield the pQSO for each subset. * * (3) Radio/X-ray associated objects, totalling 53065 without any other * attribution. The likelihoods of being a true QSO are calculated as * described in the QORG paper (Flesch & Hardcastle, 2004A&A...427..387F, * Cat. J/A+A/427/387). * * 69609 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 1957384 type-I Milliquas objects will yield 1901076 * actual type-I quasars/AGNs, so making this a true million quasars catalog. *Note (11): 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 3) shows a "2", * otherwise this is an additional radio or X-ray detection. * * 4th column: a radio lobe if the description (see note 3) 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/zope/xassist * 2XMM/2XMMi: XMM-Newton DR3, http://xmm.esac.esa.int/xsa/versions.shtml * 3XMM: XMM-Newton DR7, 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, xassist.pha.jhu.edu/zope/xassist * 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 Loading file 'catalog.dat' ! The Milliquas catalogue, V5.2 (5 August 2017) C Format for file interpretation 1 format( + F11.7,1X,F11.7,2X,A25,1X,A1,A3,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,1998464 read(1,'(A192)')ar__ read(ar__,1) + RAdeg,DEdeg,Name,Cl,Assoc,Rmag,Bmag,Comment,R,B,z,r_Name,r_z, + Qpct,XName,RName,Lobe1,Lobe2 if(ar__(77:82) .EQ. '') z = rNULL__ if(ar__(98:100) .EQ. '') Qpct = iNULL__ c ..............Just test output........... write(6,1) + RAdeg,DEdeg,Name,Cl,Assoc,Rmag,Bmag,Comment,R,B,z,r_Name,r_z, + Qpct,XName,RName,Lobe1,Lobe2 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end