Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/230/6 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/ApJS/230/6 HST -COS & -STIS absorption-line spectroscopy. II. (Keeney+, 2017) *================================================================================ *Characterizing the circumgalactic medium of nearby galaxies with HST/COS and *HST/STIS absorption-line spectroscopy. *II. Methods and models. * Keeney B.A., Stocke J.T., Danforth C.W., Shull J.M., Pratt C.T., * Froning C.S., Green J.C., Penton S.V., Savage B.D. * <Astrophys. J. Suppl. Ser., 230, 6-6 (2017)> * =2017ApJS..230....6K (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Targeted sample summary integer*4 nr__ parameter (nr__=13) ! Number of records character*135 ar__ ! Full-size record character*17 Line (nr__) ! QSO sight line identifier character*24 Gal (nr__) ! Associated galaxy identifier integer*4 cz1 (nr__) ! (km/s) [654/52914] First value of absorber velocity * (cz_abs_) integer*4 cz2 (nr__) ! (km/s) ? Second value of cz_abs_ integer*4 cz3 (nr__) ! (km/s) ? Third value of cz_abs_ integer*4 cz4 (nr__) ! (km/s) ? Fourth value of cz_abs_ integer*4 cz5 (nr__) ! (km/s) ? Fifth value of cz_abs_ integer*4 cz6 (nr__) ! (km/s) ? Sixth value of cz_abs_ character*10 Inst (nr__) ! Instrument(s) used for the absorber detections integer*4 czgal (nr__) ! (km/s) [649/52873] Galaxy's recession velocity (cz_gal_) real*4 Lum (nr__) ! (Lsun) [0.003/3.5] Rest-frame g-band galaxy luminosity (1) character*1 u_Lum (nr__) ! Uncertainty flag on Lum integer*4 rho (nr__) ! (kpc) [7/505] Impact parameter ({rho}) integer*4 Rvir (nr__) ! (kpc) [47/275] Virial radius (R_vir_) character*1 u_Rvir (nr__) ! Uncertainty flag on Rvir real*4 r_R (nr__) ! [0.1/2.1] Ratio of impact parameter to R_vir_ character*1 u_r_R (nr__) ! Uncertainty flag on r/R character*1 l_eta (nr__) ! Limit flag on eta real*4 eta (nr__) ! [0.02/4.3] {eta}_{rho}_ ratio (2) real*4 eVir (nr__) ! [0.02/4]? {eta}_vir_ ratio (3) character*1 u_eVir (nr__) ! Uncertainty flag on eVir real*4 eDelv (nr__) ! [0.01/6.8]? {eta}_{Delta}v_ ratio (4) character*1 u_eDelv (nr__) ! Uncertainty flag on eDelv *Note (1): The rest-frame g-band galaxy luminosity in L^*^ units (we adopt * M_g_^*^=-20.3, the value of Montero-Dorta & Prada 2009MNRAS.399.1106M * at h=0.7). *Note (2): The ratio of the impact parameter (in kpc) of the tabulated galaxy to * that of the next-nearest galaxy: * {eta}_{rho}_={rho}_ng_/{rho}_nng_. *Note (3): The ratio of the normalized impact parameter of the tabulated galaxy * to that of the next-nearest galaxy: * {eta}_vir_=({rho}/R_vir_)_ng_/({rho}/R_vir_)_nng_. *Note (4): The ratio of the normalized absorber-galaxy velocity difference of * the tabulated galaxy to that of the next-nearest galaxy: * {eta}_{Delta}v_=(|{Delta}v|/v_esc_)_ng_/(|{Delta}v|/v_esc_)_nng_. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Serendipitous sample summary integer*4 nr__1 parameter (nr__1=35) ! Number of records character*135 ar__1 ! Full-size record character*17 Line_1 (nr__1) ! QSO sight line identifier character*24 Gal_1 (nr__1) ! Associated galaxy identifier integer*4 cz1_1 (nr__1) ! (km/s) [654/52914] First value of absorber velocity * (cz_abs_) integer*4 cz2_1 (nr__1) ! (km/s) ? Second value of cz_abs_ integer*4 cz3_1 (nr__1) ! (km/s) ? Third value of cz_abs_ integer*4 cz4_1 (nr__1) ! (km/s) ? Fourth value of cz_abs_ integer*4 cz5_1 (nr__1) ! (km/s) ? Fifth value of cz_abs_ integer*4 cz6_1 (nr__1) ! (km/s) ? Sixth value of cz_abs_ character*10 Inst_1 (nr__1) ! Instrument(s) used for the absorber detections integer*4 czgal_1 (nr__1) ! (km/s) [649/52873] Galaxy's recession velocity (cz_gal_) real*4 Lum_1 (nr__1) ! (Lsun) [0.003/3.5] Rest-frame g-band galaxy luminosity (1) character*1 u_Lum_1 (nr__1) ! Uncertainty flag on Lum integer*4 rho_1 (nr__1) ! (kpc) [7/505] Impact parameter ({rho}) integer*4 Rvir_1 (nr__1) ! (kpc) [47/275] Virial radius (R_vir_) character*1 u_Rvir_1 (nr__1) ! Uncertainty flag on Rvir real*4 r_R_1 (nr__1) ! [0.1/2.1] Ratio of impact parameter to R_vir_ character*1 u_r_R_1 (nr__1) ! Uncertainty flag on r/R character*1 l_eta_1 (nr__1) ! Limit flag on eta real*4 eta_1 (nr__1) ! [0.02/4.3] {eta}_{rho}_ ratio (2) real*4 eVir_1 (nr__1) ! [0.02/4]? {eta}_vir_ ratio (3) character*1 u_eVir_1 (nr__1) ! Uncertainty flag on eVir real*4 eDelv_1 (nr__1) ! [0.01/6.8]? {eta}_{Delta}v_ ratio (4) character*1 u_eDelv_1 (nr__1) ! Uncertainty flag on eDelv *Note (1): The rest-frame g-band galaxy luminosity in L^*^ units (we adopt * M_g_^*^=-20.3, the value of Montero-Dorta & Prada 2009MNRAS.399.1106M * at h=0.7). *Note (2): The ratio of the impact parameter (in kpc) of the tabulated galaxy to * that of the next-nearest galaxy: * {eta}_{rho}_={rho}_ng_/{rho}_nng_. *Note (3): The ratio of the normalized impact parameter of the tabulated galaxy * to that of the next-nearest galaxy: * {eta}_vir_=({rho}/R_vir_)_ng_/({rho}/R_vir_)_nng_. *Note (4): The ratio of the normalized absorber-galaxy velocity difference of * the tabulated galaxy to that of the next-nearest galaxy: * {eta}_{Delta}v_=(|{Delta}v|/v_esc_)_ng_/(|{Delta}v|/v_esc_)_nng_. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Multi-component absorption line fits integer*4 nr__2 parameter (nr__2=467) ! Number of records character*122 ar__2 ! Full-size record integer*4 ID (nr__2) ! [1/187] System identifier character*17 Line_2 (nr__2) ! Sight line identifier character*6 Ion (nr__2) ! Ion identifier integer*4 cz (nr__2) ! (km/s) [643/52945]? Redshift from absorption line real*4 b (nr__2) ! (km/s) [4/128]? Doppler b-value character*1 u_b (nr__2) ! Uncertainty flag on b real*4 e_b (nr__2) ! (km/s) [0/57]? Uncertainty in b character*1 l_logN (nr__2) ! Limit flag on logN real*4 logN (nr__2) ! ([cm-2]) [11.6/18.4]? Log number abundance character*1 u_logN (nr__2) ! Uncertainty flag in logN real*4 e_logN (nr__2) ! ([cm-2]) [0.01/3]? Uncertainty in logN character*61 Comm (nr__2) ! Additional comments C============================================================================= C Loading file 'table1.dat' ! Targeted sample summary C Format for file interpretation 1 format( + A17,1X,A24,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,A10,1X,I5, + 1X,F5.3,A1,1X,I3,1X,I3,A1,1X,F4.2,A1,1X,A1,F4.2,1X,F4.2,A1,1X, + F4.2,A1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,13 read(1,'(A135)')ar__ read(ar__,1) + Line(i__),Gal(i__),cz1(i__),cz2(i__),cz3(i__),cz4(i__), + cz5(i__),cz6(i__),Inst(i__),czgal(i__),Lum(i__),u_Lum(i__), + rho(i__),Rvir(i__),u_Rvir(i__),r_R(i__),u_r_R(i__),l_eta(i__), + eta(i__),eVir(i__),u_eVir(i__),eDelv(i__),u_eDelv(i__) if(ar__(50:54) .EQ. '') cz2(i__) = iNULL__ if(ar__(56:60) .EQ. '') cz3(i__) = iNULL__ if(ar__(62:66) .EQ. '') cz4(i__) = iNULL__ if(ar__(68:72) .EQ. '') cz5(i__) = iNULL__ if(ar__(74:78) .EQ. '') cz6(i__) = iNULL__ if(ar__(125:128) .EQ. '') eVir(i__) = rNULL__ if(ar__(131:134) .EQ. '') eDelv(i__) = rNULL__ c ..............Just test output........... write(6,1) + Line(i__),Gal(i__),cz1(i__),cz2(i__),cz3(i__),cz4(i__), + cz5(i__),cz6(i__),Inst(i__),czgal(i__),Lum(i__),u_Lum(i__), + rho(i__),Rvir(i__),u_Rvir(i__),r_R(i__),u_r_R(i__),l_eta(i__), + eta(i__),eVir(i__),u_eVir(i__),eDelv(i__),u_eDelv(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Serendipitous sample summary C Format for file interpretation 2 format( + A17,1X,A24,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,A10,1X,I5, + 1X,F5.3,A1,1X,I3,1X,I3,A1,1X,F4.2,A1,1X,A1,F4.2,1X,F4.2,A1,1X, + F4.2,A1) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,35 read(1,'(A135)')ar__1 read(ar__1,2) + Line_1(i__),Gal_1(i__),cz1_1(i__),cz2_1(i__),cz3_1(i__), + cz4_1(i__),cz5_1(i__),cz6_1(i__),Inst_1(i__),czgal_1(i__), + Lum_1(i__),u_Lum_1(i__),rho_1(i__),Rvir_1(i__),u_Rvir_1(i__), + r_R_1(i__),u_r_R_1(i__),l_eta_1(i__),eta_1(i__),eVir_1(i__), + u_eVir_1(i__),eDelv_1(i__),u_eDelv_1(i__) if(ar__1(50:54) .EQ. '') cz2_1(i__) = iNULL__ if(ar__1(56:60) .EQ. '') cz3_1(i__) = iNULL__ if(ar__1(62:66) .EQ. '') cz4_1(i__) = iNULL__ if(ar__1(68:72) .EQ. '') cz5_1(i__) = iNULL__ if(ar__1(74:78) .EQ. '') cz6_1(i__) = iNULL__ if(ar__1(125:128) .EQ. '') eVir_1(i__) = rNULL__ if(ar__1(131:134) .EQ. '') eDelv_1(i__) = rNULL__ c ..............Just test output........... write(6,2) + Line_1(i__),Gal_1(i__),cz1_1(i__),cz2_1(i__),cz3_1(i__), + cz4_1(i__),cz5_1(i__),cz6_1(i__),Inst_1(i__),czgal_1(i__), + Lum_1(i__),u_Lum_1(i__),rho_1(i__),Rvir_1(i__),u_Rvir_1(i__), + r_R_1(i__),u_r_R_1(i__),l_eta_1(i__),eta_1(i__),eVir_1(i__), + u_eVir_1(i__),eDelv_1(i__),u_eDelv_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Multi-component absorption line fits C Format for file interpretation 3 format( + I3,1X,A17,1X,A6,1X,I5,1X,F5.1,A1,1X,F4.1,1X,A1,1X,F5.2,A1,1X, + F4.2,1X,A61) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,467 read(1,'(A122)')ar__2 read(ar__2,3) + ID(i__),Line_2(i__),Ion(i__),cz(i__),b(i__),u_b(i__),e_b(i__), + l_logN(i__),logN(i__),u_logN(i__),e_logN(i__),Comm(i__) if(ar__2(30:34) .EQ. '') cz(i__) = iNULL__ if(ar__2(36:40) .EQ. '') b(i__) = rNULL__ if(ar__2(43:46) .EQ. '') e_b(i__) = rNULL__ if(ar__2(50:54) .EQ. '') logN(i__) = rNULL__ if(ar__2(57:60) .EQ. '') e_logN(i__) = rNULL__ c ..............Just test output........... write(6,3) + ID(i__),Line_2(i__),Ion(i__),cz(i__),b(i__),u_b(i__),e_b(i__), + l_logN(i__),logN(i__),u_logN(i__),e_logN(i__),Comm(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end