Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/655/1046 into FORTRAN code for reading data files line by line.
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/ApJ/655/1046 Abundances of {phi} Her (Zavala+, 2007)
*================================================================================
*The mercury-manganese binary star {phi} Herculis: detection and properties
*of the secondary and revision of the elemental abundances of the primary.
* Zavala R.T., Adelman S.J., Hummel C.A., Gulliver A.F., Caliskan H.,
* Armstrong J.T., Hutter D.J., Johnston K.J., Pauls T.A.
* <Astrophys. J., 655, 1046-1057 (2007)>
* =2007ApJ...655.1046Z
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Astrometric solutions for {phi} Her
integer*4 nr__
parameter (nr__=25) ! Number of records
character*70 ar__ ! Full-size record
character*3 ObsM ! UT Month of NPOI observation
integer*4 ObsD ! (d) UT Day of NPOI observation
real*8 ObsY ! (yr) Julian Year of NPOI observation
real*4 rho ! (mas) Fitted binary separation
real*4 theta ! (deg) Fitted binary position angle
real*4 errMaj ! (mas) Semimajor axis of error ellipse
real*4 errMin ! (mas) Semiminor axis of error ellipse
real*4 errPA ! (deg) Position angle of error ellipse
real*4 Crho ! (mas) Calculated binary separation
real*4 Ctheta ! (deg) Calculated binary position angle
real*4 v_O_C_rho ! (mas) Binary separation (O-C) value
real*4 v_O_C_theta ! (deg) Binary position angle (O-C) value
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Abundances of {phi} Her including new lines
integer*4 nr__1
parameter (nr__1=196) ! Number of records
character*36 ar__1 ! Full-size record
character*6 Ion ! Species name (1)
character*3 Mult ! Muliplet identification (2)
real*8 lambda ! (0.1nm) ? Laboratory wavelength in Angstroms
real*4 log_gf ! ([-]) ? Log of the gf value
character*2 r_log_gf ! Reference for log(gf) (3)
integer*4 EW ! (0.1pm) ? Observed equivalent width in milliAngstroms
real*4 logN_NT ! ([-]) ? Log of the derived number abundance
*Note (1): The adopted abundance from HgII{lambda}3984 is from
* Woolf & Lambert (1999ApJ...521..414W).
*Note (2): The multiplet numbers are from Moore (1945, A Multiplet Table of
* Astrophysical Interest (Princeton: Princeton Univ. Obs.) except that I
* indicates Mn II lines from Iglesias & Velasco (1964, Publ. Inst. Opt.
* Madrid, No. 23) and J indicates Fe II lines from Johansson (1978,
* Phys. Scr., 18, 217 ).
*Note (3): The gf value references (including those for older lines):
* BG = Biemont et al. (1989A&A...209..391B for V II;
* Biemont (1981ApJ...248..867B) for Zr II
* DW = Dworetsky (1980A&A....84..350D)
* GB = Grevesse et al. (1981, in Proc. 23rd Liege Astrophys. Colloq.,
* Upper Main Sequence Stars , 211)
* HL = Hannaford et al. (1982ApJ...261..736H)
* KX = Kurucz & Bell (1995, CD-ROM 23, Atomic Data for Opacity Calculations)
* LA = Lanz & Artru (1985, Phys. Scr., 32, 115)
* LD = Lawler & Dakin (1989, J. Opt. Soc. Am. B, 6, 1457)
* MF = Fuhr et al. (1988) and Martin et al. (1988), Cat. <VI/72>
* MC = Magazzu & Cowley (1986ApJ...308..254M)
* N4 = Fuhr & Wiese (2006, J. Phys. Chem. Ref. Data, 35, 1669)
* WF = Wiese, Fuhr & Deters (1996, J. Phys. Chem. Ref. Data Monogr. Vol. 7)
* WM = Wiese & Martin (1980, Wavelengths and Transition Probabilities for
* Atoms and Atomic Ions, Part 2)
* WS = Wiese, Smith & Glennon (1966, Atomic Transition Probabilities,
* Hydrogen through Neon (NSRDS-NBS 4) and Wiese, Smith & Miles
* (1969, Atomic Transition Probabilities, Vol. 2, Sodium through
* Calcium, NSRDS-NBS 22).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Comparison of derived and solar abundances (log N/H)
integer*4 nr__2
parameter (nr__2=33) ! Number of records
character*39 ar__2 ! Full-size record
character*6 Ion_1 ! Species name (1)
real*4 AAbd ! ([-]) ? Adelman et al. (2001, Cat. <J/A+A/367/597>)
* {phi} Her A abundance
real*4 e_AAbd ! ([-]) ? Uncertainty in AAbd
real*4 ATP ! ([-]) This paper's {phi} Her A abundance
real*4 e_ATP ! ([-]) ? Uncertainty in ATP (2)
integer*4 n_ATP ! ? Number of lines used in ATP
real*4 Sun ! ([-]) Solar abundance (3)
*Note (1): The Hg II value adopted for this study is that of Woolf & Lambert
* (1999ApJ...521..414W), the Nd and Pr abundances are from Dolk et al.
* (2002A&A...385..111D).
*Note (2): When more than one line is detected the error is the standard
* deviation of the mean. No error is quoted when only one line is present.
*Note (3): From Grevesse, Noels & Sauval (1996, in ASP Conf. Ser. 99, 117).
C=============================================================================
C Loading file 'table2.dat' ! Astrometric solutions for {phi} Her
C Format for file interpretation
1 format(
+ A3,1X,I2,1X,F9.4,1X,F5.2,1X,F6.2,1X,F4.2,1X,F4.2,1X,F5.1,1X,
+ F5.2,1X,F6.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,25
read(1,'(A70)')ar__
read(ar__,1)
+ ObsM,ObsD,ObsY,rho,theta,errMaj,errMin,errPA,Crho,Ctheta,
+ v_O_C_rho,v_O_C_theta
c ..............Just test output...........
write(6,1)
+ ObsM,ObsD,ObsY,rho,theta,errMaj,errMin,errPA,Crho,Ctheta,
+ v_O_C_rho,v_O_C_theta
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Abundances of {phi} Her including new lines
C Format for file interpretation
2 format(A6,A3,1X,F8.3,F5.2,1X,A2,1X,I3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,196
read(1,'(A36)')ar__1
read(ar__1,2)Ion,Mult,lambda,log_gf,r_log_gf,EW,logN_NT
if(ar__1(11:18) .EQ. '') lambda = rNULL__
if(ar__1(19:23) .EQ. '') log_gf = rNULL__
if(ar__1(28:30) .EQ. '') EW = iNULL__
if(ar__1(32:36) .EQ. '') logN_NT = rNULL__
c ..............Just test output...........
write(6,2)Ion,Mult,lambda,log_gf,r_log_gf,EW,logN_NT
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Comparison of derived and solar abundances (log N/H)
C Format for file interpretation
3 format(A6,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,I3,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,33
read(1,'(A39)')ar__2
read(ar__2,3)Ion_1,AAbd,e_AAbd,ATP,e_ATP,n_ATP,Sun
if(ar__2(8:12) .EQ. '') AAbd = rNULL__
if(ar__2(14:17) .EQ. '') e_AAbd = rNULL__
if(ar__2(25:28) .EQ. '') e_ATP = rNULL__
if(ar__2(30:32) .EQ. '') n_ATP = iNULL__
c ..............Just test output...........
write(6,3)Ion_1,AAbd,e_AAbd,ATP,e_ATP,n_ATP,Sun
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end