Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/700/1299 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-16
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/700/1299 Gas-phase element depletions in the ISM (Jenkins, 2009)
*================================================================================
*A unified representation of gas-phase element depletions in the interstellar
*medium.
* Jenkins E.B.
* <Astrophys. J., 700, 1299-1348 (2009)>
* =2009ApJ...700.1299J
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Stellar data and line of sight information
integer*4 nr__
parameter (nr__=276) ! Number of records
character*150 ar__ ! Full-size record
C Position composed of: GLON GLAT
integer*4 HD ! ? HD number
integer*4 Vcomp ! (km/s) ? Specific Heliocentric velocity component
character*16 Name ! Name (1)
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
real*4 Vmag ! (mag) The Simbad V band magnitude (not used for
* calculating the reddenings or Dist)
character*16 SpType ! MK spectral type
character*6 r_SpType ! Reference for SpType; see refs.dat file
real*4 E_B_V ! (mag) The (B-V) extinction
real*4 Dist ! (kpc) Heliocentric distance (1)
real*4 Z ! (kpc) Distance from Galactic plane
real*4 b_logNHI ! ([cm-2]) ? Lower limit on logNHI
real*4 logNHI ! ([cm-2]) ? Log best observed HI column density
character*1 n_logNHI ! Flag on logNHI (see notes.dat file)
real*4 B_logNHI_1 ! ([cm-2]) ? Upper limit on logNHI
character*6 r_logNHI ! Reference for logNHI; see refs.dat file
real*4 b_logNH2 ! ([cm-2]) ? Lower limit on logNH2
real*4 logNH2 ! ([cm-2]) ? Log best observed H_2_ column density
real*4 B_logNH2_1 ! ([cm-2]) ? Upper limit on logNH2
character*6 r_logNH2 ! Reference for logNH2; see refs.dat file
*Note (1): Stars for which neither the Hipparcos parallaxes (Cat. I/311)
* nor the determinations by Bowen et al. 2008, Cat. J/ApJS/176/59 were
* applicable had their spectroscopic parallaxes determined with the
* standard formula for the distance d (in kpc):
* logd=(m_V_-A_V_-M_V_-10)/5, where m_V_ is the apparent V magnitude of
* the star, A_V_ is the extinction by dust in the visible (assumed to be
* 3.1E(B-V)), and M_V_ is the absolute V magnitude of the star.
* See section 4.5.1 for further details.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Element depletion parameters
integer*4 nr__1
parameter (nr__1=16) ! Number of records
character*86 ar__1 ! Full-size record
character*2 El ! Element identification "X" (1)
real*4 v_X_H_ ! ([-]) Adopted reference abundance relative to H (2)
real*4 e__X_H_ ! ([-]) Uncertainty in X/H (2)
real*4 Ax ! The A_X_ parameter; defined in Eq. 10 (3)
real*4 e_Ax ! Uncertainty in Ax
real*4 Bx ! The B_X_ parameter; defined in Eq. 10 (3)
real*4 e_Bx ! Uncertainty in Bx (4)
real*4 zx ! The z_X_ parameter; defined in Eq. 11 (3)
real*4 Xgas_H0 ! ([-]) The [X_gas_/H]_0_ parameter; defined in Eq. 13
real*4 e_Xgas_H0 ! ([-]) Uncertainty in Xgas/H0; defined in Eq. 14 (3)(4)
real*4 Xgas_H1 ! ([-]) The [X_gas_/H]_1_ parameter; defined in Eq. 15
real*4 e_Xgas_H1 ! ([-]) Uncertainty in Xgas/H1; defined in Eq. 16 (4)
real*4 chi2 ! The {chi}^2^ value
integer*4 dof ! Degrees of freedom
real*4 Prob ! Probability of worse fit
*Note (1): Coefficients for S do not appear in this table because a
* nonstandard approach was required. The coefficients are given
* in the text of Section 9.
*Note (2): On a logarithmic scale with H=12. Values and their errors
* taken from the recommended solar abundances of Lodders K.
* (2003ApJ...591.1220L).
*Note (3): The slope of a best-fit line for the observed depletions
* [X_gas_/H]_obs_=A_X_(F*-z_X_)+B_X_ for any given element X against
* F_*_ (the line-of-sight depletion strength factor). The parameter A_X_
* represents the propensity of that element to increase (the absolute
* value of) its particular depletion level as F_*_ becomes larger. For
* even the smallest observed values of F_*_, most elements still show
* some depletion (see Figure 1). We denote this level of depletion as
* [X_gas_/H]_0_. See section 3 for further details.
*Note (4): Unlike the convention for listing errors in the fit outcomes in
* Table 7, the uncertainties with the terms listed here include both the
* formal errors of the fit coefficients and the error in the adopted
* value of (X/H)_{sun}_, added together in quadrature.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Observed and synthetic log N(H) and F_*_
integer*4 nr__2
parameter (nr__2=239) ! Number of records
character*130 ar__2 ! Full-size record
integer*4 HD_1 ! ? HD number
integer*4 Vcomp_1 ! (km/s) ? Specific Heliocentric velocity component (1)
character*16 Name_1 ! Name
real*4 b_logNHo ! ([cm-2]) ? Lower limit on logNHOb
real*4 logNHo ! ([cm-2]) ? Log observed H column density
real*4 B_logNHo_1 ! ([cm-2]) ? Upper limit on logNHOb
real*4 F_o ! ? Observed line-of-sight depletion strength
* factor
real*4 e_F_o ! ? Uncertainty in F*o
real*4 logNHs ! ([cm-2]) ? Log synthetic H column density
real*4 e_logNHs ! ([cm-2]) ? Uncertainty in logNHs
real*4 F_s ! ? Synthetic line-of-sight depletion strength
* factor
real*4 e_F_s ! ? Uncertainty in F*s
real*4 Prob_1 ! ? Probability of a worse fit
character*41 El_1 ! Elements considered in fit
*Note (1): Separate velocity components, explicitly identified and not
* grouped together; see section 4.3.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Synthetic log N(H) and F_*_ for stars in the
Local Bubble
integer*4 nr__3
parameter (nr__3=29) ! Number of records
character*84 ar__3 ! Full-size record
C Position composed of: GLON GLAT
character*8 WD ! White Dwarf identification (HHMM+DDd; B1950)
character*15 AName ! Alternate name
real*4 GLON_1 ! (deg) Galactic longitude
real*4 GLAT_1 ! (deg) Galactic latitude
real*4 Vmag_1 ! (mag) The Simbad V band magnitude
integer*4 Dist_1 ! (pc) Heliocentric distance (1)
real*4 logNHs_1 ! ([-]) Log synthetic H column density
real*4 e_logNHs_1 ! ([-]) Uncertainty in logNHs
real*4 F_s_1 ! Synthetic line-of-sight depletion strength factor
real*4 e_F_s_1 ! Uncertainty in F*s
real*4 chi2_1 ! ? The {chi}^2^ value
character*9 El_2 ! Element(s) considered (2)
*Note (1): Taken from the papers that described the FUSE observations.
*Note (2): All column density data from Lehner et al. (2003ApJ...595..858L),
* except for HZ 21 (Oliveira et al. 2003ApJ...587..235O).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Observations and fits for Carbon - Krypton
integer*4 nr__4
parameter (nr__4=1475) ! Number of records
character*108 ar__4 ! Full-size record
character*2 El_3 ! Element "X"
integer*4 HD_2 ! ? HD number
integer*4 Vcomp_2 ! (km/s) ? Specific Heliocentric velocity component
character*16 Name_2 ! Name (1)
real*4 b_logNx ! ([cm-2]) ? Lower limit on logNx
real*4 logNx ! ([cm-2]) ? Log best observed column density (2)
character*1 f_logNx ! [if] infinite or fixed value of logNx (3)
real*4 B_logNx_1 ! ([cm-2]) ? Upper limit on logNx
character*1 n_logNx ! [if] infinite or fixed value of B_logNx (3)
character*6 r_logNx ! Reference codes for logNx (see refs.dat file)
real*4 F_ ! ? Line-of-sight depletion strength factor
real*4 e_F_ ! ? Uncertainty in F*
real*4 v_X_H__1 ! ([-]) ? Log expected depletion [X_gas_/H] (4)
real*4 e__X_H__1 ! ([-]) ? Uncertainty in Elgas/H
real*4 Res ! ? Observed [X/H] minus computed from Eq. 10
real*4 sigma ! ? Standard deviation (5)
*Note (1): In table 7, {omicron}^2^ Sco is a misprint for {omicron}^1^ Sco;
* corrected at CDS.
*Note (2): Column densities for CII, NI, OI, MgII, SiII, PII, SII, ClII, TiII,
* CrII, MnII, FeII, NiII, CuII, ZnII, GeII, and KrI.
*Note (3): Flag as follows:
* i = indicates value is infinite;
* f = The entry 13.44 for the adopted value of N (NiII) stated by CLMS06
* (Cartledge et al. 2006ApJ...641..327C) is incorrect; this number should
* be 13.94, as confirmed by Cartledge (private communication). The value
* listed here includes a correction for a change in the transition
* f -value by -0.106dex - see Table 1.
*Note (4): Computed using Eq. 10. The listed errors do not include an overall
* systematic uncertainty of
* C/H = 0.04,
* N/H = 0.11,
* O/H = 0.05,
* Mg/H = 0.02,
* Si/H = 0.02,
* P/H = 0.04,
* S/H = 0.04,
* Cl/H = 0.06,
* Ti/H = 0.03,
* Cr/H = 0.05,
* Mn/H = 0.03,
* Fe/H = 0.03,
* Ni/H = 0.03,
* Cu/H = 0.06,
* Zn/H = 0.04,
* Ge/H = 0.05, and
* Kr/H = 0.08
* in the solar abundance in order to show just the formal error that
* arises from the uncertainties in location of line of best fit and
* the value of F*.
*Note (5): The difference shown in the previous column divided by the
* expected difference due to the uncertainties in both the measured
* column density and the coefficients that appear in Eq. 10.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__5
parameter (nr__5=138) ! Number of records
character*69 ar__5 ! Full-size record
character*6 Ref ! Code of reference
character*19 BibCode ! Bibcode
character*21 Auth ! Author's name(s)
character*19 Comm ! Comment
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Notes of table 2
integer*4 nr__6
parameter (nr__6=44) ! Number of records
character*76 ar__6 ! Full-size record
character*1 F ! [d-r] Flag code
character*74 Note ! Text of note
C=============================================================================
C Loading file 'table2.dat' ! Stellar data and line of sight information
C Format for file interpretation
1 format(
+ I6,4X,I3,2X,A16,1X,F6.2,1X,F6.2,1X,F5.2,1X,A16,1X,A6,1X,F5.2,
+ 1X,F5.3,1X,F5.2,1X,F5.2,3X,F5.2,1X,A1,1X,F5.2,3X,A6,2X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,276
read(1,'(A150)')ar__
read(ar__,1)
+ HD,Vcomp,Name,GLON,GLAT,Vmag,SpType,r_SpType,E_B_V,Dist,Z,
+ b_logNHI,logNHI,n_logNHI,B_logNHI_1,r_logNHI,b_logNH2,logNH2,
+ B_logNH2_1,r_logNH2
if(ar__(1:6) .EQ. '') HD = iNULL__
if(ar__(11:13) .EQ. '') Vcomp = iNULL__
if(ar__(95:99) .EQ. '') b_logNHI = rNULL__
if(ar__(103:107) .EQ. '') logNHI = rNULL__
if(ar__(111:115) .EQ. '') B_logNHI_1 = rNULL__
if(ar__(127:131) .EQ. '') b_logNH2 = rNULL__
if(ar__(133:137) .EQ. '') logNH2 = rNULL__
if(ar__(139:143) .EQ. '') B_logNH2_1 = rNULL__
c ..............Just test output...........
write(6,1)
+ HD,Vcomp,Name,GLON,GLAT,Vmag,SpType,r_SpType,E_B_V,Dist,Z,
+ b_logNHI,logNHI,n_logNHI,B_logNHI_1,r_logNHI,b_logNH2,logNH2,
+ B_logNH2_1,r_logNH2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Element depletion parameters
C Format for file interpretation
2 format(
+ A2,1X,F4.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.1,1X,I3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,16
read(1,'(A86)')ar__1
read(ar__1,2)
+ El,v_X_H_,e__X_H_,Ax,e_Ax,Bx,e_Bx,zx,Xgas_H0,e_Xgas_H0,
+ Xgas_H1,e_Xgas_H1,chi2,dof,Prob
c ..............Just test output...........
write(6,2)
+ El,v_X_H_,e__X_H_,Ax,e_Ax,Bx,e_Bx,zx,Xgas_H0,e_Xgas_H0,
+ Xgas_H1,e_Xgas_H1,chi2,dof,Prob
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Observed and synthetic log N(H) and F_*_
C Format for file interpretation
3 format(
+ I6,4X,I3,2X,A16,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.3,1X,A41)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,239
read(1,'(A130)')ar__2
read(ar__2,3)
+ HD_1,Vcomp_1,Name_1,b_logNHo,logNHo,B_logNHo_1,F_o,e_F_o,
+ logNHs,e_logNHs,F_s,e_F_s,Prob_1,El_1
if(ar__2(1:6) .EQ. '') HD_1 = iNULL__
if(ar__2(11:13) .EQ. '') Vcomp_1 = iNULL__
if(ar__2(33:37) .EQ. '') b_logNHo = rNULL__
if(ar__2(39:43) .EQ. '') logNHo = rNULL__
if(ar__2(45:49) .EQ. '') B_logNHo_1 = rNULL__
if(ar__2(51:55) .EQ. '') F_o = rNULL__
if(ar__2(57:60) .EQ. '') e_F_o = rNULL__
if(ar__2(62:66) .EQ. '') logNHs = rNULL__
if(ar__2(68:71) .EQ. '') e_logNHs = rNULL__
if(ar__2(73:77) .EQ. '') F_s = rNULL__
if(ar__2(79:82) .EQ. '') e_F_s = rNULL__
if(ar__2(84:88) .EQ. '') Prob_1 = rNULL__
c ..............Just test output...........
write(6,3)
+ HD_1,Vcomp_1,Name_1,b_logNHo,logNHo,B_logNHo_1,F_o,e_F_o,
+ logNHs,e_logNHs,F_s,e_F_s,Prob_1,El_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Synthetic log N(H) and F_*_ for stars in the
* Local Bubble
C Format for file interpretation
4 format(
+ A8,1X,A15,1X,F6.2,1X,F6.2,1X,F4.1,1X,I3,1X,F5.2,1X,F4.2,1X,
+ F5.2,1X,F4.2,1X,F4.2,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,29
read(1,'(A84)')ar__3
read(ar__3,4)
+ WD,AName,GLON_1,GLAT_1,Vmag_1,Dist_1,logNHs_1,e_logNHs_1,
+ F_s_1,e_F_s_1,chi2_1,El_2
if(ar__3(71:74) .EQ. '') chi2_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ WD,AName,GLON_1,GLAT_1,Vmag_1,Dist_1,logNHs_1,e_logNHs_1,
+ F_s_1,e_F_s_1,chi2_1,El_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Observations and fits for Carbon - Krypton
C Format for file interpretation
5 format(
+ A2,1X,I6,2X,I3,2X,A16,1X,F5.2,3X,F5.2,1X,A1,1X,F5.2,1X,A1,1X,
+ A6,1X,F5.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,1475
read(1,'(A108)')ar__4
read(ar__4,5)
+ El_3,HD_2,Vcomp_2,Name_2,b_logNx,logNx,f_logNx,B_logNx_1,
+ n_logNx,r_logNx,F_,e_F_,v_X_H__1,e__X_H__1,Res,sigma
if(ar__4(4:9) .EQ. '') HD_2 = iNULL__
if(ar__4(12:14) .EQ. '') Vcomp_2 = iNULL__
if(ar__4(34:38) .EQ. '') b_logNx = rNULL__
if(ar__4(42:46) .EQ. '') logNx = rNULL__
if(ar__4(50:54) .EQ. '') B_logNx_1 = rNULL__
if(ar__4(65:69) .EQ. '') F_ = rNULL__
if(ar__4(71:74) .EQ. '') e_F_ = rNULL__
if(ar__4(76:81) .EQ. '') v_X_H__1 = rNULL__
if(ar__4(83:87) .EQ. '') e__X_H__1 = rNULL__
if(ar__4(89:94) .EQ. '') Res = rNULL__
if(ar__4(96:100) .EQ. '') sigma = rNULL__
c ..............Just test output...........
write(6,5)
+ El_3,HD_2,Vcomp_2,Name_2,b_logNx,logNx,f_logNx,B_logNx_1,
+ n_logNx,r_logNx,F_,e_F_,v_X_H__1,e__X_H__1,Res,sigma
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
6 format(A6,1X,A19,1X,A21,2X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,138
read(1,'(A69)')ar__5
read(ar__5,6)Ref,BibCode,Auth,Comm
c ..............Just test output...........
write(6,6)Ref,BibCode,Auth,Comm
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Notes of table 2
C Format for file interpretation
7 format(A1,1X,A74)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,44
read(1,'(A76)')ar__6
read(ar__6,7)F,Note
c ..............Just test output...........
write(6,7)F,Note
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end