Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/780/92 into FORTRAN code for reading data files line by line.
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 ! Source ID
real*8 RAdeg ! (deg) Right ascension (J2000.0)
real*8 DEdeg ! (deg) Declination (J2000.0)
real*8 BJD ! Barycentric Julian Date (BJD-2400000)
real*4 Jmag ! (mag) 2MASS J magnitude
real*4 e_Jmag ! (mag) J magnitude error
character*1 q_Jmag ! [A] J magnitude quality flag (1)
real*4 Hmag ! (mag) 2MASS H magnitude
real*4 e_Hmag ! (mag) H magnitude error
character*1 q_Hmag ! [AB] H magnitude quality flag (1)
real*4 Kmag ! (mag) 2MASS Ks magnitude
real*4 e_Kmag ! (mag) K magnitude error
character*1 q_Kmag ! [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 ! Modified Julian Date
real*4 mag ! (mag) Unfiltered CSS magnitude
real*4 e_mag ! (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 ! Modified Julian Date
real*4 mag_1 ! (mag) Unfiltered LINEAR magnitude
real*4 e_mag_1 ! (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 ! Modified Julian date of the u observation
real*4 umag ! (mag) ?=99.999 SDSS u magnitude
real*4 e_umag ! (mag) Error of the u magnitude
real*8 MJDg ! Modified Julian date of the g observation
real*4 gmag ! (mag) ?=99.999 SDSS g magnitude
real*4 e_gmag ! (mag) Error of the g magnitude
real*8 MJDr ! Modified Julian date of the r observation
real*4 rmag ! (mag) ?=99.999 SDSS r magnitude
real*4 e_rmag ! (mag) Error of the r magnitude
real*8 MJDi ! Modified Julian date of the i observation
real*4 imag ! (mag) ?=99.999 SDSS i magnitude
real*4 e_imag ! (mag) Error of the i magnitude
real*8 MJDz ! Modified Julian date of the z observation
real*4 zmag ! (mag) ?=99.999 SDSS z magnitude
real*4 e_zmag ! (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 ! Modified Julian date
real*4 W1mag ! (mag) WISE W1 magnitude
real*4 e_W1mag ! (mag) W1 magnitude error
real*4 W2mag ! (mag) WISE W2 magnitude
real*4 e_W2mag ! (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 ! [ugriz JHK W1W2] Photometric band
real*4 lam_eff ! (um) [0.35/4.5] Effective wavelength (micrometer)
real*4 magMin ! (mag) Magnitude at minimum not corrected for ISM
* dust extinction (1)
real*4 e_magMin ! (mag) Magnitude error at minimum
real*4 magMax ! (mag) Magnitude at maximum not corrected for ISM
* dust extinction (1)
real*4 e_magMax ! (mag) Magnitude error at maximum
real*4 Ampl ! (mag) Amplitude
real*4 mag_0_5 ! (mag) Magnitude at 0.5 phase
real*4 Ext ! (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,RAdeg,DEdeg,BJD,Jmag,e_Jmag,q_Jmag,Hmag,e_Hmag,q_Hmag,Kmag,
+ e_Kmag,q_Kmag
c ..............Just test output...........
write(6,1)
+ ID,RAdeg,DEdeg,BJD,Jmag,e_Jmag,q_Jmag,Hmag,e_Hmag,q_Hmag,Kmag,
+ e_Kmag,q_Kmag
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,mag,e_mag
c ..............Just test output...........
write(6,2)MJD,mag,e_mag
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,mag_1,e_mag_1
c ..............Just test output...........
write(6,3)MJD_1,mag_1,e_mag_1
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,umag,e_umag,MJDg,gmag,e_gmag,MJDr,rmag,e_rmag,MJDi,imag,
+ e_imag,MJDz,zmag,e_zmag
c ..............Just test output...........
write(6,4)
+ MJDu,umag,e_umag,MJDg,gmag,e_gmag,MJDr,rmag,e_rmag,MJDi,imag,
+ e_imag,MJDz,zmag,e_zmag
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,W1mag,e_W1mag,W2mag,e_W2mag
c ..............Just test output...........
write(6,5)MJD_2,W1mag,e_W1mag,W2mag,e_W2mag
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,lam_eff,magMin,e_magMin,magMax,e_magMax,Ampl,mag_0_5,Ext
c ..............Just test output...........
write(6,6)
+ Band,lam_eff,magMin,e_magMin,magMax,e_magMax,Ampl,mag_0_5,Ext
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end