FORTRAN Generation
(/./ftp/cats/J/ApJ/804/76)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/804/76 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-15
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/ApJ/804/76        ZnII lines and collision strengths        (Kisielius+, 2015)
*================================================================================
*Atomic data for Zn II: improving spectral diagnostics of chemical evolution
*in high-redshift galaxies.
*    Kisielius R., Kulkarni V.P., Ferland G.J., Bogdanovich P., Som D.,
*    Lykins M.L.
*   <Astrophys. J., 804, 76 (2015)>
*   =2015ApJ...804...76K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Comparison of calculated ZnII level energies
                              and their percentage deviations with experimental
                              data from the NIST database

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

      integer*4     N          (nr__) ! [1/27] Level index
      character*19  State      (nr__) ! Configuration
      integer*4     v2J_1      (nr__) ! [2/10] 2J+1 value
      real*8        E_NIST     (nr__) ! (cm-1) [0/114834] Experimental energy from NIST
*                                   database
      integer*4     E0         (nr__) ! (cm-1) [0/114759] our HF data with a complete
*                                   CI expansion using TRO (CI_TRO_) (1)
      real*4        dE0        (nr__) ! (%) [-1.1/1.9]? Percentage deviations of HF data
*                                   from the experimental energy
      integer*4     E1         (nr__) ! (cm-1) [0/115085] our HF data from the reduced CI
*                                   expansion calculation (CI^red1^_TRO_)
      real*4        dE1        (nr__) ! (%) [-2.8/2.5]? Percentage deviations of the
*                                   reduced CI data from the observed energy
      integer*4     E2         (nr__) ! (cm-1) [0/117659] our HF data from the reduced CI
*                                   expansion calculation (CI^red2^_TRO_)
      real*4        dE2        (nr__) ! (%) [-0.7/6.2]? Percentage deviations of the
*                                   reduced CI data from the observed energy
*Note (1): Using non-orthogonal spline Configuration Interaction (CI),
*          multiconfiguration Hartree-Fock (HF) method, and
*          Transformed Radial Orbital (TRO). Please see Section 2 in the paper
*          for all the details about the CI_HF+TRO_ approximation.

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

C  Declarations for 'table2.dat'	! Transition line strengths S for ZnII determined
                              in the CI_TRO_ approximation

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

      character*2   Type       (nr__1) ! Transition line type: E1/E2/E3 or M1/M2
      integer*4     N0         (nr__1) ! [1/26] Lower level index
      integer*4     N1         (nr__1) ! [2/27] Upper level index
      real*4        S          (nr__1) ! Transition line strength in atomic units

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

C  Declarations for 'table4.dat'	! Effective collision strengths for the
                              electron-impact excitation of ZnII at 14
                              selected temperatures determined in the
                              plane-wave Born approximation

      integer*4 nr__2
      parameter (nr__2=345)	! Number of records
      character*147 ar__2  	! Full-size record

      integer*4     N0_1       (nr__2) ! [1/26] Lower level index
      integer*4     N1_1       (nr__2) ! [2/27] Upper level index
      real*4        CS400      (nr__2) ! Effective collision strength at T=4.0x10^2^K
      real*4        CS800      (nr__2) ! Effective collision strength at T=8.0x10^2^K
      real*4        CS2k       (nr__2) ! Effective collision strength at T=2.0x10^3^K
      real*4        CS4k       (nr__2) ! Effective collision strength at T=4.0x10^3^K
      real*4        CS8k       (nr__2) ! Effective collision strength at T=8.0x10^3^K
      real*4        CS20k      (nr__2) ! Effective collision strength at T=2.0x10^4^K
      real*4        CS40k      (nr__2) ! Effective collision strength at T=4.0x10^4^K
      real*4        CS80k      (nr__2) ! Effective collision strength at T=8.0x10^4^K
      real*4        CS200k     (nr__2) ! Effective collision strength at T=2.0x10^5^K
      real*4        CS400k     (nr__2) ! Effective collision strength at T=4.0x10^5^K
      real*4        CS800k     (nr__2) ! Effective collision strength at T=8.0x10^5^K
      real*4        CS2M       (nr__2) ! Effective collision strength at T=2.0x10^6^K
      real*4        CS4M       (nr__2) ! Effective collision strength at T=4.0x10^6^K
      real*4        CS8M       (nr__2) ! Effective collision strength at T=8.0x10^6^K

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

C  Loading file 'table1.dat'	! Comparison of calculated ZnII level energies
*                              and their percentage deviations with experimental
*                              data from the NIST database

C  Format for file interpretation

    1 format(
     +  I2,1X,A19,1X,I2,1X,F9.2,1X,I6,1X,F5.2,1X,I6,1X,F5.2,1X,I6,1X,
     +  F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,27
        read(1,'(A74)')ar__
        read(ar__,1)
     +  N(i__),State(i__),v2J_1(i__),E_NIST(i__),E0(i__),dE0(i__),
     +  E1(i__),dE1(i__),E2(i__),dE2(i__)
        if(ar__(44:48) .EQ. '') dE0(i__) = rNULL__
        if(ar__(57:61) .EQ. '') dE1(i__) = rNULL__
        if(ar__(70:74) .EQ. '') dE2(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  N(i__),State(i__),v2J_1(i__),E_NIST(i__),E0(i__),dE0(i__),
     +  E1(i__),dE1(i__),E2(i__),dE2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Transition line strengths S for ZnII determined
*                              in the CI_TRO_ approximation

C  Format for file interpretation

    2 format(2X,A2,1X,I2,2X,I2,1X,E9.4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,606
        read(1,'(A21)')ar__1
        read(ar__1,2)Type(i__),N0(i__),N1(i__),S(i__)
c    ..............Just test output...........
        write(6,2)Type(i__),N0(i__),N1(i__),S(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Effective collision strengths for the
*                              electron-impact excitation of ZnII at 14
*                              selected temperatures determined in the
*                              plane-wave Born approximation

C  Format for file interpretation

    3 format(
     +  I3,I3,1X,E10.3,E10.3,E10.3,E10.3,E10.3,E10.3,E10.3,E10.3,
     +  E10.3,E10.3,E10.3,E10.3,E10.3,E10.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,345
        read(1,'(A147)')ar__2
        read(ar__2,3)
     +  N0_1(i__),N1_1(i__),CS400(i__),CS800(i__),CS2k(i__),CS4k(i__),
     +  CS8k(i__),CS20k(i__),CS40k(i__),CS80k(i__),CS200k(i__),
     +  CS400k(i__),CS800k(i__),CS2M(i__),CS4M(i__),CS8M(i__)
c    ..............Just test output...........
        write(6,3)
     +  N0_1(i__),N1_1(i__),CS400(i__),CS800(i__),CS2k(i__),CS4k(i__),
     +  CS8k(i__),CS20k(i__),CS40k(i__),CS80k(i__),CS200k(i__),
     +  CS400k(i__),CS800k(i__),CS2M(i__),CS4M(i__),CS8M(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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