Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/780/92 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 2026-Sep-10
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/780/92 Light curves of the RR Lyr SDSS J015450.17+001500.7 (Szabo+, 2014)
*================================================================================
*High-precision 2MASS JHK_s_ light curves and other data for RR Lyrae star
*SDSS J015450+001501: strong constraints for nonlinear pulsation models.
* Szabo R., Ivezic Z., Kiss L.L., Kollath Z., Jones L., Sesar B.,
* Becker A.C., Davenport J.R.A., Cutri R.M.
* <Astrophys. J., 780, 92 (2014)>
* =2014ApJ...780...92S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '2mass.dat' ! 2MASS JHKs photometry (table5)
integer*4 nr__
parameter (nr__=2969) ! Number of records
character*87 ar__ ! Full-size record
C J2000.0 position composed of: RAdeg DEdeg
integer*4 ID (nr__) ! Source ID
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000.0)
real*8 DEdeg (nr__) ! (deg) Declination (J2000.0)
real*8 BJD (nr__) ! Barycentric Julian Date (BJD-2400000)
real*4 Jmag (nr__) ! (mag) 2MASS J magnitude
real*4 e_Jmag (nr__) ! (mag) J magnitude error
character*1 q_Jmag (nr__) ! [A] J magnitude quality flag (1)
real*4 Hmag (nr__) ! (mag) 2MASS H magnitude
real*4 e_Hmag (nr__) ! (mag) H magnitude error
character*1 q_Hmag (nr__) ! [AB] H magnitude quality flag (1)
real*4 Kmag (nr__) ! (mag) 2MASS Ks magnitude
real*4 e_Kmag (nr__) ! (mag) K magnitude error
character*1 q_Kmag (nr__) ! [AB] K magnitude quality flag (1)
*Note (1): Single character flag that is related to the probability (P) that the
* extraction is a valid detection of a near-infrared source on the sky at the
* time of the observation. "A" means P>90%, "B" means 80%<P<90%.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'css.dat' ! *Catalina Sky Survey photometry (table4)
integer*4 nr__1
parameter (nr__1=270) ! Number of records
character*22 ar__1 ! Full-size record
real*8 MJD (nr__1) ! Modified Julian Date
real*4 mag (nr__1) ! (mag) Unfiltered CSS magnitude
real*4 e_mag (nr__1) ! (mag) Magnitude error
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'linear.dat' ! *LINEAR photometry (table3)
integer*4 nr__2
parameter (nr__2=225) ! Number of records
character*25 ar__2 ! Full-size record
real*8 MJD_1 (nr__2) ! Modified Julian Date
real*4 mag_1 (nr__2) ! (mag) Unfiltered LINEAR magnitude
real*4 e_mag_1 (nr__2) ! (mag) Magnitude error
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sdss.dat' ! *SDSS photometry (table2)
integer*4 nr__3
parameter (nr__3=64) ! Number of records
character*129 ar__3 ! Full-size record
real*8 MJDu (nr__3) ! Modified Julian date of the u observation
real*4 umag (nr__3) ! (mag) ?=99.999 SDSS u magnitude
real*4 e_umag (nr__3) ! (mag) Error of the u magnitude
real*8 MJDg (nr__3) ! Modified Julian date of the g observation
real*4 gmag (nr__3) ! (mag) ?=99.999 SDSS g magnitude
real*4 e_gmag (nr__3) ! (mag) Error of the g magnitude
real*8 MJDr (nr__3) ! Modified Julian date of the r observation
real*4 rmag (nr__3) ! (mag) ?=99.999 SDSS r magnitude
real*4 e_rmag (nr__3) ! (mag) Error of the r magnitude
real*8 MJDi (nr__3) ! Modified Julian date of the i observation
real*4 imag (nr__3) ! (mag) ?=99.999 SDSS i magnitude
real*4 e_imag (nr__3) ! (mag) Error of the i magnitude
real*8 MJDz (nr__3) ! Modified Julian date of the z observation
real*4 zmag (nr__3) ! (mag) ?=99.999 SDSS z magnitude
real*4 e_zmag (nr__3) ! (mag) Error of the z magnitude
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'wise.dat' ! WISE photometry extracted by Roc Cutri (IPAC)
(table6)
integer*4 nr__4
parameter (nr__4=43) ! Number of records
character*38 ar__4 ! Full-size record
real*8 MJD_2 (nr__4) ! Modified Julian date
real*4 W1mag (nr__4) ! (mag) WISE W1 magnitude
real*4 e_W1mag (nr__4) ! (mag) W1 magnitude error
real*4 W2mag (nr__4) ! (mag) WISE W2 magnitude
real*4 e_W2mag (nr__4) ! (mag) W2 magnitude error
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sed.dat' ! Spectral energy distribution (table7)
integer*4 nr__5
parameter (nr__5=10) ! Number of records
character*48 ar__5 ! Full-size record
character*2 Band (nr__5) ! [ugriz JHK W1W2] Photometric band
real*4 lam_eff (nr__5) ! (um) [0.35/4.5] Effective wavelength (micrometer)
real*4 magMin (nr__5) ! (mag) Magnitude at minimum not corrected for ISM
* dust extinction (1)
real*4 e_magMin (nr__5) ! (mag) Magnitude error at minimum
real*4 magMax (nr__5) ! (mag) Magnitude at maximum not corrected for ISM
* dust extinction (1)
real*4 e_magMax (nr__5) ! (mag) Magnitude error at maximum
real*4 Ampl (nr__5) ! (mag) Amplitude
real*4 mag_0_5 (nr__5) ! (mag) Magnitude at 0.5 phase
real*4 Ext (nr__5) ! (mag) Extinction (2)
*Note (1): The table lists 2MASS and WISE magnitudes on AB scale.
* The following Vega-to-AB magnitude offsets (mAB=mVega+offset) were used:
* 0.89, 1.37, 1.84 in JHK, and 2.6, 3.3 in W1 and W2.
*Note (2): Plausible extinction corrections are derived using the SFD value for
* the r-band extinction quoted by SDSS for this star: Ar=0.099, and coefficients
* listed in the first row of Table 1 from Berry et al. (2012ApJ...757..166B)
* for the ugrizJHK bands. For the WISE bands, the coefficients taken from Yuan,
* Liu & Xiang (2013MNRAS.430.2188Y) (AW1/Ar=0.084, AW2/Ar=0.074).
C=============================================================================
C Loading file '2mass.dat' ! 2MASS JHKs photometry (table5)
C Format for file interpretation
1 format(
+ I9,1X,F9.6,1X,F8.6,1X,F13.7,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,
+ F5.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'2mass.dat')
write(6,*) '....Loading file: 2mass.dat'
do i__=1,2969
read(1,'(A87)')ar__
read(ar__,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),BJD(i__),Jmag(i__),e_Jmag(i__),
+ q_Jmag(i__),Hmag(i__),e_Hmag(i__),q_Hmag(i__),Kmag(i__),
+ e_Kmag(i__),q_Kmag(i__)
c ..............Just test output...........
write(6,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),BJD(i__),Jmag(i__),e_Jmag(i__),
+ q_Jmag(i__),Hmag(i__),e_Hmag(i__),q_Hmag(i__),Kmag(i__),
+ e_Kmag(i__),q_Kmag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'css.dat' ! *Catalina Sky Survey photometry (table4)
C Format for file interpretation
2 format(F11.5,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'css.dat')
write(6,*) '....Loading file: css.dat'
do i__=1,270
read(1,'(A22)')ar__1
read(ar__1,2)MJD(i__),mag(i__),e_mag(i__)
c ..............Just test output...........
write(6,2)MJD(i__),mag(i__),e_mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'linear.dat' ! *LINEAR photometry (table3)
C Format for file interpretation
3 format(F12.6,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'linear.dat')
write(6,*) '....Loading file: linear.dat'
do i__=1,225
read(1,'(A25)')ar__2
read(ar__2,3)MJD_1(i__),mag_1(i__),e_mag_1(i__)
c ..............Just test output...........
write(6,3)MJD_1(i__),mag_1(i__),e_mag_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'sdss.dat' ! *SDSS photometry (table2)
C Format for file interpretation
4 format(
+ F12.6,1X,F6.3,1X,F5.3,1X,F12.6,1X,F6.3,1X,F5.3,1X,F12.6,1X,
+ F6.3,1X,F5.3,1X,F12.6,1X,F6.3,1X,F5.3,1X,F12.6,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'sdss.dat')
write(6,*) '....Loading file: sdss.dat'
do i__=1,64
read(1,'(A129)')ar__3
read(ar__3,4)
+ MJDu(i__),umag(i__),e_umag(i__),MJDg(i__),gmag(i__),
+ e_gmag(i__),MJDr(i__),rmag(i__),e_rmag(i__),MJDi(i__),
+ imag(i__),e_imag(i__),MJDz(i__),zmag(i__),e_zmag(i__)
c ..............Just test output...........
write(6,4)
+ MJDu(i__),umag(i__),e_umag(i__),MJDg(i__),gmag(i__),
+ e_gmag(i__),MJDr(i__),rmag(i__),e_rmag(i__),MJDi(i__),
+ imag(i__),e_imag(i__),MJDz(i__),zmag(i__),e_zmag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'wise.dat' ! WISE photometry extracted by Roc Cutri (IPAC)
* (table6)
C Format for file interpretation
5 format(F12.6,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'wise.dat')
write(6,*) '....Loading file: wise.dat'
do i__=1,43
read(1,'(A38)')ar__4
read(ar__4,5)
+ MJD_2(i__),W1mag(i__),e_W1mag(i__),W2mag(i__),e_W2mag(i__)
c ..............Just test output...........
write(6,5)
+ MJD_2(i__),W1mag(i__),e_W1mag(i__),W2mag(i__),e_W2mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'sed.dat' ! Spectral energy distribution (table7)
C Format for file interpretation
6 format(
+ A2,1X,F6.4,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'sed.dat')
write(6,*) '....Loading file: sed.dat'
do i__=1,10
read(1,'(A48)')ar__5
read(ar__5,6)
+ Band(i__),lam_eff(i__),magMin(i__),e_magMin(i__),magMax(i__),
+ e_magMax(i__),Ampl(i__),mag_0_5(i__),Ext(i__)
c ..............Just test output...........
write(6,6)
+ Band(i__),lam_eff(i__),magMin(i__),e_magMin(i__),magMax(i__),
+ e_magMax(i__),Ampl(i__),mag_0_5(i__),Ext(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end