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