FORTRAN Generation
(/./ftp/cats/J/A_A/363/141)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/363/141 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-29
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/363/141       Radio-loud galaxies in the ROSAT survey. IV. (Reich+, 2000)
*================================================================================
*Radio-loud active galaxies in the northern ROSAT All-Sky Survey.
*IV. Properties of sources stronger than 100 mJy at 5 GHz.
*    Reich W., Fuerst E., Reich P., Kothes R., Brinkmann W., Siebert J.
*   <Astron. Astrophys. 363, 141 (2000)>
*   =2000A&A...363..141R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Data for 741 source

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

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

      character*12  RGB         ! Green Bank catalog (<VIII/52>) name
      character*1   n_RGB       ! [*] *: source with repeated observations
      character*9   Bname       ! Name based on 1950 position (HHMM+DDMM)
      character*5   r_Pos       ! Position references (2)
      integer*4     F11cm       ! (mJy) ? Integrated 11cm (2.7GHz) flux density
      integer*4     F6cm        ! (mJy) ? Integrated 6cm (5GHz) flux density
      integer*4     F2_8cm      ! (mJy) ? Integrated 2.8cm (10.7GHz) flux density
      real*4        SI11_6      ! ? Spectral index 11cm/6cm (1)
      real*4        e_SI11_6    ! ? Error of spectral index 11cm/6cm
      real*4        F11_2_8     ! (mW/m2) Integrated flux density (11cm to 2.8cm)
      integer*4     q_F11cm     ! [1/4]? Error class see Note (3)
      integer*4     q_F6cm      ! [1/4]? Error class see Note (3)
      integer*4     q_F2_8cm    ! [1/4]? Error class see Note (3)
      real*4        SI6_2_8     ! ? Spectral index 6cm/2.8cm (1)
      real*4        e_SI6_2_8   ! ? Error of spectral index 6cm/2.8cm
      real*4        Fx          ! (mW/m2) Integrated flux density (0.1-2.4 KeV)
      real*4        z           ! ? Measured redshift
      real*4        Vmap        ! (mag) ? Optical magnitude (V at most)
      integer*4     D11a        ! (arcsec) ? RA source size 11cm
      character*1   n_D11a      ! [u] u: Unresolved (point-like sources)
      integer*4     D11d        ! (arcsec) ? DE source size 11cm
      character*1   n_D11d      ! [u] u: Unresolved (point-like sources)
      integer*4     D6a         ! (arcsec) ? RA source size 6cm
      character*1   n_D6a       ! [u] u: Unresolved (point-like sources)
      integer*4     D6d         ! (arcsec) ? DE source size 6cm
      character*1   n_D6d       ! [u] u: Unresolved (point-like sources)
      integer*4     D2_8a       ! (arcsec) ? RA source size 2.8cm
      character*1   n_D2_8a     ! [u] u: Unresolved (point-like sources)
      integer*4     D2_8d       ! (arcsec) ? DE source size 2.8cm
      character*1   n_D2_8d     ! [u] u: Unresolved (point-like sources)
      real*4        SI11_       ! ? Spectral index fit 11cm/6cm/2.8cm (1)
      real*4        e_SI11_     ! ? Error of spectral index fit
      real*4        SIx         ! ? X-ray power law photon index (1)
      real*4        e_SIx       ! ? Error of x-ray power law photon index
      integer*4     e_Fx        ! (%) Error of X-ray flux density
      character*3   OptID       ! Optical identification see Note (4)
      character*10  Name        ! Common source name
      character*1   l_Pol11     ! Limit flag on Pol11cm
      integer*4     Pol11       ! (%) ? Percentage polarization at 11cm
      integer*4     PA11        ! (deg) ? Polarisation angle at 11cm
      character*1   l_Pol6      ! Limit flag on Pol6cm
      integer*4     Pol6        ! (%) ? Percentage polarization at 6cm
      integer*4     PA6         ! (deg) ? Polarisation angle at 6cm
      character*1   l_Pol2_8    ! Limit flag on Pol2.8cm
      integer*4     Pol2_8      ! (%) ? Percentage polarization at 2.8cm
      integer*4     PA2_8       ! (deg) ? Polarisation angle at 2.8cm
      real*4        Ep11cm      ! (yr) ? Date of 11cm measurement (1900+)
      real*4        Ep6cm       ! (yr) ? Date of 6cm measurement (1900+)
      real*4        Ep2_8cm     ! (yr) ? Date of 2.8cm measurement (1900+)
      real*4        NH          ! (10+21cm-2) Column density of neutral hydrogen
      integer*4     RAh         ! (h) Right Ascension J2000 (hours) (5)
      integer*4     RAm         ! (min) Right Ascension J2000 (minutes) (5)
      real*4        RAs         ! (s) Right Ascension J2000 (seconds) (5)
      character*1   fPos        ! [*] indicates a position accurate to 1' (5)
      character*1   DE_         ! Declination J2000 (sign) (5)
      integer*4     DEd         ! (deg) Declination J2000 (degrees) (5)
      integer*4     DEm         ! (arcmin) Declination J2000 (minutes) (5)
      integer*4     DEs         ! (arcsec) ? Declination J2000 (seconds) (5)
*Note (1): Spectral index {alpha} as: S {prop.to} {nu}^{alpha}^
*Note (2): The following abbreviations are used:
*      N = Neumann et al. (1994A&AS..106..303N)
*      K = Kuehr et al. (1979 MPIfR preprint No.55, 1981A&AS...45..367K) and
*          Tabara & Inoue (1980A&AS...39..379T)
*      O = Ott et al. (1994A&A...284..331O) and
*          Tabara & Inoue (1980A&AS...39..379T)
*     GR = Reich et al., Table 1 (1998A&AS..131...11R)
*    GR1 = Reich et al., Table 3 (1998A&AS..131...11R)
*Note (3): Error class (quality) takes the following values:
*      1 = error <=  9%
*      2 = error <= 14%
*      3 = error <= 23%
*      4 = error > 23%
*Note (4): Optical identifications from Brinkmann et al. (1997A&A...323..739B)
*          and Laurent-Muehleisen et al. (1997A&AS..122..235L) as follows:
*      G = galaxy
*    QSO = quasar
*    BLL = BL Lac type object
*    GCL = galaxy cluster
*      * = star
*Note (5): the position was added in March 2007, from the log of the
*     observations at Bonn observatory. However 534 positions only could
*     be retrieved; the other positions, derived from the RGB name, are
*     accurate to 1' only and are flagged with a '*' in the fPos column.

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

C  Declarations for 'table2.dat'	! Repeated observations

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

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

      character*12  RGB_1       ! Green Bank catalog (<VIII/52>) name
      character*1   n_RGB_1     ! [*] *: source with repeated observations
      character*9   Bname_1     ! Name based on 1950 position (HHMM+DDMM)
      character*5   r_Pos_1     ! Position references (2)
      integer*4     F11cm_1     ! (mJy) ? Integrated 11cm (2.7GHz) flux density
      integer*4     F6cm_1      ! (mJy) ? Integrated 6cm (5GHz) flux density
      integer*4     F2_8cm_1    ! (mJy) ? Integrated 2.8cm (10.7GHz) flux density
      real*4        SI11_6_1    ! ? Spectral index 11cm/6cm (1)
      real*4        e_SI11_6_1  ! ? Error of spectral index 11cm/6cm
      real*4        F11_2_8_1   ! (mW/m2) Integrated flux density (11cm to 2.8cm)
      integer*4     q_F11cm_1   ! [1/4]? Error class see Note (3)
      integer*4     q_F6cm_1    ! [1/4]? Error class see Note (3)
      integer*4     q_F2_8cm_1  ! [1/4]? Error class see Note (3)
      real*4        SI6_2_8_1   ! ? Spectral index 6cm/2.8cm (1)
      real*4        e_SI6_2_8_1 ! ? Error of spectral index 6cm/2.8cm
      real*4        Fx_1        ! (mW/m2) Integrated flux density (0.1-2.4 KeV)
      real*4        z_1         ! ? Measured redshift
      real*4        Vmap_1      ! (mag) ? Optical magnitude (V at most)
      integer*4     D11a_1      ! (arcsec) ? RA source size 11cm
      character*1   n_D11a_1    ! [u] u: Unresolved (point-like sources)
      integer*4     D11d_1      ! (arcsec) ? DE source size 11cm
      character*1   n_D11d_1    ! [u] u: Unresolved (point-like sources)
      integer*4     D6a_1       ! (arcsec) ? RA source size 6cm
      character*1   n_D6a_1     ! [u] u: Unresolved (point-like sources)
      integer*4     D6d_1       ! (arcsec) ? DE source size 6cm
      character*1   n_D6d_1     ! [u] u: Unresolved (point-like sources)
      integer*4     D2_8a_1     ! (arcsec) ? RA source size 2.8cm
      character*1   n_D2_8a_1   ! [u] u: Unresolved (point-like sources)
      integer*4     D2_8d_1     ! (arcsec) ? DE source size 2.8cm
      character*1   n_D2_8d_1   ! [u] u: Unresolved (point-like sources)
      real*4        SI11__1     ! ? Spectral index fit 11cm/6cm/2.8cm (1)
      real*4        e_SI11__1   ! ? Error of spectral index fit
      real*4        SIx_1       ! ? X-ray power law photon index (1)
      real*4        e_SIx_1     ! ? Error of x-ray power law photon index
      integer*4     e_Fx_1      ! (%) Error of X-ray flux density
      character*3   OptID_1     ! Optical identification see Note (4)
      character*10  Name_1      ! Common source name
      character*1   l_Pol11_1   ! Limit flag on Pol11cm
      integer*4     Pol11_1     ! (%) ? Percentage polarization at 11cm
      integer*4     PA11_1      ! (deg) ? Polarisation angle at 11cm
      character*1   l_Pol6_1    ! Limit flag on Pol6cm
      integer*4     Pol6_1      ! (%) ? Percentage polarization at 6cm
      integer*4     PA6_1       ! (deg) ? Polarisation angle at 6cm
      character*1   l_Pol2_8_1  ! Limit flag on Pol2.8cm
      integer*4     Pol2_8_1    ! (%) ? Percentage polarization at 2.8cm
      integer*4     PA2_8_1     ! (deg) ? Polarisation angle at 2.8cm
      real*4        Ep11cm_1    ! (yr) ? Date of 11cm measurement (1900+)
      real*4        Ep6cm_1     ! (yr) ? Date of 6cm measurement (1900+)
      real*4        Ep2_8cm_1   ! (yr) ? Date of 2.8cm measurement (1900+)
      real*4        NH_1        ! (10+21cm-2) Column density of neutral hydrogen
      integer*4     RAh_1       ! (h) Right Ascension J2000 (hours) (5)
      integer*4     RAm_1       ! (min) Right Ascension J2000 (minutes) (5)
      real*4        RAs_1       ! (s) Right Ascension J2000 (seconds) (5)
      character*1   fPos_1      ! [*] indicates a position accurate to 1' (5)
      character*1   DE__1       ! Declination J2000 (sign) (5)
      integer*4     DEd_1       ! (deg) Declination J2000 (degrees) (5)
      integer*4     DEm_1       ! (arcmin) Declination J2000 (minutes) (5)
      integer*4     DEs_1       ! (arcsec) ? Declination J2000 (seconds) (5)
*Note (1): Spectral index {alpha} as: S {prop.to} {nu}^{alpha}^
*Note (2): The following abbreviations are used:
*      N = Neumann et al. (1994A&AS..106..303N)
*      K = Kuehr et al. (1979 MPIfR preprint No.55, 1981A&AS...45..367K) and
*          Tabara & Inoue (1980A&AS...39..379T)
*      O = Ott et al. (1994A&A...284..331O) and
*          Tabara & Inoue (1980A&AS...39..379T)
*     GR = Reich et al., Table 1 (1998A&AS..131...11R)
*    GR1 = Reich et al., Table 3 (1998A&AS..131...11R)
*Note (3): Error class (quality) takes the following values:
*      1 = error <=  9%
*      2 = error <= 14%
*      3 = error <= 23%
*      4 = error > 23%
*Note (4): Optical identifications from Brinkmann et al. (1997A&A...323..739B)
*          and Laurent-Muehleisen et al. (1997A&AS..122..235L) as follows:
*      G = galaxy
*    QSO = quasar
*    BLL = BL Lac type object
*    GCL = galaxy cluster
*      * = star
*Note (5): the position was added in March 2007, from the log of the
*     observations at Bonn observatory. However 534 positions only could
*     be retrieved; the other positions, derived from the RGB name, are
*     accurate to 1' only and are flagged with a '*' in the fPos column.

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

C  Loading file 'table1.dat'	! Data for 741 source

C  Format for file interpretation

    1 format(
     +  3X,A12,A1,1X,A9,1X,A5,1X,I6,1X,I5,1X,I5,1X,F5.2,1X,F4.2,1X,
     +  E10.4,1X,I1,1X,I1,1X,I1,1X,F5.2,1X,F4.2,1X,E10.4,1X,F6.4,1X,
     +  F5.2,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,
     +  F5.2,1X,F4.2,1X,F3.1,1X,F3.1,1X,I2,1X,A3,1X,A10,A1,I2,1X,I3,
     +  1X,A1,I2,1X,I3,1X,A1,I2,1X,I3,2X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
     +  1X,I2,1X,I2,1X,F4.1,A1,A1,I2,1X,I2,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,741
        read(1,'(A250)')ar__
        read(ar__,1)
     +  RGB,n_RGB,Bname,r_Pos,F11cm,F6cm,F2_8cm,SI11_6,e_SI11_6,
     +  F11_2_8,q_F11cm,q_F6cm,q_F2_8cm,SI6_2_8,e_SI6_2_8,Fx,z,Vmap,
     +  D11a,n_D11a,D11d,n_D11d,D6a,n_D6a,D6d,n_D6d,D2_8a,n_D2_8a,
     +  D2_8d,n_D2_8d,SI11_,e_SI11_,SIx,e_SIx,e_Fx,OptID,Name,l_Pol11,
     +  Pol11,PA11,l_Pol6,Pol6,PA6,l_Pol2_8,Pol2_8,PA2_8,Ep11cm,Ep6cm,
     +  Ep2_8cm,NH,RAh,RAm,RAs,fPos,DE_,DEd,DEm,DEs
        if(ar__(34:39) .EQ. '') F11cm = iNULL__
        if(ar__(41:45) .EQ. '') F6cm = iNULL__
        if(ar__(47:51) .EQ. '') F2_8cm = iNULL__
        if(ar__(53:57) .EQ. '') SI11_6 = rNULL__
        if(ar__(59:62) .EQ. '') e_SI11_6 = rNULL__
        if(ar__(75:75) .EQ. '') q_F11cm = iNULL__
        if(ar__(77:77) .EQ. '') q_F6cm = iNULL__
        if(ar__(79:79) .EQ. '') q_F2_8cm = iNULL__
        if(ar__(81:85) .EQ. '') SI6_2_8 = rNULL__
        if(ar__(87:90) .EQ. '') e_SI6_2_8 = rNULL__
        if(ar__(103:108) .EQ. '') z = rNULL__
        if(ar__(110:114) .EQ. '') Vmap = rNULL__
        if(ar__(116:118) .EQ. '') D11a = iNULL__
        if(ar__(121:123) .EQ. '') D11d = iNULL__
        if(ar__(126:128) .EQ. '') D6a = iNULL__
        if(ar__(131:133) .EQ. '') D6d = iNULL__
        if(ar__(136:138) .EQ. '') D2_8a = iNULL__
        if(ar__(141:143) .EQ. '') D2_8d = iNULL__
        if(ar__(146:150) .EQ. '') SI11_ = rNULL__
        if(ar__(152:155) .EQ. '') e_SI11_ = rNULL__
        if(ar__(157:159) .EQ. '') SIx = rNULL__
        if(ar__(161:163) .EQ. '') e_SIx = rNULL__
        if(ar__(183:184) .EQ. '') Pol11 = iNULL__
        if(ar__(186:188) .EQ. '') PA11 = iNULL__
        if(ar__(191:192) .EQ. '') Pol6 = iNULL__
        if(ar__(194:196) .EQ. '') PA6 = iNULL__
        if(ar__(199:200) .EQ. '') Pol2_8 = iNULL__
        if(ar__(202:204) .EQ. '') PA2_8 = iNULL__
        if(ar__(207:211) .EQ. '') Ep11cm = rNULL__
        if(ar__(213:217) .EQ. '') Ep6cm = rNULL__
        if(ar__(219:223) .EQ. '') Ep2_8cm = rNULL__
        if(ar__(249:250) .EQ. '') DEs = 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)
     +  RGB,n_RGB,Bname,r_Pos,F11cm,F6cm,F2_8cm,SI11_6,e_SI11_6,
     +  F11_2_8,q_F11cm,q_F6cm,q_F2_8cm,SI6_2_8,e_SI6_2_8,Fx,z,Vmap,
     +  D11a,n_D11a,D11d,n_D11d,D6a,n_D6a,D6d,n_D6d,D2_8a,n_D2_8a,
     +  D2_8d,n_D2_8d,SI11_,e_SI11_,SIx,e_SIx,e_Fx,OptID,Name,l_Pol11,
     +  Pol11,PA11,l_Pol6,Pol6,PA6,l_Pol2_8,Pol2_8,PA2_8,Ep11cm,Ep6cm,
     +  Ep2_8cm,NH,RAh,RAm,RAs,fPos,DE_,DEd,DEm,DEs
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Repeated observations

C  Format for file interpretation

    2 format(
     +  3X,A12,A1,1X,A9,1X,A5,1X,I6,1X,I5,1X,I5,1X,F5.2,1X,F4.2,1X,
     +  E10.4,1X,I1,1X,I1,1X,I1,1X,F5.2,1X,F4.2,1X,E10.4,1X,F6.4,1X,
     +  F5.2,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,I3,A1,1X,
     +  F5.2,1X,F4.2,1X,F3.1,1X,F3.1,1X,I2,1X,A3,1X,A10,A1,I2,1X,I3,
     +  1X,A1,I2,1X,I3,1X,A1,I2,1X,I3,2X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
     +  1X,I2,1X,I2,1X,F4.1,A1,A1,I2,1X,I2,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,57
        read(1,'(A250)')ar__1
        read(ar__1,2)
     +  RGB_1,n_RGB_1,Bname_1,r_Pos_1,F11cm_1,F6cm_1,F2_8cm_1,
     +  SI11_6_1,e_SI11_6_1,F11_2_8_1,q_F11cm_1,q_F6cm_1,q_F2_8cm_1,
     +  SI6_2_8_1,e_SI6_2_8_1,Fx_1,z_1,Vmap_1,D11a_1,n_D11a_1,D11d_1,
     +  n_D11d_1,D6a_1,n_D6a_1,D6d_1,n_D6d_1,D2_8a_1,n_D2_8a_1,
     +  D2_8d_1,n_D2_8d_1,SI11__1,e_SI11__1,SIx_1,e_SIx_1,e_Fx_1,
     +  OptID_1,Name_1,l_Pol11_1,Pol11_1,PA11_1,l_Pol6_1,Pol6_1,PA6_1,
     +  l_Pol2_8_1,Pol2_8_1,PA2_8_1,Ep11cm_1,Ep6cm_1,Ep2_8cm_1,NH_1,
     +  RAh_1,RAm_1,RAs_1,fPos_1,DE__1,DEd_1,DEm_1,DEs_1
        if(ar__1(34:39) .EQ. '') F11cm_1 = iNULL__
        if(ar__1(41:45) .EQ. '') F6cm_1 = iNULL__
        if(ar__1(47:51) .EQ. '') F2_8cm_1 = iNULL__
        if(ar__1(53:57) .EQ. '') SI11_6_1 = rNULL__
        if(ar__1(59:62) .EQ. '') e_SI11_6_1 = rNULL__
        if(ar__1(75:75) .EQ. '') q_F11cm_1 = iNULL__
        if(ar__1(77:77) .EQ. '') q_F6cm_1 = iNULL__
        if(ar__1(79:79) .EQ. '') q_F2_8cm_1 = iNULL__
        if(ar__1(81:85) .EQ. '') SI6_2_8_1 = rNULL__
        if(ar__1(87:90) .EQ. '') e_SI6_2_8_1 = rNULL__
        if(ar__1(103:108) .EQ. '') z_1 = rNULL__
        if(ar__1(110:114) .EQ. '') Vmap_1 = rNULL__
        if(ar__1(116:118) .EQ. '') D11a_1 = iNULL__
        if(ar__1(121:123) .EQ. '') D11d_1 = iNULL__
        if(ar__1(126:128) .EQ. '') D6a_1 = iNULL__
        if(ar__1(131:133) .EQ. '') D6d_1 = iNULL__
        if(ar__1(136:138) .EQ. '') D2_8a_1 = iNULL__
        if(ar__1(141:143) .EQ. '') D2_8d_1 = iNULL__
        if(ar__1(146:150) .EQ. '') SI11__1 = rNULL__
        if(ar__1(152:155) .EQ. '') e_SI11__1 = rNULL__
        if(ar__1(157:159) .EQ. '') SIx_1 = rNULL__
        if(ar__1(161:163) .EQ. '') e_SIx_1 = rNULL__
        if(ar__1(183:184) .EQ. '') Pol11_1 = iNULL__
        if(ar__1(186:188) .EQ. '') PA11_1 = iNULL__
        if(ar__1(191:192) .EQ. '') Pol6_1 = iNULL__
        if(ar__1(194:196) .EQ. '') PA6_1 = iNULL__
        if(ar__1(199:200) .EQ. '') Pol2_8_1 = iNULL__
        if(ar__1(202:204) .EQ. '') PA2_8_1 = iNULL__
        if(ar__1(207:211) .EQ. '') Ep11cm_1 = rNULL__
        if(ar__1(213:217) .EQ. '') Ep6cm_1 = rNULL__
        if(ar__1(219:223) .EQ. '') Ep2_8cm_1 = rNULL__
        if(ar__1(249:250) .EQ. '') DEs_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)
     +  RGB_1,n_RGB_1,Bname_1,r_Pos_1,F11cm_1,F6cm_1,F2_8cm_1,
     +  SI11_6_1,e_SI11_6_1,F11_2_8_1,q_F11cm_1,q_F6cm_1,q_F2_8cm_1,
     +  SI6_2_8_1,e_SI6_2_8_1,Fx_1,z_1,Vmap_1,D11a_1,n_D11a_1,D11d_1,
     +  n_D11d_1,D6a_1,n_D6a_1,D6d_1,n_D6d_1,D2_8a_1,n_D2_8a_1,
     +  D2_8d_1,n_D2_8d_1,SI11__1,e_SI11__1,SIx_1,e_SIx_1,e_Fx_1,
     +  OptID_1,Name_1,l_Pol11_1,Pol11_1,PA11_1,l_Pol6_1,Pol6_1,PA6_1,
     +  l_Pol2_8_1,Pol2_8_1,PA2_8_1,Ep11cm_1,Ep6cm_1,Ep2_8cm_1,NH_1,
     +  RAh_1,RAm_1,RAs_1,fPos_1,DE__1,DEd_1,DEm_1,DEs_1
        write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c    .......End.of.Just test output...........
      end do
      close(1)

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