Conversion of standardized ReadMe file for
file /./ftp/cats/J/AZh/84/714 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/AZh/84/714 Optical spectral variability of BD+48 1220 (Klochkova+, 2007)
*================================================================================
*Evolutionary stage of the spectral variable BD+48 1220=IRAS 05040+4820.
* Klochkova V.G., Chentsov E.L., Tavolzhanskaya M.S., Panchuk V.E.
* <Astron. Zh. 84, 714 (2007)>
* =2007AZh....84..714K
* =2007ARep...51..642K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Elemental abundances in the atmosphere of
BD+48 1220 (Teff=7900K, logg=0.0, xi_t_6.0)
integer*4 nr__
parameter (nr__=22) ! Number of records
character*38 ar__ ! Full-size record
character*2 E (nr__) ! Element
real*4 logeps_E (nr__) ! ([-]) Solar abundance for E element (1)
character*4 Ion (nr__) ! Ion identification
real*4 logeps_X (nr__) ! ([-]) Star abundance for X element (1)
real*4 e_logeps_X (nr__) ! ([-]) ? rms uncertainty on logeps(X)
integer*4 o_logeps_X (nr__) ! Number of lines used
real*4 v_X_Fe_ (nr__) ! ([Sun]) Solar [X/Fe] abundance value (2)
*Note (1): For logeps(H)=12.0
*Note (2): Solar abundances are taken from Lodder, 2003ApJ...591.1220
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Residual line intensities in fractions of
continuum and HRV in spectra of BD+48 1220
obtained on various dates
integer*4 nr__1
parameter (nr__1=550) ! Number of records
character*77 ar__1 ! Full-size record
character*10 El (nr__1) ! Element with multiplet number (1)
real*8 Lambda (nr__1) ! (0.1nm) ? Wavelength expressed in Angstroem
character*1 u_Lambda (nr__1) ! [:] Uncertainty flag on Lambda
real*4 Res1 (nr__1) ! ? Residual line intensity for observation
* on 2002-12-02
character*1 u_Res1 (nr__1) ! [:] Uncertainty flag on Res1
real*4 HRV1 (nr__1) ! (km/s) ? Heliocentric radial velocity for observation
* on 2002-12-02
character*1 u_HRV1 (nr__1) ! [:] Uncertainty flag on HRV1
real*4 Res2 (nr__1) ! ? Residual line intensity for observation
* on 2003-09-09
character*1 u_Res2 (nr__1) ! [:] Uncertainty flag on Res2
real*4 HRV2 (nr__1) ! (km/s) ? Heliocentric radial velocity for observation
* on 2003-09-09
character*1 u_HRV2 (nr__1) ! [:] Uncertainty flag on HRV2
real*4 Res3 (nr__1) ! ? Residual line intensity for observation
* on 2004-01-10
character*1 u_Res3 (nr__1) ! [:] Uncertainty flag on Res3
real*4 HRV3 (nr__1) ! (km/s) ? Heliocentric radial velocity for observation
* on 2004-01-10
character*1 u_HRV3 (nr__1) ! [:] Uncertainty flag on HRV3
real*4 Res4 (nr__1) ! ? Residual line intensity for observation
* on 2004-03-08
character*1 u_Res4 (nr__1) ! [:] Uncertainty flag on Res4
real*4 HRV4 (nr__1) ! (km/s) ? Heliocentric radial velocity for observation
* on 2004-03-08
character*1 u_HRV4 (nr__1) ! [:] Uncertainty flag on HRV4
character*2 n_HRV4 (nr__1) ! [em] emission
*Note (1): DIB = diffuse interstellar band
* Table 2: Diffuse bands in the spectrum of BD+48 1220
* -------------------
* lambda W_lambda_
* {AA} m{AA}
* -------------------
* 5747 42
* 5780 227
* 5797 80
* 6195 25
* 6203 64
* 6270 64
* 6379 54
* 6425 13
* 6613 104
* 6699 17
* -------------------
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Atomic parameters of lines in spectrum of
BD+48 1220 obtained on 2004-04-10, and
corresponding atmospheric elemental abundances
integer*4 nr__2
parameter (nr__2=199) ! Number of records
character*46 ar__2 ! Full-size record
character*4 Ion_1 (nr__2) ! Element
real*8 Lambda_1 (nr__2) ! (0.1nm) Wavelength expressed in Angstroem
real*8 EP (nr__2) ! (eV) Excitation potentials of the lower level
real*4 loggf (nr__2) ! ([-]) Oscillator strength (log gf)
integer*4 EW (nr__2) ! (0.1nm) Equivalent width
real*4 log_eps (nr__2) ! ([-]) Abundance of the element (1)
*Note (1): For log(eps(H)=12 in atmosphere of BD+48 1220 obtained with
* parameters Teff=7900K,logg=0.5, {xi}_t_=6.0. --------------------------------------------------------------------------------
C=============================================================================
C Loading file 'table4.dat' ! Elemental abundances in the atmosphere of
* BD+48 1220 (Teff=7900K, logg=0.0, xi_t_6.0)
C Format for file interpretation
1 format(A2,1X,F5.2,1X,A4,3X,F5.2,2X,F4.2,2X,I2,2X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,22
read(1,'(A38)')ar__
read(ar__,1)
+ E(i__),logeps_E(i__),Ion(i__),logeps_X(i__),e_logeps_X(i__),
+ o_logeps_X(i__),v_X_Fe_(i__)
if(ar__(24:27) .EQ. '') e_logeps_X(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ E(i__),logeps_E(i__),Ion(i__),logeps_X(i__),e_logeps_X(i__),
+ o_logeps_X(i__),v_X_Fe_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Residual line intensities in fractions of
* continuum and HRV in spectra of BD+48 1220
* obtained on various dates
C Format for file interpretation
2 format(
+ A10,1X,F7.2,A1,2X,F4.2,A1,1X,F5.1,A1,2X,F4.2,A1,1X,F5.1,A1,2X,
+ F4.2,A1,1X,F5.1,A1,2X,F4.2,A1,1X,F5.1,A1,A2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,550
read(1,'(A77)')ar__1
read(ar__1,2)
+ El(i__),Lambda(i__),u_Lambda(i__),Res1(i__),u_Res1(i__),
+ HRV1(i__),u_HRV1(i__),Res2(i__),u_Res2(i__),HRV2(i__),
+ u_HRV2(i__),Res3(i__),u_Res3(i__),HRV3(i__),u_HRV3(i__),
+ Res4(i__),u_Res4(i__),HRV4(i__),u_HRV4(i__),n_HRV4(i__)
if(ar__1(12:18) .EQ. '') Lambda(i__) = rNULL__
if(ar__1(22:25) .EQ. '') Res1(i__) = rNULL__
if(ar__1(28:32) .EQ. '') HRV1(i__) = rNULL__
if(ar__1(36:39) .EQ. '') Res2(i__) = rNULL__
if(ar__1(42:46) .EQ. '') HRV2(i__) = rNULL__
if(ar__1(50:53) .EQ. '') Res3(i__) = rNULL__
if(ar__1(56:60) .EQ. '') HRV3(i__) = rNULL__
if(ar__1(64:67) .EQ. '') Res4(i__) = rNULL__
if(ar__1(70:74) .EQ. '') HRV4(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ El(i__),Lambda(i__),u_Lambda(i__),Res1(i__),u_Res1(i__),
+ HRV1(i__),u_HRV1(i__),Res2(i__),u_Res2(i__),HRV2(i__),
+ u_HRV2(i__),Res3(i__),u_Res3(i__),HRV3(i__),u_HRV3(i__),
+ Res4(i__),u_Res4(i__),HRV4(i__),u_HRV4(i__),n_HRV4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Atomic parameters of lines in spectrum of
* BD+48 1220 obtained on 2004-04-10, and
* corresponding atmospheric elemental abundances
C Format for file interpretation
3 format(A4,3X,F8.3,2X,F7.4,3X,F6.3,2X,I3,2X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,199
read(1,'(A46)')ar__2
read(ar__2,3)
+ Ion_1(i__),Lambda_1(i__),EP(i__),loggf(i__),EW(i__),
+ log_eps(i__)
c ..............Just test output...........
write(6,3)
+ Ion_1(i__),Lambda_1(i__),EP(i__),loggf(i__),EW(i__),
+ log_eps(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end