FORTRAN Generation
(/./ftp/cats/J/A_A/684/A202)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/684/A202 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 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        ! IERS B1950 source name (HHMM+DDd)
      real*8        DV          ! (mas) Gaia-ICRF3 vector length
      real*4        E_DV        ! (mas) Gaia-ICRF3 vector length error
      integer*4     thetaV      ! (deg) Gaia-ICRF3 vector position-angle
      integer*4     E_thetaV    ! (deg) Gaia-ICRF3 vector position-angle error
      real*8        DG          ! (mas) Gaia-MOJAVE core vector length
      real*4        E_DG        ! (mas) Gaia-MOJAVE core vector length error
      integer*4     thetaG      ! (deg) Gaia-MOJAVE core vector position-angle
      integer*4     E_thetaG    ! (deg) Gaia-MOJAVE core vector position-angle error
      integer*4     MOJCOMP     ! [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          ! (mas) MCCG-MOJAVE core vector length
      character*3   thetaM      ! (deg) MCCG-MOJAVE core position-angle
      real*8        DGM         ! (mas) MCCG-Gaia vector length
      real*4        E_DGM       ! (mas) MCCG-Gaia vector length error
      integer*4     thetaGM     ! (deg) MCCG-Gaia vector position-angle
      integer*4     E_thetaGM   ! (deg) MCCG-Gaia vector position-angle error
      real*8        NormDG      ! Gaia-MOJAVE core normalized separation
      real*8        NormDGM     ! MCCG-Gaia normalized separation
      real*4        VG          ! (mas/yr) ?=- Gaia proper motion
      real*4        E_VG        ! (mas/yr) ?=- Gaia proper motion error
      integer*4     phiG        ! (deg) ?=- Gaia proper motion position-angle
      integer*4     E_phiG      ! (deg) ?=- Gaia proper motion position-angle error
      real*4        VM          ! (mas/yr) ?=- MCCG proper motion
      real*4        E_VM        ! (mas/yr) ?=- MCCG proper motion error
      integer*4     phiM        ! (deg) ?=- MCCG proper motion position-angle
      integer*4     E_phiM      ! (deg) ?=- MCCG proper motion position-angle error
      real*8        MOJFLUX     ! (Jy) ?=- MCCG 15 GHz flux
      integer*4     Type        ! [0/7]?=- Activity Type (NED Homogenized
*                                    Classification) (1)
      real*4        z           ! ?=- Redshift from NED
      real*4        Gmag        ! (mag) Gaia G magnitude
      real*4        OPol        ! ?=- Optical polarization from
*                                    Hutsemekers et al. (2018A&A...620A..68H,
*                                    Cat. J/A+A/620/A68)
      real*4        RPol        ! ?=- Maximum of fractional radio polarization
*                                    from Lister et al. (2018ApJS..234...12L,
*                                    Cat. J/ApJS/234/12)
      real*4        BP_RP       ! (mag) ?=- Gaia BP-RP colour index
      real*4        RUWE        ! ?=- Gaia RUWE
      real*4        VLBIStab    ! (mas) ?=- VLBI stability from geodetic time series
      integer*4     Cat         ! [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,DV,E_DV,thetaV,E_thetaV,DG,E_DG,thetaG,E_thetaG,MOJCOMP,
     +  DM,thetaM,DGM,E_DGM,thetaGM,E_thetaGM,NormDG,NormDGM,VG,E_VG,
     +  phiG,E_phiG,VM,E_VM,phiM,E_phiM,MOJFLUX,Type,z,Gmag,OPol,RPol,
     +  BP_RP,RUWE,VLBIStab,Cat
        if (MOJCOMP .EQ. 45) MOJCOMP =  iNULL__
        if (idig(ar__(110:114)).EQ.0) VG =  rNULL__
        if (idig(ar__(116:120)).EQ.0) E_VG =  rNULL__
        if (phiG .EQ. 45) phiG =  iNULL__
        if (E_phiG .EQ. 45) E_phiG =  iNULL__
        if (idig(ar__(130:134)).EQ.0) VM =  rNULL__
        if (idig(ar__(136:140)).EQ.0) E_VM =  rNULL__
        if (phiM .EQ. 45) phiM =  iNULL__
        if (E_phiM .EQ. 45) E_phiM =  iNULL__
        if (idig(ar__(150:156)).EQ.0) MOJFLUX =  rNULL__
        if (Type .EQ. 45) Type =  iNULL__
        if (idig(ar__(160:163)).EQ.0) z =  rNULL__
        if (idig(ar__(170:174)).EQ.0) OPol =  rNULL__
        if (idig(ar__(176:180)).EQ.0) RPol =  rNULL__
        if (idig(ar__(182:186)).EQ.0) BP_RP =  rNULL__
        if (idig(ar__(188:190)).EQ.0) RUWE =  rNULL__
        if (idig(ar__(192:196)).EQ.0) VLBIStab =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  IERS,DV,E_DV,thetaV,E_thetaV,DG,E_DG,thetaG,E_thetaG,MOJCOMP,
     +  DM,thetaM,DGM,E_DGM,thetaGM,E_thetaGM,NormDG,NormDGM,VG,E_VG,
     +  phiG,E_phiG,VM,E_VM,phiM,E_phiM,MOJFLUX,Type,z,Gmag,OPol,RPol,
     +  BP_RP,RUWE,VLBIStab,Cat
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