Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/151/88 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 2024-Apr-19
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/AJ/151/88 LMC NIR Synoptic Survey. II. Wesenheit relations (Bhardwaj+, 2016)
*================================================================================
*Large Magellanic Cloud near-infrared Synoptic Survey.
*II. The Wesenheit relations and their application to the distance scale.
* Bhardwaj A., Kanbur S.M., Macri L.M., Singh H.P., Ngeow C.-C.,
* Wagner-Kaiser R., Sarajedini A.
* <Astron. J., 151, 88 (2016)>
* =2016AJ....151...88B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Wesenheit magnitudes for Cepheids in the LMC
integer*4 nr__
parameter (nr__=1264) ! Number of records
character*167 ar__ ! Full-size record
character*8 ID ! Star ID
character*2 Type ! Type, Fundamental (FU) or First Overtone (FO)
real*4 logP ! ([d]) [-0.141/1.997] Log of the Period
real*4 WJHmag ! (mag) [8.908/15.478] Wesenheit W_J,H_ index
* (W_J,H_=H-1.63(J-H))
real*4 WJKsmag ! (mag) [9.065/15.598] Wesenheit W_J,Ks_ index
* (W_J,Ks_=K_S_-0.69(J-K_S_))
real*4 WHKsmag ! (mag) [9.181/15.866] Wesenheit W_H,Ks_ index
* (W_H,Ks_=K_S_-1.92(H-K_S_))
real*8 WVJmag ! (mag) [-21.718/15.999]? Wesenheit W_V,J_ index
* (W_V,J_=J-0.41(V-J))
real*4 WVHmag ! (mag) [-5.388/16.075]? Wesenheit W_V,H_ index
* (W_V,H_=H-0.22(V-H))
real*4 WVKsmag ! (mag) [1.990/15.736] Wesenheit W_V,Ks_ index
* (W_V,Ks_=K_S_-0.13(V-K_S_))
real*8 WIJmag ! (mag) [-65.743/16.048]? Wesenheit W_I,J_ index
* (W_I,J_=J-0.92(I-J))
real*8 WIHmag ! (mag) [-23.044/15.671]? Wesenheit W_I,H_ index
* (W_I,H_=H-0.42(I-H))
real*4 WIKsmag ! (mag) [-7.551/15.720]? Wesenheit W_I,Ks_ index
* (W_I,Ks_=K_S_-0.24(I-K_S_))
real*8 WHVImag ! (mag) [-20.931/15.785]? Wesenheit W^H^_V,I_ index
* (W^H^_V,I_=H-0.41(V-I))
real*8 e_WJHmag ! (mag) [0.039/127.098] Uncertainty in WJHmag (1)
real*8 e_WJKsmag ! (mag) [0.016/224.996] Uncertainty in WJKsmag (1)
real*8 e_WHKsmag ! (mag) [0.047/408.082] Uncertainty in WHKsmag (1)
real*4 e_WVJmag ! (mag) [0.023/54.039]? Uncertainty in WVJmag (1)
real*4 e_WVHmag ! (mag) [0.015/51.243]? Uncertainty in WVHmag (1)
real*8 e_WVKsmag ! (mag) [0.012/184.540]Uncertainty in WVKsmag (1)
real*4 e_WIJmag ! (mag) [0.048/67.941]? Uncertainty in WIJmag (1)
real*4 e_WIHmag ! (mag) [0.026/55.316]? Uncertainty in WIHmag (1)
real*8 e_WIKsmag ! (mag) [0.016/188.197]? Uncertainty in WIKsmag (1)
real*4 e_WHVImag ! (mag) [0.032/51.000]? Uncertainty in WHVImag (1)
*Note (1): The uncertainties were calculated by propagating the errors in mean
* magnitudes.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Fourier-fitted mean magnitudes
integer*4 nr__1
parameter (nr__1=113) ! Number of records
character*62 ar__1 ! Full-size record
character*9 ID_1 ! Star ID
character*3 Source ! Source of Light Curve (1)
real*4 Per ! (d) [1.949/68.464] Period
real*4 Jmag ! (mag) [1.773/8.843] Fourier-fitted mean J magnitude,
* 2MASS
real*4 Hmag ! (mag) [1.217/8.235] Fourier-fitted mean H magnitude,
* 2MASS
real*4 Ksmag ! (mag) [1.098/8.021] Fourier-fitted mean K_S_
* magnitude, 2MASS
real*4 e_Jmag ! (mag) [0.022/0.049] Uncertainty in Jmag (2)
real*4 e_Hmag ! (mag) [0.021/0.037] Uncertainty in Hmag (2)
real*4 e_Ksmag ! (mag) [0.022/0.068] Uncertainty in Ksmag (2)
real*4 E_B_V ! (mag) [0.039/1.236] Adopted E(B-V) colour excess (3)
*Note (1): Source as follows:
* MP = Monson & Pierce (2011ApJS..193...12M);
* BTG = Barnes et al. (1997PASP..109..645B);
* LS = Laney & Stobie (1992A&AS...93...93L);
* W = Welch et al. (1984ApJS...54..547W).
*Note (2): The error estimate includes the uncertainties from the Fourier fit
* and the photometry.
*Note (3): The color excess E(B-V) values are taken from Tammann et al.
* (2003, J/A+A/404/423).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Galactic Cepheid Distance Moduli
integer*4 nr__2
parameter (nr__2=116) ! Number of records
character*64 ar__2 ! Full-size record
character*9 ID_2 ! Star ID
real*4 IRSB ! (mag) ?=99.99 [6.68/13.41] Infrared Surface Brightness
* distance modulus (1)
real*4 e_IRSB ! (mag) ?=9.99 [0.02/9.03] Uncertainty in IRSB
real*4 MS ! (mag) ?=99.99 [6.76/12.74] Main-Sequence fitting
* distance modulus (Turner 2010Ap&SS.326..219T)
real*4 e_MS ! (mag) ?=9.99 [0.04/0.31] Uncertainty in MS
real*4 BW ! (mag) ?=99.99 [6.99/13.62] Baade-Wesselink distance
* modulus (Groenewegen 2013, J/A+A/550/A70)
real*4 e_BW ! (mag) ?=9.99 [0.02/9.56] Uncertainty in BW
real*4 Plx ! (mag) ?=99.99 [7.19/11.84] HST Parallax distance modulus
* (2)
real*4 e_Plx ! (mag) ?=9.99 [0.09/0.30] Uncertainty in Plx
real*4 WM ! (mag) [7.03/13.62] Weighted Mean distance modulus
real*4 e_WM ! (mag) [0.09/0.63] Uncertainty in WM
*Note (1): Infrared Surface Brightness (IRSB) method (Fouque et al. 2007,
* J/A+A/476/73; Storm et al. 2011, J/A+A/534/A95).
*Note (2): Hubble Space Telescope parallaxes (HST-{pi}) (Benedict et al.
* 2007AJ....133.1810B; Monson et al. 2012, J/ApJ/759/146; Riess et al.
* 2014ApJ...785..161R)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Calibrated Magnitudes for Fundamental-Mode
Galactic Cepheids
integer*4 nr__3
parameter (nr__3=103) ! Number of records
character*97 ar__3 ! Full-size record
character*9 ID_3 ! Star ID
real*4 logP_1 ! ([d]) [0.432/1.835] Log of the Period
real*4 JMag_1 ! (mag) [-8.082/-3.654] Absolute J magnitude, 2MASS
real*4 HMag_1 ! (mag) [-8.422/-3.919] Absolute H magnitude, 2MASS
real*4 KMag ! (mag) [-8.543/-3.927] Absolute K magnitude, 2MASS
real*4 e_JMag_1 ! (mag) [0.091/0.630] Uncertainty in JMag (1)
real*4 e_HMag_1 ! (mag) [0.091/0.630] Uncertainty in HMag (1)
real*4 e_KMag ! (mag) [0.091/0.630] Uncertainty in KMag (1)
character*1 F1 ! [YN] Flag if Cepheid is used in final
* period-luminosity (P-L) fit (Y=Yes, N=No)
real*4 WJHmag_1 ! (mag) [-8.975/-4.350] Wesenheit W_J,H_ index
* (W_J,H_=H-1.63(J-H))
real*4 WJKmag ! (mag) [-8.861/-4.116] Wesenheit W_J,K_ index
real*4 WHKmag ! (mag) [-8.776/-3.943] Wesenheit W_H,K_ index
real*4 e_WJHmag_1 ! (mag) [0.096/0.631] Uncertainty in WJHmag (1)
real*4 e_WJKmag ! (mag) [0.092/0.630] Uncertainty in WJKmag (1)
real*4 e_WHKmag ! (mag) [0.097/0.631] Uncertainty in WHKmag (1)
character*1 F2 ! [YN] Flag if Cepheid is used in final
* period-Wesenheit (P-W) fit (Y=Yes, N=No)
*Note (1): The uncertainties in absolute magnitudes include the errors in mean
* magnitudes and distance moduli from Table 4 and 5, errors from
* transformations to 2MASS system and reddening corrections. These
* errors are propagated to estimate uncertainty for Wesenheit
* magnitudes.
C=============================================================================
C Loading file 'table2.dat' ! Wesenheit magnitudes for Cepheids in the LMC
C Format for file interpretation
1 format(
+ A8,1X,A2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F6.3,1X,
+ F6.3,1X,F7.3,1X,F7.3,1X,F6.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F7.3,
+ 1X,F6.3,1X,F6.3,1X,F7.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,1264
read(1,'(A167)')ar__
read(ar__,1)
+ ID,Type,logP,WJHmag,WJKsmag,WHKsmag,WVJmag,WVHmag,WVKsmag,
+ WIJmag,WIHmag,WIKsmag,WHVImag,e_WJHmag,e_WJKsmag,e_WHKsmag,
+ e_WVJmag,e_WVHmag,e_WVKsmag,e_WIJmag,e_WIHmag,e_WIKsmag,
+ e_WHVImag
if(ar__(41:47) .EQ. '') WVJmag = rNULL__
if(ar__(49:54) .EQ. '') WVHmag = rNULL__
if(ar__(63:69) .EQ. '') WIJmag = rNULL__
if(ar__(71:77) .EQ. '') WIHmag = rNULL__
if(ar__(79:84) .EQ. '') WIKsmag = rNULL__
if(ar__(86:92) .EQ. '') WHVImag = rNULL__
if(ar__(118:123) .EQ. '') e_WVJmag = rNULL__
if(ar__(125:130) .EQ. '') e_WVHmag = rNULL__
if(ar__(140:145) .EQ. '') e_WIJmag = rNULL__
if(ar__(147:152) .EQ. '') e_WIHmag = rNULL__
if(ar__(154:160) .EQ. '') e_WIKsmag = rNULL__
if(ar__(162:167) .EQ. '') e_WHVImag = rNULL__
c ..............Just test output...........
write(6,1)
+ ID,Type,logP,WJHmag,WJKsmag,WHKsmag,WVJmag,WVHmag,WVKsmag,
+ WIJmag,WIHmag,WIKsmag,WHVImag,e_WJHmag,e_WJKsmag,e_WHKsmag,
+ e_WVJmag,e_WVHmag,e_WVKsmag,e_WIJmag,e_WIHmag,e_WIKsmag,
+ e_WHVImag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Fourier-fitted mean magnitudes
C Format for file interpretation
2 format(
+ A9,1X,A3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,113
read(1,'(A62)')ar__1
read(ar__1,2)
+ ID_1,Source,Per,Jmag,Hmag,Ksmag,e_Jmag,e_Hmag,e_Ksmag,E_B_V
c ..............Just test output...........
write(6,2)
+ ID_1,Source,Per,Jmag,Hmag,Ksmag,e_Jmag,e_Hmag,e_Ksmag,E_B_V
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Galactic Cepheid Distance Moduli
C Format for file interpretation
3 format(
+ A9,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,116
read(1,'(A64)')ar__2
read(ar__2,3)
+ ID_2,IRSB,e_IRSB,MS,e_MS,BW,e_BW,Plx,e_Plx,WM,e_WM
c ..............Just test output...........
write(6,3)
+ ID_2,IRSB,e_IRSB,MS,e_MS,BW,e_BW,Plx,e_Plx,WM,e_WM
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Calibrated Magnitudes for Fundamental-Mode
* Galactic Cepheids
C Format for file interpretation
4 format(
+ A9,1X,F5.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ A1,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,103
read(1,'(A97)')ar__3
read(ar__3,4)
+ ID_3,logP_1,JMag_1,HMag_1,KMag,e_JMag_1,e_HMag_1,e_KMag,F1,
+ WJHmag_1,WJKmag,WHKmag,e_WJHmag_1,e_WJKmag,e_WHKmag,F2
c ..............Just test output...........
write(6,4)
+ ID_3,logP_1,JMag_1,HMag_1,KMag,e_JMag_1,e_HMag_1,e_KMag,F1,
+ WJHmag_1,WJKmag,WHKmag,e_WJHmag_1,e_WJKmag,e_WHKmag,F2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end