Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/321/293 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/A+A/321/293 SiO production in interstellar shocks (Schilke+ 1997)
*================================================================================
*SiO production in interstellar shocks.
* Schilke P., Walmsley C.M., Pineau des Forets G., Flower D.R.
* <Astron. Astrophys. 321, 293 (1997)>
* =1997A&A...321..293S (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1a' ! *Mantle and dissociative sputtering parameters
integer*4 nr__
parameter (nr__=30) ! Number of records
character*72 ar__ ! Full-size record
character*1 Type (nr__) ! [ab] a: mantle sputtering,
* b: dissociative sputtering
character*13 R1 (nr__) ! Beginning species of the reaction (1)
character*21 R2 (nr__) ! Resulting species of the reaction
real*4 S (nr__) ! Sputtering yield factor
real*4 Eth (nr__) ! (eV) Sputtering threshold energy
*Note (1): Asterisks denote species in grain mantles
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1b' ! *Dissociation on impact parameters
integer*4 nr__1
parameter (nr__1=12) ! Number of records
character*72 ar__1 ! Full-size record
character*1 Type_1 (nr__1) ! [c] c: dissociation on impact
character*13 R1_1 (nr__1) ! Beginning species of the reaction
character*21 R2_1 (nr__1) ! Resulting species of the reaction
real*4 Y0 (nr__1) ! Multiplicative constant (1)
real*4 Ediss (nr__1) ! Dissociation energy (1)
*Note (1): The dissociative probability has an Arrhenius form, Y0exp(-Ediss/E),
* where E is the impact energy, Ediss the dissociation energy and
* Y0 the multiplicative constant
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2' ! Rate coefficient parameters, for gas-phase
reactions involving Si-bearing species
integer*4 nr__2
parameter (nr__2=100) ! Number of records
character*80 ar__2 ! Full-size record
character*14 R1_2 (nr__2) ! Beginning species of the reaction
character*17 R2_2 (nr__2) ! Resulting species of the reaction
real*4 gamma (nr__2) ! Gamma coefficient (1)
real*4 alpha (nr__2) ! ? Alpha coefficient (1)
real*8 beta (nr__2) ! (K) ? Beta coefficient (1)
integer*4 Note (nr__2) ! ? Note (2)
*Note (1): Rate coefficient : gamma.(T/300)^alpha^.exp(-beta/T) cm^3^/s
*Note (2): 1: Langer & Glassgold, 1990, ApJ 352, 123
* 2: Supposed to be the same as CH (Gredel, 1990)
* 3: MacKay, 1995MNRAS.274..694M
* 4: UMIST ratefile (Millar et al., 1996, A&AS, in press)
* 5: Herbst (1995, private communication)
* 6: Elkin & Armentrout, 1984, J. Phys. Chem. 88, 5454
* 7: See the discussion of Flower et al., 1996MNRAS.280..447F
C=============================================================================
C Loading file 'table1a' ! *Mantle and dissociative sputtering parameters
C Format for file interpretation
1 format(A1,1X,A13,1X,A21,1X,E7.2,2X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1a')
write(6,*) '....Loading file: table1a'
do i__=1,30
read(1,'(A72)')ar__
read(ar__,1)Type(i__),R1(i__),R2(i__),S(i__),Eth(i__)
c ..............Just test output...........
write(6,1)Type(i__),R1(i__),R2(i__),S(i__),Eth(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1b' ! *Dissociation on impact parameters
C Format for file interpretation
2 format(A1,1X,A13,1X,A21,1X,E7.2,2X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1b')
write(6,*) '....Loading file: table1b'
do i__=1,12
read(1,'(A72)')ar__1
read(ar__1,2)
+ Type_1(i__),R1_1(i__),R2_1(i__),Y0(i__),Ediss(i__)
c ..............Just test output...........
write(6,2)
+ Type_1(i__),R1_1(i__),R2_1(i__),Y0(i__),Ediss(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2' ! Rate coefficient parameters, for gas-phase
* reactions involving Si-bearing species
C Format for file interpretation
3 format(A14,1X,A17,1X,E8.2,2X,F5.2,1X,F7.1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table2')
write(6,*) '....Loading file: table2'
do i__=1,100
read(1,'(A80)')ar__2
read(ar__2,3)
+ R1_2(i__),R2_2(i__),gamma(i__),alpha(i__),beta(i__),Note(i__)
if(ar__2(44:48) .EQ. '') alpha(i__) = rNULL__
if(ar__2(50:56) .EQ. '') beta(i__) = rNULL__
if(ar__2(58:58) .EQ. '') Note(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ R1_2(i__),R2_2(i__),gamma(i__),alpha(i__),beta(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end