Conversion of standardized ReadMe file for
file /./ftp/cats/J/AZh/73/835 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-Apr-16
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/AZh/73/835 Identifications of objects from RC Catalog (Fletcher+, 1996)
*================================================================================
*Identification of the Objects from the RC catalog with accuracy exceeding 1"
*using the Palomar Prints and VLA maps from the MIT archive.
* Fletcher A., Conner S., Crawford F., Cartwright J., Burke B.,
* Pariiskii Yu.N., Soboleva N.S., Kopylov A.I.,Temirova A.V.,
* Verkhodanov O.V., Naugol'naya M.N.
* <Astron. Zh. 73, 835 (1996)>
* =1996AZh....73..835F
* =1996ARep...40..759F
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table.dat' ! Optical identifications
integer*4 nr__
parameter (nr__=102) ! Number of records
character*146 ar__ ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension 1950
real*8 DEdeg ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq ! Ordered source number
character*10 RC ! RC catalogue source number (J omitted)
integer*4 RAh ! (h) Right ascension (1950) (from VLA map)
integer*4 RAm ! (min) Right ascension (1950) (from VLA map)
real*4 RAs ! (s) Right ascension (1950) (from VLA map)
character*1 DE_ ! Declination sign (1950) (from VLA map)
integer*4 DEd ! (deg) Declination (1950) (from VLA map)
integer*4 DEm ! (arcmin) Declination (1950) (from VLA map)
real*4 DEs ! (arcsec) Declination (1950) (from VLA map)
real*4 Sp_Index ! ? Spectral index, alpha (-dlogS/dlog{nu})
character*1 u_Sp_Index ! [:] Uncertainty flag on spectral index
character*1 n_Sp_Index ! [s ] s - indicates that the spectrum
* becomes steeper at higher frequencies
integer*4 RAoh ! (h) ? Right ascension (1950) (optical)
integer*4 RAom ! (min) ? Right ascension (1950) (optical)
real*4 RAos ! (s) ? Right ascension (1950) (optical)
character*1 DEo_ ! ? Declination sign (1950) (optical)
integer*4 DEod ! (deg) ? Declination (1950) (optical)
integer*4 DEom ! (arcmin) ? Declination (1950) (optical)
real*4 DEos ! (arcsec) ? Declination (1950) (optical)
character*3 Notes ! [EF ?] Notes: EF - empty field on the plate
character*3 Mo ! ? Optical morphology (star/galaxy)
character*1 l_Rmag ! Limit flag on Rmag
real*4 Rmag ! (mag) ? Optical magnitude in R band
real*4 dr ! (arcsec) ? The difference between the optical and
* radio coordinates
integer*4 S1_4GHz ! (mJy) ? Flux density at 1.4GHz
character*1 u_S1_4GHz ! [:] Uncertainty flag on S1.4GHz
real*4 Rlim ! (mag) ? Limiting magnitude in R band for which
* radio luminosity is equal to optical one
real*4 dm ! (mag) ? Difference between the limiting and
* observed magnitudes
character*5 RType ! Radio morphology (1)
character*20 Comments ! Comments (2)
*Note (1): Radio morphology:
* P - point,
* D - double,
* T - triple with core component,
* M - multiple components,
* Pj- point with a jet,
* P+J - point with a distant jet,
* ext - single extended object,
* P(el) - single elliptical component
*Note (2): Comments:
* Law - object identified earlier by Lawrence et al., 1986ApJS...61..105L
* * - objects identified earlier in [7-14]
C=============================================================================
C Loading file 'table.dat' ! Optical identifications
C Format for file interpretation
1 format(
+ I2,2X,A10,2X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,2X,F5.2,
+ A1,A1,2X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,A3,3X,A3,
+ 1X,A1,F5.2,2X,F4.2,2X,I4,A1,2X,F4.1,1X,F5.2,2X,A5,2X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'table.dat')
write(6,*) '....Loading file: table.dat'
do i__=1,102
read(1,'(A146)')ar__
read(ar__,1)
+ Seq,RC,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Sp_Index,u_Sp_Index,
+ n_Sp_Index,RAoh,RAom,RAos,DEo_,DEod,DEom,DEos,Notes,Mo,l_Rmag,
+ Rmag,dr,S1_4GHz,u_S1_4GHz,Rlim,dm,RType,Comments
if(ar__(44:48) .EQ. '') Sp_Index = rNULL__
if(ar__(53:54) .EQ. '') RAoh = iNULL__
if(ar__(56:57) .EQ. '') RAom = iNULL__
if(ar__(59:63) .EQ. '') RAos = rNULL__
if(ar__(66:67) .EQ. '') DEod = iNULL__
if(ar__(69:70) .EQ. '') DEom = iNULL__
if(ar__(72:75) .EQ. '') DEos = rNULL__
if(ar__(88:92) .EQ. '') Rmag = rNULL__
if(ar__(95:98) .EQ. '') dr = rNULL__
if(ar__(101:104) .EQ. '') S1_4GHz = iNULL__
if(ar__(108:111) .EQ. '') Rlim = rNULL__
if(ar__(113:117) .EQ. '') dm = rNULL__
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)
+ Seq,RC,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Sp_Index,u_Sp_Index,
+ n_Sp_Index,RAoh,RAom,RAos,DEo_,DEod,DEom,DEos,Notes,Mo,l_Rmag,
+ Rmag,dr,S1_4GHz,u_S1_4GHz,Rlim,dm,RType,Comments
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end