Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/700/1299 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-Mar-28 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 (nr__) ! ? HD number integer*4 Vcomp (nr__) ! (km/s) ? Specific Heliocentric velocity component character*16 Name (nr__) ! Name (1) real*4 GLON (nr__) ! (deg) Galactic longitude real*4 GLAT (nr__) ! (deg) Galactic latitude real*4 Vmag (nr__) ! (mag) The Simbad V band magnitude (not used for * calculating the reddenings or Dist) character*16 SpType (nr__) ! MK spectral type character*6 r_SpType (nr__) ! Reference for SpType; see refs.dat file real*4 E_B_V (nr__) ! (mag) The (B-V) extinction real*4 Dist (nr__) ! (kpc) Heliocentric distance (1) real*4 Z (nr__) ! (kpc) Distance from Galactic plane real*4 b_logNHI (nr__) ! ([cm-2]) ? Lower limit on logNHI real*4 logNHI (nr__) ! ([cm-2]) ? Log best observed HI column density character*1 n_logNHI (nr__) ! Flag on logNHI (see notes.dat file) real*4 B_logNHI_1 (nr__) ! ([cm-2]) ? Upper limit on logNHI character*6 r_logNHI (nr__) ! Reference for logNHI; see refs.dat file real*4 b_logNH2 (nr__) ! ([cm-2]) ? Lower limit on logNH2 real*4 logNH2 (nr__) ! ([cm-2]) ? Log best observed H_2_ column density real*4 B_logNH2_1 (nr__) ! ([cm-2]) ? Upper limit on logNH2 character*6 r_logNH2 (nr__) ! 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 (nr__1) ! Element identification "X" (1) real*4 v_X_H_ (nr__1) ! ([-]) Adopted reference abundance relative to H (2) real*4 e__X_H_ (nr__1) ! ([-]) Uncertainty in X/H (2) real*4 Ax (nr__1) ! The A_X_ parameter; defined in Eq. 10 (3) real*4 e_Ax (nr__1) ! Uncertainty in Ax real*4 Bx (nr__1) ! The B_X_ parameter; defined in Eq. 10 (3) real*4 e_Bx (nr__1) ! Uncertainty in Bx (4) real*4 zx (nr__1) ! The z_X_ parameter; defined in Eq. 11 (3) real*4 Xgas_H0 (nr__1) ! ([-]) The [X_gas_/H]_0_ parameter; defined in Eq. 13 real*4 e_Xgas_H0 (nr__1) ! ([-]) Uncertainty in Xgas/H0; defined in Eq. 14 (3)(4) real*4 Xgas_H1 (nr__1) ! ([-]) The [X_gas_/H]_1_ parameter; defined in Eq. 15 real*4 e_Xgas_H1 (nr__1) ! ([-]) Uncertainty in Xgas/H1; defined in Eq. 16 (4) real*4 chi2 (nr__1) ! The {chi}^2^ value integer*4 dof (nr__1) ! Degrees of freedom real*4 Prob (nr__1) ! 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 (nr__2) ! ? HD number integer*4 Vcomp_1 (nr__2) ! (km/s) ? Specific Heliocentric velocity component (1) character*16 Name_1 (nr__2) ! Name real*4 b_logNHo (nr__2) ! ([cm-2]) ? Lower limit on logNHOb real*4 logNHo (nr__2) ! ([cm-2]) ? Log observed H column density real*4 B_logNHo_1 (nr__2) ! ([cm-2]) ? Upper limit on logNHOb real*4 F_o (nr__2) ! ? Observed line-of-sight depletion strength * factor real*4 e_F_o (nr__2) ! ? Uncertainty in F*o real*4 logNHs (nr__2) ! ([cm-2]) ? Log synthetic H column density real*4 e_logNHs (nr__2) ! ([cm-2]) ? Uncertainty in logNHs real*4 F_s (nr__2) ! ? Synthetic line-of-sight depletion strength * factor real*4 e_F_s (nr__2) ! ? Uncertainty in F*s real*4 Prob_1 (nr__2) ! ? Probability of a worse fit character*41 El_1 (nr__2) ! 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 (nr__3) ! White Dwarf identification (HHMM+DDd; B1950) character*15 AName (nr__3) ! Alternate name real*4 GLON_1 (nr__3) ! (deg) Galactic longitude real*4 GLAT_1 (nr__3) ! (deg) Galactic latitude real*4 Vmag_1 (nr__3) ! (mag) The Simbad V band magnitude integer*4 Dist_1 (nr__3) ! (pc) Heliocentric distance (1) real*4 logNHs_1 (nr__3) ! ([-]) Log synthetic H column density real*4 e_logNHs_1 (nr__3) ! ([-]) Uncertainty in logNHs real*4 F_s_1 (nr__3) ! Synthetic line-of-sight depletion strength factor real*4 e_F_s_1 (nr__3) ! Uncertainty in F*s real*4 chi2_1 (nr__3) ! ? The {chi}^2^ value character*9 El_2 (nr__3) ! 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 (nr__4) ! Element "X" integer*4 HD_2 (nr__4) ! ? HD number integer*4 Vcomp_2 (nr__4) ! (km/s) ? Specific Heliocentric velocity component character*16 Name_2 (nr__4) ! Name (1) real*4 b_logNx (nr__4) ! ([cm-2]) ? Lower limit on logNx real*4 logNx (nr__4) ! ([cm-2]) ? Log best observed column density (2) character*1 f_logNx (nr__4) ! [if] infinite or fixed value of logNx (3) real*4 B_logNx_1 (nr__4) ! ([cm-2]) ? Upper limit on logNx character*1 n_logNx (nr__4) ! [if] infinite or fixed value of B_logNx (3) character*6 r_logNx (nr__4) ! Reference codes for logNx (see refs.dat file) real*4 F_ (nr__4) ! ? Line-of-sight depletion strength factor real*4 e_F_ (nr__4) ! ? Uncertainty in F* real*4 v_X_H__1 (nr__4) ! ([-]) ? Log expected depletion [X_gas_/H] (4) real*4 e__X_H__1 (nr__4) ! ([-]) ? Uncertainty in Elgas/H real*4 Res (nr__4) ! ? Observed [X/H] minus computed from Eq. 10 real*4 sigma (nr__4) ! ? 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 (nr__5) ! Code of reference character*19 BibCode (nr__5) ! Bibcode character*21 Auth (nr__5) ! Author's name(s) character*19 Comm (nr__5) ! 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 (nr__6) ! [d-r] Flag code character*74 Note (nr__6) ! 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(i__),Vcomp(i__),Name(i__),GLON(i__),GLAT(i__),Vmag(i__), + SpType(i__),r_SpType(i__),E_B_V(i__),Dist(i__),Z(i__), + b_logNHI(i__),logNHI(i__),n_logNHI(i__),B_logNHI_1(i__), + r_logNHI(i__),b_logNH2(i__),logNH2(i__),B_logNH2_1(i__), + r_logNH2(i__) if(ar__(1:6) .EQ. '') HD(i__) = iNULL__ if(ar__(11:13) .EQ. '') Vcomp(i__) = iNULL__ if(ar__(95:99) .EQ. '') b_logNHI(i__) = rNULL__ if(ar__(103:107) .EQ. '') logNHI(i__) = rNULL__ if(ar__(111:115) .EQ. '') B_logNHI_1(i__) = rNULL__ if(ar__(127:131) .EQ. '') b_logNH2(i__) = rNULL__ if(ar__(133:137) .EQ. '') logNH2(i__) = rNULL__ if(ar__(139:143) .EQ. '') B_logNH2_1(i__) = rNULL__ c ..............Just test output........... write(6,1) + HD(i__),Vcomp(i__),Name(i__),GLON(i__),GLAT(i__),Vmag(i__), + SpType(i__),r_SpType(i__),E_B_V(i__),Dist(i__),Z(i__), + b_logNHI(i__),logNHI(i__),n_logNHI(i__),B_logNHI_1(i__), + r_logNHI(i__),b_logNH2(i__),logNH2(i__),B_logNH2_1(i__), + r_logNH2(i__) 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(i__),v_X_H_(i__),e__X_H_(i__),Ax(i__),e_Ax(i__),Bx(i__), + e_Bx(i__),zx(i__),Xgas_H0(i__),e_Xgas_H0(i__),Xgas_H1(i__), + e_Xgas_H1(i__),chi2(i__),dof(i__),Prob(i__) c ..............Just test output........... write(6,2) + El(i__),v_X_H_(i__),e__X_H_(i__),Ax(i__),e_Ax(i__),Bx(i__), + e_Bx(i__),zx(i__),Xgas_H0(i__),e_Xgas_H0(i__),Xgas_H1(i__), + e_Xgas_H1(i__),chi2(i__),dof(i__),Prob(i__) 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(i__),Vcomp_1(i__),Name_1(i__),b_logNHo(i__),logNHo(i__), + B_logNHo_1(i__),F_o(i__),e_F_o(i__),logNHs(i__),e_logNHs(i__), + F_s(i__),e_F_s(i__),Prob_1(i__),El_1(i__) if(ar__2(1:6) .EQ. '') HD_1(i__) = iNULL__ if(ar__2(11:13) .EQ. '') Vcomp_1(i__) = iNULL__ if(ar__2(33:37) .EQ. '') b_logNHo(i__) = rNULL__ if(ar__2(39:43) .EQ. '') logNHo(i__) = rNULL__ if(ar__2(45:49) .EQ. '') B_logNHo_1(i__) = rNULL__ if(ar__2(51:55) .EQ. '') F_o(i__) = rNULL__ if(ar__2(57:60) .EQ. '') e_F_o(i__) = rNULL__ if(ar__2(62:66) .EQ. '') logNHs(i__) = rNULL__ if(ar__2(68:71) .EQ. '') e_logNHs(i__) = rNULL__ if(ar__2(73:77) .EQ. '') F_s(i__) = rNULL__ if(ar__2(79:82) .EQ. '') e_F_s(i__) = rNULL__ if(ar__2(84:88) .EQ. '') Prob_1(i__) = rNULL__ c ..............Just test output........... write(6,3) + HD_1(i__),Vcomp_1(i__),Name_1(i__),b_logNHo(i__),logNHo(i__), + B_logNHo_1(i__),F_o(i__),e_F_o(i__),logNHs(i__),e_logNHs(i__), + F_s(i__),e_F_s(i__),Prob_1(i__),El_1(i__) 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(i__),AName(i__),GLON_1(i__),GLAT_1(i__),Vmag_1(i__), + Dist_1(i__),logNHs_1(i__),e_logNHs_1(i__),F_s_1(i__), + e_F_s_1(i__),chi2_1(i__),El_2(i__) if(ar__3(71:74) .EQ. '') chi2_1(i__) = rNULL__ c ..............Just test output........... write(6,4) + WD(i__),AName(i__),GLON_1(i__),GLAT_1(i__),Vmag_1(i__), + Dist_1(i__),logNHs_1(i__),e_logNHs_1(i__),F_s_1(i__), + e_F_s_1(i__),chi2_1(i__),El_2(i__) 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(i__),HD_2(i__),Vcomp_2(i__),Name_2(i__),b_logNx(i__), + logNx(i__),f_logNx(i__),B_logNx_1(i__),n_logNx(i__), + r_logNx(i__),F_(i__),e_F_(i__),v_X_H__1(i__),e__X_H__1(i__), + Res(i__),sigma(i__) if(ar__4(4:9) .EQ. '') HD_2(i__) = iNULL__ if(ar__4(12:14) .EQ. '') Vcomp_2(i__) = iNULL__ if(ar__4(34:38) .EQ. '') b_logNx(i__) = rNULL__ if(ar__4(42:46) .EQ. '') logNx(i__) = rNULL__ if(ar__4(50:54) .EQ. '') B_logNx_1(i__) = rNULL__ if(ar__4(65:69) .EQ. '') F_(i__) = rNULL__ if(ar__4(71:74) .EQ. '') e_F_(i__) = rNULL__ if(ar__4(76:81) .EQ. '') v_X_H__1(i__) = rNULL__ if(ar__4(83:87) .EQ. '') e__X_H__1(i__) = rNULL__ if(ar__4(89:94) .EQ. '') Res(i__) = rNULL__ if(ar__4(96:100) .EQ. '') sigma(i__) = rNULL__ c ..............Just test output........... write(6,5) + El_3(i__),HD_2(i__),Vcomp_2(i__),Name_2(i__),b_logNx(i__), + logNx(i__),f_logNx(i__),B_logNx_1(i__),n_logNx(i__), + r_logNx(i__),F_(i__),e_F_(i__),v_X_H__1(i__),e__X_H__1(i__), + Res(i__),sigma(i__) 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(i__),BibCode(i__),Auth(i__),Comm(i__) c ..............Just test output........... write(6,6)Ref(i__),BibCode(i__),Auth(i__),Comm(i__) 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(i__),Note(i__) c ..............Just test output........... write(6,7)F(i__),Note(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end