Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/684/A202 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-14
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/684/A202 Optical emission with radio structures (Lambert+, 2024)
*================================================================================
*Identification of the optical emission detected by Gaia with radio structures in
*parsec-scale AGN jets.
* Lambert S., Sol H., Pierron A.
* <Astron. Astrophys. 684, A202 (2024)>
* =2024A&A...684A.202L (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table.dat' ! Numerical quantities that allowed to built up
the analyses and conclusions of the article
integer*4 nr__
parameter (nr__=422) ! Number of records
character*198 ar__ ! Full-size record
character*8 IERS (nr__) ! IERS B1950 source name (HHMM+DDd)
real*8 DV (nr__) ! (mas) Gaia-ICRF3 vector length
real*4 E_DV (nr__) ! (mas) Gaia-ICRF3 vector length error
integer*4 thetaV (nr__) ! (deg) Gaia-ICRF3 vector position-angle
integer*4 E_thetaV (nr__) ! (deg) Gaia-ICRF3 vector position-angle error
real*8 DG (nr__) ! (mas) Gaia-MOJAVE core vector length
real*4 E_DG (nr__) ! (mas) Gaia-MOJAVE core vector length error
integer*4 thetaG (nr__) ! (deg) Gaia-MOJAVE core vector position-angle
integer*4 E_thetaG (nr__) ! (deg) Gaia-MOJAVE core vector position-angle error
integer*4 MOJCOMP (nr__) ! [0/57]?=- Identifier of MOJAVE component
* closest to Gaia (MCCG) following the
* numbering of Lister et al.
* (2021ApJ...923...30L, Cat. J/ApJ/923/30)
real*4 DM (nr__) ! (mas) MCCG-MOJAVE core vector length
character*3 thetaM (nr__) ! (deg) MCCG-MOJAVE core position-angle
real*8 DGM (nr__) ! (mas) MCCG-Gaia vector length
real*4 E_DGM (nr__) ! (mas) MCCG-Gaia vector length error
integer*4 thetaGM (nr__) ! (deg) MCCG-Gaia vector position-angle
integer*4 E_thetaGM (nr__) ! (deg) MCCG-Gaia vector position-angle error
real*8 NormDG (nr__) ! Gaia-MOJAVE core normalized separation
real*8 NormDGM (nr__) ! MCCG-Gaia normalized separation
real*4 VG (nr__) ! (mas/yr) ?=- Gaia proper motion
real*4 E_VG (nr__) ! (mas/yr) ?=- Gaia proper motion error
integer*4 phiG (nr__) ! (deg) ?=- Gaia proper motion position-angle
integer*4 E_phiG (nr__) ! (deg) ?=- Gaia proper motion position-angle error
real*4 VM (nr__) ! (mas/yr) ?=- MCCG proper motion
real*4 E_VM (nr__) ! (mas/yr) ?=- MCCG proper motion error
integer*4 phiM (nr__) ! (deg) ?=- MCCG proper motion position-angle
integer*4 E_phiM (nr__) ! (deg) ?=- MCCG proper motion position-angle error
real*8 MOJFLUX (nr__) ! (Jy) ?=- MCCG 15 GHz flux
integer*4 Type (nr__) ! [0/7]?=- Activity Type (NED Homogenized
* Classification) (1)
real*4 z (nr__) ! ?=- Redshift from NED
real*4 Gmag (nr__) ! (mag) Gaia G magnitude
real*4 OPol (nr__) ! ?=- Optical polarization from
* Hutsemekers et al. (2018A&A...620A..68H,
* Cat. J/A+A/620/A68)
real*4 RPol (nr__) ! ?=- Maximum of fractional radio polarization
* from Lister et al. (2018ApJS..234...12L,
* Cat. J/ApJS/234/12)
real*4 BP_RP (nr__) ! (mag) ?=- Gaia BP-RP colour index
real*4 RUWE (nr__) ! ?=- Gaia RUWE
real*4 VLBIStab (nr__) ! (mas) ?=- VLBI stability from geodetic time series
integer*4 Cat (nr__) ! [1/4] Category (1: C, 2: B, 3: J, 4: O)
*Note (1): Activity Type as follows:
* 0 = FSRQ
* 1 = QSO
* 2 = BL Lac
* 3 = Seyfert 1
* 4 = Seyfert 2
* 5 = GPS
* 6 = CSS
* 7 = other
C=============================================================================
C Loading file 'table.dat' ! Numerical quantities that allowed to built up
* the analyses and conclusions of the article
C Format for file interpretation
1 format(
+ A8,1X,F7.3,1X,F6.3,1X,I3,1X,I3,1X,F7.3,1X,F6.3,1X,I3,1X,I3,1X,
+ I2,1X,F6.3,1X,A3,1X,F7.3,1X,F6.3,1X,I3,1X,I4,1X,F7.2,1X,F7.2,
+ 1X,F5.3,1X,F5.3,1X,I3,1X,I3,1X,F5.3,1X,F5.3,1X,I3,1X,I3,1X,
+ F7.2,1X,I1,1X,F4.1,1X,F4.1,1X,F5.2,1X,F5.2,1X,F5.3,1X,F3.1,1X,
+ F5.3,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table.dat')
write(6,*) '....Loading file: table.dat'
do i__=1,422
read(1,'(A198)')ar__
read(ar__,1)
+ IERS(i__),DV(i__),E_DV(i__),thetaV(i__),E_thetaV(i__),DG(i__),
+ E_DG(i__),thetaG(i__),E_thetaG(i__),MOJCOMP(i__),DM(i__),
+ thetaM(i__),DGM(i__),E_DGM(i__),thetaGM(i__),E_thetaGM(i__),
+ NormDG(i__),NormDGM(i__),VG(i__),E_VG(i__),phiG(i__),
+ E_phiG(i__),VM(i__),E_VM(i__),phiM(i__),E_phiM(i__),
+ MOJFLUX(i__),Type(i__),z(i__),Gmag(i__),OPol(i__),RPol(i__),
+ BP_RP(i__),RUWE(i__),VLBIStab(i__),Cat(i__)
if (MOJCOMP(i__) .EQ. 45) MOJCOMP(i__) = iNULL__
if (idig(ar__(110:114)).EQ.0) VG(i__) = rNULL__
if (idig(ar__(116:120)).EQ.0) E_VG(i__) = rNULL__
if (phiG(i__) .EQ. 45) phiG(i__) = iNULL__
if (E_phiG(i__) .EQ. 45) E_phiG(i__) = iNULL__
if (idig(ar__(130:134)).EQ.0) VM(i__) = rNULL__
if (idig(ar__(136:140)).EQ.0) E_VM(i__) = rNULL__
if (phiM(i__) .EQ. 45) phiM(i__) = iNULL__
if (E_phiM(i__) .EQ. 45) E_phiM(i__) = iNULL__
if (idig(ar__(150:156)).EQ.0) MOJFLUX(i__) = rNULL__
if (Type(i__) .EQ. 45) Type(i__) = iNULL__
if (idig(ar__(160:163)).EQ.0) z(i__) = rNULL__
if (idig(ar__(170:174)).EQ.0) OPol(i__) = rNULL__
if (idig(ar__(176:180)).EQ.0) RPol(i__) = rNULL__
if (idig(ar__(182:186)).EQ.0) BP_RP(i__) = rNULL__
if (idig(ar__(188:190)).EQ.0) RUWE(i__) = rNULL__
if (idig(ar__(192:196)).EQ.0) VLBIStab(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ IERS(i__),DV(i__),E_DV(i__),thetaV(i__),E_thetaV(i__),DG(i__),
+ E_DG(i__),thetaG(i__),E_thetaG(i__),MOJCOMP(i__),DM(i__),
+ thetaM(i__),DGM(i__),E_DGM(i__),thetaGM(i__),E_thetaGM(i__),
+ NormDG(i__),NormDGM(i__),VG(i__),E_VG(i__),phiG(i__),
+ E_phiG(i__),VM(i__),E_VM(i__),phiM(i__),E_phiM(i__),
+ MOJFLUX(i__),Type(i__),z(i__),Gmag(i__),OPol(i__),RPol(i__),
+ BP_RP(i__),RUWE(i__),VLBIStab(i__),Cat(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end