Conversion of standardized ReadMe file for
file /./ftp/cats/J/PASJ/58/389 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-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/PASJ/58/389 O/Fe abundances of 15 RR Lyrae (Takeda+, 2006)
*================================================================================
*On the spectroscopic determination of atmospheric parameters and O/Fe
*abundances of RR Lyrae stars.
* Takeda Y., Honda S., Aoki W., Takada-hidai M., Zhao G., Chen Y.-Q.,
* Shi J.-R.
* <Publ. Astron. Soc. Jap., 58, 389-406 (2006)>
* =2006PASJ...58..389T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ewidlist.dat' ! Fe line data and the measured equivalent widths
integer*4 nr__
parameter (nr__=220) ! Number of records
character*139 ar__ ! Full-size record
integer*4 Line (nr__) ! Line number (7000- for Fe I 8000- for Fe II)
real*4 Code (nr__) ! Species code (26.00: Fe I, 26.01: Fe II)
real*8 Lambda (nr__) ! (0.1nm) Wavelength
real*4 Ep (nr__) ! (eV) Lower excitation potential
real*4 loggf (nr__) ! log gf value
real*4 Gammar (nr__) ! Radiation damping parameter
real*4 Gammas (nr__) ! Quadratic Stark effect damping parameter
real*4 Gammaw (nr__) ! van der Waals effect damping parameter
real*4 EW_rrlyr1 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum rrlyr-1
real*4 EW_rrlyr2 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum rrlyr-2
real*4 EW_rrlyr3 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum rrlyr-3
real*4 EW_rrlyr4 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum rrlyr-4
real*4 EW_rrlyr5 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum rrlyr-5
real*4 EW_dxdel1 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum dxdel-1
real*4 EW_dxdel2 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum dxdel-2
real*4 EW_dxdel3 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum dxdel-3
real*4 EW_dhpeg1 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum dhpeg-1
real*4 EW_dhpeg2 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum dhpeg-2
real*4 EW_dhpeg3 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum dhpeg-3
real*4 EW_vyser1 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum vyser-1
real*4 EW_vyser2 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum vyser-2
real*4 EW_vyser3 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum vyser-3
real*4 EW_vyser4 (nr__) ! (0.1pm) ?=0.0 Equivalent width for spectrum vyser-4
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabe2/*' ! EW data for TT Lyr, DX Del, DH Peg, VY Ser and Sun
integer*4 nr__1
parameter (nr__1=16) ! Number of records
character*55 ar__1 ! Full-size record
integer*4 Line_1 (nr__1) ! Line number
real*4 Code_1 (nr__1) ! FeI: 26.00 FeII: 26.01
real*8 Lambda_1 (nr__1) ! (0.1nm) Wavelength
real*8 Ep_1 (nr__1) ! (eV) Lower excitation potential (chi)
real*8 loggf_1 (nr__1) ! Adopted log(gf) value
real*4 EW (nr__1) ! (0.1pm) Equivalent width
real*8 Abund1 (nr__1) ! Fe abundance (sign inversed for Fe II lines)
real*8 Abund2 (nr__1) ! Fe abundance (sign inversed for Fe I lines)
C=============================================================================
C Loading file 'ewidlist.dat' ! Fe line data and the measured equivalent widths
C Format for file interpretation
1 format(
+ I4,F6.2,F9.3,F6.2,F6.2,F6.2,F6.2,F6.2,F6.1,F6.1,F6.1,F6.1,
+ F6.1,F6.1,F6.1,F6.1,F6.1,F6.1,F6.1,F6.1,F6.1,F6.1,F6.1)
C Effective file loading
open(unit=1,status='old',file=
+'ewidlist.dat')
write(6,*) '....Loading file: ewidlist.dat'
do i__=1,220
read(1,'(A139)')ar__
read(ar__,1)
+ Line(i__),Code(i__),Lambda(i__),Ep(i__),loggf(i__),
+ Gammar(i__),Gammas(i__),Gammaw(i__),EW_rrlyr1(i__),
+ EW_rrlyr2(i__),EW_rrlyr3(i__),EW_rrlyr4(i__),EW_rrlyr5(i__),
+ EW_dxdel1(i__),EW_dxdel2(i__),EW_dxdel3(i__),EW_dhpeg1(i__),
+ EW_dhpeg2(i__),EW_dhpeg3(i__),EW_vyser1(i__),EW_vyser2(i__),
+ EW_vyser3(i__),EW_vyser4(i__)
c ..............Just test output...........
write(6,1)
+ Line(i__),Code(i__),Lambda(i__),Ep(i__),loggf(i__),
+ Gammar(i__),Gammas(i__),Gammaw(i__),EW_rrlyr1(i__),
+ EW_rrlyr2(i__),EW_rrlyr3(i__),EW_rrlyr4(i__),EW_rrlyr5(i__),
+ EW_dxdel1(i__),EW_dxdel2(i__),EW_dxdel3(i__),EW_dhpeg1(i__),
+ EW_dhpeg2(i__),EW_dhpeg3(i__),EW_vyser1(i__),EW_vyser2(i__),
+ EW_vyser3(i__),EW_vyser4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabe2/*' ! EW data for TT Lyr, DX Del, DH Peg, VY Ser and Sun
C Format for file interpretation
2 format(I4,F6.2,F9.3,F7.3,F7.3,F6.1,F8.3,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'tabe2/*')
write(6,*) '....Loading file: tabe2/*'
do i__=1,16
read(1,'(A55)')ar__1
read(ar__1,2)
+ Line_1(i__),Code_1(i__),Lambda_1(i__),Ep_1(i__),loggf_1(i__),
+ EW(i__),Abund1(i__),Abund2(i__)
c ..............Just test output...........
write(6,2)
+ Line_1(i__),Code_1(i__),Lambda_1(i__),Ep_1(i__),loggf_1(i__),
+ EW(i__),Abund1(i__),Abund2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end