Conversion of standardized ReadMe file for
file /./ftp/cats/VII/166 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-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 (nr__) ! (deg) Right Ascension 1950
real*8 DEdeg (nr__) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*1 FLAG (nr__) ! [*] A flag which is blank if the object has
* been detected as a radio source and which
* is '*' if not.
character*17 Name (nr__) ! Most common name of the object.
integer*4 RAh (nr__) ! (h) Right Ascension 1950 (hours)
character*1 n_RAh (nr__) ! 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 (nr__) ! (min) RA1950 (minutes)
real*4 RAs (nr__) ! (s) RA1950 (seconds)
character*1 DE_ (nr__) ! Declination 1950 (Sign)
integer*4 DEd (nr__) ! (deg) DE1950 (degrees)
integer*4 DEm (nr__) ! (arcmin) DE1950 (')
integer*4 DEs (nr__) ! (arcsec) DE1950 (")
character*1 n_Vmag (nr__) ! [/R] Note on the "V" magnitude
* '/' indicates photographic magnitude
* and colours when the colours are given.
* 'R' indicates a red magnitude
real*4 Vmag (nr__) ! (mag) ? V magnitude, or photographic or red
* magnitude depending on n_Vmag
real*4 B_V (nr__) ! (mag) ? Colour index
real*4 U_B (nr__) ! (mag) ? Colour index
character*1 l_z (nr__) ! [/>] 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 (nr__) ! ? Value of the redshift
real*4 F6cm (nr__) ! (Jy) ? 6cm flux (Jy)
integer*4 r_F6cm (nr__) ! [1,1486]?=0 Reference of 6cm Flux
real*4 F11cm (nr__) ! (Jy) ? 11cm flux (Jy)
integer*4 r_F11cm (nr__) ! [1,1486]?=0 Reference of 11cm flux
character*1 Xray (nr__) ! [X] 'X' when the object was detected
* in the X-ray range
character*2 Class (nr__) ! [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 (nr__) ! [1,1486]? Reference of the Finding Chart
integer*4 r_Vmag (nr__) ! [1,1486]? Reference of the Photometry
integer*4 r_z (nr__) ! [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 (nr__1) ! (deg) Right Ascension 1950
real*8 DEdeg_1 (nr__1) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*1 FLAG_1 (nr__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 (nr__1) ! Most common name of the object.
integer*4 RAh_1 (nr__1) ! (h) Right Ascension 1950 (hours)
character*1 n_RAh_1 (nr__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 (nr__1) ! (min) RA1950 (minutes)
real*4 RAs_1 (nr__1) ! (s) RA1950 (seconds)
character*1 DE__1 (nr__1) ! Declination 1950 (Sign)
integer*4 DEd_1 (nr__1) ! (deg) DE1950 (degrees)
integer*4 DEm_1 (nr__1) ! (arcmin) DE1950 (')
integer*4 DEs_1 (nr__1) ! (arcsec) DE1950 (")
character*1 n_Vmag_1 (nr__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 (nr__1) ! (mag) ? V magnitude, or photographic or red
* magnitude depending on n_Vmag
real*4 B_V_1 (nr__1) ! (mag) ? Colour index
real*4 U_B_1 (nr__1) ! (mag) ? Colour index
character*1 l_z_1 (nr__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 (nr__1) ! ? Value of the redshift
real*4 F6cm_1 (nr__1) ! (Jy) ? 6cm flux (Jy)
integer*4 r_F6cm_1 (nr__1) ! [1,1486]?=0 Reference of 6cm Flux
real*4 F11cm_1 (nr__1) ! (Jy) ? 11cm flux (Jy)
integer*4 r_F11cm_1 (nr__1) ! [1,1486]?=0 Reference of 11cm flux
character*1 Xray_1 (nr__1) ! [X] 'X' when the object was detected
* in the X-ray range
character*2 Class_1 (nr__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 (nr__1) ! [1,1486]? Reference of the Finding Chart
integer*4 r_Vmag_1 (nr__1) ! [1,1486]? Reference of the Photometry
integer*4 r_z_1 (nr__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 (nr__2) ! (deg) Right Ascension 1950
real*8 DEdeg_2 (nr__2) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*1 FLAG_2 (nr__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 (nr__2) ! Most common name of the object.
integer*4 RAh_2 (nr__2) ! (h) Right Ascension 1950 (hours)
character*1 n_RAh_2 (nr__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 (nr__2) ! (min) RA1950 (minutes)
real*4 RAs_2 (nr__2) ! (s) RA1950 (seconds)
character*1 DE__2 (nr__2) ! Declination 1950 (Sign)
integer*4 DEd_2 (nr__2) ! (deg) DE1950 (degrees)
integer*4 DEm_2 (nr__2) ! (arcmin) DE1950 (')
integer*4 DEs_2 (nr__2) ! (arcsec) DE1950 (")
character*1 n_Vmag_2 (nr__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 (nr__2) ! (mag) ? V magnitude, or photographic or red
* magnitude depending on n_Vmag
real*4 B_V_2 (nr__2) ! (mag) ? Colour index
real*4 U_B_2 (nr__2) ! (mag) ? Colour index
character*1 l_z_2 (nr__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 (nr__2) ! ? Value of the redshift
real*4 F6cm_2 (nr__2) ! (Jy) ? 6cm flux (Jy)
integer*4 r_F6cm_2 (nr__2) ! [1,1486]?=0 Reference of 6cm Flux
real*4 F11cm_2 (nr__2) ! (Jy) ? 11cm flux (Jy)
integer*4 r_F11cm_2 (nr__2) ! [1,1486]?=0 Reference of 11cm flux
character*1 Xray_2 (nr__2) ! [X] 'X' when the object was detected
* in the X-ray range
character*2 Class_2 (nr__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 (nr__2) ! [1,1486]? Reference of the Finding Chart
integer*4 r_Vmag_2 (nr__2) ! [1,1486]? Reference of the Photometry
integer*4 r_z_2 (nr__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 (nr__3) ! (deg) Right Ascension 1950
real*8 DEdeg_3 (nr__3) ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
character*1 FLAG_3 (nr__3) ! [*] Flag similar to table1
character*11 Name_3 (nr__3) ! Most common name of the object.
character*6 v4C (nr__3) ! Its number in 4C radio source catalogue
integer*4 RAh_3 (nr__3) ! (h) Right Ascension 1950 (hours)
character*1 n_RAh_3 (nr__3) ! Note on the position, as in table1
integer*4 RAm_3 (nr__3) ! (min) RA1950 (minutes)
real*4 RAs_3 (nr__3) ! (s) RA1950 (seconds)
character*1 DE__3 (nr__3) ! Declination 1950 (Sign)
integer*4 DEd_3 (nr__3) ! (deg) DE1950 (degrees)
integer*4 DEm_3 (nr__3) ! (arcmin) DE1950 (')
integer*4 DEs_3 (nr__3) ! (arcsec) DE1950 (")
character*1 n_Vmag_3 (nr__3) ! [/R] Note on Vmag, as in table1
real*4 Vmag_3 (nr__3) ! (mag) ? V magnitude, as in table1
real*4 B_V_3 (nr__3) ! (mag) ? Colour index
real*4 U_B_3 (nr__3) ! (mag) ? Colour index
character*1 l_z_3 (nr__3) ! [/>] Note on the redshift, as in table1
real*4 z_3 (nr__3) ! ? Value of the redshift
character*2 Class0 (nr__3) ! Original classification, as in table1
character*2 Class1 (nr__3) ! [*G? ] New classification, with symbols
* 'G' = galaxy
* '*' = Star
integer*4 FC_3 (nr__3) ! [1,1486]? Reference of the Finding Chart
integer*4 r_Vmag_3 (nr__3) ! [1,1486]? Reference of the Photometry
integer*4 r_z_3 (nr__3) ! [1,1486]? Reference of the redshift
integer*4 r_Class1 (nr__3) ! [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 (nr__4) ! Reference number
character*76 Text (nr__4) ! 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(i__),Name(i__),RAh(i__),n_RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),n_Vmag(i__),Vmag(i__),
+ B_V(i__),U_B(i__),l_z(i__),z(i__),F6cm(i__),r_F6cm(i__),
+ F11cm(i__),r_F11cm(i__),Xray(i__),Class(i__),FC(i__),
+ r_Vmag(i__),r_z(i__)
if(ar__(40:44) .EQ. '') Vmag(i__) = rNULL__
if(ar__(46:50) .EQ. '') B_V(i__) = rNULL__
if(ar__(51:55) .EQ. '') U_B(i__) = rNULL__
if(ar__(57:61) .EQ. '') z(i__) = rNULL__
if(ar__(62:67) .EQ. '') F6cm(i__) = rNULL__
if(ar__(72:76) .EQ. '') F11cm(i__) = rNULL__
if(ar__(87:90) .EQ. '') FC(i__) = iNULL__
if(ar__(92:95) .EQ. '') r_Vmag(i__) = iNULL__
if(ar__(97:100) .EQ. '') r_z(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)
+ FLAG(i__),Name(i__),RAh(i__),n_RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),n_Vmag(i__),Vmag(i__),
+ B_V(i__),U_B(i__),l_z(i__),z(i__),F6cm(i__),r_F6cm(i__),
+ F11cm(i__),r_F11cm(i__),Xray(i__),Class(i__),FC(i__),
+ r_Vmag(i__),r_z(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' ! 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(i__),Name_1(i__),RAh_1(i__),n_RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
+ n_Vmag_1(i__),Vmag_1(i__),B_V_1(i__),U_B_1(i__),l_z_1(i__),
+ z_1(i__),F6cm_1(i__),r_F6cm_1(i__),F11cm_1(i__),
+ r_F11cm_1(i__),Xray_1(i__),Class_1(i__),FC_1(i__),
+ r_Vmag_1(i__),r_z_1(i__)
if(ar__1(40:44) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__1(46:50) .EQ. '') B_V_1(i__) = rNULL__
if(ar__1(51:55) .EQ. '') U_B_1(i__) = rNULL__
if(ar__1(57:61) .EQ. '') z_1(i__) = rNULL__
if(ar__1(62:67) .EQ. '') F6cm_1(i__) = rNULL__
if(ar__1(72:76) .EQ. '') F11cm_1(i__) = rNULL__
if(ar__1(87:90) .EQ. '') FC_1(i__) = iNULL__
if(ar__1(92:95) .EQ. '') r_Vmag_1(i__) = iNULL__
if(ar__1(97:100) .EQ. '') r_z_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)
+ FLAG_1(i__),Name_1(i__),RAh_1(i__),n_RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
+ n_Vmag_1(i__),Vmag_1(i__),B_V_1(i__),U_B_1(i__),l_z_1(i__),
+ z_1(i__),F6cm_1(i__),r_F6cm_1(i__),F11cm_1(i__),
+ r_F11cm_1(i__),Xray_1(i__),Class_1(i__),FC_1(i__),
+ r_Vmag_1(i__),r_z_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=============================================================================
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(i__),Name_2(i__),RAh_2(i__),n_RAh_2(i__),RAm_2(i__),
+ RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
+ n_Vmag_2(i__),Vmag_2(i__),B_V_2(i__),U_B_2(i__),l_z_2(i__),
+ z_2(i__),F6cm_2(i__),r_F6cm_2(i__),F11cm_2(i__),
+ r_F11cm_2(i__),Xray_2(i__),Class_2(i__),FC_2(i__),
+ r_Vmag_2(i__),r_z_2(i__)
if(ar__2(40:44) .EQ. '') Vmag_2(i__) = rNULL__
if(ar__2(46:50) .EQ. '') B_V_2(i__) = rNULL__
if(ar__2(51:55) .EQ. '') U_B_2(i__) = rNULL__
if(ar__2(57:61) .EQ. '') z_2(i__) = rNULL__
if(ar__2(62:67) .EQ. '') F6cm_2(i__) = rNULL__
if(ar__2(72:76) .EQ. '') F11cm_2(i__) = rNULL__
if(ar__2(87:90) .EQ. '') FC_2(i__) = iNULL__
if(ar__2(92:95) .EQ. '') r_Vmag_2(i__) = iNULL__
if(ar__2(97:100) .EQ. '') r_z_2(i__) = iNULL__
RAdeg_2(i__) = rNULL__
DEdeg_2(i__) = 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(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c ..............Just test output...........
write(6,3)
+ FLAG_2(i__),Name_2(i__),RAh_2(i__),n_RAh_2(i__),RAm_2(i__),
+ RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
+ n_Vmag_2(i__),Vmag_2(i__),B_V_2(i__),U_B_2(i__),l_z_2(i__),
+ z_2(i__),F6cm_2(i__),r_F6cm_2(i__),F11cm_2(i__),
+ r_F11cm_2(i__),Xray_2(i__),Class_2(i__),FC_2(i__),
+ r_Vmag_2(i__),r_z_2(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
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(i__),Name_3(i__),v4C(i__),RAh_3(i__),n_RAh_3(i__),
+ RAm_3(i__),RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),
+ DEs_3(i__),n_Vmag_3(i__),Vmag_3(i__),B_V_3(i__),U_B_3(i__),
+ l_z_3(i__),z_3(i__),Class0(i__),Class1(i__),FC_3(i__),
+ r_Vmag_3(i__),r_z_3(i__),r_Class1(i__)
if(ar__3(40:44) .EQ. '') Vmag_3(i__) = rNULL__
if(ar__3(46:50) .EQ. '') B_V_3(i__) = rNULL__
if(ar__3(51:55) .EQ. '') U_B_3(i__) = rNULL__
if(ar__3(57:61) .EQ. '') z_3(i__) = rNULL__
if(ar__3(68:71) .EQ. '') FC_3(i__) = iNULL__
if(ar__3(73:76) .EQ. '') r_Vmag_3(i__) = iNULL__
if(ar__3(78:81) .EQ. '') r_z_3(i__) = iNULL__
if(ar__3(83:86) .EQ. '') r_Class1(i__) = iNULL__
RAdeg_3(i__) = rNULL__
DEdeg_3(i__) = 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(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c ..............Just test output...........
write(6,4)
+ FLAG_3(i__),Name_3(i__),v4C(i__),RAh_3(i__),n_RAh_3(i__),
+ RAm_3(i__),RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),
+ DEs_3(i__),n_Vmag_3(i__),Vmag_3(i__),B_V_3(i__),U_B_3(i__),
+ l_z_3(i__),z_3(i__),Class0(i__),Class1(i__),FC_3(i__),
+ r_Vmag_3(i__),r_z_3(i__),r_Class1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
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(i__),Text(i__)
c ..............Just test output...........
write(6,5)Ref(i__),Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end