Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/449/855 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-18
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/449/855 Linelists for ammonia NH3, NH2D, ND2H, and ND3 (Coudert+, 2006)
*================================================================================
*Linelists for NH3, NH2D, ND2H, and ND3 with quadrupole coupling hyperfine
*components.
* Coudert L.H., Roueff E.
* <Astron. Astrophys. 449, 855 (2006)>
* =2006A&A...449..855C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Linelists for NH_3_
integer*4 nr__
parameter (nr__=632) ! Number of records
character*75 ar__ ! Full-size record
real*8 Freq (nr__) ! (MHz) Calculated or observed line frequency
real*8 e_Freq (nr__) ! (MHz) Frequency uncertainty
real*8 logInt (nr__) ! (nm+2.MHz) Base 10 logarithm of the integrated intensity
* in units of nm^2^MHz at 300K
integer*4 DOF (nr__) ! [3] Number of degrees of freedom
real*8 Elow (nr__) ! (cm-1) Lower state energy
integer*4 Gup (nr__) ! Upper state degeneracy
integer*4 Species (nr__) ! Species code (G1)
integer*4 QNFMT (nr__) ! [1304] JPL-type number describing the quantum
* numbers format (1304 for NH3 and ND3)
integer*4 J_ (nr__) ! Upper level rotational quantum number J
integer*4 K_ (nr__) ! Upper level rotational quantum number K
integer*4 v_ (nr__) ! Upper level vibrational quantum number v
* 0=>0s, 1=>0a
integer*4 F_ (nr__) ! Upper level total angular momentum F, F=I+J
integer*4 J__1 (nr__) ! Lower level rotational quantum number J
integer*4 K__1 (nr__) ! Lower level rotational quantum number K
integer*4 v__1 (nr__) ! Lower level vibrational quantum number v
* 0=>0s, 1=>0a
integer*4 F__1 (nr__) ! Lower level total angular momentum F, F=I+J
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Linelists for NH_2_D
integer*4 nr__1
parameter (nr__1=5910) ! Number of records
character*77 ar__1 ! Full-size record
real*8 Freq_1 (nr__1) ! (MHz) Calculated or observed line frequency
real*8 e_Freq_1 (nr__1) ! (MHz) Frequency Uncertainty
real*8 logInt_1 (nr__1) ! (nm+2.MHz) Base 10 logarithm of the integrated
* intensity in units of nm**2 MHz at 300 K
integer*4 DOF_1 (nr__1) ! [3] Number of degrees of freedom.
real*8 Elow_1 (nr__1) ! (cm-1) Lower state energy in cm-1
integer*4 Gup_1 (nr__1) ! Upper state degeneracy
integer*4 Species_1 (nr__1) ! Species code (G1)
integer*4 QNFMT_1 (nr__1) ! [1405] JPL-type number describing the quantum
* numbers format (14005 for NH2D and ND2H)
integer*4 J__2 (nr__1) ! Upper level rotational quantum number J
integer*4 Ka_ (nr__1) ! Upper level rotational quantum number Ka
integer*4 Kc_ (nr__1) ! Upper level rotational quantum number Kc
integer*4 v__2 (nr__1) ! Upper level vibrational quantum number
* 0=>0s, 1=>0a
integer*4 F__2 (nr__1) ! Upper level total angular momentum F, F=I+J
integer*4 J__3 (nr__1) ! Lower level rotational quantum number J
integer*4 Ka__1 (nr__1) ! Lower level rotational quantum number Ka
integer*4 Kc__1 (nr__1) ! Lower level rotational quantum number Kc
integer*4 v__3 (nr__1) ! Lower level vibrational quantum number
* 0=>0s, 1=>0a
integer*4 F__3 (nr__1) ! Lower level total angular momentum F, F=I+J
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! Linelists for ND_2_H
integer*4 nr__2
parameter (nr__2=6148) ! Number of records
character*77 ar__2 ! Full-size record
real*8 Freq_2 (nr__2) ! (MHz) Calculated or observed line frequency
real*8 e_Freq_2 (nr__2) ! (MHz) Frequency Uncertainty
real*8 logInt_2 (nr__2) ! (nm+2.MHz) Base 10 logarithm of the integrated
* intensity in units of nm**2 MHz at 300 K
integer*4 DOF_2 (nr__2) ! [3] Number of degrees of freedom.
real*8 Elow_2 (nr__2) ! (cm-1) Lower state energy in cm-1
integer*4 Gup_2 (nr__2) ! Upper state degeneracy
integer*4 Species_2 (nr__2) ! Species code (G1)
integer*4 QNFMT_2 (nr__2) ! [1405] JPL-type number describing the quantum
* numbers format (14005 for NH2D and ND2H)
integer*4 J__4 (nr__2) ! Upper level rotational quantum number J
integer*4 Ka__2 (nr__2) ! Upper level rotational quantum number Ka
integer*4 Kc__2 (nr__2) ! Upper level rotational quantum number Kc
integer*4 v__4 (nr__2) ! Upper level vibrational quantum number
* 0=>0s, 1=>0a
integer*4 F__4 (nr__2) ! Upper level total angular momentum F, F=I+J
integer*4 J__5 (nr__2) ! Lower level rotational quantum number J
integer*4 Ka__3 (nr__2) ! Lower level rotational quantum number Ka
integer*4 Kc__3 (nr__2) ! Lower level rotational quantum number Kc
integer*4 v__5 (nr__2) ! Lower level vibrational quantum number
* 0=>0s, 1=>0a
integer*4 F__5 (nr__2) ! Lower level total angular momentum F, F=I+J
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Linelists for ND_3_
integer*4 nr__3
parameter (nr__3=824) ! Number of records
character*75 ar__3 ! Full-size record
real*8 Freq_3 (nr__3) ! (MHz) Calculated or observed line frequency
real*8 e_Freq_3 (nr__3) ! (MHz) Frequency uncertainty
real*8 logInt_3 (nr__3) ! (nm+2.MHz) Base 10 logarithm of the integrated intensity
* in units of nm^2^MHz at 300K
integer*4 DOF_3 (nr__3) ! [3] Number of degrees of freedom
real*8 Elow_3 (nr__3) ! (cm-1) Lower state energy
integer*4 Gup_3 (nr__3) ! Upper state degeneracy
integer*4 Species_3 (nr__3) ! Species code (G1)
integer*4 QNFMT_3 (nr__3) ! [1304] JPL-type number describing the quantum
* numbers format (1304 for NH3 and ND3)
integer*4 J__6 (nr__3) ! Upper level rotational quantum number J
integer*4 K__2 (nr__3) ! Upper level rotational quantum number K
integer*4 v__6 (nr__3) ! Upper level vibrational quantum number v
* 0=>0s, 1=>0a
integer*4 F__6 (nr__3) ! Upper level total angular momentum F, F=I+J
integer*4 J__7 (nr__3) ! Lower level rotational quantum number J
integer*4 K__3 (nr__3) ! Lower level rotational quantum number K
integer*4 v__7 (nr__3) ! Lower level vibrational quantum number v
* 0=>0s, 1=>0a
integer*4 F__7 (nr__3) ! Lower level total angular momentum F, F=I+J
C=============================================================================
C Loading file 'table8.dat' ! Linelists for NH_3_
C Format for file interpretation
1 format(
+ F13.4,1X,F7.4,F8.4,1X,I1,1X,F9.4,I3,1X,I6,I4,1X,I1,I2,1X,I1,
+ 1X,I1,4X,I2,I2,1X,I1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,632
read(1,'(A75)')ar__
read(ar__,1)
+ Freq(i__),e_Freq(i__),logInt(i__),DOF(i__),Elow(i__),Gup(i__),
+ Species(i__),QNFMT(i__),J_(i__),K_(i__),v_(i__),F_(i__),
+ J__1(i__),K__1(i__),v__1(i__),F__1(i__)
c ..............Just test output...........
write(6,1)
+ Freq(i__),e_Freq(i__),logInt(i__),DOF(i__),Elow(i__),Gup(i__),
+ Species(i__),QNFMT(i__),J_(i__),K_(i__),v_(i__),F_(i__),
+ J__1(i__),K__1(i__),v__1(i__),F__1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Linelists for NH_2_D
C Format for file interpretation
2 format(
+ F13.4,F8.4,F8.4,1X,I1,F10.4,I3,I7,I4,I2,I2,I2,I2,I2,2X,I2,I2,
+ I2,I2,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,5910
read(1,'(A77)')ar__1
read(ar__1,2)
+ Freq_1(i__),e_Freq_1(i__),logInt_1(i__),DOF_1(i__),
+ Elow_1(i__),Gup_1(i__),Species_1(i__),QNFMT_1(i__),J__2(i__),
+ Ka_(i__),Kc_(i__),v__2(i__),F__2(i__),J__3(i__),Ka__1(i__),
+ Kc__1(i__),v__3(i__),F__3(i__)
c ..............Just test output...........
write(6,2)
+ Freq_1(i__),e_Freq_1(i__),logInt_1(i__),DOF_1(i__),
+ Elow_1(i__),Gup_1(i__),Species_1(i__),QNFMT_1(i__),J__2(i__),
+ Ka_(i__),Kc_(i__),v__2(i__),F__2(i__),J__3(i__),Ka__1(i__),
+ Kc__1(i__),v__3(i__),F__3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! Linelists for ND_2_H
C Format for file interpretation
3 format(
+ F13.4,F8.4,F8.4,1X,I1,F10.4,I3,I7,I4,I2,I2,I2,I2,I2,2X,I2,I2,
+ I2,I2,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,6148
read(1,'(A77)')ar__2
read(ar__2,3)
+ Freq_2(i__),e_Freq_2(i__),logInt_2(i__),DOF_2(i__),
+ Elow_2(i__),Gup_2(i__),Species_2(i__),QNFMT_2(i__),J__4(i__),
+ Ka__2(i__),Kc__2(i__),v__4(i__),F__4(i__),J__5(i__),
+ Ka__3(i__),Kc__3(i__),v__5(i__),F__5(i__)
c ..............Just test output...........
write(6,3)
+ Freq_2(i__),e_Freq_2(i__),logInt_2(i__),DOF_2(i__),
+ Elow_2(i__),Gup_2(i__),Species_2(i__),QNFMT_2(i__),J__4(i__),
+ Ka__2(i__),Kc__2(i__),v__4(i__),F__4(i__),J__5(i__),
+ Ka__3(i__),Kc__3(i__),v__5(i__),F__5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Linelists for ND_3_
C Format for file interpretation
4 format(
+ F13.4,1X,F7.4,F8.4,1X,I1,1X,F9.4,I3,1X,I6,I4,1X,I1,I2,1X,I1,
+ 1X,I1,4X,I2,I2,1X,I1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,824
read(1,'(A75)')ar__3
read(ar__3,4)
+ Freq_3(i__),e_Freq_3(i__),logInt_3(i__),DOF_3(i__),
+ Elow_3(i__),Gup_3(i__),Species_3(i__),QNFMT_3(i__),J__6(i__),
+ K__2(i__),v__6(i__),F__6(i__),J__7(i__),K__3(i__),v__7(i__),
+ F__7(i__)
c ..............Just test output...........
write(6,4)
+ Freq_3(i__),e_Freq_3(i__),logInt_3(i__),DOF_3(i__),
+ Elow_3(i__),Gup_3(i__),Species_3(i__),QNFMT_3(i__),J__6(i__),
+ K__2(i__),v__6(i__),F__6(i__),J__7(i__),K__3(i__),v__7(i__),
+ F__7(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end