FORTRAN Generation
(/./ftp/cats/J/ApJS/232/19)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/232/19 into FORTRAN code for reading data files line by line.

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-28
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/ApJS/232/19     H_2_, D_2_, and HD c^3^{Pi}_u_^-^(v,N) levels     (Liu+, 2017)
*================================================================================
*H_2_X ^1^{Sigma}^+^_g_ - c^3^{Pi}_u_ excitation by electron impact: energies,
*spectra, emission yields, cross-sections, and H(1s) kinetic energy
*distributions.
*    Liu X., Shemansky D.E., Yoshii J., Liu M.J., Johnson P.V., Malone C.P.,
*    Khakoo M.A.
*    <Astrophys. J. Suppl. Ser., 232, 19-19 (2017)>
*    =2017ApJS..232...19L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table4.dat'	! Non-adiabatic transition energies and vibrational
                             overlap integrals of the
                             X^1^{Sigma}_g_^+^(v_i_,N_i_)-c^3^{Pi}_u_(v_j_,N_j_)

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

      integer*4     vj          ! [0/21] The X vibrational state
      integer*4     Nj          ! [1/17] The X rotational state
      integer*4     vi          ! [0] The c vibrational state
      integer*4     Ni          ! [0/15] The c rotational state
      real*8        Eij         ! (cm-1) [87724/118377] Transition energy
      real*4        VOI         ! [-0.5/0.5] Vibrational Overlap Integral (1)
*Note (1): The rotationally dependent FCF equals to the square of vibrational
*          overlap integral, |<vi,Ni|vj,Nj>|^2^.

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

C  Declarations for 'table5.dat'	! Energies, transition frequencies, transition
                             probabilities and Franck-Condon Factors of the
                             H_2_ a^3^{Sigma}_g_^+^-c^3^{Pi}_u_^-^ band systems

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

      integer*4     N1          ! [1/15] Upper rotational level
      integer*4     N0          ! [1/15] Lower rotational level
      integer*4     v1          ! [0/21] Upper vibrational state
      integer*4     v0          ! [0/20] Lower vibrational state
      real*8        Elow        ! (cm-1) [94941/118377] Lower state energy
      real*8        Ea_Ec       ! (cm-1) [-23233.5/23432.2] Transition frequency (1)
      real*4        A           ! (s-1) [/94299] Transition probability (2)
      real*4        FCF         ! [/1] Franck-Condon factor (FCF=|<v1,N1|v0,N0>|^2^)
*Note (1): The transition frequency is always defined as the energy difference
*          between the a^3^{Sigma}_g_^+^ and c^3^{Pi}_u_^-^ states.
*Note (2): A is positive even when the transition frequency is negative.

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

C  Declarations for 'table7.dat'	! Predissociation rates, kinetic energy release, and
                             FCFs of the c^3^{Pi}_u_^+^(v,N) levels

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

      integer*4     v           ! [0/21] The v vibrational state
      integer*4     N           ! [1/15] The N rotational level
      real*4        Width       ! (cm-1) [0.0004/3.5] Predissociation  width
      real*4        Rate        ! (s-1) Predissociation rate
      real*8        Ek          ! (eV) [7.2/10.2] Kenetic energy
      real*4        FCF_1       ! Franck-Condon factor in units of 1/hartree;
*                                FCF=|<c,v,N|b,Ek,N>|^2^

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

C  Declarations for 'table8.dat'	! Predissociation rates of the c^3^{Pi}_u_^-^ (v,N,J)
                             levels

      integer*4 nr__3
      parameter (nr__3=105)	! Number of records
      character*36 ar__3  	! Full-size record

      integer*4     v_1         ! [0/6] The v vibrational state
      integer*4     N_1         ! [1/15] The N rotational level
      real*4        F1          ! (s-1) [3/2447] F_1_ fine structure predissociation rate
      real*8        F2          ! (s-1) [243/270505] F_2_ fine structure predissociation
*                               rate
      real*4        F3          ! (s-1) [2/1073]? F_3_ fine structure predissociation rate
      real*8        Avg         ! (s-1) Average predissociation rate (1)
*Note (1): Average predissociation rate of F_1_, F_2_, and F_3_ fine-structure
*          components based on (2J+1) degeneracy.

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

C  Declarations for 'table9.dat'	! Calculated energies of the v=0-10 and N=1-15 levels
                             for the H_2_, HD, and D_2_ c^3^{Pi}_u_^-^ state

      integer*4 nr__4
      parameter (nr__4=33)	! Number of records
      character*155 ar__4  	! Full-size record

      character*2   ID          ! Calculation identifier (1)
      integer*4     v_2         ! [0/10] The v vibrational state
      real*8        Eng_N1      ! (cm-1) [94941/113080] Calculated energy
*                                for N=1 rotational level
      real*8        Eng_N2      ! (cm-1) Calculated energy for N=2 rotational level
      real*8        Eng_N3      ! (cm-1) Calculated energy for N=3 rotational level
      real*8        Eng_N4      ! (cm-1) Calculated energy for N=4 rotational level
      real*8        Eng_N5      ! (cm-1) Calculated energy for N=5 rotational level
      real*8        Eng_N6      ! (cm-1) Calculated energy for N=6 rotational level
      real*8        Eng_N7      ! (cm-1) Calculated energy for N=7 rotational level
      real*8        Eng_N8      ! (cm-1) Calculated energy for N=8 rotational level
      real*8        Eng_N9      ! (cm-1) Calculated energy for N=9 rotational level
      real*8        Eng_N10     ! (cm-1) Calculated energy for N=10 rotational level
      real*8        Eng_N11     ! (cm-1) Calculated energy for N=11 rotational level
      real*8        Eng_N12     ! (cm-1) Calculated energy for N=12 rotational level
      real*8        Eng_N13     ! (cm-1) Calculated energy for N=13 rotational level
      real*8        Eng_N14     ! (cm-1) Calculated energy for N=14 rotational level
      real*8        Eng_N15     ! (cm-1) [98600/116384] Calculated energy
*                                for N=15 rotational level
*Note (1): Identifier as follows:
*    H2 = Calculated H_2_ triplet c-(v,N) energy with {beta}=-0.055.
*    HD = Calculated HD triplet c-(v,N) energy with {beta}=-0.055.
*    D2 = Calculated D_2_ triplet c-(v,N) energy with {beta}=-0.055.

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

C  Loading file 'table4.dat'	! Non-adiabatic transition energies and vibrational
*                             overlap integrals of the
*                             X^1^{Sigma}_g_^+^(v_i_,N_i_)-c^3^{Pi}_u_(v_j_,N_j_)

C  Format for file interpretation

    1 format(I2,1X,I2,1X,I1,1X,I2,1X,F9.2,1X,E12.5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,1333
        read(1,'(A33)')ar__
        read(ar__,1)vj,Nj,vi,Ni,Eij,VOI
c    ..............Just test output...........
        write(6,1)vj,Nj,vi,Ni,Eij,VOI
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Energies, transition frequencies, transition
*                             probabilities and Franck-Condon Factors of the
*                             H_2_ a^3^{Sigma}_g_^+^-c^3^{Pi}_u_^-^ band systems

C  Format for file interpretation

    2 format(I2,1X,I2,1X,I2,1X,I2,1X,F9.2,1X,F9.2,1X,E10.4,1X,E10.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,4649
        read(1,'(A53)')ar__1
        read(ar__1,2)N1,N0,v1,v0,Elow,Ea_Ec,A,FCF
c    ..............Just test output...........
        write(6,2)N1,N0,v1,v0,Elow,Ea_Ec,A,FCF
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Predissociation rates, kinetic energy release, and
*                             FCFs of the c^3^{Pi}_u_^+^(v,N) levels

C  Format for file interpretation

    3 format(I2,1X,I2,1X,E9.3,1X,E9.3,1X,F8.5,1X,E9.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,257
        read(1,'(A44)')ar__2
        read(ar__2,3)v,N,Width,Rate,Ek,FCF_1
c    ..............Just test output...........
        write(6,3)v,N,Width,Rate,Ek,FCF_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Predissociation rates of the c^3^{Pi}_u_^-^ (v,N,J)
*                             levels

C  Format for file interpretation

    4 format(I1,1X,I2,1X,F6.1,1X,F8.1,1X,F6.1,1X,F8.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,105
        read(1,'(A36)')ar__3
        read(ar__3,4)v_1,N_1,F1,F2,F3,Avg
        if(ar__3(22:27) .EQ. '') F3 = rNULL__
c    ..............Just test output...........
        write(6,4)v_1,N_1,F1,F2,F3,Avg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Calculated energies of the v=0-10 and N=1-15 levels
*                             for the H_2_, HD, and D_2_ c^3^{Pi}_u_^-^ state

C  Format for file interpretation

    5 format(
     +  A2,1X,I2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,
     +  F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,1X,F9.2,
     +  1X,F9.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,33
        read(1,'(A155)')ar__4
        read(ar__4,5)
     +  ID,v_2,Eng_N1,Eng_N2,Eng_N3,Eng_N4,Eng_N5,Eng_N6,Eng_N7,
     +  Eng_N8,Eng_N9,Eng_N10,Eng_N11,Eng_N12,Eng_N13,Eng_N14,Eng_N15
c    ..............Just test output...........
        write(6,5)
     +  ID,v_2,Eng_N1,Eng_N2,Eng_N3,Eng_N4,Eng_N5,Eng_N6,Eng_N7,
     +  Eng_N8,Eng_N9,Eng_N10,Eng_N11,Eng_N12,Eng_N13,Eng_N14,Eng_N15
c    .......End.of.Just test output...........
      end do
      close(1)

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