Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/484/119 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-Mar-29 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/484/119 Multi-epoch VLBI survey of CJF sources (Britzen+, 2008) *================================================================================ *A multi-epoch VLBI Survey of the kinematics of Caltech-Jodrell Bank *flat-spectrum sources. * Britzen S., Vermeulen R.C., Campbell R.M., Taylor G.B., Pearson T.J., * Readhead A.C.S., Xu W., Browne I.W., Henstock D.R., Wilkinson P. * <Astron. Astrophys. 484, 119 (2008)> * =2008A&A...484..119B C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'sources.dat' ! Source list integer*4 nr__ parameter (nr__=266) ! Number of records character*72 ar__ ! Full-size record character*8 Name (nr__) ! IAU-name of the source HHMM+DDd (B1950), * CJF BHHMM+DDd in Simbad character*19 Name2000 (nr__) ! Name with J2000 position character*13 FileNameR (nr__) ! Image name of pmtot, pmr and pmphi as function * of radial distance from the core, * in subdirectory fig character*13 FileNameT (nr__) ! Image name x and y of each component as a * function of time, in subdirectory fig (1) character*14 FileNameXY (nr__) ! Image name of positions of all components at * all epochs, in subdirectory fig (2) *Note (1): For each component, a "stationary" dotted line is plotted for * reference. *Note (2): The components are color-coded and the epochs are represented by * different plotting symbols. The modeled positions at each epoch are * shown in small black symbols. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablec1.dat' ! The kinematic properties of 779 AGN-components integer*4 nr__1 parameter (nr__1=779) ! Number of records character*188 ar__1 ! Full-size record character*8 Name_1 (nr__1) ! IAU-name of the source HHMM+DDd (B1950), * CJF BHHMM+DDd in Simbad character*10 Class (nr__1) ! Optical classification (1) character*2 l_z (nr__1) ! [>=] Limit flag on z real*4 z (nr__1) ! ? Redshift, see details in paper character*1 u_z (nr__1) ! [?] Uncertainty flag on z character*3 Comp (nr__1) ! Jet-component identification (2) integer*4 q (nr__1) ! [1/3] Quality class (3) integer*4 NDet (nr__1) ! Number of times this component has been * detected real*4 Span (nr__1) ! (yr) Time span from the first to the last detection real*8 Epoch (nr__1) ! (d) Reference epoch -2400000 real*8 X0 (nr__1) ! (mas) Component position at the reference epoch in * rectangular X-coordinate real*4 e_X0 (nr__1) ! (mas) ? rms uncertainty on X0 (4) real*8 Y0 (nr__1) ! (mas) Component position at the reference epoch in * rectangular Y-coordinate real*4 e_Y0 (nr__1) ! (mas) ? rms uncertainty on Y0 (4) real*4 pmX (nr__1) ! (mas/yr) Proper motion components in X real*4 e_pmX (nr__1) ! (mas/yr) ? rms uncertainty on pmX (4) real*4 pmY (nr__1) ! (mas/yr) Proper motion components in Y real*4 e_pmY (nr__1) ! (mas/yr) ? rms uncertainty on pmY (4) real*4 pmr (nr__1) ! (mas/yr) Radial proper motion component real*4 e_pmr (nr__1) ! (mas/yr) ? rms uncertainty on pmr (4) real*4 pmphi (nr__1) ! (mas/yr) Orthogonal proper motion component real*4 e_pmphi (nr__1) ! (mas/yr) ? rms uncertainty on pmphi (4) real*4 pmtot (nr__1) ! (mas/yr) Total proper motion, calculated from pmx, pmy real*4 e_pmtot (nr__1) ! (mas/yr) ? rms uncertainty on pmtot (4) real*8 betar (nr__1) ! (c) ? Radial apparent velocity real*4 e_betar (nr__1) ! (c) ? rms uncertainty on betar (4) real*8 betaphi (nr__1) ! (c) ? Orthogonal apparent velocity real*4 e_betaphi (nr__1) ! (c) ? rms uncertainty on betaphi (4) real*8 betatot (nr__1) ! (c) ? Total apparent velocity real*4 e_betatot (nr__1) ! (c) ? rms uncertainty on betatot (4) *Note (1): Optical classification as follows: * Q = quasar * B = BL Lac object * G = galaxy * U = unclassified *Note (2): Jet-component identification as: * C plus a number increasing with ordinal separation from the core * CC in case of counter-jet denotes a counter-jet component *Note (3): Quality flag is 1 for best to 3 for suffers from some problems * in identification *Note (4): Components that have observations at only two epochs do not have * an associated uncertainty ; with no degrees of freedom in the fit for * the kinematic model (four constraints and four unknowns), the scaling * to ~chi=1 is ill-defined. C============================================================================= C Loading file 'sources.dat' ! Source list C Format for file interpretation 1 format(A8,1X,A19,2X,A13,1X,A13,1X,A14) C Effective file loading open(unit=1,status='old',file= +'sources.dat') write(6,*) '....Loading file: sources.dat' do i__=1,266 read(1,'(A72)')ar__ read(ar__,1) + Name(i__),Name2000(i__),FileNameR(i__),FileNameT(i__), + FileNameXY(i__) c ..............Just test output........... write(6,1) + Name(i__),Name2000(i__),FileNameR(i__),FileNameT(i__), + FileNameXY(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablec1.dat' ! The kinematic properties of 779 AGN-components C Format for file interpretation 2 format( + A8,1X,A10,A2,F5.3,A1,1X,A3,1X,I1,1X,I1,1X,F4.2,1X,F10.4,1X, + F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3, + 1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F7.3,1X, + F6.3,1X,F7.3,1X,F6.3,1X,F7.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'tablec1.dat') write(6,*) '....Loading file: tablec1.dat' do i__=1,779 read(1,'(A188)')ar__1 read(ar__1,2) + Name_1(i__),Class(i__),l_z(i__),z(i__),u_z(i__),Comp(i__), + q(i__),NDet(i__),Span(i__),Epoch(i__),X0(i__),e_X0(i__), + Y0(i__),e_Y0(i__),pmX(i__),e_pmX(i__),pmY(i__),e_pmY(i__), + pmr(i__),e_pmr(i__),pmphi(i__),e_pmphi(i__),pmtot(i__), + e_pmtot(i__),betar(i__),e_betar(i__),betaphi(i__), + e_betaphi(i__),betatot(i__),e_betatot(i__) if(ar__1(22:26) .EQ. '') z(i__) = rNULL__ if(ar__1(61:65) .EQ. '') e_X0(i__) = rNULL__ if(ar__1(75:79) .EQ. '') e_Y0(i__) = rNULL__ if(ar__1(88:92) .EQ. '') e_pmX(i__) = rNULL__ if(ar__1(101:105) .EQ. '') e_pmY(i__) = rNULL__ if(ar__1(114:118) .EQ. '') e_pmr(i__) = rNULL__ if(ar__1(127:131) .EQ. '') e_pmphi(i__) = rNULL__ if(ar__1(139:143) .EQ. '') e_pmtot(i__) = rNULL__ if(ar__1(145:151) .EQ. '') betar(i__) = rNULL__ if(ar__1(153:158) .EQ. '') e_betar(i__) = rNULL__ if(ar__1(160:166) .EQ. '') betaphi(i__) = rNULL__ if(ar__1(168:173) .EQ. '') e_betaphi(i__) = rNULL__ if(ar__1(175:181) .EQ. '') betatot(i__) = rNULL__ if(ar__1(183:188) .EQ. '') e_betatot(i__) = rNULL__ c ..............Just test output........... write(6,2) + Name_1(i__),Class(i__),l_z(i__),z(i__),u_z(i__),Comp(i__), + q(i__),NDet(i__),Span(i__),Epoch(i__),X0(i__),e_X0(i__), + Y0(i__),e_Y0(i__),pmX(i__),e_pmX(i__),pmY(i__),e_pmY(i__), + pmr(i__),e_pmr(i__),pmphi(i__),e_pmphi(i__),pmtot(i__), + e_pmtot(i__),betar(i__),e_betar(i__),betaphi(i__), + e_betaphi(i__),betatot(i__),e_betatot(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end