FORTRAN Generation
(/./ftp/cats/J/ApJ/459/686)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/459/686 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-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/459/686 Energy Deposition in Interstellar Dust Grains   (Dwek+ 1996)
*================================================================================
*Energy Deposition and Photoelectric Emission from the Interaction of 10 eV
*to 1 MeV Photons with Interstellar Dust Particles
*     Dwek E., Smith R.K.
*    <Astrophys. J. 459, 686 (1996)>
*    =1996ApJ...459..686D
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'silicate.dat'	! Deposition in silicate

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

      real*8        Ephot      (nr__) ! (eV) Energy of incoming photon
      real*4        E5nm       (nr__) ! (eV) *Energy deposition*Q_abs for grain size=5nm
      real*4        E7nm       (nr__) ! (eV) *Energy deposition*Q_abs for grain size=7nm
      real*4        E10nm      (nr__) ! (eV) *Energy deposition*Q_abs for grain size=10nm
      real*4        E15nm      (nr__) ! (eV) *Energy deposition*Q_abs for grain size=15nm
      real*4        E20nm      (nr__) ! (eV) *Energy deposition*Q_abs for grain size=20nm
      real*4        E30nm      (nr__) ! (eV) *Energy deposition*Q_abs for grain size=30nm
      real*4        E50nm      (nr__) ! (eV) *Energy deposition*Q_abs for grain size=50nm
      real*4        E70nm      (nr__) ! (eV) *Energy deposition*Q_abs for grain size=70nm
      real*4        E100nm     (nr__) ! (eV) *Energy deposition*Q_abs for grain size=100nm
      real*4        E150nm     (nr__) ! (eV) *Energy deposition*Q_abs for grain size=150nm
      real*4        E200nm     (nr__) ! (eV) *Energy deposition*Q_abs for grain size
      real*4        E300nm     (nr__) ! (eV) *Energy deposition*Q_abs for grain size=300nm
      real*4        E500nm     (nr__) ! (eV) *Energy deposition*Q_abs for grain size=500nm
      real*4        E700nm     (nr__) ! (eV) *Energy deposition*Q_abs for grain size=700nm
      real*4        E1um       (nr__) ! (eV) *Energy deposition*Q_abs for grain size=1micron
*Note on E5nm, E7nm, E10nm, E15nm, E20nm, E30nm, E50nm, E70nm, E100nm, E150nm:
*Note on E200nm, E300nm, E500nm, E700nm, E1um:
*   These fields contain the energy deposition multiplied by Q_abs. Where
*   Q_abs is the dust absorption efficiency calculated from Mie theory.
*   The grain size indicated is the dust particle radius.
*   The total energy deposited in a dust particle of radius a is therefore
*   flux * pi a^2 * dt * table_value; the flux has units of photons per
*   cm^2 per second, pi a^2 is the geometric area of the dust, dt is the
*   length of time considered, and the table value contains both Q_abs
*   (which when multiplied by the geometric cross section gives the cross
*   section for absorption) and the energy deposited when that photon is
*   absorbed, which may be the entire photon's energy or some lesser
*   value.

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

C  Declarations for 'graphite.dat'	! Deposition in graphite

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

      real*8        Ephot_1    (nr__1) ! (eV) Energy of incoming photon
      real*4        E5nm_1     (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=5nm
      real*4        E7nm_1     (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=7nm
      real*4        E10nm_1    (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=10nm
      real*4        E15nm_1    (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=15nm
      real*4        E20nm_1    (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=20nm
      real*4        E30nm_1    (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=30nm
      real*4        E50nm_1    (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=50nm
      real*4        E70nm_1    (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=70nm
      real*4        E100nm_1   (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=100nm
      real*4        E150nm_1   (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=150nm
      real*4        E200nm_1   (nr__1) ! (eV) *Energy deposition*Q_abs for grain size
      real*4        E300nm_1   (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=300nm
      real*4        E500nm_1   (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=500nm
      real*4        E700nm_1   (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=700nm
      real*4        E1um_1     (nr__1) ! (eV) *Energy deposition*Q_abs for grain size=1micron
*Note on E5nm, E7nm, E10nm, E15nm, E20nm, E30nm, E50nm, E70nm, E100nm, E150nm:
*Note on E200nm, E300nm, E500nm, E700nm, E1um:
*   These fields contain the energy deposition multiplied by Q_abs. Where
*   Q_abs is the dust absorption efficiency calculated from Mie theory.
*   The grain size indicated is the dust particle radius.
*   The total energy deposited in a dust particle of radius a is therefore
*   flux * pi a^2 * dt * table_value; the flux has units of photons per
*   cm^2 per second, pi a^2 is the geometric area of the dust, dt is the
*   length of time considered, and the table value contains both Q_abs
*   (which when multiplied by the geometric cross section gives the cross
*   section for absorption) and the energy deposited when that photon is
*   absorbed, which may be the entire photon's energy or some lesser
*   value.

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

C  Loading file 'silicate.dat'	! Deposition in silicate

C  Format for file interpretation

    1 format(
     +  F10.2,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
     +  2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
     +  E10.4,2X,E10.4,2X,E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'silicate.dat')
      write(6,*) '....Loading file: silicate.dat'
      do i__=1,60
        read(1,'(A191)')ar__
        read(ar__,1)
     +  Ephot(i__),E5nm(i__),E7nm(i__),E10nm(i__),E15nm(i__),
     +  E20nm(i__),E30nm(i__),E50nm(i__),E70nm(i__),E100nm(i__),
     +  E150nm(i__),E200nm(i__),E300nm(i__),E500nm(i__),E700nm(i__),
     +  E1um(i__)
c    ..............Just test output...........
        write(6,1)
     +  Ephot(i__),E5nm(i__),E7nm(i__),E10nm(i__),E15nm(i__),
     +  E20nm(i__),E30nm(i__),E50nm(i__),E70nm(i__),E100nm(i__),
     +  E150nm(i__),E200nm(i__),E300nm(i__),E500nm(i__),E700nm(i__),
     +  E1um(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'graphite.dat'	! Deposition in graphite

C  Format for file interpretation

    2 format(
     +  F10.2,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
     +  2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
     +  E10.4,2X,E10.4,2X,E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'graphite.dat')
      write(6,*) '....Loading file: graphite.dat'
      do i__=1,60
        read(1,'(A193)')ar__1
        read(ar__1,2)
     +  Ephot_1(i__),E5nm_1(i__),E7nm_1(i__),E10nm_1(i__),
     +  E15nm_1(i__),E20nm_1(i__),E30nm_1(i__),E50nm_1(i__),
     +  E70nm_1(i__),E100nm_1(i__),E150nm_1(i__),E200nm_1(i__),
     +  E300nm_1(i__),E500nm_1(i__),E700nm_1(i__),E1um_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Ephot_1(i__),E5nm_1(i__),E7nm_1(i__),E10nm_1(i__),
     +  E15nm_1(i__),E20nm_1(i__),E30nm_1(i__),E50nm_1(i__),
     +  E70nm_1(i__),E100nm_1(i__),E150nm_1(i__),E200nm_1(i__),
     +  E300nm_1(i__),E500nm_1(i__),E700nm_1(i__),E1um_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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