Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/482/483 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 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/A+A/482/483 GHz-peaked spectrum (GPS) sources (Torniainen+, 2008)
*================================================================================
*Cluster analyses of gigahertz-peaked spectrum sources with self-organising maps.
* Torniainen I., Tornikoski M., Turunen M., Lainela M., Lahteenmaki A.,
* Hovatta T., Mingaliev M.G., Aller M.F., Aller H.D.
* <Astron. Astrophys. 482, 483 (2008)>
* =2008A&A...482..483T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table.dat' ! Data for each source
integer*4 nr__
parameter (nr__=206) ! Number of records
character*368 ar__ ! Full-size record
character*10 SRC (nr__) ! Source name (HHMM+DDd/d)
character*6 ID (nr__) ! Optical ID (1)
character*7 r_ID (nr__) ! Reference for optical ID
character*1 u_ID (nr__) ! [)] Uncertainty flag on r_ID
character*5 OClass (nr__) ! Original classification, as GPS or HFP (High-
* Frequency Peaker)
character*7 r_OClass (nr__) ! Reference for original classification
character*5 Class (nr__) ! Auxiliary classification (3)
real*4 z (nr__) ! ?=9.99 Redshift
character*6 r_z (nr__) ! Reference for redshift
character*1 f_z (nr__) ! [*] * for photometric redshift
real*4 Szkpc (nr__) ! (kpc) ?=99.99 Size of the source
character*6 r_Szkpc (nr__) ! Reference for size
character*1 f_Szkpc (nr__) ! [*] * for recalculated size (2)
real*4 P5GHz (nr__) ! ([W/Hz]) ?=99.99 log(10) of source power at 5GHz (5)
character*4 r_P5GHz (nr__) ! [calc] Reference for P5GHz
real*4 polOpt (nr__) ! (%) ?=99.99 Optical polarization
character*8 r_polOpt (nr__) ! Reference for optical polarization
real*4 polRad (nr__) ! (%) ?=99.99 Radio polarization
character*7 r_polRad (nr__) ! Reference for radio polarization
character*1 f_polRad (nr__) ! [*] * for zuk99 and step95 reference
character*6 n_polRad (nr__) ! Note on polRad
real*4 Bmag (nr__) ! (mag) ?=99.99 B magnitude
character*5 r_Bmag (nr__) ! Reference for B magnitude
real*4 Vmag (nr__) ! (mag) ?=99.99 V magnitude
character*6 r_Vmag (nr__) ! Reference for V magnitude
real*4 Rmag (nr__) ! (mag) ?=99.99 R magnitude
character*7 r_Rmag (nr__) ! Reference for R magnitude
character*1 f_Rmag (nr__) ! [*] * when calculated from r magnitude
real*4 Imag (nr__) ! (mag) ?=99.99 I magnitude
character*6 r_Imag (nr__) ! Reference for I magnitude
character*1 f_Imag (nr__) ! [*] * when calculated from i magnitude
real*4 V_R (nr__) ! (mag) ?=99.99 V-R colour index
character*5 r_V_R (nr__) ! Reference for V-R colour
real*4 B_V (nr__) ! (mag) ?=99.99 B-V colour index
character*4 r_B_V (nr__) ! Reference for B-V colour
real*4 U_B (nr__) ! (mag) ?=99.99 U-B colour index
character*4 r_U_B (nr__) ! [vv03] Reference for U-B colour
real*4 R_I (nr__) ! (mag) ?=9.99 R-I colour index
character*6 r_R_I (nr__) ! Reference for R-I colour
real*4 Freq (nr__) ! (GHz) ?=99.99 Rest-frame frequency of the turnover
* of the spectrum
real*4 ab (nr__) ! ?=99.99 Spectral index below the turnover
real*4 aa (nr__) ! ?=999.99 Spectral index above the turnover
real*4 symm (nr__) ! ?=999.99 Symmetry of the spectrum
real*4 curv (nr__) ! ?=99.99 Curvature of the spectrum
real*4 FWHM (nr__) ! ?=99.99 Full width at half maximum of the
* spectrum
real*4 Mid (nr__) ! ?=999.99 Mid-frequency of the bin used in
* calculations
real*4 VI (nr__) ! ?=99.99 Variability index
integer*4 o_VI (nr__) ! ?=9999 Number of data points used for VI
real*8 NHX (nr__) ! (10+22cm-2) ?=9999.99 Hydrogen column density in X-rays
character*5 r_NHX (nr__) ! Reference for N_H_X
real*4 Gamma (nr__) ! ?=9.99 Power-law slope in X-rays
character*5 r_Gamma (nr__) ! Reference for power-law slope in X-rays
real*4 NH21 (nr__) ! (10+20cm-2) ?=99.99 Neutral hydrogen column density
character*6 r_NH21 (nr__) ! Reference for NH21
real*4 Emag (nr__) ! (mag) ?=99.99 E (photographic red) magnitude
character*7 r_Emag (nr__) ! reference for E magnitude
real*4 Omag (nr__) ! (mag) ?=99.99 O (photographic blue) magnitude
character*7 r_Omag (nr__) ! Reference for O magnitude
real*4 O_E (nr__) ! (mag) ?=99.99 O-E colour index
character*7 r_O_E (nr__) ! Reference for O-E colour
character*5 Morph (nr__) ! VLBI morphology (4)
character*6 r_Morph (nr__) ! Reference for VLBI morphology
*Note (1): Optical identification as follows:
* BLO = BL Lac object
* BLO/G = BL Lac object/galaxy
* EF = empty field
* GAL = galaxy
* GA/QS = galaxy/QSO
* HPQ = high polarization quasar
* LPQ = low polarization quasar
* LPQ&bl = low polarization quasar and bl
* S = stellar
* QSO = quasar
*Note (2): Unless otherwise stated * means that the linear size has
* been recalculated from the angular size with H_0_=71km/s/Mpc,
* Omega_M_=0.27, Omega_vac_=0.73 with Edward L. Wright's Cosmology
* Calculator (http://www.astro.ucla.edu/~wright/CosmoCalc.html).
*Note (3): Auxiliary classification as follows:
* n = not enough data for GPS identification
* s = steep spectrum
* f = flat spectrum
* f/s = flat at low frequencies, steep at high frequencies
* c = convex spectrum
* idb = inverted during bursts
* v = variability of Var_DeltaS>3
* gps = gigahertz peaked spectrum
*Note (4): VLBI morphology as follows:
* cd = compact double
* cj = core-jet object
* cj/cd = core-jet object/compact double
* cso = compact symmetric object
* ct = compact triple
* cx = complex
* gl = gravitational lens
* ln = linear
* s = stellar
* unres = unresolved
* mr = marginally resolved
* p = point-like
*Note (5): Unit [W/Hz] was not specified in the paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__1
parameter (nr__1=64) ! Number of records
character*74 ar__1 ! Full-size record
character*7 Ref (nr__1) ! Reference code
character*19 BibCode (nr__1) ! BibCode
character*22 Aut (nr__1) ! Author's name
character*23 Com (nr__1) ! Comments
C=============================================================================
C Loading file 'table.dat' ! Data for each source
C Format for file interpretation
1 format(
+ A10,1X,A6,1X,A7,A1,A5,1X,A7,1X,A5,1X,F4.2,1X,A6,A1,1X,F5.2,1X,
+ A6,A1,1X,F5.2,1X,A4,1X,F5.2,1X,A8,1X,F5.2,1X,A7,A1,A6,1X,F5.2,
+ 1X,A5,1X,F5.2,1X,A6,1X,F5.2,1X,A7,A1,1X,F5.2,1X,A6,A1,1X,F5.2,
+ 1X,A5,1X,F5.2,1X,A4,1X,F5.2,1X,A4,1X,F4.2,1X,A6,1X,F5.2,1X,
+ F5.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,I4,1X,
+ F7.2,1X,A5,1X,F4.2,1X,A5,1X,F5.2,1X,A6,1X,F5.2,1X,A7,1X,F5.2,
+ 1X,A7,1X,F5.2,1X,A7,1X,A5,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table.dat')
write(6,*) '....Loading file: table.dat'
do i__=1,206
read(1,'(A368)')ar__
read(ar__,1)
+ SRC(i__),ID(i__),r_ID(i__),u_ID(i__),OClass(i__),
+ r_OClass(i__),Class(i__),z(i__),r_z(i__),f_z(i__),Szkpc(i__),
+ r_Szkpc(i__),f_Szkpc(i__),P5GHz(i__),r_P5GHz(i__),polOpt(i__),
+ r_polOpt(i__),polRad(i__),r_polRad(i__),f_polRad(i__),
+ n_polRad(i__),Bmag(i__),r_Bmag(i__),Vmag(i__),r_Vmag(i__),
+ Rmag(i__),r_Rmag(i__),f_Rmag(i__),Imag(i__),r_Imag(i__),
+ f_Imag(i__),V_R(i__),r_V_R(i__),B_V(i__),r_B_V(i__),U_B(i__),
+ r_U_B(i__),R_I(i__),r_R_I(i__),Freq(i__),ab(i__),aa(i__),
+ symm(i__),curv(i__),FWHM(i__),Mid(i__),VI(i__),o_VI(i__),
+ NHX(i__),r_NHX(i__),Gamma(i__),r_Gamma(i__),NH21(i__),
+ r_NH21(i__),Emag(i__),r_Emag(i__),Omag(i__),r_Omag(i__),
+ O_E(i__),r_O_E(i__),Morph(i__),r_Morph(i__)
c ..............Just test output...........
write(6,1)
+ SRC(i__),ID(i__),r_ID(i__),u_ID(i__),OClass(i__),
+ r_OClass(i__),Class(i__),z(i__),r_z(i__),f_z(i__),Szkpc(i__),
+ r_Szkpc(i__),f_Szkpc(i__),P5GHz(i__),r_P5GHz(i__),polOpt(i__),
+ r_polOpt(i__),polRad(i__),r_polRad(i__),f_polRad(i__),
+ n_polRad(i__),Bmag(i__),r_Bmag(i__),Vmag(i__),r_Vmag(i__),
+ Rmag(i__),r_Rmag(i__),f_Rmag(i__),Imag(i__),r_Imag(i__),
+ f_Imag(i__),V_R(i__),r_V_R(i__),B_V(i__),r_B_V(i__),U_B(i__),
+ r_U_B(i__),R_I(i__),r_R_I(i__),Freq(i__),ab(i__),aa(i__),
+ symm(i__),curv(i__),FWHM(i__),Mid(i__),VI(i__),o_VI(i__),
+ NHX(i__),r_NHX(i__),Gamma(i__),r_Gamma(i__),NH21(i__),
+ r_NH21(i__),Emag(i__),r_Emag(i__),Omag(i__),r_Omag(i__),
+ O_E(i__),r_O_E(i__),Morph(i__),r_Morph(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
2 format(A7,1X,A19,1X,A22,1X,A23)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,64
read(1,'(A74)')ar__1
read(ar__1,2)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,2)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end