FORTRAN Generation
(/./ftp/cats/J/A_A/494/591)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/494/591 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.5, on 2013-May-24
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__=-1.e37)     	! NULL real number
      parameter  (iNULL__=-2147483647)	! NULL int  number

C=============================================================================
Cat. J/A+A/494/591        Abundances in planetary nebulae         (Chiappini+, 2009)
*================================================================================
*Abundances in the Galactic bulge: results from planetary nebulae and
*giant stars.
*    Chiappini C., Gorny S.K., Stasinska G., Barbuy B.
*   <Astron. Astrophys., 494, 591-610 (2009)>
*   =2009A&A...494..591C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1a.dat'	! Plasma diagnostics and abundances of
                             Galactic bulge planetary nebulae

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

      character*11  PNG        (nr__) ! Planetary Nebula PNG Name (see Cat. <V/84>)
*                                   (for tables 1a and 1b only)
      character*11  Name       (nr__) ! Usual name
      real*4        Ne_SII_    (nr__) ! (cm-3) ? Electron density from [SII] 6731/6717
      character*1   q_Ne_SII_  (nr__) ! Quality flag on Ne[SII] (1)
      integer*4     Te_OIII_   (nr__) ! (K) ? Electron temperature from [OIII] 4363/5007
      character*1   q_Te_OIII_ (nr__) ! Quality flag on Te[OIII] (1)
      integer*4     Te_NII_    (nr__) ! (K) ? Electron temperature from [NII] 5755/6584
      character*1   q_Te_NII_  (nr__) ! Quality flag on Te[NII] (1)
      real*4        He_H       (nr__) ! ? He/H ratio
      character*1   q_He_H     (nr__) ! Quality flag on He/H (1)
      real*4        N_H        (nr__) ! ? N/H ratio
      character*1   q_N_H      (nr__) ! Quality flag on N/H (1)
      real*4        O_H        (nr__) ! ? O/H ratio
      character*1   q_O_H      (nr__) ! Quality flag on O/H (1)
      real*4        Ne_H       (nr__) ! ? Ne/H ratio
      character*1   q_Ne_H     (nr__) ! Quality flag on Ne/H (1)
      real*4        S_H        (nr__) ! ? S/H ratio
      character*1   q_S_H      (nr__) ! Quality flag on S/H (1)
      real*4        Ar_H       (nr__) ! ? Ar/H ratio
      character*1   q_Ar_H     (nr__) ! Quality flag on Ar/H (1)
      real*4        Cl_H       (nr__) ! ? Cl/H ratio
      character*1   q_Cl_H     (nr__) ! Quality flag on Cl/H (1)
      integer*4     Note       (nr__) ! [1]? 1 = O[III] temperature used for all ions
      real*4        E_Ne_SII_  (nr__) ! (cm-3) ? Error on Ne[SII] (upper limit)
      integer*4     E_Te_OIII_ (nr__) ! (K) ? Error on Te[OIII] (upper limit)
      integer*4     E_Te_NII_  (nr__) ! (K) ? Error on Te[NII] (upper limit)
      real*4        E_He_H     (nr__) ! ? Error on He/H (upper limit)
      real*4        E_N_H      (nr__) ! ? Error on N/H (upper limit)
      real*4        E_O_H      (nr__) ! ? Error on O/H (upper limit)
      real*4        E_Ne_H     (nr__) ! ? Error on Ne/H (upper limit)
      real*4        E_S_H      (nr__) ! ? Error on S/H (upper limit)
      real*4        E_Ar_H     (nr__) ! ? Error on Ar/H (upper limit)
      real*4        E_Cl_H     (nr__) ! ? Error on Cl/H (upper limit)
      real*4        e_Ne_SII__1(nr__) ! (cm-3) ? Error on Ne[SII] (lower limit)
      integer*4     e_Te_OIII__1(nr__) ! (K) ? Error on Te[OIII] (lower limit)
      integer*4     e_Te_NII__1(nr__) ! (K) ? Error on Te[NII] (lower limit)
      real*4        e_He_H_1   (nr__) ! ? Error on He/H (lower limit)
      real*4        e_N_H_1    (nr__) ! ? Error on N/H (lower limit)
      real*4        e_O_H_1    (nr__) ! ? Error on O/H (lower limit)
      real*4        e_Ne_H_1   (nr__) ! ? Error on Ne/H (lower limit)
      real*4        e_S_H_1    (nr__) ! ? Error on S/H (lower limit)
      real*4        e_Ar_H_1   (nr__) ! ? Error on Ar/H (lower limit)
      real*4        e_Cl_H_1   (nr__) ! ? Error on Cl/H (lower limit)
*Note (1): Quality flags on derived parameters are defined as follows:
*      + = error is smaller than the computed median error of the parameter
*          for the entire PNe sample
*      : = larger uncertainty but error not exceeding the 0.3dex limit
*      ; = data rejected due to error superior to 0.3dex or large
*          uncertainties arising from the ICFs
*      - = parameter could not be derived

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

C  Declarations for 'table1b.dat'	! Plasma diagnostics and abundances of
                             Galactic inner disk planetary nebulae

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

      character*11  PNG_1      (nr__1) ! Planetary Nebula PNG Name (see Cat. <V/84>)
*                                   (for tables 1a and 1b only)
      character*11  Name_1     (nr__1) ! Usual name
      real*4        Ne_SII__1  (nr__1) ! (cm-3) ? Electron density from [SII] 6731/6717
      character*1   q_Ne_SII__1(nr__1) ! Quality flag on Ne[SII] (1)
      integer*4     Te_OIII__1 (nr__1) ! (K) ? Electron temperature from [OIII] 4363/5007
      character*1   q_Te_OIII__1(nr__1) ! Quality flag on Te[OIII] (1)
      integer*4     Te_NII__1  (nr__1) ! (K) ? Electron temperature from [NII] 5755/6584
      character*1   q_Te_NII__1(nr__1) ! Quality flag on Te[NII] (1)
      real*4        He_H_1     (nr__1) ! ? He/H ratio
      character*1   q_He_H_1   (nr__1) ! Quality flag on He/H (1)
      real*4        N_H_1      (nr__1) ! ? N/H ratio
      character*1   q_N_H_1    (nr__1) ! Quality flag on N/H (1)
      real*4        O_H_1      (nr__1) ! ? O/H ratio
      character*1   q_O_H_1    (nr__1) ! Quality flag on O/H (1)
      real*4        Ne_H_1     (nr__1) ! ? Ne/H ratio
      character*1   q_Ne_H_1   (nr__1) ! Quality flag on Ne/H (1)
      real*4        S_H_1      (nr__1) ! ? S/H ratio
      character*1   q_S_H_1    (nr__1) ! Quality flag on S/H (1)
      real*4        Ar_H_1     (nr__1) ! ? Ar/H ratio
      character*1   q_Ar_H_1   (nr__1) ! Quality flag on Ar/H (1)
      real*4        Cl_H_1     (nr__1) ! ? Cl/H ratio
      character*1   q_Cl_H_1   (nr__1) ! Quality flag on Cl/H (1)
      integer*4     Note_1     (nr__1) ! [1]? 1 = O[III] temperature used for all ions
      real*4        E_Ne_SII__2(nr__1) ! (cm-3) ? Error on Ne[SII] (upper limit)
      integer*4     E_Te_OIII__2(nr__1) ! (K) ? Error on Te[OIII] (upper limit)
      integer*4     E_Te_NII__2(nr__1) ! (K) ? Error on Te[NII] (upper limit)
      real*4        E_He_H_2   (nr__1) ! ? Error on He/H (upper limit)
      real*4        E_N_H_2    (nr__1) ! ? Error on N/H (upper limit)
      real*4        E_O_H_2    (nr__1) ! ? Error on O/H (upper limit)
      real*4        E_Ne_H_2   (nr__1) ! ? Error on Ne/H (upper limit)
      real*4        E_S_H_2    (nr__1) ! ? Error on S/H (upper limit)
      real*4        E_Ar_H_2   (nr__1) ! ? Error on Ar/H (upper limit)
      real*4        E_Cl_H_2   (nr__1) ! ? Error on Cl/H (upper limit)
      real*4        e_Ne_SII__3(nr__1) ! (cm-3) ? Error on Ne[SII] (lower limit)
      integer*4     e_Te_OIII__3(nr__1) ! (K) ? Error on Te[OIII] (lower limit)
      integer*4     e_Te_NII__3(nr__1) ! (K) ? Error on Te[NII] (lower limit)
      real*4        e_He_H_3   (nr__1) ! ? Error on He/H (lower limit)
      real*4        e_N_H_3    (nr__1) ! ? Error on N/H (lower limit)
      real*4        e_O_H_3    (nr__1) ! ? Error on O/H (lower limit)
      real*4        e_Ne_H_3   (nr__1) ! ? Error on Ne/H (lower limit)
      real*4        e_S_H_3    (nr__1) ! ? Error on S/H (lower limit)
      real*4        e_Ar_H_3   (nr__1) ! ? Error on Ar/H (lower limit)
      real*4        e_Cl_H_3   (nr__1) ! ? Error on Cl/H (lower limit)
*Note (1): Quality flags on derived parameters are defined as follows:
*      + = error is smaller than the computed median error of the parameter
*          for the entire PNe sample
*      : = larger uncertainty but error not exceeding the 0.3dex limit
*      ; = data rejected due to error superior to 0.3dex or large
*          uncertainties arising from the ICFs
*      - = parameter could not be derived

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

C  Declarations for 'table1c.dat'	! Plasma diagnostics and abundances of
                             LMC planetary nebulae

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

      character*11  PNG_2      (nr__2) ! Planetary Nebula PNG Name (see Cat. <V/84>)
*                                   (for tables 1a and 1b only)
      character*11  Name_2     (nr__2) ! Usual name
      real*4        Ne_SII__2  (nr__2) ! (cm-3) ? Electron density from [SII] 6731/6717
      character*1   q_Ne_SII__2(nr__2) ! Quality flag on Ne[SII] (1)
      integer*4     Te_OIII__2 (nr__2) ! (K) ? Electron temperature from [OIII] 4363/5007
      character*1   q_Te_OIII__2(nr__2) ! Quality flag on Te[OIII] (1)
      integer*4     Te_NII__2  (nr__2) ! (K) ? Electron temperature from [NII] 5755/6584
      character*1   q_Te_NII__2(nr__2) ! Quality flag on Te[NII] (1)
      real*4        He_H_2     (nr__2) ! ? He/H ratio
      character*1   q_He_H_2   (nr__2) ! Quality flag on He/H (1)
      real*4        N_H_2      (nr__2) ! ? N/H ratio
      character*1   q_N_H_2    (nr__2) ! Quality flag on N/H (1)
      real*4        O_H_2      (nr__2) ! ? O/H ratio
      character*1   q_O_H_2    (nr__2) ! Quality flag on O/H (1)
      real*4        Ne_H_2     (nr__2) ! ? Ne/H ratio
      character*1   q_Ne_H_2   (nr__2) ! Quality flag on Ne/H (1)
      real*4        S_H_2      (nr__2) ! ? S/H ratio
      character*1   q_S_H_2    (nr__2) ! Quality flag on S/H (1)
      real*4        Ar_H_2     (nr__2) ! ? Ar/H ratio
      character*1   q_Ar_H_2   (nr__2) ! Quality flag on Ar/H (1)
      real*4        Cl_H_2     (nr__2) ! ? Cl/H ratio
      character*1   q_Cl_H_2   (nr__2) ! Quality flag on Cl/H (1)
      integer*4     Note_2     (nr__2) ! [1]? 1 = O[III] temperature used for all ions
      real*4        E_Ne_SII__4(nr__2) ! (cm-3) ? Error on Ne[SII] (upper limit)
      integer*4     E_Te_OIII__4(nr__2) ! (K) ? Error on Te[OIII] (upper limit)
      integer*4     E_Te_NII__4(nr__2) ! (K) ? Error on Te[NII] (upper limit)
      real*4        E_He_H_4   (nr__2) ! ? Error on He/H (upper limit)
      real*4        E_N_H_4    (nr__2) ! ? Error on N/H (upper limit)
      real*4        E_O_H_4    (nr__2) ! ? Error on O/H (upper limit)
      real*4        E_Ne_H_4   (nr__2) ! ? Error on Ne/H (upper limit)
      real*4        E_S_H_4    (nr__2) ! ? Error on S/H (upper limit)
      real*4        E_Ar_H_4   (nr__2) ! ? Error on Ar/H (upper limit)
      real*4        E_Cl_H_4   (nr__2) ! ? Error on Cl/H (upper limit)
      real*4        e_Ne_SII__5(nr__2) ! (cm-3) ? Error on Ne[SII] (lower limit)
      integer*4     e_Te_OIII__5(nr__2) ! (K) ? Error on Te[OIII] (lower limit)
      integer*4     e_Te_NII__5(nr__2) ! (K) ? Error on Te[NII] (lower limit)
      real*4        e_He_H_5   (nr__2) ! ? Error on He/H (lower limit)
      real*4        e_N_H_5    (nr__2) ! ? Error on N/H (lower limit)
      real*4        e_O_H_5    (nr__2) ! ? Error on O/H (lower limit)
      real*4        e_Ne_H_5   (nr__2) ! ? Error on Ne/H (lower limit)
      real*4        e_S_H_5    (nr__2) ! ? Error on S/H (lower limit)
      real*4        e_Ar_H_5   (nr__2) ! ? Error on Ar/H (lower limit)
      real*4        e_Cl_H_5   (nr__2) ! ? Error on Cl/H (lower limit)
*Note (1): Quality flags on derived parameters are defined as follows:
*      + = error is smaller than the computed median error of the parameter
*          for the entire PNe sample
*      : = larger uncertainty but error not exceeding the 0.3dex limit
*      ; = data rejected due to error superior to 0.3dex or large
*          uncertainties arising from the ICFs
*      - = parameter could not be derived

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

C  Loading file 'table1a.dat'	! Plasma diagnostics and abundances of
*                             Galactic bulge planetary nebulae

C  Format for file interpretation

    1 format(
     +  A11,1X,A11,1X,E8.3,A1,1X,I5,A1,1X,I5,A1,1X,E8.3,A1,1X,E8.3,A1,
     +  1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,I1,
     +  1X,E8.3,1X,I5,1X,I5,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,I5,1X,I5,1X,E8.3,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)

C  Effective file loading

      open(unit=1,file='table1a.dat', status='old')
      write(6,*) '....Loading file: table1a.dat'
      do i__=1,168
        read(1,'(A287)')ar__
        read(ar__,1)
     +  PNG(i__),Name(i__),Ne_SII_(i__),q_Ne_SII_(i__),Te_OIII_(i__),
     +  q_Te_OIII_(i__),Te_NII_(i__),q_Te_NII_(i__),He_H(i__),
     +  q_He_H(i__),N_H(i__),q_N_H(i__),O_H(i__),q_O_H(i__),Ne_H(i__),
     +  q_Ne_H(i__),S_H(i__),q_S_H(i__),Ar_H(i__),q_Ar_H(i__),
     +  Cl_H(i__),q_Cl_H(i__),Note(i__),E_Ne_SII_(i__),
     +  E_Te_OIII_(i__),E_Te_NII_(i__),E_He_H(i__),E_N_H(i__),
     +  E_O_H(i__),E_Ne_H(i__),E_S_H(i__),E_Ar_H(i__),E_Cl_H(i__),
     +  e_Ne_SII__1(i__),e_Te_OIII__1(i__),e_Te_NII__1(i__),
     +  e_He_H_1(i__),e_N_H_1(i__),e_O_H_1(i__),e_Ne_H_1(i__),
     +  e_S_H_1(i__),e_Ar_H_1(i__),e_Cl_H_1(i__)
        if(ar__(25:32) .EQ. '') Ne_SII_(i__) = rNULL__
        if(ar__(35:39) .EQ. '') Te_OIII_(i__) = iNULL__
        if(ar__(42:46) .EQ. '') Te_NII_(i__) = iNULL__
        if(ar__(49:56) .EQ. '') He_H(i__) = rNULL__
        if(ar__(59:66) .EQ. '') N_H(i__) = rNULL__
        if(ar__(69:76) .EQ. '') O_H(i__) = rNULL__
        if(ar__(79:86) .EQ. '') Ne_H(i__) = rNULL__
        if(ar__(89:96) .EQ. '') S_H(i__) = rNULL__
        if(ar__(99:106) .EQ. '') Ar_H(i__) = rNULL__
        if(ar__(109:116) .EQ. '') Cl_H(i__) = rNULL__
        if(ar__(119:119) .EQ. '') Note(i__) = iNULL__
        if(ar__(121:128) .EQ. '') E_Ne_SII_(i__) = rNULL__
        if(ar__(130:134) .EQ. '') E_Te_OIII_(i__) = iNULL__
        if(ar__(136:140) .EQ. '') E_Te_NII_(i__) = iNULL__
        if(ar__(142:149) .EQ. '') E_He_H(i__) = rNULL__
        if(ar__(151:158) .EQ. '') E_N_H(i__) = rNULL__
        if(ar__(160:167) .EQ. '') E_O_H(i__) = rNULL__
        if(ar__(169:176) .EQ. '') E_Ne_H(i__) = rNULL__
        if(ar__(178:185) .EQ. '') E_S_H(i__) = rNULL__
        if(ar__(187:194) .EQ. '') E_Ar_H(i__) = rNULL__
        if(ar__(196:203) .EQ. '') E_Cl_H(i__) = rNULL__
        if(ar__(205:212) .EQ. '') e_Ne_SII__1(i__) = rNULL__
        if(ar__(214:218) .EQ. '') e_Te_OIII__1(i__) = iNULL__
        if(ar__(220:224) .EQ. '') e_Te_NII__1(i__) = iNULL__
        if(ar__(226:233) .EQ. '') e_He_H_1(i__) = rNULL__
        if(ar__(235:242) .EQ. '') e_N_H_1(i__) = rNULL__
        if(ar__(244:251) .EQ. '') e_O_H_1(i__) = rNULL__
        if(ar__(253:260) .EQ. '') e_Ne_H_1(i__) = rNULL__
        if(ar__(262:269) .EQ. '') e_S_H_1(i__) = rNULL__
        if(ar__(271:278) .EQ. '') e_Ar_H_1(i__) = rNULL__
        if(ar__(280:287) .EQ. '') e_Cl_H_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  PNG(i__),Name(i__),Ne_SII_(i__),q_Ne_SII_(i__),Te_OIII_(i__),
     +  q_Te_OIII_(i__),Te_NII_(i__),q_Te_NII_(i__),He_H(i__),
     +  q_He_H(i__),N_H(i__),q_N_H(i__),O_H(i__),q_O_H(i__),Ne_H(i__),
     +  q_Ne_H(i__),S_H(i__),q_S_H(i__),Ar_H(i__),q_Ar_H(i__),
     +  Cl_H(i__),q_Cl_H(i__),Note(i__),E_Ne_SII_(i__),
     +  E_Te_OIII_(i__),E_Te_NII_(i__),E_He_H(i__),E_N_H(i__),
     +  E_O_H(i__),E_Ne_H(i__),E_S_H(i__),E_Ar_H(i__),E_Cl_H(i__),
     +  e_Ne_SII__1(i__),e_Te_OIII__1(i__),e_Te_NII__1(i__),
     +  e_He_H_1(i__),e_N_H_1(i__),e_O_H_1(i__),e_Ne_H_1(i__),
     +  e_S_H_1(i__),e_Ar_H_1(i__),e_Cl_H_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table1b.dat'	! Plasma diagnostics and abundances of
*                             Galactic inner disk planetary nebulae

C  Format for file interpretation

    2 format(
     +  A11,1X,A11,1X,E8.3,A1,1X,I5,A1,1X,I5,A1,1X,E8.3,A1,1X,E8.3,A1,
     +  1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,I1,
     +  1X,E8.3,1X,I5,1X,I5,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,I5,1X,I5,1X,E8.3,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)

C  Effective file loading

      open(unit=1,file='table1b.dat', status='old')
      write(6,*) '....Loading file: table1b.dat'
      do i__=1,61
        read(1,'(A287)')ar__1
        read(ar__1,2)
     +  PNG_1(i__),Name_1(i__),Ne_SII__1(i__),q_Ne_SII__1(i__),
     +  Te_OIII__1(i__),q_Te_OIII__1(i__),Te_NII__1(i__),
     +  q_Te_NII__1(i__),He_H_1(i__),q_He_H_1(i__),N_H_1(i__),
     +  q_N_H_1(i__),O_H_1(i__),q_O_H_1(i__),Ne_H_1(i__),
     +  q_Ne_H_1(i__),S_H_1(i__),q_S_H_1(i__),Ar_H_1(i__),
     +  q_Ar_H_1(i__),Cl_H_1(i__),q_Cl_H_1(i__),Note_1(i__),
     +  E_Ne_SII__2(i__),E_Te_OIII__2(i__),E_Te_NII__2(i__),
     +  E_He_H_2(i__),E_N_H_2(i__),E_O_H_2(i__),E_Ne_H_2(i__),
     +  E_S_H_2(i__),E_Ar_H_2(i__),E_Cl_H_2(i__),e_Ne_SII__3(i__),
     +  e_Te_OIII__3(i__),e_Te_NII__3(i__),e_He_H_3(i__),e_N_H_3(i__),
     +  e_O_H_3(i__),e_Ne_H_3(i__),e_S_H_3(i__),e_Ar_H_3(i__),
     +  e_Cl_H_3(i__)
        if(ar__1(25:32) .EQ. '') Ne_SII__1(i__) = rNULL__
        if(ar__1(35:39) .EQ. '') Te_OIII__1(i__) = iNULL__
        if(ar__1(42:46) .EQ. '') Te_NII__1(i__) = iNULL__
        if(ar__1(49:56) .EQ. '') He_H_1(i__) = rNULL__
        if(ar__1(59:66) .EQ. '') N_H_1(i__) = rNULL__
        if(ar__1(69:76) .EQ. '') O_H_1(i__) = rNULL__
        if(ar__1(79:86) .EQ. '') Ne_H_1(i__) = rNULL__
        if(ar__1(89:96) .EQ. '') S_H_1(i__) = rNULL__
        if(ar__1(99:106) .EQ. '') Ar_H_1(i__) = rNULL__
        if(ar__1(109:116) .EQ. '') Cl_H_1(i__) = rNULL__
        if(ar__1(119:119) .EQ. '') Note_1(i__) = iNULL__
        if(ar__1(121:128) .EQ. '') E_Ne_SII__2(i__) = rNULL__
        if(ar__1(130:134) .EQ. '') E_Te_OIII__2(i__) = iNULL__
        if(ar__1(136:140) .EQ. '') E_Te_NII__2(i__) = iNULL__
        if(ar__1(142:149) .EQ. '') E_He_H_2(i__) = rNULL__
        if(ar__1(151:158) .EQ. '') E_N_H_2(i__) = rNULL__
        if(ar__1(160:167) .EQ. '') E_O_H_2(i__) = rNULL__
        if(ar__1(169:176) .EQ. '') E_Ne_H_2(i__) = rNULL__
        if(ar__1(178:185) .EQ. '') E_S_H_2(i__) = rNULL__
        if(ar__1(187:194) .EQ. '') E_Ar_H_2(i__) = rNULL__
        if(ar__1(196:203) .EQ. '') E_Cl_H_2(i__) = rNULL__
        if(ar__1(205:212) .EQ. '') e_Ne_SII__3(i__) = rNULL__
        if(ar__1(214:218) .EQ. '') e_Te_OIII__3(i__) = iNULL__
        if(ar__1(220:224) .EQ. '') e_Te_NII__3(i__) = iNULL__
        if(ar__1(226:233) .EQ. '') e_He_H_3(i__) = rNULL__
        if(ar__1(235:242) .EQ. '') e_N_H_3(i__) = rNULL__
        if(ar__1(244:251) .EQ. '') e_O_H_3(i__) = rNULL__
        if(ar__1(253:260) .EQ. '') e_Ne_H_3(i__) = rNULL__
        if(ar__1(262:269) .EQ. '') e_S_H_3(i__) = rNULL__
        if(ar__1(271:278) .EQ. '') e_Ar_H_3(i__) = rNULL__
        if(ar__1(280:287) .EQ. '') e_Cl_H_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  PNG_1(i__),Name_1(i__),Ne_SII__1(i__),q_Ne_SII__1(i__),
     +  Te_OIII__1(i__),q_Te_OIII__1(i__),Te_NII__1(i__),
     +  q_Te_NII__1(i__),He_H_1(i__),q_He_H_1(i__),N_H_1(i__),
     +  q_N_H_1(i__),O_H_1(i__),q_O_H_1(i__),Ne_H_1(i__),
     +  q_Ne_H_1(i__),S_H_1(i__),q_S_H_1(i__),Ar_H_1(i__),
     +  q_Ar_H_1(i__),Cl_H_1(i__),q_Cl_H_1(i__),Note_1(i__),
     +  E_Ne_SII__2(i__),E_Te_OIII__2(i__),E_Te_NII__2(i__),
     +  E_He_H_2(i__),E_N_H_2(i__),E_O_H_2(i__),E_Ne_H_2(i__),
     +  E_S_H_2(i__),E_Ar_H_2(i__),E_Cl_H_2(i__),e_Ne_SII__3(i__),
     +  e_Te_OIII__3(i__),e_Te_NII__3(i__),e_He_H_3(i__),e_N_H_3(i__),
     +  e_O_H_3(i__),e_Ne_H_3(i__),e_S_H_3(i__),e_Ar_H_3(i__),
     +  e_Cl_H_3(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table1c.dat'	! Plasma diagnostics and abundances of
*                             LMC planetary nebulae

C  Format for file interpretation

    3 format(
     +  A11,1X,A11,1X,E8.3,A1,1X,I5,A1,1X,I5,A1,1X,E8.3,A1,1X,E8.3,A1,
     +  1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,E8.3,A1,1X,I1,
     +  1X,E8.3,1X,I5,1X,I5,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,I5,1X,I5,1X,E8.3,1X,E8.3,1X,E8.3,
     +  1X,E8.3,1X,E8.3,1X,E8.3,1X,E8.3)

C  Effective file loading

      open(unit=1,file='table1c.dat', status='old')
      write(6,*) '....Loading file: table1c.dat'
      do i__=1,110
        read(1,'(A287)')ar__2
        read(ar__2,3)
     +  PNG_2(i__),Name_2(i__),Ne_SII__2(i__),q_Ne_SII__2(i__),
     +  Te_OIII__2(i__),q_Te_OIII__2(i__),Te_NII__2(i__),
     +  q_Te_NII__2(i__),He_H_2(i__),q_He_H_2(i__),N_H_2(i__),
     +  q_N_H_2(i__),O_H_2(i__),q_O_H_2(i__),Ne_H_2(i__),
     +  q_Ne_H_2(i__),S_H_2(i__),q_S_H_2(i__),Ar_H_2(i__),
     +  q_Ar_H_2(i__),Cl_H_2(i__),q_Cl_H_2(i__),Note_2(i__),
     +  E_Ne_SII__4(i__),E_Te_OIII__4(i__),E_Te_NII__4(i__),
     +  E_He_H_4(i__),E_N_H_4(i__),E_O_H_4(i__),E_Ne_H_4(i__),
     +  E_S_H_4(i__),E_Ar_H_4(i__),E_Cl_H_4(i__),e_Ne_SII__5(i__),
     +  e_Te_OIII__5(i__),e_Te_NII__5(i__),e_He_H_5(i__),e_N_H_5(i__),
     +  e_O_H_5(i__),e_Ne_H_5(i__),e_S_H_5(i__),e_Ar_H_5(i__),
     +  e_Cl_H_5(i__)
        if(ar__2(25:32) .EQ. '') Ne_SII__2(i__) = rNULL__
        if(ar__2(35:39) .EQ. '') Te_OIII__2(i__) = iNULL__
        if(ar__2(42:46) .EQ. '') Te_NII__2(i__) = iNULL__
        if(ar__2(49:56) .EQ. '') He_H_2(i__) = rNULL__
        if(ar__2(59:66) .EQ. '') N_H_2(i__) = rNULL__
        if(ar__2(69:76) .EQ. '') O_H_2(i__) = rNULL__
        if(ar__2(79:86) .EQ. '') Ne_H_2(i__) = rNULL__
        if(ar__2(89:96) .EQ. '') S_H_2(i__) = rNULL__
        if(ar__2(99:106) .EQ. '') Ar_H_2(i__) = rNULL__
        if(ar__2(109:116) .EQ. '') Cl_H_2(i__) = rNULL__
        if(ar__2(119:119) .EQ. '') Note_2(i__) = iNULL__
        if(ar__2(121:128) .EQ. '') E_Ne_SII__4(i__) = rNULL__
        if(ar__2(130:134) .EQ. '') E_Te_OIII__4(i__) = iNULL__
        if(ar__2(136:140) .EQ. '') E_Te_NII__4(i__) = iNULL__
        if(ar__2(142:149) .EQ. '') E_He_H_4(i__) = rNULL__
        if(ar__2(151:158) .EQ. '') E_N_H_4(i__) = rNULL__
        if(ar__2(160:167) .EQ. '') E_O_H_4(i__) = rNULL__
        if(ar__2(169:176) .EQ. '') E_Ne_H_4(i__) = rNULL__
        if(ar__2(178:185) .EQ. '') E_S_H_4(i__) = rNULL__
        if(ar__2(187:194) .EQ. '') E_Ar_H_4(i__) = rNULL__
        if(ar__2(196:203) .EQ. '') E_Cl_H_4(i__) = rNULL__
        if(ar__2(205:212) .EQ. '') e_Ne_SII__5(i__) = rNULL__
        if(ar__2(214:218) .EQ. '') e_Te_OIII__5(i__) = iNULL__
        if(ar__2(220:224) .EQ. '') e_Te_NII__5(i__) = iNULL__
        if(ar__2(226:233) .EQ. '') e_He_H_5(i__) = rNULL__
        if(ar__2(235:242) .EQ. '') e_N_H_5(i__) = rNULL__
        if(ar__2(244:251) .EQ. '') e_O_H_5(i__) = rNULL__
        if(ar__2(253:260) .EQ. '') e_Ne_H_5(i__) = rNULL__
        if(ar__2(262:269) .EQ. '') e_S_H_5(i__) = rNULL__
        if(ar__2(271:278) .EQ. '') e_Ar_H_5(i__) = rNULL__
        if(ar__2(280:287) .EQ. '') e_Cl_H_5(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  PNG_2(i__),Name_2(i__),Ne_SII__2(i__),q_Ne_SII__2(i__),
     +  Te_OIII__2(i__),q_Te_OIII__2(i__),Te_NII__2(i__),
     +  q_Te_NII__2(i__),He_H_2(i__),q_He_H_2(i__),N_H_2(i__),
     +  q_N_H_2(i__),O_H_2(i__),q_O_H_2(i__),Ne_H_2(i__),
     +  q_Ne_H_2(i__),S_H_2(i__),q_S_H_2(i__),Ar_H_2(i__),
     +  q_Ar_H_2(i__),Cl_H_2(i__),q_Cl_H_2(i__),Note_2(i__),
     +  E_Ne_SII__4(i__),E_Te_OIII__4(i__),E_Te_NII__4(i__),
     +  E_He_H_4(i__),E_N_H_4(i__),E_O_H_4(i__),E_Ne_H_4(i__),
     +  E_S_H_4(i__),E_Ar_H_4(i__),E_Cl_H_4(i__),e_Ne_SII__5(i__),
     +  e_Te_OIII__5(i__),e_Te_NII__5(i__),e_He_H_5(i__),e_N_H_5(i__),
     +  e_O_H_5(i__),e_Ne_H_5(i__),e_S_H_5(i__),e_Ar_H_5(i__),
     +  e_Cl_H_5(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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