Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/AJ/129/433 Abundances and equivalent widths of 14 Cepheids (Kovtyukh+, 2005)
*================================================================================
*Phase-dependent variation of the fundamental parameters of cepheids.
*II. Periods longer than 10 days.
* Kovtyukh V.V., Andrievsky S.M., Belik S.I., Luck R.E.
* <Astron. J., 129, 433-453 (2005)>
* =2005AJ....129..433K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Observing log
integer*4 nr__
parameter (nr__=186) ! Number of records
character*64 ar__ ! Full-size record
character*6 Name (nr__) ! Cepheid name
real*8 Per (nr__) ! (d) Period
character*7 Tel (nr__) ! Telescope
character*10 ObsDate (nr__) ! ("YYYY-MM-DD") Date of the observation (YYYY-MM-DD)
real*8 JD (nr__) ! (d) Julian Date of observation
real*4 Phase (nr__) ! Phase of Cepheid
integer*4 ExpTime (nr__) ! (s) ? Exposure time
integer*4 S_N (nr__) ! Mean signal-to-noise ratio
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Line equivalent widths
integer*4 nr__1
parameter (nr__1=122016) ! Number of records
character*55 ar__1 ! Full-size record
character*6 Name_1 (nr__1) ! Cepheid name
character*6 Ion (nr__1) ! Ion
real*8 lambda (nr__1) ! (0.1nm) Wavelength
real*4 Z (nr__1) ! Atomic number
real*4 EP (nr__1) ! (eV) Lower excitation potential
real*8 loggf (nr__1) ! Oscillator strength
real*4 Phase_1 (nr__1) ! Phase in Cepheid
integer*4 EW (nr__1) ! (pm) ? Line equivalent width at Phase
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Per-phase abundances
integer*4 nr__2
parameter (nr__2=6138) ! Number of records
character*33 ar__2 ! Full-size record
real*4 Z_1 (nr__2) ! Atomic number (1)
character*6 Name_2 (nr__2) ! Cepheid name
real*4 Phase_2 (nr__2) ! Phase
real*4 Abund (nr__2) ! ([Sun]) ?=-9.99 Abundance of element Z
real*4 e_Abund (nr__2) ! ([Sun]) ?=0 Standard deviation of Abundance
integer*4 o_Abund (nr__2) ! ?=0 Number of lines for the species
*Note (1): The first decimal is 0 for the neutral atom, and 1 for the ionized
* atom (e.g. 26.0 = Fe I and 26.1 = Fe II).
C=============================================================================
C Loading file 'table5.dat' ! Observing log
C Format for file interpretation
1 format(A6,1X,F11.8,1X,A7,1X,A10,1X,F11.3,1X,F5.3,1X,I4,1X,I3)
C Effective file loading
open(unit=1,file='table5.dat', status='old')
write(6,*) '....Loading file: table5.dat'
do i__=1,186
read(1,'(A64)')ar__
read(ar__,1)
+ Name(i__),Per(i__),Tel(i__),ObsDate(i__),JD(i__),Phase(i__),
+ ExpTime(i__),S_N(i__)
if(ar__(57:60) .EQ. '') ExpTime(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),Per(i__),Tel(i__),ObsDate(i__),JD(i__),Phase(i__),
+ ExpTime(i__),S_N(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Line equivalent widths
C Format for file interpretation
2 format(A6,1X,A6,1X,F8.3,1X,F4.1,3X,F6.3,1X,F7.3,1X,F5.3,1X,I4)
C Effective file loading
open(unit=1,file='table6.dat', status='old')
write(6,*) '....Loading file: table6.dat'
do i__=1,122016
read(1,'(A55)')ar__1
read(ar__1,2)
+ Name_1(i__),Ion(i__),lambda(i__),Z(i__),EP(i__),loggf(i__),
+ Phase_1(i__),EW(i__)
if(ar__1(52:55) .EQ. '') EW(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Ion(i__),lambda(i__),Z(i__),EP(i__),loggf(i__),
+ Phase_1(i__),EW(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Per-phase abundances
C Format for file interpretation
3 format(F4.1,1X,A6,2X,F5.3,1X,F5.2,1X,F4.2,1X,I3)
C Effective file loading
open(unit=1,file='table7.dat', status='old')
write(6,*) '....Loading file: table7.dat'
do i__=1,6138
read(1,'(A33)')ar__2
read(ar__2,3)
+ Z_1(i__),Name_2(i__),Phase_2(i__),Abund(i__),e_Abund(i__),
+ o_Abund(i__)
c ..............Just test output...........
write(6,3)
+ Z_1(i__),Name_2(i__),Phase_2(i__),Abund(i__),e_Abund(i__),
+ o_Abund(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end