Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/733/91 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-Mar-29 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/733/91 Diffuse interstellar band equivalent widths (Xiang+, 2011) *================================================================================ *A tale of two mysteries in interstellar astrophysics: the 2175 {AA} extinction *bump and diffuse interstellar bands. * Xiang F.Y., Li A., Zhong J.X. * <Astrophys. J., 733, 91 (2011)> * =2011ApJ...733...91X C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! The 2175{AA} extinction bump strengths (W_2175_) and DIB equivalent widths (WDIB) for 84 interstellar sightlines integer*4 nr__ parameter (nr__=84) ! Number of records character*97 ar__ ! Full-size record character*7 Star ! HD or BD name real*4 E_B_V ! (mag) Interstellar reddening E(B-V) (1) real*4 Rv ! ? Total-to-selective extinction ratio * (R_V_=A_V_/E(B-V); A_V_= V-band extinction) real*4 e_Rv ! ? Rv uncertainty integer*4 r_Rv ! ? Rv reference (2) real*4 c3 ! (um-2) Parameter c_3_ that relates to the strength * of the bump integer*4 r_c3 ! c3 reference (2) real*4 gamma ! (um-1) Band width {gamma} (extinction parameter) integer*4 r_gamma ! Gamma reference (2) real*4 W2175 ! (um-1) W_2175_/E(B-V) (3) real*4 W5707 ! (0.1pm) ? 5707 equivalent width in milli{AA} (1) real*4 W5780 ! (0.1pm) 5780 equivalent width in milli{AA} (1) real*4 W5797 ! (0.1pm) ? 5797 equivalent width in milli{AA} (1) real*4 W6195 ! (0.1pm) ? 6195/6196 equivalent width in milli{AA} (1) real*4 W6203 ! (0.1pm) ? 6203 equivalent width in milli{AA} (1) real*4 W6269 ! (0.1pm) ? 6269/6270 equivalent width in milli{AA} (1) real*4 W6284 ! (0.1pm) ? 6284 equivalent width in milli{AA} (1) real*4 W6376 ! (0.1pm) ? 6376/6379 equivalent width in milli{AA} (1) real*4 W6614 ! (0.1pm) ? 6614 equivalent width in milli{AA} (1) *Note (1): Sources for E(B-V) and W_DIB_ are * * Dorschner et al. (1977A&A....58..201D), * * Herbig (1975ApJ...196..129H, 1993ApJ...407..142H, 2000ApJ...542..334H), * * Seab & Snow (1984ApJ...277..200S), * * Benvenuti & Porceddu (1989A&A...223..329B), * * Desert et al. (1995A&A...303..223D), * * Sonnentrucker al. (1997A&A...327.1215S), * * Snow et al. (2002ApJ...573..670S), * * Wszolek & Godlowski (2003MNRAS.338..990W), * * Thorburn+ (2003ApJ...584..339T), * * Galazutdinov et al. (2004MNRAS.355..169G), * * Megier et al. (2005, Cat. J/MNRAS/358/563), * * Cox et al. (2007A&A...465..899C), * * Hobbs+ (2008ApJ...680.1256H, 2009, Cat. J/ApJ/705/32), * * McCall et al. (2010, Cat. J/ApJ/708/1628), * * Friedman et al. (2011ApJ...727...33F). * See Appendix and Tables 3 and 4 for details on how we obtain the DIB * equivalent widths adopted here. *Note (2): Sources as follows: * 1 = Wegner (2003, Cat. J/AN/324/219), * 2 = Geminale (2006, PhD thesis, Universita degli Studi di Padova); * 3 = Jenniskens & Greenberg (1993A&A...274..439J), * 4 = Fitzpatrick & Massa (2005AJ....130.1127F), * 5 = Whittet et al. (2004ApJ...602..291W), * 6 = Fitzpatrick & Massa (1990ApJS...72..163F), * 7 = Lewis et al. (2005ApJ...619..357L), * 8 = Sofia et al. (2005ApJ...625..167S). *Note (3): We use the black shaded area (W_2175_/E(B-V)={pi}c_3_/2{gamma}) * above the "linear background extinction" in the Fitzpatrick & Massa * (FM90; 1990ApJS...72..163F) decomposition scheme as a measure of the * strength of the 2175{AA} extinction bump (normalized by E(B-V)). * See section 1. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Sources for the equivalent widths of the 5707, 5780, 5797, 6203, 6284, and 6614 DIBs integer*4 nr__1 parameter (nr__1=797) ! Number of records character*87 ar__1 ! Full-size record integer*4 DIB ! (0.1nm) [5707/6614] DIB wavelength identification; in * Angstroms character*7 Star_1 ! [2905/+631964] Stellar HD or BD ID real*4 E_B_V_1 ! (mag) [0.05/1.28] E(B-V) real*4 W2175_1 ! (um-1) The 2175 equivalent width divided by E(B-V) real*4 W ! (0.1pm) ? Equivalent width of DIB (in milli{AA}) real*4 e_W ! (0.1pm) ? Uncertainty in W (in milli{AA}) character*1 q_W ! [AX] Adopt/Reject code (G1) character*22 Aut ! Author's name for Ref character*19 Ref ! Bibcode of source reference real*4 Wado ! (0.1pm) [5/1888]? Adopted equivalent width of DIB (in * milli-Angstroms) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Sources for the equivalent widths of the 6195/6996, 6269/6270, and 6376/6379 DIBs integer*4 nr__2 parameter (nr__2=385) ! Number of records character*104 ar__2 ! Full-size record integer*4 DIB1 ! (0.1nm) DIBs pair wavelength identification; in {AA} integer*4 DIB2 ! (0.1nm) DIBs pair wavelength identification; in {AA} character*7 Star_2 ! [2905/+631964] Stellar HD or BD ID real*4 E_B_V_1 ! (mag) [0.05/1.28] E(B-V) real*4 W2175_2 ! (um-1) The 2175 equivalent width divided by E(B-V) character*1 l_W1 ! Limit flag on W1 real*4 W1 ! (0.1pm) ? Equivalent width of 1st DIB in pair (in m{AA}) real*4 e_W1 ! (0.1pm) ? Uncertainty in W1 (in milli{AA}) character*1 q_W1 ! [AX] Adopt/Reject code for 1st DIB in pair (G1) real*4 W2 ! (0.1pm) ? Equivalent width of 2nd DIB in pair (in m{AA}) real*4 e_W2 ! (0.1pm) ? Uncertainty in W2 (in milli{AA}) character*1 q_W2 ! [AX] Adopt/Reject code for 2nd DIB in pair (G1) character*22 Aut_1 ! Author's name for Ref character*19 Ref_1 ! Bibcode of source reference real*4 Wado_1 ! (0.1pm) [2.5/246.9]? Adopted DIB equivalent * widths (in milli-Angstroms) C============================================================================= C Loading file 'table1.dat' ! The 2175{AA} extinction bump strengths (W_2175_) * and DIB equivalent widths (WDIB) for 84 * interstellar sightlines C Format for file interpretation 1 format( + A7,1X,F4.2,1X,F4.2,1X,F4.2,1X,I1,1X,F4.2,1X,I1,1X,F4.2,1X,I1, + 1X,F4.2,1X,F5.1,1X,F5.1,1X,F5.1,1X,F4.1,1X,F5.1,1X,F5.1,1X, + F6.1,1X,F5.1,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,84 read(1,'(A97)')ar__ read(ar__,1) + Star,E_B_V,Rv,e_Rv,r_Rv,c3,r_c3,gamma,r_gamma,W2175,W5707, + W5780,W5797,W6195,W6203,W6269,W6284,W6376,W6614 if(ar__(14:17) .EQ. '') Rv = rNULL__ if(ar__(19:22) .EQ. '') e_Rv = rNULL__ if(ar__(24:24) .EQ. '') r_Rv = iNULL__ if(ar__(45:49) .EQ. '') W5707 = rNULL__ if(ar__(57:61) .EQ. '') W5797 = rNULL__ if(ar__(63:66) .EQ. '') W6195 = rNULL__ if(ar__(68:72) .EQ. '') W6203 = rNULL__ if(ar__(74:78) .EQ. '') W6269 = rNULL__ if(ar__(80:85) .EQ. '') W6284 = rNULL__ if(ar__(87:91) .EQ. '') W6376 = rNULL__ if(ar__(93:97) .EQ. '') W6614 = rNULL__ c ..............Just test output........... write(6,1) + Star,E_B_V,Rv,e_Rv,r_Rv,c3,r_c3,gamma,r_gamma,W2175,W5707, + W5780,W5797,W6195,W6203,W6269,W6284,W6376,W6614 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Sources for the equivalent widths of the 5707, * 5780, 5797, 6203, 6284, and 6614 DIBs C Format for file interpretation 2 format( + I4,1X,A7,1X,F4.2,1X,F4.2,1X,F6.1,1X,F5.1,1X,A1,1X,A22,1X,A19, + 1X,F6.1) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,797 read(1,'(A87)')ar__1 read(ar__1,2)DIB,Star_1,E_B_V_1,W2175_1,W,e_W,q_W,Aut,Ref,Wado if(ar__1(24:29) .EQ. '') W = rNULL__ if(ar__1(31:35) .EQ. '') e_W = rNULL__ if(ar__1(82:87) .EQ. '') Wado = rNULL__ c ..............Just test output........... write(6,2)DIB,Star_1,E_B_V_1,W2175_1,W,e_W,q_W,Aut,Ref,Wado c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Sources for the equivalent widths of the * 6195/6996, 6269/6270, and 6376/6379 DIBs C Format for file interpretation 3 format( + I4,1X,I4,1X,A7,1X,F4.2,1X,F4.2,1X,A1,F5.1,1X,F4.1,1X,A1,1X, + F6.1,1X,F4.1,1X,A1,1X,A22,1X,A19,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,385 read(1,'(A104)')ar__2 read(ar__2,3) + DIB1,DIB2,Star_2,E_B_V_1,W2175_2,l_W1,W1,e_W1,q_W1,W2,e_W2, + q_W2,Aut_1,Ref_1,Wado_1 if(ar__2(30:34) .EQ. '') W1 = rNULL__ if(ar__2(36:39) .EQ. '') e_W1 = rNULL__ if(ar__2(43:48) .EQ. '') W2 = rNULL__ if(ar__2(50:53) .EQ. '') e_W2 = rNULL__ if(ar__2(100:104) .EQ. '') Wado_1 = rNULL__ c ..............Just test output........... write(6,3) + DIB1,DIB2,Star_2,E_B_V_1,W2175_2,l_W1,W1,e_W1,q_W1,W2,e_W2, + q_W2,Aut_1,Ref_1,Wado_1 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end