Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/849/150 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/849/150 RR Lyrae candidates in the MW outer halo (Cohen+, 2017) *================================================================================ *The outer halo of the Milky Way as probed by RR Lyr variables from the Palomar *Transient Facility. * Cohen J.G., Sesar B., Bahnolzer S., He K., Kulkarni S.R., Prince T.A., * Bellm E., Laher R.R. * <Astrophys. J., 849, 150 (2017)> * =2017ApJ...849..150C C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Light curve parameters for RR Lyr candidates integer*4 nr__ parameter (nr__=116) ! Number of records character*63 ar__ ! Full-size record C J2000 position composed of: RAdeg DEdeg character*1 Sgr (nr__) ! Sgr, a=Probable Sgr Stream real*8 RAdeg (nr__) ! (deg) Right ascension, decimal degrees (J2000) real*8 DEdeg (nr__) ! (deg) Declination, decimal degrees (J2000) real*8 Per (nr__) ! (d) [0.4/0.9] Period real*8 HJD (nr__) ! (d) [55007.9/57538.7] Epoch of maximum light in * heliocentric Julian date; HJD-2400000 (1) real*4 AmpR (nr__) ! (mag) [0.4/1.2] Amplitude of variation in PTF-R band real*4 Rmag (nr__) ! (mag) [0/20.7]? Mean PTF-R band AB magnitude integer*4 N_R (nr__) ! [0/723]? Number of PTF-R band epochs star detected * in late 2016 integer*4 N_g (nr__) ! [0/385]? Number of PTF-g band epochs star detected * in late 2016 integer*4 Qual (nr__) ! [0/3]? Quality (2) *Note (1): Epoch of maximum light in heliocentric Julian date - 2400000 days. * This choice, made for ease of computations, requires 7 digits in the * period for accurate phasing at the present epoch *Note (2): Observed R light curve resembles that of a RR Lyr; code as follows: * 1 = excellent; * 2 = probable; * 3 = uncertain. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Radial velocities for RR Lyr candidates integer*4 nr__1 parameter (nr__1=116) ! Number of records character*54 ar__1 ! Full-size record C J2000 position composed of: RAdeg DEdeg character*1 Sgr_1 (nr__1) ! Sgr, a=Probable Sgr Stream real*8 RAdeg_1 (nr__1) ! (deg) Right ascension, decimal degrees (J2000) real*8 DEdeg_1 (nr__1) ! (deg) Declination, decimal degrees (J2000) real*4 Dist (nr__1) ! (kpc) [50/109.1] Heliocentric distance integer*4 Nsp (nr__1) ! [1/2] Number of spectra real*4 RV_GSR (nr__1) ! (km/s) [-256.2/224] Radial velocity, Galactic system * of rest (3) real*4 e_RV_GSR (nr__1) ! (km/s) [12/29] Uncertainty in RVel character*8 Date (nr__1) ! Approximate RV date of observation MM/YYYY or * MM/DD/YYYY; "dup" = 2 spectra taken character*1 f_Nsp (nr__1) ! Flag on spectra taken (4) *Note (3): RV-GSR is corrected to the systemic velocity, then to the * heliocentric velocity, and then to the Galactic system of rest. *Note (4): Flag on spectra taken -- * d = A DBSP spectrum taken with the Hale Telescope at Palomar Observatory. C============================================================================= C Loading file 'table1.dat' ! Light curve parameters for RR Lyr candidates C Format for file interpretation 1 format( + A1,1X,F10.6,1X,F9.6,1X,F9.7,1X,F8.2,1X,F4.2,1X,F5.2,1X,I3,1X, + I3,2X,I1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,116 read(1,'(A63)')ar__ read(ar__,1) + Sgr(i__),RAdeg(i__),DEdeg(i__),Per(i__),HJD(i__),AmpR(i__), + Rmag(i__),N_R(i__),N_g(i__),Qual(i__) if(ar__(48:52) .EQ. '') Rmag(i__) = rNULL__ if(ar__(54:56) .EQ. '') N_R(i__) = iNULL__ if(ar__(58:60) .EQ. '') N_g(i__) = iNULL__ if(ar__(63:63) .EQ. '') Qual(i__) = iNULL__ c ..............Just test output........... write(6,1) + Sgr(i__),RAdeg(i__),DEdeg(i__),Per(i__),HJD(i__),AmpR(i__), + Rmag(i__),N_R(i__),N_g(i__),Qual(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Radial velocities for RR Lyr candidates C Format for file interpretation 2 format( + A1,1X,F10.6,1X,F9.6,1X,F5.1,1X,I1,1X,F6.1,1X,F4.1,1X,A8,2X,A1) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,116 read(1,'(A54)')ar__1 read(ar__1,2) + Sgr_1(i__),RAdeg_1(i__),DEdeg_1(i__),Dist(i__),Nsp(i__), + RV_GSR(i__),e_RV_GSR(i__),Date(i__),f_Nsp(i__) c ..............Just test output........... write(6,2) + Sgr_1(i__),RAdeg_1(i__),DEdeg_1(i__),Dist(i__),Nsp(i__), + RV_GSR(i__),e_RV_GSR(i__),Date(i__),f_Nsp(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end