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