Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/95/535 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 2026-Jun-13
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/95/535 Infrared CO line list for X 1{Sigma}+ state (Goorvitch, 1994)
*================================================================================
*Infrared CO line list for the X ^1^{Sigma}^+^ state
* Goorvitch D.
* <Astrophys. J. Suppl. Ser. 95, 535 (1994)>
* =1994ApJS...95..535G (SIMBAD/NED Reference)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Transition frequencies, gf values, and J values
for the bandheads of the R-branch overtone
bands of 12C16O
integer*4 nr__
parameter (nr__=21) ! Number of records
character*74 ar__ ! Full-size record
integer*4 v (nr__) ! Lower state vibrational quantum number
real*8 sigma1 (nr__) ! (cm-1) Frequency, fundamental
real*4 log_gf1 (nr__) ! gf-value, fundamental
integer*4 J1 (nr__) ! J-value, fundamental
real*8 sigma2 (nr__) ! (cm-1) Frequency, first overtone
real*4 log_gf2 (nr__) ! gf-value, first overtone
integer*4 J2 (nr__) ! J-value, first overtone
real*8 sigma3 (nr__) ! (cm-1) Frequency, second overtone
real*4 log_gf3 (nr__) ! gf-value, second overtone
integer*4 J3 (nr__) ! J-value, second overtone
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Rotationless dipole moments, M(0), for the
pure rotational transitions (D)
integer*4 nr__1
parameter (nr__1=6) ! Number of records
character*86 ar__1 ! Full-size record
integer*4 v_1 (nr__1) ! Vibrational quantum number
real*4 M26 (nr__1) ! (D) Rotationless dipole moment, 12C16O
real*4 M36 (nr__1) ! (D) Rotationless dipole moment, 13C16O
real*4 M28 (nr__1) ! (D) Rotationless dipole moment, 12C18O
real*4 M27 (nr__1) ! (D) Rotationless dipole moment, 12C17O
real*4 M38 (nr__1) ! (D) Rotationless dipole moment, 13C18O
real*4 M46 (nr__1) ! (D) Rotationless dipole moment, 14C16O
real*4 M37 (nr__1) ! (D) Rotationless dipole moment, 13C17O
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Polynomial fits to [F(m)]^{1/2} for the pure
rotational transitions
integer*4 nr__2
parameter (nr__2=42) ! Number of records
character*64 ar__2 ! Full-size record
character*2 Isotope (nr__2) ! Code for the isotope, e.g. '26' for 12C16O
integer*4 v_2 (nr__2) ! Vibrational quantum number
real*4 b0 (nr__2) ! b_0 coefficient
real*4 b1 (nr__2) ! b_1 coefficient
real*4 b2 (nr__2) ! b_2 coefficient
real*4 b3 (nr__2) ! b_3 coefficient
real*4 b4 (nr__2) ! b_4 coefficient
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Line list of 12C16O Delta v = +1 transitions
integer*4 nr__3
parameter (nr__3=147) ! Number of records
character*72 ar__3 ! Full-size record
real*8 sigma (nr__3) ! (cm-1) Transition frequency
real*4 R2 (nr__3) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A (nr__3) ! (Hz) Einstein transition probability (A-value)
real*8 E (nr__3) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf (nr__3) ! gf-value
real*4 S (nr__3) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu (nr__3) ! Upper state vibrational quantum number
integer*4 vl (nr__3) ! Lower state vibrational quantum number
character*1 Type (nr__3) ! [RP] Transition type
integer*4 Jl (nr__3) ! Lower state angular momentum number
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! Line list of 12C16O Delta v = +2 transitions
integer*4 nr__4
parameter (nr__4=147) ! Number of records
character*72 ar__4 ! Full-size record
real*8 sigma_1 (nr__4) ! (cm-1) Transition frequency
real*4 R2_1 (nr__4) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A_1 (nr__4) ! (Hz) Einstein transition probability (A-value)
real*8 E_1 (nr__4) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf_1 (nr__4) ! gf-value
real*4 S_1 (nr__4) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu_1 (nr__4) ! Upper state vibrational quantum number
integer*4 vl_1 (nr__4) ! Lower state vibrational quantum number
character*1 Type_1 (nr__4) ! [RP] Transition type
integer*4 Jl_1 (nr__4) ! Lower state angular momentum number
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Line list of 12C16O Delta v = +3 transitions
integer*4 nr__5
parameter (nr__5=147) ! Number of records
character*72 ar__5 ! Full-size record
real*8 sigma_2 (nr__5) ! (cm-1) Transition frequency
real*4 R2_2 (nr__5) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A_2 (nr__5) ! (Hz) Einstein transition probability (A-value)
real*8 E_2 (nr__5) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf_1 (nr__5) ! gf-value
real*4 S_2 (nr__5) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu_2 (nr__5) ! Upper state vibrational quantum number
integer*4 vl_2 (nr__5) ! Lower state vibrational quantum number
character*1 Type_2 (nr__5) ! [RP] Transition type
integer*4 Jl_2 (nr__5) ! Lower state angular momentum number
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Complete line list of the pure rotational
transitions
integer*4 nr__6
parameter (nr__6=2562) ! Number of records
character*75 ar__6 ! Full-size record
real*8 sigma_3 (nr__6) ! (cm-1) Transition frequency
real*4 R2_3 (nr__6) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A_3 (nr__6) ! (Hz) Einstein transition probability (A-value)
real*8 E_3 (nr__6) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf_1 (nr__6) ! gf-value
real*4 S_3 (nr__6) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu_3 (nr__6) ! Upper state vibrational quantum number
integer*4 vl_3 (nr__6) ! Lower state vibrational quantum number
character*1 Type_3 (nr__6) ! [RP] Transition type
integer*4 Jl_3 (nr__6) ! Lower state angular momentum number
character*2 Isotope_1 (nr__6) ! Code for the isotope, e.g. '26' for 12C16O
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! Complete line list of the fundamental
transitions
integer*4 nr__7
parameter (nr__7=43953) ! Number of records
character*75 ar__7 ! Full-size record
real*8 sigma_4 (nr__7) ! (cm-1) Transition frequency
real*4 R2_4 (nr__7) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A_4 (nr__7) ! (Hz) Einstein transition probability (A-value)
real*8 E_4 (nr__7) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf_1 (nr__7) ! gf-value
real*4 S_4 (nr__7) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu_4 (nr__7) ! Upper state vibrational quantum number
integer*4 vl_4 (nr__7) ! Lower state vibrational quantum number
character*1 Type_4 (nr__7) ! [RP] Transition type
integer*4 Jl_4 (nr__7) ! Lower state angular momentum number
character*2 Isotope_2 (nr__7) ! Code for the isotope, e.g. '26' for 12C16O
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table14.dat' ! Complete line list of the first overtone
transitions
integer*4 nr__8
parameter (nr__8=43953) ! Number of records
character*75 ar__8 ! Full-size record
real*8 sigma_5 (nr__8) ! (cm-1) Transition frequency
real*4 R2_5 (nr__8) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A_5 (nr__8) ! (Hz) Einstein transition probability (A-value)
real*8 E_5 (nr__8) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf_1 (nr__8) ! gf-value
real*4 S_5 (nr__8) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu_5 (nr__8) ! Upper state vibrational quantum number
integer*4 vl_5 (nr__8) ! Lower state vibrational quantum number
character*1 Type_5 (nr__8) ! [RP] Transition type
integer*4 Jl_5 (nr__8) ! Lower state angular momentum number
character*2 Isotope_3 (nr__8) ! Code for the isotope, e.g. '26' for 12C16O
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table15.dat' ! Complete line list of the second overtone
transitions
integer*4 nr__9
parameter (nr__9=43953) ! Number of records
character*75 ar__9 ! Full-size record
real*8 sigma_6 (nr__9) ! (cm-1) Transition frequency
real*4 R2_6 (nr__9) ! (D2) Expectation value, effective electric dipole
* matrix operator squared
real*4 A_6 (nr__9) ! (Hz) Einstein transition probability (A-value)
real*8 E_6 (nr__9) ! (cm-1) Lower state term energy, lowest level = 0
real*4 log_gf_1 (nr__9) ! gf-value
real*4 S_6 (nr__9) ! (cm/mol) Strength evaluated at T = 3000 K
integer*4 vu_6 (nr__9) ! Upper state vibrational quantum number
integer*4 vl_6 (nr__9) ! Lower state vibrational quantum number
character*1 Type_6 (nr__9) ! [RP] Transition type
integer*4 Jl_6 (nr__9) ! Lower state angular momentum number
character*2 Isotope_4 (nr__9) ! Code for the isotope, e.g. '26' for 12C16O
C=============================================================================
C Loading file 'table6.dat' ! Transition frequencies, gf values, and J values
* for the bandheads of the R-branch overtone
* bands of 12C16O
C Format for file interpretation
1 format(
+ I2,1X,F10.4,1X,E9.3,1X,I2,1X,F10.4,1X,E9.3,1X,I2,1X,F10.4,1X,
+ E9.3,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,21
read(1,'(A74)')ar__
read(ar__,1)
+ v(i__),sigma1(i__),log_gf1(i__),J1(i__),sigma2(i__),
+ log_gf2(i__),J2(i__),sigma3(i__),log_gf3(i__),J3(i__)
c ..............Just test output...........
write(6,1)
+ v(i__),sigma1(i__),log_gf1(i__),J1(i__),sigma2(i__),
+ log_gf2(i__),J2(i__),sigma3(i__),log_gf3(i__),J3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Rotationless dipole moments, M(0), for the
* pure rotational transitions (D)
C Format for file interpretation
2 format(
+ 1X,I1,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,
+ 1X,E11.4)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,6
read(1,'(A86)')ar__1
read(ar__1,2)
+ v_1(i__),M26(i__),M36(i__),M28(i__),M27(i__),M38(i__),
+ M46(i__),M37(i__)
c ..............Just test output...........
write(6,2)
+ v_1(i__),M26(i__),M36(i__),M28(i__),M27(i__),M38(i__),
+ M46(i__),M37(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Polynomial fits to [F(m)]^{1/2} for the pure
* rotational transitions
C Format for file interpretation
3 format(A2,1X,I1,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,42
read(1,'(A64)')ar__2
read(ar__2,3)
+ Isotope(i__),v_2(i__),b0(i__),b1(i__),b2(i__),b3(i__),b4(i__)
c ..............Just test output...........
write(6,3)
+ Isotope(i__),v_2(i__),b0(i__),b1(i__),b2(i__),b3(i__),b4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Line list of 12C16O Delta v = +1 transitions
C Format for file interpretation
4 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,147
read(1,'(A72)')ar__3
read(ar__3,4)
+ sigma(i__),R2(i__),A(i__),E(i__),log_gf(i__),S(i__),vu(i__),
+ vl(i__),Type(i__),Jl(i__)
c ..............Just test output...........
write(6,4)
+ sigma(i__),R2(i__),A(i__),E(i__),log_gf(i__),S(i__),vu(i__),
+ vl(i__),Type(i__),Jl(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! Line list of 12C16O Delta v = +2 transitions
C Format for file interpretation
5 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,147
read(1,'(A72)')ar__4
read(ar__4,5)
+ sigma_1(i__),R2_1(i__),A_1(i__),E_1(i__),log_gf_1(i__),
+ S_1(i__),vu_1(i__),vl_1(i__),Type_1(i__),Jl_1(i__)
c ..............Just test output...........
write(6,5)
+ sigma_1(i__),R2_1(i__),A_1(i__),E_1(i__),log_gf_1(i__),
+ S_1(i__),vu_1(i__),vl_1(i__),Type_1(i__),Jl_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Line list of 12C16O Delta v = +3 transitions
C Format for file interpretation
6 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,147
read(1,'(A72)')ar__5
read(ar__5,6)
+ sigma_2(i__),R2_2(i__),A_2(i__),E_2(i__),log_gf_1(i__),
+ S_2(i__),vu_2(i__),vl_2(i__),Type_2(i__),Jl_2(i__)
c ..............Just test output...........
write(6,6)
+ sigma_2(i__),R2_2(i__),A_2(i__),E_2(i__),log_gf_1(i__),
+ S_2(i__),vu_2(i__),vl_2(i__),Type_2(i__),Jl_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! Complete line list of the pure rotational
* transitions
C Format for file interpretation
7 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,I4,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,2562
read(1,'(A75)')ar__6
read(ar__6,7)
+ sigma_3(i__),R2_3(i__),A_3(i__),E_3(i__),log_gf_1(i__),
+ S_3(i__),vu_3(i__),vl_3(i__),Type_3(i__),Jl_3(i__),
+ Isotope_1(i__)
c ..............Just test output...........
write(6,7)
+ sigma_3(i__),R2_3(i__),A_3(i__),E_3(i__),log_gf_1(i__),
+ S_3(i__),vu_3(i__),vl_3(i__),Type_3(i__),Jl_3(i__),
+ Isotope_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! Complete line list of the fundamental
* transitions
C Format for file interpretation
8 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,I4,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,43953
read(1,'(A75)')ar__7
read(ar__7,8)
+ sigma_4(i__),R2_4(i__),A_4(i__),E_4(i__),log_gf_1(i__),
+ S_4(i__),vu_4(i__),vl_4(i__),Type_4(i__),Jl_4(i__),
+ Isotope_2(i__)
c ..............Just test output...........
write(6,8)
+ sigma_4(i__),R2_4(i__),A_4(i__),E_4(i__),log_gf_1(i__),
+ S_4(i__),vu_4(i__),vl_4(i__),Type_4(i__),Jl_4(i__),
+ Isotope_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table14.dat' ! Complete line list of the first overtone
* transitions
C Format for file interpretation
9 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,I4,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table14.dat')
write(6,*) '....Loading file: table14.dat'
do i__=1,43953
read(1,'(A75)')ar__8
read(ar__8,9)
+ sigma_5(i__),R2_5(i__),A_5(i__),E_5(i__),log_gf_1(i__),
+ S_5(i__),vu_5(i__),vl_5(i__),Type_5(i__),Jl_5(i__),
+ Isotope_3(i__)
c ..............Just test output...........
write(6,9)
+ sigma_5(i__),R2_5(i__),A_5(i__),E_5(i__),log_gf_1(i__),
+ S_5(i__),vu_5(i__),vl_5(i__),Type_5(i__),Jl_5(i__),
+ Isotope_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table15.dat' ! Complete line list of the second overtone
* transitions
C Format for file interpretation
10 format(
+ F9.4,1X,E9.3,1X,E9.3,1X,F10.4,1X,E9.3,1X,E9.3,1X,I2,1X,I2,1X,
+ A1,I4,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'table15.dat')
write(6,*) '....Loading file: table15.dat'
do i__=1,43953
read(1,'(A75)')ar__9
read(ar__9,10)
+ sigma_6(i__),R2_6(i__),A_6(i__),E_6(i__),log_gf_1(i__),
+ S_6(i__),vu_6(i__),vl_6(i__),Type_6(i__),Jl_6(i__),
+ Isotope_4(i__)
c ..............Just test output...........
write(6,10)
+ sigma_6(i__),R2_6(i__),A_6(i__),E_6(i__),log_gf_1(i__),
+ S_6(i__),vu_6(i__),vl_6(i__),Type_6(i__),Jl_6(i__),
+ Isotope_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end