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

Conversion of standardized ReadMe file for file /./ftp/cats/VII/166 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-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. VII/166  Quasars and Active Galactic Nuclei (6th Ed.)     (Veron+ 1993)
*================================================================================
*Quasars and Active Galactic Nuclei (6th Ed.)
*     Veron-Cetty M.P., Veron P.
*     <ESO Scientific Report 13 (1993)>
*     =1993ESOSR..13....1V            (SIMBAD/NED Reference)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Quasars (defined as
                                    brighter than absolute B magnitude -23)

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

C  1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension 1950
      real*8        DEdeg       ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      character*1   FLAG        ! [*] A flag which is blank if the object has
*                                      been detected as a radio source and which
*                                      is '*' if not.
      character*17  Name        ! Most common name of the object.
      integer*4     RAh         ! (h) Right Ascension 1950 (hours)
      character*1   n_RAh       ! Note on the position:
*                                    'A' means low-accuracy position
*                                    'O' means optical position
*                                    'R' means radio position measured
*                                    both 'O' and 'R" positions measured with
*                                    an accuracy better than 1 arcsec
      integer*4     RAm         ! (min) RA1950   (minutes)
      real*4        RAs         ! (s) RA1950   (seconds)
      character*1   DE_         ! Declination 1950 (Sign)
      integer*4     DEd         ! (deg) DE1950   (degrees)
      integer*4     DEm         ! (arcmin) DE1950   (')
      integer*4     DEs         ! (arcsec) DE1950   (")
      character*1   n_Vmag      ! [/R] Note on the "V" magnitude
*                                    '/' indicates photographic magnitude
*                                        and colours when the colours are given.
*                                    'R' indicates a red magnitude
      real*4        Vmag        ! (mag) ? V magnitude, or photographic or red
*                                        magnitude depending on n_Vmag
      real*4        B_V         ! (mag) ? Colour index
      real*4        U_B         ! (mag) ? Colour index
      character*1   l_z         ! [/>] Note on the redshift:
*                                    '/' means that the redshift was estimated
*                                        from a low dispersion slitless spectrum.
*                                    '>' indicates that z is a lower limit
      real*4        z           ! ? Value of the redshift
      real*4        F6cm        ! (Jy) ? 6cm flux (Jy)
      integer*4     r_F6cm      ! [1,1486]?=0 Reference of 6cm Flux
      real*4        F11cm       ! (Jy) ? 11cm flux (Jy)
      integer*4     r_F11cm     ! [1,1486]?=0 Reference of 11cm flux
      character*1   Xray        ! [X] 'X' when the object was detected
*                                      in the X-ray range
      character*2   Class       ! [S123?HPBL ] Classification of the object:
*                                    S1: Seyfert 1 spectrum
*                                    S2: Seyfert 2 spectrum
*                                    S3: Seyfert 3 or liner
*                                    S?: possibly a Seyfert
*                                    H2: nuclear HII region
*                                    HP: high optical polarization (>3%)
*                                    BL: confirmed BL Lac object
      integer*4     FC          ! [1,1486]? Reference of the Finding Chart
      integer*4     r_Vmag      ! [1,1486]? Reference of the Photometry
      integer*4     r_z         ! [1,1486]? Reference of the redshift

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

C  Declarations for 'table2.dat'	! BL Lac objects

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

C  1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_1     ! (deg) Right Ascension 1950
      real*8        DEdeg_1     ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      character*1   FLAG_1      ! [*] A flag which is blank if the object has
*                                      been detected as a radio source and which
*                                      is '*' if not.
      character*17  Name_1      ! Most common name of the object.
      integer*4     RAh_1       ! (h) Right Ascension 1950 (hours)
      character*1   n_RAh_1     ! Note on the position:
*                                    'A' means low-accuracy position
*                                    'O' means optical position
*                                    'R' means radio position measured
*                                    both 'O' and 'R" positions measured with
*                                    an accuracy better than 1 arcsec
      integer*4     RAm_1       ! (min) RA1950   (minutes)
      real*4        RAs_1       ! (s) RA1950   (seconds)
      character*1   DE__1       ! Declination 1950 (Sign)
      integer*4     DEd_1       ! (deg) DE1950   (degrees)
      integer*4     DEm_1       ! (arcmin) DE1950   (')
      integer*4     DEs_1       ! (arcsec) DE1950   (")
      character*1   n_Vmag_1    ! [/R] Note on the "V" magnitude
*                                    '/' indicates photographic magnitude
*                                        and colours when the colours are given.
*                                    'R' indicates a red magnitude
      real*4        Vmag_1      ! (mag) ? V magnitude, or photographic or red
*                                        magnitude depending on n_Vmag
      real*4        B_V_1       ! (mag) ? Colour index
      real*4        U_B_1       ! (mag) ? Colour index
      character*1   l_z_1       ! [/>] Note on the redshift:
*                                    '/' means that the redshift was estimated
*                                        from a low dispersion slitless spectrum.
*                                    '>' indicates that z is a lower limit
      real*4        z_1         ! ? Value of the redshift
      real*4        F6cm_1      ! (Jy) ? 6cm flux (Jy)
      integer*4     r_F6cm_1    ! [1,1486]?=0 Reference of 6cm Flux
      real*4        F11cm_1     ! (Jy) ? 11cm flux (Jy)
      integer*4     r_F11cm_1   ! [1,1486]?=0 Reference of 11cm flux
      character*1   Xray_1      ! [X] 'X' when the object was detected
*                                      in the X-ray range
      character*2   Class_1     ! [S123?HPBL ] Classification of the object:
*                                    S1: Seyfert 1 spectrum
*                                    S2: Seyfert 2 spectrum
*                                    S3: Seyfert 3 or liner
*                                    S?: possibly a Seyfert
*                                    H2: nuclear HII region
*                                    HP: high optical polarization (>3%)
*                                    BL: confirmed BL Lac object
      integer*4     FC_1        ! [1,1486]? Reference of the Finding Chart
      integer*4     r_Vmag_1    ! [1,1486]? Reference of the Photometry
      integer*4     r_z_1       ! [1,1486]? Reference of the redshift

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

C  Declarations for 'table3.dat'	! Active Galaxy Nuclei (defined as
                                    fainter than  absolute B magnitude -23)

      integer*4 nr__2
      parameter (nr__2=1855)	! Number of records
      character*100 ar__2  	! Full-size record

C  1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_2     ! (deg) Right Ascension 1950
      real*8        DEdeg_2     ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      character*1   FLAG_2      ! [*] A flag which is blank if the object has
*                                      been detected as a radio source and which
*                                      is '*' if not.
      character*17  Name_2      ! Most common name of the object.
      integer*4     RAh_2       ! (h) Right Ascension 1950 (hours)
      character*1   n_RAh_2     ! Note on the position:
*                                    'A' means low-accuracy position
*                                    'O' means optical position
*                                    'R' means radio position measured
*                                    both 'O' and 'R" positions measured with
*                                    an accuracy better than 1 arcsec
      integer*4     RAm_2       ! (min) RA1950   (minutes)
      real*4        RAs_2       ! (s) RA1950   (seconds)
      character*1   DE__2       ! Declination 1950 (Sign)
      integer*4     DEd_2       ! (deg) DE1950   (degrees)
      integer*4     DEm_2       ! (arcmin) DE1950   (')
      integer*4     DEs_2       ! (arcsec) DE1950   (")
      character*1   n_Vmag_2    ! [/R] Note on the "V" magnitude
*                                    '/' indicates photographic magnitude
*                                        and colours when the colours are given.
*                                    'R' indicates a red magnitude
      real*4        Vmag_2      ! (mag) ? V magnitude, or photographic or red
*                                        magnitude depending on n_Vmag
      real*4        B_V_2       ! (mag) ? Colour index
      real*4        U_B_2       ! (mag) ? Colour index
      character*1   l_z_2       ! [/>] Note on the redshift:
*                                    '/' means that the redshift was estimated
*                                        from a low dispersion slitless spectrum.
*                                    '>' indicates that z is a lower limit
      real*4        z_2         ! ? Value of the redshift
      real*4        F6cm_2      ! (Jy) ? 6cm flux (Jy)
      integer*4     r_F6cm_2    ! [1,1486]?=0 Reference of 6cm Flux
      real*4        F11cm_2     ! (Jy) ? 11cm flux (Jy)
      integer*4     r_F11cm_2   ! [1,1486]?=0 Reference of 11cm flux
      character*1   Xray_2      ! [X] 'X' when the object was detected
*                                      in the X-ray range
      character*2   Class_2     ! [S123?HPBL ] Classification of the object:
*                                    S1: Seyfert 1 spectrum
*                                    S2: Seyfert 2 spectrum
*                                    S3: Seyfert 3 or liner
*                                    S?: possibly a Seyfert
*                                    H2: nuclear HII region
*                                    HP: high optical polarization (>3%)
*                                    BL: confirmed BL Lac object
      integer*4     FC_2        ! [1,1486]? Reference of the Finding Chart
      integer*4     r_Vmag_2    ! [1,1486]? Reference of the Photometry
      integer*4     r_z_2       ! [1,1486]? Reference of the redshift

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

C  Declarations for 'table4.dat'	! Objects with a published redshift, but
                                         later shown not to be a quasar.

      integer*4 nr__3
      parameter (nr__3=52)	! Number of records
      character*86 ar__3  	! Full-size record

C  1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg_3     ! (deg) Right Ascension 1950
      real*8        DEdeg_3     ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      character*1   FLAG_3      ! [*] Flag similar to table1
      character*11  Name_3      ! Most common name of the object.
      character*6   v4C         ! Its number in 4C radio source catalogue
      integer*4     RAh_3       ! (h) Right Ascension 1950 (hours)
      character*1   n_RAh_3     ! Note on the position, as in table1
      integer*4     RAm_3       ! (min) RA1950   (minutes)
      real*4        RAs_3       ! (s) RA1950   (seconds)
      character*1   DE__3       ! Declination 1950 (Sign)
      integer*4     DEd_3       ! (deg) DE1950   (degrees)
      integer*4     DEm_3       ! (arcmin) DE1950   (')
      integer*4     DEs_3       ! (arcsec) DE1950   (")
      character*1   n_Vmag_3    ! [/R] Note on Vmag, as in table1
      real*4        Vmag_3      ! (mag) ? V magnitude, as in table1
      real*4        B_V_3       ! (mag) ? Colour index
      real*4        U_B_3       ! (mag) ? Colour index
      character*1   l_z_3       ! [/>] Note on the redshift, as in table1
      real*4        z_3         ! ? Value of the redshift
      character*2   Class0      ! Original classification, as in table1
      character*2   Class1      ! [*G? ] New classification, with symbols
*                                  'G' = galaxy
*                                  '*' = Star
      integer*4     FC_3        ! [1,1486]? Reference of the Finding Chart
      integer*4     r_Vmag_3    ! [1,1486]? Reference of the Photometry
      integer*4     r_z_3       ! [1,1486]? Reference of the redshift
      integer*4     r_Class1    ! [1,1486]? Reference of New classification

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

C  Declarations for 'refs.dat'	! References

      integer*4 nr__4
      parameter (nr__4=1834)	! Number of records
      character*81 ar__4  	! Full-size record

      integer*4     Ref         ! Reference number
      character*76  Text        ! Text of the reference. A blank in column 6
*                                       indicates a continuation line

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

C  Loading file 'table1.dat'	! Quasars (defined as
*                                    brighter than absolute B magnitude -23)

C  Format for file interpretation

    1 format(
     +  A1,A17,1X,I2,A1,I2,1X,F4.1,A1,I2,1X,I2,1X,I2,A1,F5.2,1X,F5.2,
     +  F5.2,A1,F5.3,F6.3,I4,F5.2,I4,3X,A1,A2,I4,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,7437
        read(1,'(A100)')ar__
        read(ar__,1)
     +  FLAG,Name,RAh,n_RAh,RAm,RAs,DE_,DEd,DEm,DEs,n_Vmag,Vmag,B_V,
     +  U_B,l_z,z,F6cm,r_F6cm,F11cm,r_F11cm,Xray,Class,FC,r_Vmag,r_z
        if(ar__(40:44) .EQ. '') Vmag = rNULL__
        if(ar__(46:50) .EQ. '') B_V = rNULL__
        if(ar__(51:55) .EQ. '') U_B = rNULL__
        if(ar__(57:61) .EQ. '') z = rNULL__
        if(ar__(62:67) .EQ. '') F6cm = rNULL__
        if(ar__(72:76) .EQ. '') F11cm = rNULL__
        if(ar__(87:90) .EQ. '') FC = iNULL__
        if(ar__(92:95) .EQ. '') r_Vmag = iNULL__
        if(ar__(97:100) .EQ. '') r_z = iNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  FLAG,Name,RAh,n_RAh,RAm,RAs,DE_,DEd,DEm,DEs,n_Vmag,Vmag,B_V,
     +  U_B,l_z,z,F6cm,r_F6cm,F11cm,r_F11cm,Xray,Class,FC,r_Vmag,r_z
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! BL Lac objects

C  Format for file interpretation

    2 format(
     +  A1,A17,1X,I2,A1,I2,1X,F4.1,A1,I2,1X,I2,1X,I2,A1,F5.2,1X,F5.2,
     +  F5.2,A1,F5.3,F6.3,I4,F5.2,I4,3X,A1,A2,I4,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,171
        read(1,'(A100)')ar__1
        read(ar__1,2)
     +  FLAG_1,Name_1,RAh_1,n_RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,
     +  DEs_1,n_Vmag_1,Vmag_1,B_V_1,U_B_1,l_z_1,z_1,F6cm_1,r_F6cm_1,
     +  F11cm_1,r_F11cm_1,Xray_1,Class_1,FC_1,r_Vmag_1,r_z_1
        if(ar__1(40:44) .EQ. '') Vmag_1 = rNULL__
        if(ar__1(46:50) .EQ. '') B_V_1 = rNULL__
        if(ar__1(51:55) .EQ. '') U_B_1 = rNULL__
        if(ar__1(57:61) .EQ. '') z_1 = rNULL__
        if(ar__1(62:67) .EQ. '') F6cm_1 = rNULL__
        if(ar__1(72:76) .EQ. '') F11cm_1 = rNULL__
        if(ar__1(87:90) .EQ. '') FC_1 = iNULL__
        if(ar__1(92:95) .EQ. '') r_Vmag_1 = iNULL__
        if(ar__1(97:100) .EQ. '') r_z_1 = iNULL__
        RAdeg_1 = rNULL__
        DEdeg_1 = rNULL__
c  Derive coordinates RAdeg_1 and DEdeg_1 from input data
c  (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
        if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
        if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
        if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
        if(DEd_1 .GE. 0) DEdeg_1=DEd_1
        if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
        if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
        if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c    ..............Just test output...........
        write(6,2)
     +  FLAG_1,Name_1,RAh_1,n_RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,
     +  DEs_1,n_Vmag_1,Vmag_1,B_V_1,U_B_1,l_z_1,z_1,F6cm_1,r_F6cm_1,
     +  F11cm_1,r_F11cm_1,Xray_1,Class_1,FC_1,r_Vmag_1,r_z_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Active Galaxy Nuclei (defined as
*                                    fainter than  absolute B magnitude -23)

C  Format for file interpretation

    3 format(
     +  A1,A17,1X,I2,A1,I2,1X,F4.1,A1,I2,1X,I2,1X,I2,A1,F5.2,1X,F5.2,
     +  F5.2,A1,F5.3,F6.3,I4,F5.2,I4,3X,A1,A2,I4,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,1855
        read(1,'(A100)')ar__2
        read(ar__2,3)
     +  FLAG_2,Name_2,RAh_2,n_RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,
     +  DEs_2,n_Vmag_2,Vmag_2,B_V_2,U_B_2,l_z_2,z_2,F6cm_2,r_F6cm_2,
     +  F11cm_2,r_F11cm_2,Xray_2,Class_2,FC_2,r_Vmag_2,r_z_2
        if(ar__2(40:44) .EQ. '') Vmag_2 = rNULL__
        if(ar__2(46:50) .EQ. '') B_V_2 = rNULL__
        if(ar__2(51:55) .EQ. '') U_B_2 = rNULL__
        if(ar__2(57:61) .EQ. '') z_2 = rNULL__
        if(ar__2(62:67) .EQ. '') F6cm_2 = rNULL__
        if(ar__2(72:76) .EQ. '') F11cm_2 = rNULL__
        if(ar__2(87:90) .EQ. '') FC_2 = iNULL__
        if(ar__2(92:95) .EQ. '') r_Vmag_2 = iNULL__
        if(ar__2(97:100) .EQ. '') r_z_2 = iNULL__
        RAdeg_2 = rNULL__
        DEdeg_2 = rNULL__
c  Derive coordinates RAdeg_2 and DEdeg_2 from input data
c  (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
        if(RAh_2 .GT. -180) RAdeg_2=RAh_2*15.
        if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4.
        if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240.
        if(DEd_2 .GE. 0) DEdeg_2=DEd_2
        if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60.
        if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600.
        if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c    ..............Just test output...........
        write(6,3)
     +  FLAG_2,Name_2,RAh_2,n_RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,
     +  DEs_2,n_Vmag_2,Vmag_2,B_V_2,U_B_2,l_z_2,z_2,F6cm_2,r_F6cm_2,
     +  F11cm_2,r_F11cm_2,Xray_2,Class_2,FC_2,r_Vmag_2,r_z_2
        write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Objects with a published redshift, but
*                                         later shown not to be a quasar.

C  Format for file interpretation

    4 format(
     +  A1,A11,A6,1X,I2,A1,I2,1X,F4.1,A1,I2,1X,I2,1X,I2,A1,F5.2,1X,
     +  F5.2,F5.2,A1,F5.3,A2,1X,A2,1X,I4,1X,I4,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,52
        read(1,'(A86)')ar__3
        read(ar__3,4)
     +  FLAG_3,Name_3,v4C,RAh_3,n_RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,
     +  DEs_3,n_Vmag_3,Vmag_3,B_V_3,U_B_3,l_z_3,z_3,Class0,Class1,
     +  FC_3,r_Vmag_3,r_z_3,r_Class1
        if(ar__3(40:44) .EQ. '') Vmag_3 = rNULL__
        if(ar__3(46:50) .EQ. '') B_V_3 = rNULL__
        if(ar__3(51:55) .EQ. '') U_B_3 = rNULL__
        if(ar__3(57:61) .EQ. '') z_3 = rNULL__
        if(ar__3(68:71) .EQ. '') FC_3 = iNULL__
        if(ar__3(73:76) .EQ. '') r_Vmag_3 = iNULL__
        if(ar__3(78:81) .EQ. '') r_z_3 = iNULL__
        if(ar__3(83:86) .EQ. '') r_Class1 = iNULL__
        RAdeg_3 = rNULL__
        DEdeg_3 = rNULL__
c  Derive coordinates RAdeg_3 and DEdeg_3 from input data
c  (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
        if(RAh_3 .GT. -180) RAdeg_3=RAh_3*15.
        if(RAm_3 .GT. -180) RAdeg_3=RAdeg_3+RAm_3/4.
        if(RAs_3 .GT. -180) RAdeg_3=RAdeg_3+RAs_3/240.
        if(DEd_3 .GE. 0) DEdeg_3=DEd_3
        if(DEm_3 .GE. 0) DEdeg_3=DEdeg_3+DEm_3/60.
        if(DEs_3 .GE. 0) DEdeg_3=DEdeg_3+DEs_3/3600.
        if(DE__3.EQ.'-'.AND.DEdeg_3.GE.0) DEdeg_3=-DEdeg_3
c    ..............Just test output...........
        write(6,4)
     +  FLAG_3,Name_3,v4C,RAh_3,n_RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,
     +  DEs_3,n_Vmag_3,Vmag_3,B_V_3,U_B_3,l_z_3,z_3,Class0,Class1,
     +  FC_3,r_Vmag_3,r_z_3,r_Class1
        write(6,'(6H Pos: 2F8.4)') RAdeg_3,DEdeg_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    5 format(I4,1X,A76)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,1834
        read(1,'(A81)')ar__4
        read(ar__4,5)Ref,Text
c    ..............Just test output...........
        write(6,5)Ref,Text
c    .......End.of.Just test output...........
      end do
      close(1)

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