Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/660/1428 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-18 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/ApJ/660/1428 Astrometry around RX J0720.4-3125 (Kaplan+, 2007) *================================================================================ *The distance to the isolated neutron star RX J0720.4-3125. * Kaplan D.L., Van Kerkwijk M.H., Anderson J. * <Astrophys. J., 660, 1428-1443 (2007)> * =2007ApJ...660.1428K C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! *Astrometry of reference sources integer*4 nr__ parameter (nr__=23) ! Number of records character*70 ar__ ! Full-size record character*3 ID (nr__) ! Object identifier (G1) integer*4 o_ID (nr__) ! [1/8] Number of epochs for which a source * could be measured real*4 mag (nr__) ! (mag) Instrumental magnitude (2) real*8 DRA (nr__) ! (arcsec) Offset in right ascension (3) integer*4 e_DRA (nr__) ! (10mas) rms uncertainty on DRA real*8 DDE (nr__) ! (arcsec) Offset in declination (3) integer*4 e_DDE (nr__) ! (10mas) rms uncertainty on DDE real*4 pmRA (nr__) ! (mas/yr) Proper motion along RA (4) real*4 e_pmRA (nr__) ! (mas/yr) rms uncertainty on pmRA real*4 pmDE (nr__) ! (mas/yr) Proper motion along DE (4) real*4 e_pmDE (nr__) ! (mas/yr) rms uncertainty on pmDE real*4 Rchi2 (nr__) ! Reduced {chi}^2^ (5) *Note (2): Instrumental magnitudes defined as -2.5log_10_a, where a is the * average amplitude of the effective PSF (ePSF) required to fit a given * star in a single exposure. *Note (3): Position offsets at MJD 52645.1, relative to star A, excluding * parallactic offsets. Reference position: 07:20:24.4837-31:25:51.786 * (J2000.0) *Note (4): Proper motion relative to that of star A. The systematic * uncertainty on the proper motion, from assuming that star A has mu=0, * is ~1mas/yr. *Note (5): Reduced {chi}^2^ values for fitting, with degrees of freedom * 2Nepoch-Npar, where the number of parameters Npar=5 for RX * J0720.4-3125, 0 for A, and 4 for all other stars (note that these * {chi}^2^_red_ values for individual sources are not robust estimates. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! ePSF position measurements for all exposures integer*4 nr__1 parameter (nr__1=1172) ! Number of records character*130 ar__1 ! Full-size record character*3 ID_1 (nr__1) ! Object identifier (G1) character*5 Exp (nr__1) ! Exposure identifier (hrcNN) real*8 DXpix (nr__1) ! (pix) Distortion-corrected X position real*8 e_DXpix (nr__1) ! (pix) Uncertainty in Xpix real*8 DYpix (nr__1) ! (pix) Distortion-corrected Y position real*8 e_DYpix (nr__1) ! (pix) Uncertainty in Ypix real*4 Num (nr__1) ! Number of bad pixel and distance to * nearest bad one (1) real*8 Rchi2_1 (nr__1) ! Reduced {chi}^2^ of the ePSF fit real*8 Xpix (nr__1) ! (pix) Measured X position on chip real*8 e_Xpix (nr__1) ! (pix) Uncertainty in Xpix real*8 Ypix (nr__1) ! (pix) Measured Y position on chip real*8 e_Ypix (nr__1) ! (pix) Uncertainty in Ypix real*8 Amp (nr__1) ! Amplitude of effective PSF (ePSF) real*4 e_Amp (nr__1) ! Uncertainty in Amp real*4 Sky (nr__1) ! Sky level real*4 e_Sky (nr__1) ! Uncertainty in Sky (0 = fixed sky value) *Note (1): This is the sum of the number of bad pixels and the distance from * the initial guess of the centroid to the nearest bad pixel, divided by * 1000. Only measurements with values of 0 (no bad pixels) were used in * our analysis. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Merged positions for each epoch integer*4 nr__2 parameter (nr__2=156) ! Number of records character*56 ar__2 ! Full-size record character*3 ID_2 (nr__2) ! Object identifier (G1) character*4 Exp_1 (nr__2) ! Exposure identifier (hrcN) real*8 DXpix_1 (nr__2) ! (pix) Distortion-corrected X position real*8 e_DXpix_1 (nr__2) ! (pix) Uncertainty in D-Xpix real*8 DYpix_1 (nr__2) ! (pix) Distortion-corrected Y position real*8 e_DYpix_1 (nr__2) ! (pix) Uncertainty in D-Ypix integer*4 Num_1 (nr__2) ! Number of measurements (2x number of exposures) real*4 Rchi2_2 (nr__2) ! Reduced {chi}^2^ (={chi}^2^/(Num-6)) C============================================================================= C Loading file 'table4.dat' ! *Astrometry of reference sources C Format for file interpretation 1 format( + A3,1X,I1,1X,F5.1,1X,F9.5,1X,I2,2X,F9.5,1X,I3,2X,F6.2,1X,F4.2, + 2X,F5.2,1X,F4.2,2X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,23 read(1,'(A70)')ar__ read(ar__,1) + ID(i__),o_ID(i__),mag(i__),DRA(i__),e_DRA(i__),DDE(i__), + e_DDE(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__), + Rchi2(i__) c ..............Just test output........... write(6,1) + ID(i__),o_ID(i__),mag(i__),DRA(i__),e_DRA(i__),DDE(i__), + e_DDE(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),e_pmDE(i__), + Rchi2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! ePSF position measurements for all exposures C Format for file interpretation 2 format( + A3,1X,A5,1X,F10.5,1X,F7.5,1X,F10.5,1X,F7.5,1X,F5.3,1X,F8.3,1X, + F10.5,1X,F7.5,1X,F10.5,1X,F7.5,1X,F9.2,1X,F6.2,1X,F6.3,1X,F5.3) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,1172 read(1,'(A130)')ar__1 read(ar__1,2) + ID_1(i__),Exp(i__),DXpix(i__),e_DXpix(i__),DYpix(i__), + e_DYpix(i__),Num(i__),Rchi2_1(i__),Xpix(i__),e_Xpix(i__), + Ypix(i__),e_Ypix(i__),Amp(i__),e_Amp(i__),Sky(i__),e_Sky(i__) c ..............Just test output........... write(6,2) + ID_1(i__),Exp(i__),DXpix(i__),e_DXpix(i__),DYpix(i__), + e_DYpix(i__),Num(i__),Rchi2_1(i__),Xpix(i__),e_Xpix(i__), + Ypix(i__),e_Ypix(i__),Amp(i__),e_Amp(i__),Sky(i__),e_Sky(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! Merged positions for each epoch C Format for file interpretation 3 format(A3,1X,A4,1X,F10.5,1X,F7.5,1X,F10.5,1X,F7.5,1X,I2,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,156 read(1,'(A56)')ar__2 read(ar__2,3) + ID_2(i__),Exp_1(i__),DXpix_1(i__),e_DXpix_1(i__),DYpix_1(i__), + e_DYpix_1(i__),Num_1(i__),Rchi2_2(i__) c ..............Just test output........... write(6,3) + ID_2(i__),Exp_1(i__),DXpix_1(i__),e_DXpix_1(i__),DYpix_1(i__), + e_DYpix_1(i__),Num_1(i__),Rchi2_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end