FORTRAN Generation
(/./ftp/cats/J/ApJ/733/91)

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