FORTRAN Generation
(/./ftp/cats/J/ApJ/649/L9)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/649/L9 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-Apr-19
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/649/L9     INTEGRAL IBIS census of the sky beyond 100keV  (Bazzano+, 2006)
*================================================================================
*INTEGRAL IBIS census of the sky beyond 100 keV.
*    Bazzano A., Stephen J.B., Fiocchi M., Bird A.J., Bassani L., Dean A.J.,
*    Malizia A., Ubertini P., Lebrun F., Walter R., Winkler C.
*   <Astrophys. J., 649, L9-L12 (2006)>
*   =2006ApJ...649L...9B
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! IBIS detections in (100-150)keV and (150-300)keV
                             energy bands

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

      character*14  Name       (nr__) ! Source name (1)
      character*11  Type       (nr__) ! Bird et al. (2006, Cat. <J/ApJ/636/765>)
*                                      classification of source
      real*8        S_20_40_I  (nr__) ! IBIS 20-40keV flux in mCrabs (1)
      real*4        e_S_20_40_I(nr__) ! Uncertainty in S(20-40)I
      real*4        S_40_100_I (nr__) ! IBIS 40-100keV flux in mCrabs (1)
      real*4        e_S_40_100_I(nr__) ! Uncertainty in S(40-100)I
      real*4        S_100_150_I(nr__) ! IBIS 100-150keV flux in mCrabs (1)
      real*4        e_S_100_150_I(nr__) ! Uncertainty in S(100-150)I
      integer*4     S_150_300_I(nr__) ! ? IBIS 150-300keV flux in mCrabs (1)
      integer*4     e_S_150_300_I(nr__) ! ? Uncertainty in IBIS S(150-300)I
      real*4        S_100_150_B(nr__) ! ? Beppo-Sax 100-150keV flux in mCrabs (1)
      real*4        E_S_100_150_B(nr__) ! ? Upper limit uncertainty in S(100-150)B
      real*4        e_S_100_150_B_1(nr__) ! ? Lower limit uncertainty in S(100-150)B
      character*1   Mo         (nr__) ! [a-c] Model flag on the Beppo-Sax data (2)
      real*4        S_150_200_B(nr__) ! ? Beppo-Sax 150-200keV flux in mCrabs (1)
      real*4        E_S_150_200_B(nr__) ! ? Upper limit uncertainty in S(150-200)B
      real*4        e_S_150_200_B_1(nr__) ! ? Lower limit uncertainty in S(150-200)B
      integer*4     r_Name     (nr__) ! ? Reference, in refs.dat file
*Note (1): Only sources detected with signal to noise >4 are reported.
*     1mCrab=3.2E-12erg/cm^2^/s (3.2fW/m^2^) in the 100-150keV range and
*     1mCrab=4.8E-12erg/cm^2^/s (4.8fW/m^2^) in the 150-300keV range.
*Note (2): Model flag as follows:
*      a = simple power law
*      b = cut-off power law
*      c = broken power law

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

C  Declarations for 'refs.dat'	! References

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

      integer*4     Ref        (nr__1) ! Reference code
      character*19  BibCode    (nr__1) ! BibCode
      character*20  Aut        (nr__1) ! Author's name
      character*141 Com        (nr__1) ! Comments

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

C  Declarations for 'types.dat'	! Source types (from Bird et al. 2007, J/ApJS/170/175)

      integer*4 nr__2
      parameter (nr__2=40)	! Number of records
      character*48 ar__2  	! Full-size record

      character*9   Type_1     (nr__2) ! Source type
      character*38  Def        (nr__2) ! Definition

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

C  Loading file 'table1.dat'	! IBIS detections in (100-150)keV and (150-300)keV
*                             energy bands

C  Format for file interpretation

    1 format(
     +  A14,1X,A11,1X,F7.2,1X,F4.2,1X,F6.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,
     +  I4,1X,I1,1X,F6.1,1X,F3.1,1X,F3.1,1X,A1,1X,F6.1,1X,F3.1,1X,
     +  F3.1,1X,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,49
        read(1,'(A103)')ar__
        read(ar__,1)
     +  Name(i__),Type(i__),S_20_40_I(i__),e_S_20_40_I(i__),
     +  S_40_100_I(i__),e_S_40_100_I(i__),S_100_150_I(i__),
     +  e_S_100_150_I(i__),S_150_300_I(i__),e_S_150_300_I(i__),
     +  S_100_150_B(i__),E_S_100_150_B(i__),e_S_100_150_B_1(i__),
     +  Mo(i__),S_150_200_B(i__),E_S_150_200_B(i__),
     +  e_S_150_200_B_1(i__),r_Name(i__)
        if(ar__(63:66) .EQ. '') S_150_300_I(i__) = iNULL__
        if(ar__(68:68) .EQ. '') e_S_150_300_I(i__) = iNULL__
        if(ar__(70:75) .EQ. '') S_100_150_B(i__) = rNULL__
        if(ar__(77:79) .EQ. '') E_S_100_150_B(i__) = rNULL__
        if(ar__(81:83) .EQ. '') e_S_100_150_B_1(i__) = rNULL__
        if(ar__(87:92) .EQ. '') S_150_200_B(i__) = rNULL__
        if(ar__(94:96) .EQ. '') E_S_150_200_B(i__) = rNULL__
        if(ar__(98:100) .EQ. '') e_S_150_200_B_1(i__) = rNULL__
        if(ar__(102:103) .EQ. '') r_Name(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),Type(i__),S_20_40_I(i__),e_S_20_40_I(i__),
     +  S_40_100_I(i__),e_S_40_100_I(i__),S_100_150_I(i__),
     +  e_S_100_150_I(i__),S_150_300_I(i__),e_S_150_300_I(i__),
     +  S_100_150_B(i__),E_S_100_150_B(i__),e_S_100_150_B_1(i__),
     +  Mo(i__),S_150_200_B(i__),E_S_150_200_B(i__),
     +  e_S_150_200_B_1(i__),r_Name(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    2 format(I2,1X,A19,1X,A20,1X,A141)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,23
        read(1,'(A185)')ar__1
        read(ar__1,2)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c    ..............Just test output...........
        write(6,2)Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'types.dat'	! Source types (from Bird et al. 2007, J/ApJS/170/175)

C  Format for file interpretation

    3 format(A9,1X,A38)

C  Effective file loading

      open(unit=1,status='old',file=
     +'types.dat')
      write(6,*) '....Loading file: types.dat'
      do i__=1,40
        read(1,'(A48)')ar__2
        read(ar__2,3)Type_1(i__),Def(i__)
c    ..............Just test output...........
        write(6,3)Type_1(i__),Def(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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