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