FORTRAN Generation
(/./ftp/cats/J/A_A/657/A24)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/657/A24 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-Aug-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/657/A24         HCCCHCN laboratory observations           (Cabezas+, 2022)
*================================================================================
*Laboratory observation and astronomical search of 1-cyano propargyl radical,
*HCCCHCN.
*    Cabezas C., Nakajima M., Chang C.H., Agundez M., Endo Y., Cernicharo J.
*    <Astron. Astrophys. 657, A24 (2022)>
*    =2022A&A...657A..24C        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'tablea2.dat'	! Predicted frequencies for HC_3_HCN at 10K

      integer*4 nr__
      parameter (nr__=60708)	! Number of records
      character*67 ar__   	! Full-size record

      real*8        Freq       (nr__) ! (MHz) Calculated transition frequency
      real*8        O_C        (nr__) ! (MHz) Uncertainty of calculated transition frequency
      real*8        logInt     (nr__) ! ([nm+2.MHz]) Base 10 logarithm of the integrated intensity
      real*8        Elo        (nr__) ! (cm-1) The energy of the lower state
      integer*4     gu         (nr__) ! Degeneracy of the upper level
      integer*4     J_         (nr__) ! Quantum number J for upper level
      integer*4     Ka_        (nr__) ! Quantum number Ka for upper level
      integer*4     Kc_        (nr__) ! Quantum number Kc for upper level
      integer*4     TSP_       (nr__) ! Aggregate spin number for upper level
      integer*4     It_        (nr__) ! Quantum number Itotal for upper level
      integer*4     J__1       (nr__) ! Quantum number J for lower level
      integer*4     Ka__1      (nr__) ! Quantum number Ka for lower level
      integer*4     Kc__1      (nr__) ! Quantum number Kc for lower level
      integer*4     TSP__1     (nr__) ! Aggregate spin number for lower level
      integer*4     It__1      (nr__) ! Quantum number Itotal for lower level

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'tablea3.dat'	! Predicted frequencies for HC_3_HCN at 300K

      integer*4 nr__1
      parameter (nr__1=43202)	! Number of records
      character*67 ar__1  	! Full-size record

      real*8        Freq_1     (nr__1) ! (MHz) Calculated transition frequency
      real*8        O_C_1      (nr__1) ! (MHz) Uncertainty of calculated transition frequency
      real*8        logInt_1   (nr__1) ! ([nm+2.MHz]) Base 10 logarithm of the integrated intensity
      real*8        Elo_1      (nr__1) ! (cm-1) The energy of the lower state
      integer*4     gu_1       (nr__1) ! Degeneracy of the upper level
      integer*4     J__2       (nr__1) ! Quantum number J for upper level
      integer*4     Ka__2      (nr__1) ! Quantum number Ka for upper level
      integer*4     Kc__2      (nr__1) ! Quantum number Kc for upper level
      integer*4     TSP__2     (nr__1) ! Aggregate spin number for upper level
      integer*4     It__2      (nr__1) ! Quantum number Itotal for upper level
      integer*4     J__3       (nr__1) ! Quantum number J for lower level
      integer*4     Ka__3      (nr__1) ! Quantum number Ka for lower level
      integer*4     Kc__3      (nr__1) ! Quantum number Kc for lower level
      integer*4     TSP__3     (nr__1) ! Aggregate spin number for lower level
      integer*4     It__3      (nr__1) ! Quantum number Itotal for lower level

C=============================================================================

C  Loading file 'tablea2.dat'	! Predicted frequencies for HC_3_HCN at 10K

C  Format for file interpretation

    1 format(
     +  2X,F11.4,1X,F7.4,1X,F8.4,2X,F8.4,1X,I2,2X,I2,I2,I2,I2,I2,1X,
     +  I3,I2,I2,I2,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea2.dat')
      write(6,*) '....Loading file: tablea2.dat'
      do i__=1,60708
        read(1,'(A67)')ar__
        read(ar__,1)
     +  Freq(i__),O_C(i__),logInt(i__),Elo(i__),gu(i__),J_(i__),
     +  Ka_(i__),Kc_(i__),TSP_(i__),It_(i__),J__1(i__),Ka__1(i__),
     +  Kc__1(i__),TSP__1(i__),It__1(i__)
c    ..............Just test output...........
        write(6,1)
     +  Freq(i__),O_C(i__),logInt(i__),Elo(i__),gu(i__),J_(i__),
     +  Ka_(i__),Kc_(i__),TSP_(i__),It_(i__),J__1(i__),Ka__1(i__),
     +  Kc__1(i__),TSP__1(i__),It__1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'tablea3.dat'	! Predicted frequencies for HC_3_HCN at 300K

C  Format for file interpretation

    2 format(
     +  2X,F11.4,1X,F7.4,1X,F8.4,2X,F8.4,1X,I2,2X,I2,I2,I2,I2,I2,1X,
     +  I3,I2,I2,I2,I2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablea3.dat')
      write(6,*) '....Loading file: tablea3.dat'
      do i__=1,43202
        read(1,'(A67)')ar__1
        read(ar__1,2)
     +  Freq_1(i__),O_C_1(i__),logInt_1(i__),Elo_1(i__),gu_1(i__),
     +  J__2(i__),Ka__2(i__),Kc__2(i__),TSP__2(i__),It__2(i__),
     +  J__3(i__),Ka__3(i__),Kc__3(i__),TSP__3(i__),It__3(i__)
c    ..............Just test output...........
        write(6,2)
     +  Freq_1(i__),O_C_1(i__),logInt_1(i__),Elo_1(i__),gu_1(i__),
     +  J__2(i__),Ka__2(i__),Kc__2(i__),TSP__2(i__),It__2(i__),
     +  J__3(i__),Ka__3(i__),Kc__3(i__),TSP__3(i__),It__3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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