FORTRAN Generation
(/./ftp/cats/VI/131)

Conversion of standardized ReadMe file for file /./ftp/cats/VI/131 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-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       ! [1,168110] Unique index number of state
      integer*4     J           ! [0,60] Rotational quantum number of lower state
      integer*4     P           ! [0,1] Parity of the lower state
      integer*4     n           ! [1,2220] Number of lower state in J-parity block
      real*8        Eai         ! (cm-1) Value of the calculated energy level
      integer*4     iso         ! [-2,1] Isomer configuration (G1)
      integer*4     v1          ! [-2,5] H-(CN) stretch approximate quantum number
      integer*4     v2          ! [-2,14] Bend approximate quantum number of state
      integer*4     l           ! [-2,13] Vibrational angular momentum approximate
*                                        quantum number of state
      integer*4     v3          ! [-2,5] C-N stretch approximate quantum number
      real*8        Elab        ! (cm-1) ? Lab/empirical energy of the state
      real*4        e_Elab      ! (cm-1) ? Estimated error on the lab empirical energy
      character*1   n_Elab      ! [ec] Label identifying lab determined (e)
*                                     or interpolated (c) data
      character*1   f_Elab      ! [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      ! (0.1nm) [7692/999897] Wavelength of the transition
      real*4        gf          ! gf in absorption
      real*8        E           ! (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        ! (cm-1) Wavenumber of the transition
      real*8        E_          ! (cm-1) Energy of the lower state
      integer*4     J_          ! [0,60] Rotational quantum number of lower state
      integer*4     P_          ! [0,1] Parity of the lower state (even=1, odd=0)
      integer*4     n_          ! [1,2195] Number of lower state in J-parity block
      integer*4     J__1        ! [0,60] Rotational quantum number of upper state
      integer*4     P__1        ! [0,1] Parity of the upper state (even=1, odd=0)
      integer*4     n__1        ! [1,2217] Number of upper state in J-parity block
      real*4        A           ! (s-1) Einstein A coefficient of the transition
      integer*4     Index_      ! [1/168090] Unique index number of lower state
      integer*4     Index__1    ! [2/168099] Unique index number of upper state
      integer*4     iso_        ! [-2/1] Isomer configuration of lower state (G1)
      integer*4     v1_         ! [-2/5] H-(CN) stretch approximate quantum number
*                                       of lower state
      integer*4     v2_         ! [-2/14] Bend approximate quantum number of
*                                        lower state
      integer*4     l_          ! [-2/13] Vibrational angular momentum approximate
*                                        quantum number of lower state
      integer*4     v3_         ! [-2/5] C-N stretch approximate quantum number of
*                                       lower state
      integer*4     iso__1      ! [-2/1] Isomer configuration of upper state (G1)
      integer*4     v1__1       ! [-2/5] H-(CN) stretch approximate quantum number
*                                       of upper state
      integer*4     v2__1       ! [-2/14] Bend approximate quantum number of
*                                        upper state
      integer*4     l__1        ! [-2/13] Vibrational angular momentum approximate
*                                        quantum number of upper state
      integer*4     v3__1       ! [-2/5] C-N stretch approximate quantum number of
*                                 upper state
      character*2   Lab         ! [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       ! [0] H-C stretch quantum number for upper state
      integer*4     v2__2       ! [0,1] Bend quantum number for upper state
      integer*4     l__2        ! [0,1] vibrational angular momentum
*                                       quantum number for upper state
      integer*4     v3__2       ! [1] C-N stretch quantum number for upper state
      integer*4     v1__3       ! [0] H-C stretch quantum number for lower state
      integer*4     v2__3       ! [0,1] Bend quantum number for lower state
      integer*4     l__3        ! [0,1] vibrational angular momentum
*                                       quantum number for lower state
      integer*4     v3__3       ! [0] C-N stretch quantum number for lower state
      real*4        VIF         ! [1] Vibrational intensity factor
      real*8        mu          ! (D) Band transition dipole
      real*8        q1          ! First order Herman-Wallace constant
      real*8        q2          ! Second order Herman-Wallace constant
      real*8        qP          ! 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,J,P,n,Eai,iso,v1,v2,l,v3,Elab,e_Elab,n_Elab,f_Elab
        if(ar__(46:57) .EQ. '') Elab = rNULL__
        if(ar__(61:68) .EQ. '') e_Elab = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Index,J,P,n,Eai,iso,v1,v2,l,v3,Elab,e_Elab,n_Elab,f_Elab
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,gf,E
c    ..............Just test output...........
        write(6,2)lambda,gf,E
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,E_,J_,P_,n_,J__1,P__1,n__1,A,Index_,Index__1,iso_,v1_,
     +  v2_,l_,v3_,iso__1,v1__1,v2__1,l__1,v3__1,Lab
c    ..............Just test output...........
        write(6,3)
     +  Freq,E_,J_,P_,n_,J__1,P__1,n__1,A,Index_,Index__1,iso_,v1_,
     +  v2_,l_,v3_,iso__1,v1__1,v2__1,l__1,v3__1,Lab
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,v2__2,l__2,v3__2,v1__3,v2__3,l__3,v3__3,VIF,mu,q1,q2,qP
c    ..............Just test output...........
        write(6,4)
     +  v1__2,v2__2,l__2,v3__2,v1__3,v2__3,l__3,v3__3,VIF,mu,q1,q2,qP
c    .......End.of.Just test output...........
      end do
      close(1)

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