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 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 ! (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