Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/121/2032 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/AJ/121/2032 Direct variables in M33A field (Mochejska+, 2001)
*================================================================================
*DIRECT distances to nearby galaxies using detached eclipsing binaries and
*Cepheids. VII. Additional variables in the field M33A discovered with
*image subtraction.
* Mochejska B.J., Kaluzny J., Stanek K.Z., Sasselov D.D., Szentgyorgyi A.H.
* <Astron. J. 121, 2032 (2001)>
* =2001AJ....121.2032M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! DIRECT eclipsing binaries in M33A
integer*4 nr__
parameter (nr__=53) ! Number of records
character*64 ar__ ! Full-size record
character*18 D33 (nr__) ! Binary name, JHHMMSS.s+DDMMSS.s (G1)
real*8 Per (nr__) ! (d) Period of the system
real*4 Vmag (nr__) ! (mag) ? Maximum V band magnitude outside of eclipse
real*4 Bmag (nr__) ! (mag) ? Maximum B band magnitude outside of eclipse
real*4 R1 (nr__) ! Estimated primary component radius in units
* of the orbital separation
real*4 R2 (nr__) ! Estimated secondary component radius in units
* of the orbital separation
real*4 Inc (nr__) ! (deg) Estimated inclination angle of the orbit
real*4 Eccen (nr__) ! Estimated eccentricity of the orbit
integer*4 Com (nr__) ! ? Comments (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! DIRECT flux eclipsing binaries in M33A
integer*4 nr__1
parameter (nr__1=10) ! Number of records
character*30 ar__1 ! Full-size record
character*18 D33_1 (nr__1) ! Name, JHHMMSS.s+DDMMSS.s (G1)
real*8 Per_1 (nr__1) ! (d) Period
integer*4 Com_1 (nr__1) ! ? Comment (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Light curves of eclipsing binaries in M33A
integer*4 nr__2
parameter (nr__2=14935) ! Number of records
character*55 ar__2 ! Full-size record
character*18 D33_2 (nr__2) ! Object name, JHHMMSS.s+DDMMSS.s (G1)
character*1 Filter (nr__2) ! Filter passband (1)
real*8 HJD (nr__2) ! (d) Julian date of observation
real*8 Data (nr__2) ! (mag) Observed magnitude or differential flux
real*8 e_Data (nr__2) ! (mag) Uncertainty in Data
*Note (1):
* B = Indicates Data is the Johnson B passband magnitude;
* V = Indicates Data is the Johnson V passband magnitude;
* b = Indicates Data is the Johnson B passband differential flux;
* v = Indicates Data is the Johnson V passband differential flux.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! DIRECT Cepheids in M33A
integer*4 nr__3
parameter (nr__3=242) ! Number of records
character*60 ar__3 ! Full-size record
character*18 D33_3 (nr__3) ! Cepheid or object name, JHHMMSS.s+DDMMSS.s (G1)
real*4 Per_2 (nr__3) ! (d) Period of the system
real*4 Vmag_1 (nr__3) ! (mag) ? Flux-weighted average V band magnitude
real*4 Bmag_1 (nr__3) ! (mag) ? Flux-weighted average B band magnitude
real*4 Vamp (nr__3) ! (mag) ? The V band amplitude
real*4 Bamp (nr__3) ! (mag) ? The B band amplitude
character*9 Com_2 (nr__3) ! Comments (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! DIRECT flux cepheids in M33A
integer*4 nr__4
parameter (nr__4=27) ! Number of records
character*30 ar__4 ! Full-size record
character*18 D33_4 (nr__4) ! Name, JHHMMSS.s+DDMMSS.s (G1)
real*8 Per_3 (nr__4) ! (d) Period
integer*4 Com_3 (nr__4) ! ? Comment (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Light curves of Cepheids in M33A
integer*4 nr__5
parameter (nr__5=71491) ! Number of records
character*55 ar__5 ! Full-size record
character*18 D33_5 (nr__5) ! Object name, JHHMMSS.s+DDMMSS.s (G1)
character*1 Filter_1 (nr__5) ! Filter passband (1)
real*8 HJD_1 (nr__5) ! (d) Julian date of observation
real*8 Data_1 (nr__5) ! (mag) Observed magnitude or differential flux
real*8 e_Data_1 (nr__5) ! (mag) Uncertainty in Data
*Note (1):
* B = Indicates Data is the Johnson B passband magnitude;
* V = Indicates Data is the Johnson V passband magnitude;
* b = Indicates Data is the Johnson B passband differential flux;
* v = Indicates Data is the Johnson V passband differential flux.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! DIRECT other periodic variables in M33A
integer*4 nr__6
parameter (nr__6=13) ! Number of records
character*52 ar__6 ! Full-size record
character*18 D33_6 (nr__6) ! Variable name, JHHMMSS.s+DDMMSS.s (G1)
real*4 Per_4 (nr__6) ! (d) Period
real*4 Vmag_2 (nr__6) ! (mag) ? Maximum V magnitude
real*4 Bmag_2 (nr__6) ! (mag) ? Maximum B magnitude
real*4 v_Vmag_ (nr__6) ! (mag) ? Flux-weighted average magnitudes
real*4 v_Bmag_ (nr__6) ! (mag) ? Flux-weighted average magnitudes
character*2 Com_4 (nr__6) ! Comments (EB = Eclipsing binary)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Light curves of other periodic variables in M33A
integer*4 nr__7
parameter (nr__7=3725) ! Number of records
character*55 ar__7 ! Full-size record
character*18 D33_7 (nr__7) ! Object name, JHHMMSS.s+DDMMSS.s (G1)
character*1 Filter_2 (nr__7) ! Filter passband (1)
real*8 HJD_2 (nr__7) ! (d) Julian date of observation
real*8 Data_2 (nr__7) ! (mag) Observed magnitude or differential flux
real*8 e_Data_2 (nr__7) ! (mag) Uncertainty in Data
*Note (1):
* B = Indicates Data is the Johnson B passband magnitude;
* V = Indicates Data is the Johnson V passband magnitude;
* b = Indicates Data is the Johnson B passband differential flux;
* v = Indicates Data is the Johnson V passband differential flux.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! DIRECT miscellaneous variables in M33A
integer*4 nr__8
parameter (nr__8=50) ! Number of records
character*45 ar__8 ! Full-size record
character*18 D33_8 (nr__8) ! Cepheid or object name, JHHMMSS.s+DDMMSS.s (G1)
real*4 Vmag_3 (nr__8) ! (mag) ? The average V band magnitude
real*4 Bmag_3 (nr__8) ! (mag) ? The average B band magnitude
real*4 Vamp_1 (nr__8) ! (mag) ? The V band amplitude
real*4 Bamp_1 (nr__8) ! (mag) ? The B band amplitude
integer*4 Com_5 (nr__8) ! ? Comment (G2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Light curves of miscellaneous variables in M33A
integer*4 nr__9
parameter (nr__9=10958) ! Number of records
character*55 ar__9 ! Full-size record
character*18 D33_9 (nr__9) ! Object name, JHHMMSS.s+DDMMSS.s (G1)
character*1 Filter_3 (nr__9) ! Filter passband (1)
real*8 HJD_3 (nr__9) ! (d) Julian date of observation
real*8 Data_3 (nr__9) ! (mag) Observed magnitude or differential flux
real*8 e_Data_3 (nr__9) ! (mag) Uncertainty in Data
*Note (1):
* B = Indicates Data is the Johnson B passband magnitude;
* V = Indicates Data is the Johnson V passband magnitude;
* b = Indicates Data is the Johnson B passband differential flux;
* v = Indicates Data is the Johnson V passband differential flux.
C=============================================================================
C Loading file 'table2.dat' ! DIRECT eclipsing binaries in M33A
C Format for file interpretation
1 format(
+ 3X,A18,1X,F7.4,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,53
read(1,'(A64)')ar__
read(ar__,1)
+ D33(i__),Per(i__),Vmag(i__),Bmag(i__),R1(i__),R2(i__),
+ Inc(i__),Eccen(i__),Com(i__)
if(ar__(31:35) .EQ. '') Vmag(i__) = rNULL__
if(ar__(37:41) .EQ. '') Bmag(i__) = rNULL__
if(ar__(64:64) .EQ. '') Com(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ D33(i__),Per(i__),Vmag(i__),Bmag(i__),R1(i__),R2(i__),
+ Inc(i__),Eccen(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! DIRECT flux eclipsing binaries in M33A
C Format for file interpretation
2 format(3X,A18,1X,F7.4,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,10
read(1,'(A30)')ar__1
read(ar__1,2)D33_1(i__),Per_1(i__),Com_1(i__)
if(ar__1(30:30) .EQ. '') Com_1(i__) = iNULL__
c ..............Just test output...........
write(6,2)D33_1(i__),Per_1(i__),Com_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Light curves of eclipsing binaries in M33A
C Format for file interpretation
3 format(3X,A18,1X,A1,1X,F12.4,1X,F10.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,14935
read(1,'(A55)')ar__2
read(ar__2,3)
+ D33_2(i__),Filter(i__),HJD(i__),Data(i__),e_Data(i__)
c ..............Just test output...........
write(6,3)
+ D33_2(i__),Filter(i__),HJD(i__),Data(i__),e_Data(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! DIRECT Cepheids in M33A
C Format for file interpretation
4 format(3X,A18,1X,F6.3,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,242
read(1,'(A60)')ar__3
read(ar__3,4)
+ D33_3(i__),Per_2(i__),Vmag_1(i__),Bmag_1(i__),Vamp(i__),
+ Bamp(i__),Com_2(i__)
if(ar__3(30:34) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__3(36:40) .EQ. '') Bmag_1(i__) = rNULL__
if(ar__3(42:45) .EQ. '') Vamp(i__) = rNULL__
if(ar__3(47:50) .EQ. '') Bamp(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ D33_3(i__),Per_2(i__),Vmag_1(i__),Bmag_1(i__),Vamp(i__),
+ Bamp(i__),Com_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! DIRECT flux cepheids in M33A
C Format for file interpretation
5 format(3X,A18,1X,F7.4,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,27
read(1,'(A30)')ar__4
read(ar__4,5)D33_4(i__),Per_3(i__),Com_3(i__)
if(ar__4(30:30) .EQ. '') Com_3(i__) = iNULL__
c ..............Just test output...........
write(6,5)D33_4(i__),Per_3(i__),Com_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Light curves of Cepheids in M33A
C Format for file interpretation
6 format(3X,A18,1X,A1,1X,F12.4,1X,F10.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,71491
read(1,'(A55)')ar__5
read(ar__5,6)
+ D33_5(i__),Filter_1(i__),HJD_1(i__),Data_1(i__),e_Data_1(i__)
c ..............Just test output...........
write(6,6)
+ D33_5(i__),Filter_1(i__),HJD_1(i__),Data_1(i__),e_Data_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! DIRECT other periodic variables in M33A
C Format for file interpretation
7 format(3X,A18,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,13
read(1,'(A52)')ar__6
read(ar__6,7)
+ D33_6(i__),Per_4(i__),Vmag_2(i__),Bmag_2(i__),v_Vmag_(i__),
+ v_Bmag_(i__),Com_4(i__)
if(ar__6(29:33) .EQ. '') Vmag_2(i__) = rNULL__
if(ar__6(35:39) .EQ. '') Bmag_2(i__) = rNULL__
if(ar__6(41:44) .EQ. '') v_Vmag_(i__) = rNULL__
if(ar__6(46:49) .EQ. '') v_Bmag_(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ D33_6(i__),Per_4(i__),Vmag_2(i__),Bmag_2(i__),v_Vmag_(i__),
+ v_Bmag_(i__),Com_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Light curves of other periodic variables in M33A
C Format for file interpretation
8 format(3X,A18,1X,A1,1X,F12.4,1X,F10.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,3725
read(1,'(A55)')ar__7
read(ar__7,8)
+ D33_7(i__),Filter_2(i__),HJD_2(i__),Data_2(i__),e_Data_2(i__)
c ..............Just test output...........
write(6,8)
+ D33_7(i__),Filter_2(i__),HJD_2(i__),Data_2(i__),e_Data_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! DIRECT miscellaneous variables in M33A
C Format for file interpretation
9 format(3X,A18,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,50
read(1,'(A45)')ar__8
read(ar__8,9)
+ D33_8(i__),Vmag_3(i__),Bmag_3(i__),Vamp_1(i__),Bamp_1(i__),
+ Com_5(i__)
if(ar__8(23:27) .EQ. '') Vmag_3(i__) = rNULL__
if(ar__8(29:33) .EQ. '') Bmag_3(i__) = rNULL__
if(ar__8(35:38) .EQ. '') Vamp_1(i__) = rNULL__
if(ar__8(40:43) .EQ. '') Bamp_1(i__) = rNULL__
if(ar__8(45:45) .EQ. '') Com_5(i__) = iNULL__
c ..............Just test output...........
write(6,9)
+ D33_8(i__),Vmag_3(i__),Bmag_3(i__),Vamp_1(i__),Bamp_1(i__),
+ Com_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Light curves of miscellaneous variables in M33A
C Format for file interpretation
10 format(3X,A18,1X,A1,1X,F12.4,1X,F10.3,1X,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,10958
read(1,'(A55)')ar__9
read(ar__9,10)
+ D33_9(i__),Filter_3(i__),HJD_3(i__),Data_3(i__),e_Data_3(i__)
c ..............Just test output...........
write(6,10)
+ D33_9(i__),Filter_3(i__),HJD_3(i__),Data_3(i__),e_Data_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end