Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/608/A100 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/A+A/608/A100 CMD and mass distribution of Ba stars (Escorza+, 2017)
*================================================================================
*Hertzsprung-Russell diagram and mass distribution of barium stars.
* Escorza A., Boffin H.M.J., Jorissen A., Van Eck S., Siess L.,
* Van Winckel H., Karinkuzhi D., Shetye S., Pourbaix D.
* <Astron. Astrophys. 608, A100 (2017)>
* =2017A&A...608A.100E (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Effective temperatures and luminosities obtained
for a sample of barium and related stars
integer*4 nr__
parameter (nr__=437) ! Number of records
character*262 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*14 Name (nr__) ! Name of the star
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
character*5 Type (nr__) ! Star type
real*8 Per (nr__) ! (d) ? Orbital period (1)
real*4 e (nr__) ! ? Orbital eccentricity (1)
real*4 v_Fe_H_ (nr__) ! ([Sun]) ? Metallicity with respect to solar
character*22 r__Fe_H_ (nr__) ! Reference for [Fe/H] (2)
real*4 plx (nr__) ! (mas) Parallax
real*4 e_plx (nr__) ! (mas) Error on the parallax
character*29 r_plx (nr__) ! Reference for the parallax (3)
integer*4 b_T_0_0 (nr__) ! (K) Lower limit of Teff for grid with [Fe/H]=0.0
integer*4 T_0_0 (nr__) ! (K) Effective temperature for grid with [Fe/H]=0.0
integer*4 B_T_0_0_1 (nr__) ! (K) Upper limit of Teff for grid with [Fe/H]=0.0
real*8 b_L_0_0 (nr__) ! (Lsun) Lower limit of luminosity for grid
* with [Fe/H]=0.0
real*8 L_0_0 (nr__) ! (Lsun) Luminosity for grid with [Fe/H]=0.0
real*8 B_L_0_0_1 (nr__) ! (Lsun) Upper limit of luminosity for grid
* with [Fe/H]=0.0
integer*4 b_T__0_25 (nr__) ! (K) Lower limit of Teff for grid with [Fe/H]=-0.25
integer*4 T__0_25 (nr__) ! (K) Effective temperature for grid
* with [Fe/H]=-0.25
integer*4 B_T__0_25_1(nr__) ! (K) Upper limit of Teff for grid with [Fe/H]=-0.25
real*8 b_L__0_25 (nr__) ! (Lsun) Lower limit of luminosity for grid
* with [Fe/H]=-0.25
real*8 L__0_25 (nr__) ! (Lsun) Luminosity for grid with [Fe/H]=-0.25
real*8 B_L__0_25_1(nr__) ! (Lsun) Upper limit of luminosity for grid
* with [Fe/H]=-0.25
integer*4 b_T__0_5 (nr__) ! (K) Lower limit of Teff for grid with [Fe/H]=-0.5
integer*4 T__0_5 (nr__) ! (K) Effective temperature for grid
* with [Fe/H]=-0.5
integer*4 B_T__0_5_1 (nr__) ! (K) Upper limit of Teff for grid with [Fe/H]=-0.5
real*8 b_L__0_5 (nr__) ! (Lsun) Lower limit of luminosity for grid
* with [Fe/H]=-0.5
real*8 L__0_5 (nr__) ! (Lsun) Luminosity for grid with [Fe/H]=-0.5
real*8 B_L__0_5_1 (nr__) ! (Lsun) Upper limit of luminosity for grid
* with [Fe/H]=-0.5
*Note (1): Orbital periods and eccentricities are determined by Jorissen et al.
* (1998A&A...332..877J) and Jorissen et al. (2017, in preparation) for the giants
* and by Escorza et al. (2017, in preparation) for the main sequence
* and subgiant stars.
*Note (2): Karinkuzhi for Karinkuzhi et al. (2017, in prep).
*Note (3): References as follows:
* Hipparcos (ESA, 1997) = HIP, 1997HIP...C......0E, Cat. I/239
* Pourbaix & Jorissen (2000) = Pourbaix & Jorissen, 2000A&AS..145..161P
* TGAS (Lindegren et al., 2016) = TGAS, 2016A&A...595A...1G, Cat. I/337
C=============================================================================
C Loading file 'tablea1.dat' ! Effective temperatures and luminosities obtained
* for a sample of barium and related stars
C Format for file interpretation
1 format(
+ A14,1X,F12.8,1X,F12.8,1X,A5,1X,F8.2,1X,F5.3,1X,F5.2,1X,A22,1X,
+ F5.1,1X,F3.1,1X,A29,1X,I4,1X,I4,1X,I4,1X,F8.3,1X,F9.3,1X,F9.3,
+ 1X,I4,1X,I4,1X,I4,1X,F8.3,1X,F9.3,1X,F9.3,1X,I4,1X,I4,1X,I4,
+ 1X,F8.3,1X,F9.3,1X,F9.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,437
read(1,'(A262)')ar__
read(ar__,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),Type(i__),Per(i__),e(i__),
+ v_Fe_H_(i__),r__Fe_H_(i__),plx(i__),e_plx(i__),r_plx(i__),
+ b_T_0_0(i__),T_0_0(i__),B_T_0_0_1(i__),b_L_0_0(i__),
+ L_0_0(i__),B_L_0_0_1(i__),b_T__0_25(i__),T__0_25(i__),
+ B_T__0_25_1(i__),b_L__0_25(i__),L__0_25(i__),B_L__0_25_1(i__),
+ b_T__0_5(i__),T__0_5(i__),B_T__0_5_1(i__),b_L__0_5(i__),
+ L__0_5(i__),B_L__0_5_1(i__)
if(ar__(48:55) .EQ. '') Per(i__) = rNULL__
if(ar__(57:61) .EQ. '') e(i__) = rNULL__
if(ar__(63:67) .EQ. '') v_Fe_H_(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),Type(i__),Per(i__),e(i__),
+ v_Fe_H_(i__),r__Fe_H_(i__),plx(i__),e_plx(i__),r_plx(i__),
+ b_T_0_0(i__),T_0_0(i__),B_T_0_0_1(i__),b_L_0_0(i__),
+ L_0_0(i__),B_L_0_0_1(i__),b_T__0_25(i__),T__0_25(i__),
+ B_T__0_25_1(i__),b_L__0_25(i__),L__0_25(i__),B_L__0_25_1(i__),
+ b_T__0_5(i__),T__0_5(i__),B_T__0_5_1(i__),b_L__0_5(i__),
+ L__0_5(i__),B_L__0_5_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end