Conversion of standardized ReadMe file for
file /./ftp/cats/VI/131 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-19
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/131 H13CN/HN13C linelist (Harris+, 2008)
*================================================================================
*A H^13^CN/HN^13^C linelist, model atmospheres and synthetic spectra for carbon
*stars.
* Harris G.J., Larner F.C., Tennyson J., Kaminsky B.M., Pavlenko Ya.V.,
* Jones H.R.A.
* <Mon. Not. R. Astron. Soc., 390, 143-148 (2008)>
* =2008MNRAS.390..143H
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'enlvs13.dat' ! Energy levels for H13CN/HN13C
integer*4 nr__
parameter (nr__=168110) ! Number of records
character*74 ar__ ! Full-size record
integer*4 Index (nr__) ! [1,168110] Unique index number of state
integer*4 J (nr__) ! [0,60] Rotational quantum number of lower state
integer*4 P (nr__) ! [0,1] Parity of the lower state
integer*4 n (nr__) ! [1,2220] Number of lower state in J-parity block
real*8 Eai (nr__) ! (cm-1) Value of the calculated energy level
integer*4 iso (nr__) ! [-2,1] Isomer configuration (G1)
integer*4 v1 (nr__) ! [-2,5] H-(CN) stretch approximate quantum number
integer*4 v2 (nr__) ! [-2,14] Bend approximate quantum number of state
integer*4 l (nr__) ! [-2,13] Vibrational angular momentum approximate
* quantum number of state
integer*4 v3 (nr__) ! [-2,5] C-N stretch approximate quantum number
real*8 Elab (nr__) ! (cm-1) ? Lab/empirical energy of the state
real*4 e_Elab (nr__) ! (cm-1) ? Estimated error on the lab empirical energy
character*1 n_Elab (nr__) ! [ec] Label identifying lab determined (e)
* or interpolated (c) data
character*1 f_Elab (nr__) ! [IC] Label identifying an Interpolated
* or ab initio Calculated energy
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h13cnlin.dat' ! Line list of H13CN/HN13C in spectroscopic format
integer*4 nr__1
parameter (nr__1=11739981) ! Number of records
character*26 ar__1 ! Full-size record
real*8 lambda (nr__1) ! (0.1nm) [7692/999897] Wavelength of the transition
real*4 gf (nr__1) ! gf in absorption
real*8 E (nr__1) ! (eV) Energy of the lower state
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h13cnlis.dat' ! Line list of H13CN/HN13C in format for
stellar modelling
integer*4 nr__2
parameter (nr__2=34179477) ! Number of records
character*104 ar__2 ! Full-size record
real*8 Freq (nr__2) ! (cm-1) Wavenumber of the transition
real*8 E_ (nr__2) ! (cm-1) Energy of the lower state
integer*4 J_ (nr__2) ! [0,60] Rotational quantum number of lower state
integer*4 P_ (nr__2) ! [0,1] Parity of the lower state (even=1, odd=0)
integer*4 n_ (nr__2) ! [1,2195] Number of lower state in J-parity block
integer*4 J__1 (nr__2) ! [0,60] Rotational quantum number of upper state
integer*4 P__1 (nr__2) ! [0,1] Parity of the upper state (even=1, odd=0)
integer*4 n__1 (nr__2) ! [1,2217] Number of upper state in J-parity block
real*4 A (nr__2) ! (s-1) Einstein A coefficient of the transition
integer*4 Index_ (nr__2) ! [1/168090] Unique index number of lower state
integer*4 Index__1 (nr__2) ! [2/168099] Unique index number of upper state
integer*4 iso_ (nr__2) ! [-2/1] Isomer configuration of lower state (G1)
integer*4 v1_ (nr__2) ! [-2/5] H-(CN) stretch approximate quantum number
* of lower state
integer*4 v2_ (nr__2) ! [-2/14] Bend approximate quantum number of
* lower state
integer*4 l_ (nr__2) ! [-2/13] Vibrational angular momentum approximate
* quantum number of lower state
integer*4 v3_ (nr__2) ! [-2/5] C-N stretch approximate quantum number of
* lower state
integer*4 iso__1 (nr__2) ! [-2/1] Isomer configuration of upper state (G1)
integer*4 v1__1 (nr__2) ! [-2/5] H-(CN) stretch approximate quantum number
* of upper state
integer*4 v2__1 (nr__2) ! [-2/14] Bend approximate quantum number of
* upper state
integer*4 l__1 (nr__2) ! [-2/13] Vibrational angular momentum approximate
* quantum number of upper state
integer*4 v3__1 (nr__2) ! [-2/5] C-N stretch approximate quantum number of
* upper state
character*2 Lab (nr__2) ! [ai lb] ab initio (ai) or lab determined (lb)
* frequency
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'dipoles.dat' ! Laboratory determined transition dipole and
Herman-Wallice constants for the CN stretch
fundamental and first hot band
integer*4 nr__3
parameter (nr__3=2) ! Number of records
character*55 ar__3 ! Full-size record
integer*4 v1__2 (nr__3) ! [0] H-C stretch quantum number for upper state
integer*4 v2__2 (nr__3) ! [0,1] Bend quantum number for upper state
integer*4 l__2 (nr__3) ! [0,1] vibrational angular momentum
* quantum number for upper state
integer*4 v3__2 (nr__3) ! [1] C-N stretch quantum number for upper state
integer*4 v1__3 (nr__3) ! [0] H-C stretch quantum number for lower state
integer*4 v2__3 (nr__3) ! [0,1] Bend quantum number for lower state
integer*4 l__3 (nr__3) ! [0,1] vibrational angular momentum
* quantum number for lower state
integer*4 v3__3 (nr__3) ! [0] C-N stretch quantum number for lower state
real*4 VIF (nr__3) ! [1] Vibrational intensity factor
real*8 mu (nr__3) ! (D) Band transition dipole
real*8 q1 (nr__3) ! First order Herman-Wallace constant
real*8 q2 (nr__3) ! Second order Herman-Wallace constant
real*8 qP (nr__3) ! Perpendicular Herman-Wallace constant
C=============================================================================
C Loading file 'enlvs13.dat' ! Energy levels for H13CN/HN13C
C Format for file interpretation
1 format(
+ I6,1X,I2,1X,I1,1X,I4,1X,F12.6,1X,I2,1X,I2,1X,I2,1X,I2,1X,I2,
+ 1X,F12.6,3X,E8.3,2X,A1,2X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'enlvs13.dat')
write(6,*) '....Loading file: enlvs13.dat'
do i__=1,168110
read(1,'(A74)')ar__
read(ar__,1)
+ Index(i__),J(i__),P(i__),n(i__),Eai(i__),iso(i__),v1(i__),
+ v2(i__),l(i__),v3(i__),Elab(i__),e_Elab(i__),n_Elab(i__),
+ f_Elab(i__)
if(ar__(46:57) .EQ. '') Elab(i__) = rNULL__
if(ar__(61:68) .EQ. '') e_Elab(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Index(i__),J(i__),P(i__),n(i__),Eai(i__),iso(i__),v1(i__),
+ v2(i__),l(i__),v3(i__),Elab(i__),e_Elab(i__),n_Elab(i__),
+ f_Elab(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h13cnlin.dat' ! Line list of H13CN/HN13C in spectroscopic format
C Format for file interpretation
2 format(F9.2,E10.3,F7.3)
C Effective file loading
open(unit=1,status='old',file=
+'h13cnlin.dat')
write(6,*) '....Loading file: h13cnlin.dat'
do i__=1,11739981
read(1,'(A26)')ar__1
read(ar__1,2)lambda(i__),gf(i__),E(i__)
c ..............Just test output...........
write(6,2)lambda(i__),gf(i__),E(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h13cnlis.dat' ! Line list of H13CN/HN13C in format for
* stellar modelling
C Format for file interpretation
3 format(
+ F12.6,F13.6,I3,I2,I5,I3,I2,I5,E10.3,I7,I7,1X,I3,I3,I3,I3,I3,
+ 1X,I3,I3,I3,I3,I3,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'h13cnlis.dat')
write(6,*) '....Loading file: h13cnlis.dat'
do i__=1,34179477
read(1,'(A104)')ar__2
read(ar__2,3)
+ Freq(i__),E_(i__),J_(i__),P_(i__),n_(i__),J__1(i__),P__1(i__),
+ n__1(i__),A(i__),Index_(i__),Index__1(i__),iso_(i__),v1_(i__),
+ v2_(i__),l_(i__),v3_(i__),iso__1(i__),v1__1(i__),v2__1(i__),
+ l__1(i__),v3__1(i__),Lab(i__)
c ..............Just test output...........
write(6,3)
+ Freq(i__),E_(i__),J_(i__),P_(i__),n_(i__),J__1(i__),P__1(i__),
+ n__1(i__),A(i__),Index_(i__),Index__1(i__),iso_(i__),v1_(i__),
+ v2_(i__),l_(i__),v3_(i__),iso__1(i__),v1__1(i__),v2__1(i__),
+ l__1(i__),v3__1(i__),Lab(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'dipoles.dat' ! Laboratory determined transition dipole and
* Herman-Wallice constants for the CN stretch
* fundamental and first hot band
C Format for file interpretation
4 format(I1,I2,I2,1X,I2,2X,I2,I2,I2,I2,F5.1,F9.5,F9.1,F7.1,F7.1)
C Effective file loading
open(unit=1,status='old',file=
+'dipoles.dat')
write(6,*) '....Loading file: dipoles.dat'
do i__=1,2
read(1,'(A55)')ar__3
read(ar__3,4)
+ v1__2(i__),v2__2(i__),l__2(i__),v3__2(i__),v1__3(i__),
+ v2__3(i__),l__3(i__),v3__3(i__),VIF(i__),mu(i__),q1(i__),
+ q2(i__),qP(i__)
c ..............Just test output...........
write(6,4)
+ v1__2(i__),v2__2(i__),l__2(i__),v3__2(i__),v1__3(i__),
+ v2__3(i__),l__3(i__),v3__3(i__),VIF(i__),mu(i__),q1(i__),
+ q2(i__),qP(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end