Conversion of standardized ReadMe file for
file /./home/cats/J/A_A/707/A30 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-17
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/707/A30 Qasar pairs with Gaia astrometric data (Chen+, 2026)
*================================================================================
*Search for quasar pairs with Gaia astrometric data. I. Method and candidates.
* Chen Q., Jing L., Zhu X., Fang Y., He Z., Deng Z., Xiang C., Wu J.
* <Astron. Astrophys. 707, A30 (2026)>
* =2026A&A...707A..30C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mgqpc.dat' ! Quasar Pair Candidates derived from the
cross-match of Million Quasars catalogue (MQC)
and Gaia
integer*4 nr__
parameter (nr__=4112) ! Number of records
character*370 ar__ ! Full-size record
C Position composed of: RAAdeg DEAdeg (Epoch=2000)
character*23 System (nr__) ! Pair name designated by the center coordinates
* at J2000 (ICRS), JHHMMSSsss+DDMMSSsss
* (SystemName)
character*25 NameA (nr__) ! The name of pair member A retrieved from the
* Million Quasars catalogue (NAME)
character*4 TypeA (nr__) ! The object type of pair member A retrieved
* from the MQCv8 (TYPE)
real*8 zA (nr__) ! Redshift of pair member A, updated with the
* MQCv8 and SIMBAD (Z_A)
real*8 seplim (nr__) ! (arcsec) Cross-matching radius calculated with Z_A and
* 100 kpc (ec seplim)
integer*8 GaiaDR3A (nr__) ! Unique Gaia DR3 identifier (source_id) of
* pair member A (DR3Name_A)
real*8 RAAdeg (nr__) ! (deg) Gaia (ICRS) at Ep=2000 right ascension of
* pair member A (RA_A)
real*8 DEAdeg (nr__) ! (deg) Gaia (ICRS) at Ep=2000 declination of pair
* member A (DEC_A)
real*8 GmagA (nr__) ! Gaia G-band magnitude pair member A (Gmag_A)
real*8 BPmagA (nr__) ! ? Gaia BP-band magnitude pair member A
* (BPmag_A)
real*8 RPmagA (nr__) ! ? Gaia RP-band magnitude pair member A
* (RPmag_A)
real*4 f0plxA (nr__) ! ? Probability density of zero parallax
* calculated for pair member A (f_0plx_A)
real*4 f0pmA (nr__) ! ? Probability density of zero proper motion
* calculated for pair member A (f_0pm_A)
integer*8 GaiaDR3B (nr__) ! Unique Gaia DR3 identifier (source_id) of
* pair member B (DR3Name_B)
real*8 RABdeg (nr__) ! (deg) Gaia (ICRS) at Ep=2000 right ascension of
* pair member B (RA_B)
real*8 DEBdeg (nr__) ! (deg) Gaia (ICRS) at Ep=2000 declination of
* pair member B (DEC_B)
real*8 GmagB (nr__) ! Gaia G-band magnitude pair member B (Gmag_B)
real*8 BPmagB (nr__) ! ? Gaia BP-band magnitude pair member B
* (BPmag_B)
real*8 RPmagB (nr__) ! ? Gaia RP-band magnitude pair member B
* (RPmag_B)
real*4 f0plxB (nr__) ! Probability density of zero parallax
* calculated for pair member B (f_0plx_B)
real*8 f0pmB (nr__) ! Probability density of zero proper motion
* calculated for pair member B (f_0pm_B)
real*8 sepAB (nr__) ! (arcsec) Angular separation between the two members in
* the pair (sep_AB)
real*8 GLONB (nr__) ! (deg) Galactic longitude of pair member B (l_B)
real*8 GLATB (nr__) ! (deg) Galactic latitude of pair member B (b_B)
character*6 SkyLab (nr__) ! Label of the sky region where pair member B
* is located (SkyLab) (1)
character*6 MTypeB (nr__) ! Morphological type of pair member B retrieved
* from DESI-LS DR10 (MorphoType_B)
character*5 VIClass (nr__) ! Pair classification assigned by visual
* Inspection (VIClass) (2)
*Note (1): e.g., "Gb2030" represents the sky region of 20 deg <= |b| < 30 deg.
*Note (2): Pair classification as follows:
* QPC = quasar pair candidate
* LQC = lensed quasar candidate
* WSLQC = wide-separation lensed quasar candidate
C=============================================================================
C Loading file 'mgqpc.dat' ! Quasar Pair Candidates derived from the
* cross-match of Million Quasars catalogue (MQC)
* and Gaia
C Format for file interpretation
1 format(
+ A23,1X,A25,1X,A4,1X,F8.6,1X,F11.8,10X,I19,1X,F13.9,1X,F13.9,
+ 1X,F9.6,1X,F10.6,1X,F9.6,1X,E11.9,1X,E22.16,10X,I19,1X,F13.9,
+ 1X,F13.9,1X,F9.6,1X,F10.6,1X,F10.6,1X,E12.9,1X,F12.9,1X,F9.6,
+ 1X,F13.9,1X,F12.8,1X,A6,1X,A6,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'mgqpc.dat')
write(6,*) '....Loading file: mgqpc.dat'
do i__=1,4112
read(1,'(A370)')ar__
read(ar__,1)
+ System(i__),NameA(i__),TypeA(i__),zA(i__),seplim(i__),
+ GaiaDR3A(i__),RAAdeg(i__),DEAdeg(i__),GmagA(i__),BPmagA(i__),
+ RPmagA(i__),f0plxA(i__),f0pmA(i__),GaiaDR3B(i__),RABdeg(i__),
+ DEBdeg(i__),GmagB(i__),BPmagB(i__),RPmagB(i__),f0plxB(i__),
+ f0pmB(i__),sepAB(i__),GLONB(i__),GLATB(i__),SkyLab(i__),
+ MTypeB(i__),VIClass(i__)
if(ar__(144:153) .EQ. '') BPmagA(i__) = rNULL__
if(ar__(155:163) .EQ. '') RPmagA(i__) = rNULL__
if(ar__(165:175) .EQ. '') f0plxA(i__) = rNULL__
if(ar__(177:198) .EQ. '') f0pmA(i__) = rNULL__
if(ar__(267:276) .EQ. '') BPmagB(i__) = rNULL__
if(ar__(278:287) .EQ. '') RPmagB(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ System(i__),NameA(i__),TypeA(i__),zA(i__),seplim(i__),
+ GaiaDR3A(i__),RAAdeg(i__),DEAdeg(i__),GmagA(i__),BPmagA(i__),
+ RPmagA(i__),f0plxA(i__),f0pmA(i__),GaiaDR3B(i__),RABdeg(i__),
+ DEBdeg(i__),GmagB(i__),BPmagB(i__),RPmagB(i__),f0plxB(i__),
+ f0pmB(i__),sepAB(i__),GLONB(i__),GLATB(i__),SkyLab(i__),
+ MTypeB(i__),VIClass(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end