Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/470/864 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-17
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/470/864 JHKL photometry of the Galactic Center (Blum+, 1996)
*================================================================================
*JHKL photometry and the K-band luminosity function at the Galactic Center.
* Blum R.D., Sellgren K., Depoy D.L.
* <Astrophys. J. 470, 864 (1996)>
* =1996ApJ...470..864B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Galactic center observed photometry
integer*4 nr__
parameter (nr__=185) ! Number of records
character*89 ar__ ! Full-size record
integer*4 v_BSD96_ (nr__) ! Sequential number
character*1 n__BSD96_ (nr__) ! [*] *: variable star, in table2.dat file
character*11 Name (nr__) ! Name (G1)
real*4 oRA (nr__) ! (arcsec) Offset in right ascension from IRS 7 (2)
real*4 oDE (nr__) ! (arcsec) Offset in declination from IRS 7 (2)
real*4 Kmag (nr__) ! (mag) K magnitude (3)
real*4 e_Kmag (nr__) ! (mag) rms uncertainty on Kmag
real*4 J_K (nr__) ! (mag) ? J-K colour index (3)
real*4 e_J_K (nr__) ! (mag) ? rms uncertainty on J-K
real*4 H_K (nr__) ! (mag) ? H-K colour index (3)
real*4 e_H_K (nr__) ! (mag) ? rms uncertainty on H-K
real*4 K_L (nr__) ! (mag) ? K-L colour index (3)
real*4 e_K_L (nr__) ! (mag) ? rms uncertainty on K-L
character*12 Notes (nr__) ! Comma-separated notes, see file notes.dat
*Note (2): IRS 7: RA = 17 45 39.987 DE = -29 00 22.24 ICRS 2000.0
* (Position from Simbad, ref: 1989IAUS..136..443Y)
*Note (3): All photometry is from our primary OSIRIS data set unless otherwise
* noted. The OSIRIS magnitudes include both measurement and calibration
* uncertainties. The DePoy & Sharp (1991AJ....101.1324D, hereafter DS91)
* data include only measurement uncertainty; calibration uncertainty is
* not included.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Notes to table1.dat
integer*4 nr__1
parameter (nr__1=43) ! Number of records
character*80 ar__1 ! Full-size record
character*4 Note (nr__1) ! Symbol used in table1
character*72 Text (nr__1) ! Text of note
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Variable stars in the Galactic Center
integer*4 nr__2
parameter (nr__2=12) ! Number of records
character*83 ar__2 ! Full-size record
character*7 Name_1 (nr__2) ! Star name (IRS 7, 9 and 12N)
real*4 Jmag (nr__2) ! (mag) ? J magnitude
real*4 e_Jmag (nr__2) ! (mag) ? rms uncertainty on Jmag
real*4 Hmag (nr__2) ! (mag) ? H magnitude
real*4 e_Hmag (nr__2) ! (mag) ? rms uncertainty on Hmag
real*4 Kmag_1 (nr__2) ! (mag) ? K magnitude
real*4 e_Kmag_1 (nr__2) ! (mag) ? rms uncertainty on Kmag
integer*4 Obs_Y (nr__2) ! (yr) Observation date (year)
character*7 Obs_M (nr__2) ! Observation date (month/month)
character*20 Set (nr__2) ! Data set (1)
character*9 Filt (nr__2) ! Filters used to derive J, H, and K magnitudes (2)
*Note (1): Data set:
* Becklin et al., 1978ApJ...220..831B
* DS91: DePoy & Sharp, 1991AJ....101.1324D
* DS91b: DS91 photometry is consistent with J, H, and K magnitudes of
* IRS 7 being the same as for Becklin et al., 1978ApJ...220..831B;
* see text.
*Note (2): J, H, and K are broad-band filters.
* KND is broad-band K with a 1% neutral density filter.
* NBK is K magnitude derived from two narrow-band filters
* ({Delta}{lambda}/{lambda}=1%) near 2.2mm; see text.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Galactic Center reddening-corrected photometry
integer*4 nr__3
parameter (nr__3=159) ! Number of records
character*81 ar__3 ! Full-size record
integer*4 v_BSD96__1 (nr__3) ! Sequential number
character*11 Name_2 (nr__3) ! Name (G1)
character*1 n_Name (nr__3) ! [*] Individual note on IRS 10EL (1)
real*4 K0mag (nr__3) ! (mag) Dereddened K magnitude (2)
real*4 e_K0mag (nr__3) ! (mag) rms uncertainty on K0mag
real*4 AK (nr__3) ! (mag) Absorption in K band (3)
real*4 e_AK (nr__3) ! (mag) ? rms uncertainty on AK
real*4 v_J_K_0 (nr__3) ! (mag) ? Dereddened J-K colour index (2)
real*4 e__J_K_0 (nr__3) ! (mag) ? rms uncertainty on (J-K)0
real*4 v_H_K_0 (nr__3) ! (mag) ? Dereddened H-K colour index (2)
real*4 e__H_K_0 (nr__3) ! (mag) ? rms uncertainty on (H-K)0
real*4 v_K_L_0 (nr__3) ! (mag) ? Dereddened K-L colour index (2)
real*4 e__K_L_0 (nr__3) ! (mag) ? rms uncertainty on (K-L)0
character*10 Notes_1 (nr__3) ! Notes (4)
*Note (1): See notes for Table 1. Since this star is variable and K and L were
* observed at different times, it is possible that the observed K-L is
* too red and A_K_ overestimated. The source is probably quite red,
* however, since it was not detected at H or J.
*Note (2): Includes photometric uncertainty and uncertainty in derived A_K_.
* K0 results from average of DS91 and OSIRIS K when these differed by
* less than 0.2mag.
*Note (3): Includes photometric uncertainty only. No uncertainty given for
* the case of mean reddening (A_K_=3.58) or A_K_ = lower limit; see
* notes c and d. For the case of two or more observed colors, only
* observed colors that fell within +0.5 mag of the reddening line were
* used to determine A_K_; see text.
*Note (4): Notes:
* c: A_K_ is an estimated lower limit resulting from H-K, where H is
* taken as the limiting H magnitude (14.25); see text.
* d: Mean A_K_ adopted from stars with measured H-K only; see text.
C=============================================================================
C Loading file 'table1.dat' ! Galactic center observed photometry
C Format for file interpretation
1 format(
+ I3,A1,1X,A11,1X,F6.2,2X,F6.2,2X,F5.2,1X,F4.2,2X,F4.2,1X,F4.2,
+ 2X,F4.2,1X,F4.2,2X,F4.2,1X,F4.2,1X,A12)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,185
read(1,'(A89)')ar__
read(ar__,1)
+ v_BSD96_(i__),n__BSD96_(i__),Name(i__),oRA(i__),oDE(i__),
+ Kmag(i__),e_Kmag(i__),J_K(i__),e_J_K(i__),H_K(i__),e_H_K(i__),
+ K_L(i__),e_K_L(i__),Notes(i__)
if(ar__(46:49) .EQ. '') J_K(i__) = rNULL__
if(ar__(51:54) .EQ. '') e_J_K(i__) = rNULL__
if(ar__(57:60) .EQ. '') H_K(i__) = rNULL__
if(ar__(62:65) .EQ. '') e_H_K(i__) = rNULL__
if(ar__(68:71) .EQ. '') K_L(i__) = rNULL__
if(ar__(73:76) .EQ. '') e_K_L(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ v_BSD96_(i__),n__BSD96_(i__),Name(i__),oRA(i__),oDE(i__),
+ Kmag(i__),e_Kmag(i__),J_K(i__),e_J_K(i__),H_K(i__),e_H_K(i__),
+ K_L(i__),e_K_L(i__),Notes(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Notes to table1.dat
C Format for file interpretation
2 format(A4,4X,A72)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,43
read(1,'(A80)')ar__1
read(ar__1,2)Note(i__),Text(i__)
c ..............Just test output...........
write(6,2)Note(i__),Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Variable stars in the Galactic Center
C Format for file interpretation
3 format(
+ A7,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I4,1X,
+ A7,1X,A20,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,12
read(1,'(A83)')ar__2
read(ar__2,3)
+ Name_1(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),
+ Kmag_1(i__),e_Kmag_1(i__),Obs_Y(i__),Obs_M(i__),Set(i__),
+ Filt(i__)
if(ar__2(9:13) .EQ. '') Jmag(i__) = rNULL__
if(ar__2(15:18) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__2(20:24) .EQ. '') Hmag(i__) = rNULL__
if(ar__2(26:29) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__2(31:34) .EQ. '') Kmag_1(i__) = rNULL__
if(ar__2(36:39) .EQ. '') e_Kmag_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_1(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),e_Hmag(i__),
+ Kmag_1(i__),e_Kmag_1(i__),Obs_Y(i__),Obs_M(i__),Set(i__),
+ Filt(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Galactic Center reddening-corrected photometry
C Format for file interpretation
4 format(
+ I3,1X,A11,A1,1X,F4.2,1X,F4.2,2X,F4.2,1X,F4.2,2X,F4.2,1X,F4.2,
+ 2X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A10)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,159
read(1,'(A81)')ar__3
read(ar__3,4)
+ v_BSD96__1(i__),Name_2(i__),n_Name(i__),K0mag(i__),
+ e_K0mag(i__),AK(i__),e_AK(i__),v_J_K_0(i__),e__J_K_0(i__),
+ v_H_K_0(i__),e__H_K_0(i__),v_K_L_0(i__),e__K_L_0(i__),
+ Notes_1(i__)
if(ar__3(34:37) .EQ. '') e_AK(i__) = rNULL__
if(ar__3(40:43) .EQ. '') v_J_K_0(i__) = rNULL__
if(ar__3(45:48) .EQ. '') e__J_K_0(i__) = rNULL__
if(ar__3(51:54) .EQ. '') v_H_K_0(i__) = rNULL__
if(ar__3(56:59) .EQ. '') e__H_K_0(i__) = rNULL__
if(ar__3(61:65) .EQ. '') v_K_L_0(i__) = rNULL__
if(ar__3(67:70) .EQ. '') e__K_L_0(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ v_BSD96__1(i__),Name_2(i__),n_Name(i__),K0mag(i__),
+ e_K0mag(i__),AK(i__),e_AK(i__),v_J_K_0(i__),e__J_K_0(i__),
+ v_H_K_0(i__),e__H_K_0(i__),v_K_L_0(i__),e__K_L_0(i__),
+ Notes_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end