FORTRAN Generation
(/./ftp/cats/J/ApJ/705/32)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/705/32 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/705/32    Diffuse interstellar bands from HD 183143   (Hobbs+, 2009)
*================================================================================
*Studies of the diffuse interstellar bands.
*III. HD 183143.
*    Hobbs L.M., York D.G., Thorburn J.A., Snow T.P., Bishof M., Friedman S.D.,
*    McCall B.J., Oka T., Rachford B., Sonnentrucker P., Welty D.E.
*   <Astrophys. J., 705, 32-45 (2009)>
*   =2009ApJ...705...32H
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! Diffuse interstellar bands (DIBs) in the
                             spectrum of HD 183143

      integer*4 nr__
      parameter (nr__=418)	! Number of records
      character*103 ar__   	! Full-size record

      real*8        JD94       (nr__) ! (0.1nm) ? Jenniskens & Desert (1994A&AS..106...39J)
*                               central wavelength
      character*1   f_JD94     (nr__) ! [b] blend (see section 3.2 for details)
      real*8        GM00       (nr__) ! (0.1nm) ? Galazutdinov et al. (2000PASP..112..648G)
*                               central wavelength
      real*8        TC00       (nr__) ! (0.1nm) ? Tuairisg et al. (2000A&AS..142..225T) central
*                               wavelength
      character*1   f_TC00     (nr__) ! [b] blend (see section 3.2 for details)
      real*8        WS00       (nr__) ! (0.1nm) ? Weselak et al. (2000A&AS..142..239W) central
*                               wavelength
      character*1   f_WS00     (nr__) ! [b] blend (see section 3.2 for details)
      real*8        Wc         (nr__) ! (0.1nm) Band central wavelength in Angstroms
      real*4        FWHM       (nr__) ! (0.1nm) ? Band Full-Width at Half-Maximum
      character*1   l_EW       (nr__) ! Limit flag on EW
      real*4        EW         (nr__) ! (10-13m) ? Band equivalent width in milli-Angstroms
      character*1   f_EW       (nr__) ! [p] p: present
      real*4        e_EW       (nr__) ! (10-13m) ? Uncertainty in EW
      character*5   n_EW       (nr__) ! Warning note on EW (1)
      real*8        RWc        (nr__) ! (0.1nm) Band central wavelength in stellar rest frame
      character*1   l_BWc      (nr__) ! Limit flag on BWc (2)
      real*8        BWc        (nr__) ! (0.1nm) ? {beta} Ori central wavelength
      character*1   f_BWc      (nr__) ! [P] Indicates BWc shows a P Cygni profile
      real*4        BEW        (nr__) ! (10-13m) ? {beta} Ori equivalent width in milli-Angstroms
      character*1   f_BEW      (nr__) ! [b] blend (3)
      integer*4     Seq        (nr__) ! ? Running index number (only for table 2)
*Note (1): Note as follows:
*    d = blend with another DIB (diffuse interstellar band);
*    s = blend with stellar line;
*    t = blend with telluric line;
*  det = detector artifact;
*   s1 = blend with stellar NI;
*   s2 = blend with stellar NII.
*Note (2): Indicates more than one (isolated) stellar line interfers,
*     typically for a broad DIB (diffuse interstellar band). The
*     corresponding sum of the equivalent widths is given in BEW.
*Note (3):
*    b = the stellar line measured is itself a blend of at least two lines.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! 71 possible diffuse interstellar bands in the
                             spectrum of HD 183143

      integer*4 nr__1
      parameter (nr__1=71)	! Number of records
      character*103 ar__1  	! Full-size record

      real*8        JD94_1     (nr__1) ! (0.1nm) ? Jenniskens & Desert (1994A&AS..106...39J)
*                               central wavelength
      character*1   f_JD94_1   (nr__1) ! [b] blend (see section 3.2 for details)
      real*8        GM00_1     (nr__1) ! (0.1nm) ? Galazutdinov et al. (2000PASP..112..648G)
*                               central wavelength
      real*8        TC00_1     (nr__1) ! (0.1nm) ? Tuairisg et al. (2000A&AS..142..225T) central
*                               wavelength
      character*1   f_TC00_1   (nr__1) ! [b] blend (see section 3.2 for details)
      real*8        WS00_1     (nr__1) ! (0.1nm) ? Weselak et al. (2000A&AS..142..239W) central
*                               wavelength
      character*1   f_WS00_1   (nr__1) ! [b] blend (see section 3.2 for details)
      real*8        Wc_1       (nr__1) ! (0.1nm) Band central wavelength in Angstroms
      real*4        FWHM_1     (nr__1) ! (0.1nm) ? Band Full-Width at Half-Maximum
      character*1   l_EW_1     (nr__1) ! Limit flag on EW
      real*4        EW_1       (nr__1) ! (10-13m) ? Band equivalent width in milli-Angstroms
      character*1   f_EW_1     (nr__1) ! [p] p: present
      real*4        e_EW_1     (nr__1) ! (10-13m) ? Uncertainty in EW
      character*5   n_EW_1     (nr__1) ! Warning note on EW (1)
      real*8        RWc_1      (nr__1) ! (0.1nm) Band central wavelength in stellar rest frame
      character*1   l_BWc_1    (nr__1) ! Limit flag on BWc (2)
      real*8        BWc_1      (nr__1) ! (0.1nm) ? {beta} Ori central wavelength
      character*1   f_BWc_1    (nr__1) ! [P] Indicates BWc shows a P Cygni profile
      real*4        BEW_1      (nr__1) ! (10-13m) ? {beta} Ori equivalent width in milli-Angstroms
      character*1   f_BEW_1    (nr__1) ! [b] blend (3)
      integer*4     Seq_1      (nr__1) ! ? Running index number (only for table 2)
*Note (1): Note as follows:
*    d = blend with another DIB (diffuse interstellar band);
*    s = blend with stellar line;
*    t = blend with telluric line;
*  det = detector artifact;
*   s1 = blend with stellar NI;
*   s2 = blend with stellar NII.
*Note (2): Indicates more than one (isolated) stellar line interfers,
*     typically for a broad DIB (diffuse interstellar band). The
*     corresponding sum of the equivalent widths is given in BEW.
*Note (3):
*    b = the stellar line measured is itself a blend of at least two lines.

C=============================================================================

C  Loading file 'table2.dat'	! Diffuse interstellar bands (DIBs) in the
*                             spectrum of HD 183143

C  Format for file interpretation

    1 format(
     +  F7.2,1X,A1,1X,F7.2,1X,F7.2,1X,A1,1X,F7.2,1X,A1,1X,F7.2,1X,
     +  F5.2,1X,A1,F6.1,1X,A1,1X,F4.1,1X,A5,1X,F7.2,1X,A1,F7.2,1X,A1,
     +  1X,F5.1,1X,A1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,418
        read(1,'(A103)')ar__
        read(ar__,1)
     +  JD94(i__),f_JD94(i__),GM00(i__),TC00(i__),f_TC00(i__),
     +  WS00(i__),f_WS00(i__),Wc(i__),FWHM(i__),l_EW(i__),EW(i__),
     +  f_EW(i__),e_EW(i__),n_EW(i__),RWc(i__),l_BWc(i__),BWc(i__),
     +  f_BWc(i__),BEW(i__),f_BEW(i__),Seq(i__)
        if(ar__(1:7) .EQ. '') JD94(i__) = rNULL__
        if(ar__(11:17) .EQ. '') GM00(i__) = rNULL__
        if(ar__(19:25) .EQ. '') TC00(i__) = rNULL__
        if(ar__(29:35) .EQ. '') WS00(i__) = rNULL__
        if(ar__(47:51) .EQ. '') FWHM(i__) = rNULL__
        if(ar__(54:59) .EQ. '') EW(i__) = rNULL__
        if(ar__(63:66) .EQ. '') e_EW(i__) = rNULL__
        if(ar__(83:89) .EQ. '') BWc(i__) = rNULL__
        if(ar__(93:97) .EQ. '') BEW(i__) = rNULL__
        if(ar__(101:103) .EQ. '') Seq(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  JD94(i__),f_JD94(i__),GM00(i__),TC00(i__),f_TC00(i__),
     +  WS00(i__),f_WS00(i__),Wc(i__),FWHM(i__),l_EW(i__),EW(i__),
     +  f_EW(i__),e_EW(i__),n_EW(i__),RWc(i__),l_BWc(i__),BWc(i__),
     +  f_BWc(i__),BEW(i__),f_BEW(i__),Seq(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! 71 possible diffuse interstellar bands in the
*                             spectrum of HD 183143

C  Format for file interpretation

    2 format(
     +  F7.2,1X,A1,1X,F7.2,1X,F7.2,1X,A1,1X,F7.2,1X,A1,1X,F7.2,1X,
     +  F5.2,1X,A1,F6.1,1X,A1,1X,F4.1,1X,A5,1X,F7.2,1X,A1,F7.2,1X,A1,
     +  1X,F5.1,1X,A1,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,71
        read(1,'(A103)')ar__1
        read(ar__1,2)
     +  JD94_1(i__),f_JD94_1(i__),GM00_1(i__),TC00_1(i__),
     +  f_TC00_1(i__),WS00_1(i__),f_WS00_1(i__),Wc_1(i__),FWHM_1(i__),
     +  l_EW_1(i__),EW_1(i__),f_EW_1(i__),e_EW_1(i__),n_EW_1(i__),
     +  RWc_1(i__),l_BWc_1(i__),BWc_1(i__),f_BWc_1(i__),BEW_1(i__),
     +  f_BEW_1(i__),Seq_1(i__)
        if(ar__1(1:7) .EQ. '') JD94_1(i__) = rNULL__
        if(ar__1(11:17) .EQ. '') GM00_1(i__) = rNULL__
        if(ar__1(19:25) .EQ. '') TC00_1(i__) = rNULL__
        if(ar__1(29:35) .EQ. '') WS00_1(i__) = rNULL__
        if(ar__1(47:51) .EQ. '') FWHM_1(i__) = rNULL__
        if(ar__1(54:59) .EQ. '') EW_1(i__) = rNULL__
        if(ar__1(63:66) .EQ. '') e_EW_1(i__) = rNULL__
        if(ar__1(83:89) .EQ. '') BWc_1(i__) = rNULL__
        if(ar__1(93:97) .EQ. '') BEW_1(i__) = rNULL__
        if(ar__1(101:103) .EQ. '') Seq_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  JD94_1(i__),f_JD94_1(i__),GM00_1(i__),TC00_1(i__),
     +  f_TC00_1(i__),WS00_1(i__),f_WS00_1(i__),Wc_1(i__),FWHM_1(i__),
     +  l_EW_1(i__),EW_1(i__),f_EW_1(i__),e_EW_1(i__),n_EW_1(i__),
     +  RWc_1(i__),l_BWc_1(i__),BWc_1(i__),f_BWc_1(i__),BEW_1(i__),
     +  f_BEW_1(i__),Seq_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end