Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/366/1265 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-Sep-06
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/MNRAS/366/1265 CENSORS infrared imaging (Brookes+, 2006)
*================================================================================
*CENSORS: A Combined EIS-NVSS Survey of Radio Sources -
*II. Infrared imaging and the K-z relation.
* Brookes M.H., Best P.N., Rengelink R., Roettgering H.J.A.
* <Mon. Not. R. Astron. Soc., 366, 1265-1288 (2006)>
* =2006MNRAS.366.1265B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Details of observing runs for K-band imaging
of CENSORS.
integer*4 nr__
parameter (nr__=6) ! Number of records
character*38 ar__ ! Full-size record
integer*4 Run (nr__) ! Run number
character*1 n_Run (nr__) ! [*] Note on Run (1)
character*21 Dates (nr__) ! Dates of observation
character*5 Tel (nr__) ! Telescope name
character*5 Inst (nr__) ! Instrument name
integer*4 Night (nr__) ! Number of photometric nights
*Note (1): *: included in this data set are two observations made in
* service mode: CENSORS 51 and 98 in 2002 November 27 and
* 2003 February 4, respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Details of K-band observations of CENSORS and
results.
integer*4 nr__1
parameter (nr__1=154) ! Number of records
character*117 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 CENSORS (nr__1) ! CENSORS sequential number
character*1 m_CENSORS (nr__1) ! [ABCDE] Muntiplicity index on CENSORS number
integer*4 EISD (nr__1) ! EISD sequential number
integer*4 Run1 (nr__1) ! ? Observation run
integer*4 Run2 (nr__1) ! ? Observation run
integer*4 Run3 (nr__1) ! ? Observation run
real*8 ExpTime1 (nr__1) ! (s) ? Exposure time
integer*4 ExpTime2 (nr__1) ! (s) ? Exposure time
integer*4 ExpTime3 (nr__1) ! (s) ? Exposure time
integer*4 RAh (nr__1) ! (h) ? Right ascension (J2000)
integer*4 RAm (nr__1) ! (min) ? Right ascension (J2000)
real*4 RAs (nr__1) ! (s) ? Right ascension (J2000)
character*1 DE_ (nr__1) ! Declination sign (J2000)
integer*4 DEd (nr__1) ! (deg) ? Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) ? Declination (J2000)
real*4 DEs (nr__1) ! (arcsec) ? Declination (J2000)
real*4 Kmag1 (nr__1) ! (mag) ? K magnitude in aperture of radius 1arcsec
real*4 e_Kmag1 (nr__1) ! (mag) ? rms uncertainty on Kmag1
real*4 Kmag1_5 (nr__1) ! (mag) ? K magnitude in aperture of radius 1.5arcsec
real*4 e_Kmag1_5 (nr__1) ! (mag) ? rms uncertainty on Kmag1.5
real*4 Kmag2_5 (nr__1) ! (mag) ? K magnitude in aperture of radius 2.5arcsec
real*4 e_Kmag2_5 (nr__1) ! (mag) ? rms uncertainty on Kmag2.5
real*4 Kmag4_5 (nr__1) ! (mag) ? K magnitude in aperture of radius 4.5arcsec
real*4 e_Kmag4_5 (nr__1) ! (mag) ? rms uncertainty on Kmag4.5
real*4 Kmago (nr__1) ! (mag) ? K magnitude in aperture of radius theta
real*4 e_Kmago (nr__1) ! (mag) ? rms uncertainty on Kmago
character*3 theta (nr__1) ! Radius for the Kmago aperture (0.5 or rec)
real*4 KmagLim (nr__1) ! (mag) ? 2{sigma} limit K magnitude, for undetected
* sources, based upon a 1" radius aperture
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Aperture-corrected K-band magnitudes for all
sources with a measured spectroscopic redshift
(presented in Paper 3, in prep.) and a K-band
detection.
integer*4 nr__2
parameter (nr__2=77) ! Number of records
character*35 ar__2 ! Full-size record
integer*4 CENSORS_1 (nr__2) ! CENSORS sequential number
integer*4 EISD_1 (nr__2) ! EISD sequential number
real*4 z (nr__2) ! Redshift
real*4 thetam (nr__2) ! (arcsec) Aperture radius used to measure the
* K-band magnitude
real*4 Kmag_thetam(nr__2) ! (mag) Measured aperture magnitude
* (in thetam radius)
real*4 Kmag (nr__2) ! (mag) Magnitude after correction to 63.9kpc
integer*4 e_Kmag (nr__2) ! (0.01mag) rms uncertainty on Kmag
character*1 Class (nr__2) ! [GQS] Classification (1)
*Note (1): Classification as follows:
* G = radio galaxies
* Q = quasars
* S = star forming galaxies
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Individual notes
integer*4 nr__3
parameter (nr__3=123) ! Number of records
character*80 ar__3 ! Full-size record
integer*4 CENSORS_2 (nr__3) ! CENSORS sequential number
character*76 Note (nr__3) ! Text of the note
C=============================================================================
C Loading file 'table1.dat' ! Details of observing runs for K-band imaging
* of CENSORS.
C Format for file interpretation
1 format(I1,A1,1X,A21,1X,A5,1X,A5,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,6
read(1,'(A38)')ar__
read(ar__,1)
+ Run(i__),n_Run(i__),Dates(i__),Tel(i__),Inst(i__),Night(i__)
c ..............Just test output...........
write(6,1)
+ Run(i__),n_Run(i__),Dates(i__),Tel(i__),Inst(i__),Night(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Details of K-band observations of CENSORS and
* results.
C Format for file interpretation
2 format(
+ I3,A1,1X,I3,1X,I1,1X,I1,1X,I1,1X,F7.2,1X,I4,1X,I4,1X,I2,1X,I2,
+ 1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2,1X,F3.2,1X,F5.2,1X,
+ F3.2,1X,F5.2,1X,F3.2,1X,F5.2,1X,F3.2,1X,F5.2,1X,F3.2,1X,A3,2X,
+ F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,154
read(1,'(A117)')ar__1
read(ar__1,2)
+ CENSORS(i__),m_CENSORS(i__),EISD(i__),Run1(i__),Run2(i__),
+ Run3(i__),ExpTime1(i__),ExpTime2(i__),ExpTime3(i__),RAh(i__),
+ RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
+ Kmag1(i__),e_Kmag1(i__),Kmag1_5(i__),e_Kmag1_5(i__),
+ Kmag2_5(i__),e_Kmag2_5(i__),Kmag4_5(i__),e_Kmag4_5(i__),
+ Kmago(i__),e_Kmago(i__),theta(i__),KmagLim(i__)
if(ar__1(10:10) .EQ. '') Run1(i__) = iNULL__
if(ar__1(12:12) .EQ. '') Run2(i__) = iNULL__
if(ar__1(14:14) .EQ. '') Run3(i__) = iNULL__
if(ar__1(16:22) .EQ. '') ExpTime1(i__) = rNULL__
if(ar__1(24:27) .EQ. '') ExpTime2(i__) = iNULL__
if(ar__1(29:32) .EQ. '') ExpTime3(i__) = iNULL__
if(ar__1(34:35) .EQ. '') RAh(i__) = iNULL__
if(ar__1(37:38) .EQ. '') RAm(i__) = iNULL__
if(ar__1(40:44) .EQ. '') RAs(i__) = rNULL__
if(ar__1(47:48) .EQ. '') DEd(i__) = iNULL__
if(ar__1(50:51) .EQ. '') DEm(i__) = iNULL__
if(ar__1(53:57) .EQ. '') DEs(i__) = rNULL__
if(ar__1(59:63) .EQ. '') Kmag1(i__) = rNULL__
if(ar__1(65:67) .EQ. '') e_Kmag1(i__) = rNULL__
if(ar__1(69:73) .EQ. '') Kmag1_5(i__) = rNULL__
if(ar__1(75:77) .EQ. '') e_Kmag1_5(i__) = rNULL__
if(ar__1(79:83) .EQ. '') Kmag2_5(i__) = rNULL__
if(ar__1(85:87) .EQ. '') e_Kmag2_5(i__) = rNULL__
if(ar__1(89:93) .EQ. '') Kmag4_5(i__) = rNULL__
if(ar__1(95:97) .EQ. '') e_Kmag4_5(i__) = rNULL__
if(ar__1(99:103) .EQ. '') Kmago(i__) = rNULL__
if(ar__1(105:107) .EQ. '') e_Kmago(i__) = rNULL__
if(ar__1(114:117) .EQ. '') KmagLim(i__) = rNULL__
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,2)
+ CENSORS(i__),m_CENSORS(i__),EISD(i__),Run1(i__),Run2(i__),
+ Run3(i__),ExpTime1(i__),ExpTime2(i__),ExpTime3(i__),RAh(i__),
+ RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
+ Kmag1(i__),e_Kmag1(i__),Kmag1_5(i__),e_Kmag1_5(i__),
+ Kmag2_5(i__),e_Kmag2_5(i__),Kmag4_5(i__),e_Kmag4_5(i__),
+ Kmago(i__),e_Kmago(i__),theta(i__),KmagLim(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 'table4.dat' ! Aperture-corrected K-band magnitudes for all
* sources with a measured spectroscopic redshift
* (presented in Paper 3, in prep.) and a K-band
* detection.
C Format for file interpretation
3 format(I3,1X,I3,1X,F6.4,1X,F3.1,1X,F5.2,1X,F5.2,1X,I2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,77
read(1,'(A35)')ar__2
read(ar__2,3)
+ CENSORS_1(i__),EISD_1(i__),z(i__),thetam(i__),
+ Kmag_thetam(i__),Kmag(i__),e_Kmag(i__),Class(i__)
c ..............Just test output...........
write(6,3)
+ CENSORS_1(i__),EISD_1(i__),z(i__),thetam(i__),
+ Kmag_thetam(i__),Kmag(i__),e_Kmag(i__),Class(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Individual notes
C Format for file interpretation
4 format(I3,1X,A76)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,123
read(1,'(A80)')ar__3
read(ar__3,4)CENSORS_2(i__),Note(i__)
c ..............Just test output...........
write(6,4)CENSORS_2(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end