FORTRAN Generation
(/./ftp/cats/VI/56)

Conversion of standardized ReadMe file for file /./ftp/cats/VI/56 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. VI/56               Tests of shock chemistry in IC 443G   (Turner+, 1992)
*================================================================================
*Tests of Shock Chemistry in IC 443G
*     Turner B.E., Chan K.-W., Green S., Lubowich D.A.
*    <Astrophys. J. 399, 114 (1992)>
*    =1992ApJ...399..114T
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'csh2rate.dat'	! CS-H2  Collisional excitation rates

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

      integer*4     Ji         (nr__) ! [0/20] Initial rotational level
      integer*4     Jf         (nr__) ! [0/20] Final rotational level
      real*4        Rate_20    (nr__) ! (cm3/s) Collisional excitation rate at 20K
      real*4        Rate_40    (nr__) ! (cm3/s) Collisional excitation rate at 40K
      real*4        Rate_70    (nr__) ! (cm3/s) Collisional excitation rate at 70K
      real*4        Rate_100   (nr__) ! (cm3/s) Collisional excitation rate at 100K
      real*4        Rate_150   (nr__) ! (cm3/s) Collisional excitation rate at 150K
      real*4        Rate_200   (nr__) ! (cm3/s) Collisional excitation rate at 200K
      real*4        Rate_250   (nr__) ! (cm3/s) Collisional excitation rate at 250K
      real*4        Rate_300   (nr__) ! (cm3/s) Collisional excitation rate at 300K

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

C  Declarations for 'sioh2rat.dat'	! SiO-H2 Collisional excitation rates

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

      integer*4     Ji_1       (nr__1) ! [0/20] Initial rotational level
      integer*4     Jf_1       (nr__1) ! [0/20] Final rotational level
      real*4        Rate_20_1  (nr__1) ! (cm3/s) Collisional excitation rate at 20K
      real*4        Rate_40_1  (nr__1) ! (cm3/s) Collisional excitation rate at 40K
      real*4        Rate_70_1  (nr__1) ! (cm3/s) Collisional excitation rate at 70K
      real*4        Rate_100_1 (nr__1) ! (cm3/s) Collisional excitation rate at 100K
      real*4        Rate_150_1 (nr__1) ! (cm3/s) Collisional excitation rate at 150K
      real*4        Rate_200_1 (nr__1) ! (cm3/s) Collisional excitation rate at 200K
      real*4        Rate_250_1 (nr__1) ! (cm3/s) Collisional excitation rate at 250K
      real*4        Rate_300_1 (nr__1) ! (cm3/s) Collisional excitation rate at 300K

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

C  Loading file 'csh2rate.dat'	! CS-H2  Collisional excitation rates

C  Format for file interpretation

    1 format(
     +  3X,I2,1X,I2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,
     +  1X,E8.2,1X,E8.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'csh2rate.dat')
      write(6,*) '....Loading file: csh2rate.dat'
      do i__=1,420
        read(1,'(A80)')ar__
        read(ar__,1)
     +  Ji(i__),Jf(i__),Rate_20(i__),Rate_40(i__),Rate_70(i__),
     +  Rate_100(i__),Rate_150(i__),Rate_200(i__),Rate_250(i__),
     +  Rate_300(i__)
c    ..............Just test output...........
        write(6,1)
     +  Ji(i__),Jf(i__),Rate_20(i__),Rate_40(i__),Rate_70(i__),
     +  Rate_100(i__),Rate_150(i__),Rate_200(i__),Rate_250(i__),
     +  Rate_300(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'sioh2rat.dat'	! SiO-H2 Collisional excitation rates

C  Format for file interpretation

    2 format(
     +  3X,I2,1X,I2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,
     +  1X,E8.2,1X,E8.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'sioh2rat.dat')
      write(6,*) '....Loading file: sioh2rat.dat'
      do i__=1,420
        read(1,'(A80)')ar__1
        read(ar__1,2)
     +  Ji_1(i__),Jf_1(i__),Rate_20_1(i__),Rate_40_1(i__),
     +  Rate_70_1(i__),Rate_100_1(i__),Rate_150_1(i__),
     +  Rate_200_1(i__),Rate_250_1(i__),Rate_300_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  Ji_1(i__),Jf_1(i__),Rate_20_1(i__),Rate_40_1(i__),
     +  Rate_70_1(i__),Rate_100_1(i__),Rate_150_1(i__),
     +  Rate_200_1(i__),Rate_250_1(i__),Rate_300_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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