Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/363/141 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-15
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 (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*12 RGB (nr__) ! Green Bank catalog (<VIII/52>) name
character*1 n_RGB (nr__) ! [*] *: source with repeated observations
character*9 Bname (nr__) ! Name based on 1950 position (HHMM+DDMM)
character*5 r_Pos (nr__) ! Position references (2)
integer*4 F11cm (nr__) ! (mJy) ? Integrated 11cm (2.7GHz) flux density
integer*4 F6cm (nr__) ! (mJy) ? Integrated 6cm (5GHz) flux density
integer*4 F2_8cm (nr__) ! (mJy) ? Integrated 2.8cm (10.7GHz) flux density
real*4 SI11_6 (nr__) ! ? Spectral index 11cm/6cm (1)
real*4 e_SI11_6 (nr__) ! ? Error of spectral index 11cm/6cm
real*4 F11_2_8 (nr__) ! (mW/m2) Integrated flux density (11cm to 2.8cm)
integer*4 q_F11cm (nr__) ! [1/4]? Error class see Note (3)
integer*4 q_F6cm (nr__) ! [1/4]? Error class see Note (3)
integer*4 q_F2_8cm (nr__) ! [1/4]? Error class see Note (3)
real*4 SI6_2_8 (nr__) ! ? Spectral index 6cm/2.8cm (1)
real*4 e_SI6_2_8 (nr__) ! ? Error of spectral index 6cm/2.8cm
real*4 Fx (nr__) ! (mW/m2) Integrated flux density (0.1-2.4 KeV)
real*4 z (nr__) ! ? Measured redshift
real*4 Vmap (nr__) ! (mag) ? Optical magnitude (V at most)
integer*4 D11a (nr__) ! (arcsec) ? RA source size 11cm
character*1 n_D11a (nr__) ! [u] u: Unresolved (point-like sources)
integer*4 D11d (nr__) ! (arcsec) ? DE source size 11cm
character*1 n_D11d (nr__) ! [u] u: Unresolved (point-like sources)
integer*4 D6a (nr__) ! (arcsec) ? RA source size 6cm
character*1 n_D6a (nr__) ! [u] u: Unresolved (point-like sources)
integer*4 D6d (nr__) ! (arcsec) ? DE source size 6cm
character*1 n_D6d (nr__) ! [u] u: Unresolved (point-like sources)
integer*4 D2_8a (nr__) ! (arcsec) ? RA source size 2.8cm
character*1 n_D2_8a (nr__) ! [u] u: Unresolved (point-like sources)
integer*4 D2_8d (nr__) ! (arcsec) ? DE source size 2.8cm
character*1 n_D2_8d (nr__) ! [u] u: Unresolved (point-like sources)
real*4 SI11_ (nr__) ! ? Spectral index fit 11cm/6cm/2.8cm (1)
real*4 e_SI11_ (nr__) ! ? Error of spectral index fit
real*4 SIx (nr__) ! ? X-ray power law photon index (1)
real*4 e_SIx (nr__) ! ? Error of x-ray power law photon index
integer*4 e_Fx (nr__) ! (%) Error of X-ray flux density
character*3 OptID (nr__) ! Optical identification see Note (4)
character*10 Name (nr__) ! Common source name
character*1 l_Pol11 (nr__) ! Limit flag on Pol11cm
integer*4 Pol11 (nr__) ! (%) ? Percentage polarization at 11cm
integer*4 PA11 (nr__) ! (deg) ? Polarisation angle at 11cm
character*1 l_Pol6 (nr__) ! Limit flag on Pol6cm
integer*4 Pol6 (nr__) ! (%) ? Percentage polarization at 6cm
integer*4 PA6 (nr__) ! (deg) ? Polarisation angle at 6cm
character*1 l_Pol2_8 (nr__) ! Limit flag on Pol2.8cm
integer*4 Pol2_8 (nr__) ! (%) ? Percentage polarization at 2.8cm
integer*4 PA2_8 (nr__) ! (deg) ? Polarisation angle at 2.8cm
real*4 Ep11cm (nr__) ! (yr) ? Date of 11cm measurement (1900+)
real*4 Ep6cm (nr__) ! (yr) ? Date of 6cm measurement (1900+)
real*4 Ep2_8cm (nr__) ! (yr) ? Date of 2.8cm measurement (1900+)
real*4 NH (nr__) ! (10+21cm-2) Column density of neutral hydrogen
integer*4 RAh (nr__) ! (h) Right Ascension J2000 (hours) (5)
integer*4 RAm (nr__) ! (min) Right Ascension J2000 (minutes) (5)
real*4 RAs (nr__) ! (s) Right Ascension J2000 (seconds) (5)
character*1 fPos (nr__) ! [*] indicates a position accurate to 1' (5)
character*1 DE_ (nr__) ! Declination J2000 (sign) (5)
integer*4 DEd (nr__) ! (deg) Declination J2000 (degrees) (5)
integer*4 DEm (nr__) ! (arcmin) Declination J2000 (minutes) (5)
integer*4 DEs (nr__) ! (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 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*12 RGB_1 (nr__1) ! Green Bank catalog (<VIII/52>) name
character*1 n_RGB_1 (nr__1) ! [*] *: source with repeated observations
character*9 Bname_1 (nr__1) ! Name based on 1950 position (HHMM+DDMM)
character*5 r_Pos_1 (nr__1) ! Position references (2)
integer*4 F11cm_1 (nr__1) ! (mJy) ? Integrated 11cm (2.7GHz) flux density
integer*4 F6cm_1 (nr__1) ! (mJy) ? Integrated 6cm (5GHz) flux density
integer*4 F2_8cm_1 (nr__1) ! (mJy) ? Integrated 2.8cm (10.7GHz) flux density
real*4 SI11_6_1 (nr__1) ! ? Spectral index 11cm/6cm (1)
real*4 e_SI11_6_1 (nr__1) ! ? Error of spectral index 11cm/6cm
real*4 F11_2_8_1 (nr__1) ! (mW/m2) Integrated flux density (11cm to 2.8cm)
integer*4 q_F11cm_1 (nr__1) ! [1/4]? Error class see Note (3)
integer*4 q_F6cm_1 (nr__1) ! [1/4]? Error class see Note (3)
integer*4 q_F2_8cm_1 (nr__1) ! [1/4]? Error class see Note (3)
real*4 SI6_2_8_1 (nr__1) ! ? Spectral index 6cm/2.8cm (1)
real*4 e_SI6_2_8_1(nr__1) ! ? Error of spectral index 6cm/2.8cm
real*4 Fx_1 (nr__1) ! (mW/m2) Integrated flux density (0.1-2.4 KeV)
real*4 z_1 (nr__1) ! ? Measured redshift
real*4 Vmap_1 (nr__1) ! (mag) ? Optical magnitude (V at most)
integer*4 D11a_1 (nr__1) ! (arcsec) ? RA source size 11cm
character*1 n_D11a_1 (nr__1) ! [u] u: Unresolved (point-like sources)
integer*4 D11d_1 (nr__1) ! (arcsec) ? DE source size 11cm
character*1 n_D11d_1 (nr__1) ! [u] u: Unresolved (point-like sources)
integer*4 D6a_1 (nr__1) ! (arcsec) ? RA source size 6cm
character*1 n_D6a_1 (nr__1) ! [u] u: Unresolved (point-like sources)
integer*4 D6d_1 (nr__1) ! (arcsec) ? DE source size 6cm
character*1 n_D6d_1 (nr__1) ! [u] u: Unresolved (point-like sources)
integer*4 D2_8a_1 (nr__1) ! (arcsec) ? RA source size 2.8cm
character*1 n_D2_8a_1 (nr__1) ! [u] u: Unresolved (point-like sources)
integer*4 D2_8d_1 (nr__1) ! (arcsec) ? DE source size 2.8cm
character*1 n_D2_8d_1 (nr__1) ! [u] u: Unresolved (point-like sources)
real*4 SI11__1 (nr__1) ! ? Spectral index fit 11cm/6cm/2.8cm (1)
real*4 e_SI11__1 (nr__1) ! ? Error of spectral index fit
real*4 SIx_1 (nr__1) ! ? X-ray power law photon index (1)
real*4 e_SIx_1 (nr__1) ! ? Error of x-ray power law photon index
integer*4 e_Fx_1 (nr__1) ! (%) Error of X-ray flux density
character*3 OptID_1 (nr__1) ! Optical identification see Note (4)
character*10 Name_1 (nr__1) ! Common source name
character*1 l_Pol11_1 (nr__1) ! Limit flag on Pol11cm
integer*4 Pol11_1 (nr__1) ! (%) ? Percentage polarization at 11cm
integer*4 PA11_1 (nr__1) ! (deg) ? Polarisation angle at 11cm
character*1 l_Pol6_1 (nr__1) ! Limit flag on Pol6cm
integer*4 Pol6_1 (nr__1) ! (%) ? Percentage polarization at 6cm
integer*4 PA6_1 (nr__1) ! (deg) ? Polarisation angle at 6cm
character*1 l_Pol2_8_1 (nr__1) ! Limit flag on Pol2.8cm
integer*4 Pol2_8_1 (nr__1) ! (%) ? Percentage polarization at 2.8cm
integer*4 PA2_8_1 (nr__1) ! (deg) ? Polarisation angle at 2.8cm
real*4 Ep11cm_1 (nr__1) ! (yr) ? Date of 11cm measurement (1900+)
real*4 Ep6cm_1 (nr__1) ! (yr) ? Date of 6cm measurement (1900+)
real*4 Ep2_8cm_1 (nr__1) ! (yr) ? Date of 2.8cm measurement (1900+)
real*4 NH_1 (nr__1) ! (10+21cm-2) Column density of neutral hydrogen
integer*4 RAh_1 (nr__1) ! (h) Right Ascension J2000 (hours) (5)
integer*4 RAm_1 (nr__1) ! (min) Right Ascension J2000 (minutes) (5)
real*4 RAs_1 (nr__1) ! (s) Right Ascension J2000 (seconds) (5)
character*1 fPos_1 (nr__1) ! [*] indicates a position accurate to 1' (5)
character*1 DE__1 (nr__1) ! Declination J2000 (sign) (5)
integer*4 DEd_1 (nr__1) ! (deg) Declination J2000 (degrees) (5)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination J2000 (minutes) (5)
integer*4 DEs_1 (nr__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(i__),n_RGB(i__),Bname(i__),r_Pos(i__),F11cm(i__),
+ F6cm(i__),F2_8cm(i__),SI11_6(i__),e_SI11_6(i__),F11_2_8(i__),
+ q_F11cm(i__),q_F6cm(i__),q_F2_8cm(i__),SI6_2_8(i__),
+ e_SI6_2_8(i__),Fx(i__),z(i__),Vmap(i__),D11a(i__),n_D11a(i__),
+ D11d(i__),n_D11d(i__),D6a(i__),n_D6a(i__),D6d(i__),n_D6d(i__),
+ D2_8a(i__),n_D2_8a(i__),D2_8d(i__),n_D2_8d(i__),SI11_(i__),
+ e_SI11_(i__),SIx(i__),e_SIx(i__),e_Fx(i__),OptID(i__),
+ Name(i__),l_Pol11(i__),Pol11(i__),PA11(i__),l_Pol6(i__),
+ Pol6(i__),PA6(i__),l_Pol2_8(i__),Pol2_8(i__),PA2_8(i__),
+ Ep11cm(i__),Ep6cm(i__),Ep2_8cm(i__),NH(i__),RAh(i__),RAm(i__),
+ RAs(i__),fPos(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__)
if(ar__(34:39) .EQ. '') F11cm(i__) = iNULL__
if(ar__(41:45) .EQ. '') F6cm(i__) = iNULL__
if(ar__(47:51) .EQ. '') F2_8cm(i__) = iNULL__
if(ar__(53:57) .EQ. '') SI11_6(i__) = rNULL__
if(ar__(59:62) .EQ. '') e_SI11_6(i__) = rNULL__
if(ar__(75:75) .EQ. '') q_F11cm(i__) = iNULL__
if(ar__(77:77) .EQ. '') q_F6cm(i__) = iNULL__
if(ar__(79:79) .EQ. '') q_F2_8cm(i__) = iNULL__
if(ar__(81:85) .EQ. '') SI6_2_8(i__) = rNULL__
if(ar__(87:90) .EQ. '') e_SI6_2_8(i__) = rNULL__
if(ar__(103:108) .EQ. '') z(i__) = rNULL__
if(ar__(110:114) .EQ. '') Vmap(i__) = rNULL__
if(ar__(116:118) .EQ. '') D11a(i__) = iNULL__
if(ar__(121:123) .EQ. '') D11d(i__) = iNULL__
if(ar__(126:128) .EQ. '') D6a(i__) = iNULL__
if(ar__(131:133) .EQ. '') D6d(i__) = iNULL__
if(ar__(136:138) .EQ. '') D2_8a(i__) = iNULL__
if(ar__(141:143) .EQ. '') D2_8d(i__) = iNULL__
if(ar__(146:150) .EQ. '') SI11_(i__) = rNULL__
if(ar__(152:155) .EQ. '') e_SI11_(i__) = rNULL__
if(ar__(157:159) .EQ. '') SIx(i__) = rNULL__
if(ar__(161:163) .EQ. '') e_SIx(i__) = rNULL__
if(ar__(183:184) .EQ. '') Pol11(i__) = iNULL__
if(ar__(186:188) .EQ. '') PA11(i__) = iNULL__
if(ar__(191:192) .EQ. '') Pol6(i__) = iNULL__
if(ar__(194:196) .EQ. '') PA6(i__) = iNULL__
if(ar__(199:200) .EQ. '') Pol2_8(i__) = iNULL__
if(ar__(202:204) .EQ. '') PA2_8(i__) = iNULL__
if(ar__(207:211) .EQ. '') Ep11cm(i__) = rNULL__
if(ar__(213:217) .EQ. '') Ep6cm(i__) = rNULL__
if(ar__(219:223) .EQ. '') Ep2_8cm(i__) = rNULL__
if(ar__(249:250) .EQ. '') DEs(i__) = iNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ RGB(i__),n_RGB(i__),Bname(i__),r_Pos(i__),F11cm(i__),
+ F6cm(i__),F2_8cm(i__),SI11_6(i__),e_SI11_6(i__),F11_2_8(i__),
+ q_F11cm(i__),q_F6cm(i__),q_F2_8cm(i__),SI6_2_8(i__),
+ e_SI6_2_8(i__),Fx(i__),z(i__),Vmap(i__),D11a(i__),n_D11a(i__),
+ D11d(i__),n_D11d(i__),D6a(i__),n_D6a(i__),D6d(i__),n_D6d(i__),
+ D2_8a(i__),n_D2_8a(i__),D2_8d(i__),n_D2_8d(i__),SI11_(i__),
+ e_SI11_(i__),SIx(i__),e_SIx(i__),e_Fx(i__),OptID(i__),
+ Name(i__),l_Pol11(i__),Pol11(i__),PA11(i__),l_Pol6(i__),
+ Pol6(i__),PA6(i__),l_Pol2_8(i__),Pol2_8(i__),PA2_8(i__),
+ Ep11cm(i__),Ep6cm(i__),Ep2_8cm(i__),NH(i__),RAh(i__),RAm(i__),
+ RAs(i__),fPos(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),n_RGB_1(i__),Bname_1(i__),r_Pos_1(i__),
+ F11cm_1(i__),F6cm_1(i__),F2_8cm_1(i__),SI11_6_1(i__),
+ e_SI11_6_1(i__),F11_2_8_1(i__),q_F11cm_1(i__),q_F6cm_1(i__),
+ q_F2_8cm_1(i__),SI6_2_8_1(i__),e_SI6_2_8_1(i__),Fx_1(i__),
+ z_1(i__),Vmap_1(i__),D11a_1(i__),n_D11a_1(i__),D11d_1(i__),
+ n_D11d_1(i__),D6a_1(i__),n_D6a_1(i__),D6d_1(i__),n_D6d_1(i__),
+ D2_8a_1(i__),n_D2_8a_1(i__),D2_8d_1(i__),n_D2_8d_1(i__),
+ SI11__1(i__),e_SI11__1(i__),SIx_1(i__),e_SIx_1(i__),
+ e_Fx_1(i__),OptID_1(i__),Name_1(i__),l_Pol11_1(i__),
+ Pol11_1(i__),PA11_1(i__),l_Pol6_1(i__),Pol6_1(i__),PA6_1(i__),
+ l_Pol2_8_1(i__),Pol2_8_1(i__),PA2_8_1(i__),Ep11cm_1(i__),
+ Ep6cm_1(i__),Ep2_8cm_1(i__),NH_1(i__),RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),fPos_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__)
if(ar__1(34:39) .EQ. '') F11cm_1(i__) = iNULL__
if(ar__1(41:45) .EQ. '') F6cm_1(i__) = iNULL__
if(ar__1(47:51) .EQ. '') F2_8cm_1(i__) = iNULL__
if(ar__1(53:57) .EQ. '') SI11_6_1(i__) = rNULL__
if(ar__1(59:62) .EQ. '') e_SI11_6_1(i__) = rNULL__
if(ar__1(75:75) .EQ. '') q_F11cm_1(i__) = iNULL__
if(ar__1(77:77) .EQ. '') q_F6cm_1(i__) = iNULL__
if(ar__1(79:79) .EQ. '') q_F2_8cm_1(i__) = iNULL__
if(ar__1(81:85) .EQ. '') SI6_2_8_1(i__) = rNULL__
if(ar__1(87:90) .EQ. '') e_SI6_2_8_1(i__) = rNULL__
if(ar__1(103:108) .EQ. '') z_1(i__) = rNULL__
if(ar__1(110:114) .EQ. '') Vmap_1(i__) = rNULL__
if(ar__1(116:118) .EQ. '') D11a_1(i__) = iNULL__
if(ar__1(121:123) .EQ. '') D11d_1(i__) = iNULL__
if(ar__1(126:128) .EQ. '') D6a_1(i__) = iNULL__
if(ar__1(131:133) .EQ. '') D6d_1(i__) = iNULL__
if(ar__1(136:138) .EQ. '') D2_8a_1(i__) = iNULL__
if(ar__1(141:143) .EQ. '') D2_8d_1(i__) = iNULL__
if(ar__1(146:150) .EQ. '') SI11__1(i__) = rNULL__
if(ar__1(152:155) .EQ. '') e_SI11__1(i__) = rNULL__
if(ar__1(157:159) .EQ. '') SIx_1(i__) = rNULL__
if(ar__1(161:163) .EQ. '') e_SIx_1(i__) = rNULL__
if(ar__1(183:184) .EQ. '') Pol11_1(i__) = iNULL__
if(ar__1(186:188) .EQ. '') PA11_1(i__) = iNULL__
if(ar__1(191:192) .EQ. '') Pol6_1(i__) = iNULL__
if(ar__1(194:196) .EQ. '') PA6_1(i__) = iNULL__
if(ar__1(199:200) .EQ. '') Pol2_8_1(i__) = iNULL__
if(ar__1(202:204) .EQ. '') PA2_8_1(i__) = iNULL__
if(ar__1(207:211) .EQ. '') Ep11cm_1(i__) = rNULL__
if(ar__1(213:217) .EQ. '') Ep6cm_1(i__) = rNULL__
if(ar__1(219:223) .EQ. '') Ep2_8cm_1(i__) = rNULL__
if(ar__1(249:250) .EQ. '') DEs_1(i__) = iNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ RGB_1(i__),n_RGB_1(i__),Bname_1(i__),r_Pos_1(i__),
+ F11cm_1(i__),F6cm_1(i__),F2_8cm_1(i__),SI11_6_1(i__),
+ e_SI11_6_1(i__),F11_2_8_1(i__),q_F11cm_1(i__),q_F6cm_1(i__),
+ q_F2_8cm_1(i__),SI6_2_8_1(i__),e_SI6_2_8_1(i__),Fx_1(i__),
+ z_1(i__),Vmap_1(i__),D11a_1(i__),n_D11a_1(i__),D11d_1(i__),
+ n_D11d_1(i__),D6a_1(i__),n_D6a_1(i__),D6d_1(i__),n_D6d_1(i__),
+ D2_8a_1(i__),n_D2_8a_1(i__),D2_8d_1(i__),n_D2_8d_1(i__),
+ SI11__1(i__),e_SI11__1(i__),SIx_1(i__),e_SIx_1(i__),
+ e_Fx_1(i__),OptID_1(i__),Name_1(i__),l_Pol11_1(i__),
+ Pol11_1(i__),PA11_1(i__),l_Pol6_1(i__),Pol6_1(i__),PA6_1(i__),
+ l_Pol2_8_1(i__),Pol2_8_1(i__),PA2_8_1(i__),Ep11cm_1(i__),
+ Ep6cm_1(i__),Ep2_8cm_1(i__),NH_1(i__),RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),fPos_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end