FORTRAN Generation
(/./ftp/cats/J/ApJ/831/L3)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/831/L3 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-Mar-28
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/831/L3       NIR spectra of 10 PNe in LMC and SMC       (Mashburn+, 2016)
*================================================================================
*Neutron-capture element abundances in Magellanic Cloud planetary nebulae.
*    Mashburn A.L., Sterling N.C., Madonna S., Dinerstein H.L., Roederer I.U.,
*    Geballe T.R.
*   <Astrophys. J., 831, L3-L3 (2016)>
*   =2016ApJ...831L...3M    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Observing log and nebular properties of
                                 the 10 planetary nebulae (PNe)

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

      character*3   Gal        (nr__) ! Galaxy name
      integer*4     SMP        (nr__) ! PN name (<SMP {S/L}MC NN> in Simbad)
      character*10  Date       (nr__) ! ("Y/M/D") Observed date
      character*5   Inst       (nr__) ! Instrument (FIRE or GNRIS,
*                                  see the "Description" section above)
      integer*4     Tint       (nr__) ! (s) [1120/8000] Integration time
      real*4        TeO3       (nr__) ! (10+3K) [11.7/16]? T_e_[OIII] nebular temperature
      real*4        ne         (nr__) ! (10+3/cm3) [2.2/31.4]? Electron density
      real*4        CHb        (nr__) ! [0/0.5]? H{beta} extinction coefficient
      real*4        m5007      (nr__) ! [14.6/16.2]? Apparent [OIII] 5007 magnitude (1)
      real*4        logC_O     (nr__) ! [-0.9/0.7]? C/O abundance
      character*7   r_logC_O   (nr__) ! Reference for C/O abundance (2)
      real*4        logN_O     (nr__) ! ? N/O abundance
      character*4   r_logN_O   (nr__) ! Reference for N/O abundance (2)
      real*4        Teff       (nr__) ! (10+3K) [38.8/150]? Central star temperature
      character*4   r_Teff     (nr__) ! Teff reference (2)
      character*7   Minit      (nr__) ! (Msun) Estimated progenitor mass M_init_ (3)
*Note (1): Apparent [OIII] 5007 magnitudes computed from absolute fluxes
*          (corrected for foreground extinction) measured with the Hubble Space
*          Telescope (Stanghellini+ 2003ApJ...596..997S;
*          Shaw+ 2006ApJS..167..201S), with the exceptions of SMP 85 and SMP 99
*          (Leisy & Dennefeld, 2006, J/A+A/456/451), using the relation
*          m_5007_=-2.5logF_5007_-13.74 (Jacoby 1989ApJ...339...39J).
*Note (2): References as follows:
*   S10  = Shaw+ (2010ApJ...717..562S)
*   T03  = Tsamis+ (2003MNRAS.345..186T)
*   MD91 = Meatheringham & Dopita (1991, J/ApJS/76/1085)
*   LD06 = Leisy & Dennefeld (2006, J/A+A/456/451)
*   A87  = Aller+ 1987ApJ...320..159A
*   D94  = Dopita+ 1994ApJ...426..150D
*   S09  = Stanghellini+ 2009ApJ...702..733S
*   V03  = Villaver+ 2003ApJ...597..298V
*   V04  = Villaver+ 2004ApJ...614..716V
*Note (3): from Villaver+ (2003ApJ...597..298V, 2004ApJ...614..716V)

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

C  Declarations for 'table2.dat'	! Line identifications and intensities

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

      character*3   Gal_1      (nr__1) ! Galaxy name (SMC or LMC)
      integer*4     SMP_1      (nr__1) ! Planetary Nebular Name
*                                   (<SMP {S/L}MC NN> in Simbad)
      character*5   Inst_1     (nr__1) ! Spectrometer used (1)
      character*12  Line       (nr__1) ! Line Identification, ion
      real*8        lamobs     (nr__1) ! (um) [0.831/2.462]? Observed wavelength, microns
      real*8        lamlab     (nr__1) ! (um) [0.830/2.455] Laboratory wavelength, microns
      character*2   l_F_F_Brg  (nr__1) ! [<= ] Limit flag on F/F(Brg)
      real*4        F_F_Brg    (nr__1) ! [0.07/1410]? Measured Line Flux,
*                                   normalized (2)
      real*4        e_F_F_Brg  (nr__1) ! [0.01/54]? Uncertainty in F/F(Brg)
      character*2   l_I_I_Brg  (nr__1) ! [<= ] Limit flag on I/I(Brg)
      real*4        I_I_Brg    (nr__1) ! [0.08/1570]? Measured Line Intensity,
*                                   normalized (2)
      real*4        e_I_I_Brg  (nr__1) ! [0.01/62]? Uncertainty in I/I(Brg)
      character*5   Comm       (nr__1) ! Comment (3)
*Note (1): Spectrometer used as follows:
*    FIRE  = Folded-Port InfraRed Echellette (FIRE), 6.5-m Baade Telescope,
*             Las Campanas Observatory (823 instances);
*    GNIRS = Gemini Near-InfraRed Spectrograph (GNIRS), 8.1-m Gemini South
*             telescope (24 instances).
*Note (2): Measured fluxes and intensities are on the scale
*          F(HI, Br{gamma}) = I(HI, Br{gamma}) = 100
*Note (3): Comment as follows:
*    :   = Marginal detection;
*    Fit = Line fit;
*    *   = Blended feature.

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

C  Loading file 'table1.dat'	! Observing log and nebular properties of
*                                 the 10 planetary nebulae (PNe)

C  Format for file interpretation

    1 format(
     +  A3,5X,I2,1X,A10,1X,A5,1X,I4,1X,F4.1,1X,F5.2,1X,F4.2,1X,F5.2,
     +  1X,F5.2,1X,A7,1X,F5.2,1X,A4,1X,F5.1,1X,A4,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,11
        read(1,'(A98)')ar__
        read(ar__,1)
     +  Gal(i__),SMP(i__),Date(i__),Inst(i__),Tint(i__),TeO3(i__),
     +  ne(i__),CHb(i__),m5007(i__),logC_O(i__),r_logC_O(i__),
     +  logN_O(i__),r_logN_O(i__),Teff(i__),r_Teff(i__),Minit(i__)
        if(ar__(34:37) .EQ. '') TeO3(i__) = rNULL__
        if(ar__(39:43) .EQ. '') ne(i__) = rNULL__
        if(ar__(45:48) .EQ. '') CHb(i__) = rNULL__
        if(ar__(50:54) .EQ. '') m5007(i__) = rNULL__
        if(ar__(56:60) .EQ. '') logC_O(i__) = rNULL__
        if(ar__(70:74) .EQ. '') logN_O(i__) = rNULL__
        if(ar__(81:85) .EQ. '') Teff(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Gal(i__),SMP(i__),Date(i__),Inst(i__),Tint(i__),TeO3(i__),
     +  ne(i__),CHb(i__),m5007(i__),logC_O(i__),r_logC_O(i__),
     +  logN_O(i__),r_logN_O(i__),Teff(i__),r_Teff(i__),Minit(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Line identifications and intensities

C  Format for file interpretation

    2 format(
     +  A3,5X,I2,2X,A5,2X,A12,1X,F7.5,F7.5,1X,A2,2X,E8.3,1X,E8.3,2X,
     +  A2,2X,E8.3,1X,E8.3,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,847
        read(1,'(A97)')ar__1
        read(ar__1,2)
     +  Gal_1(i__),SMP_1(i__),Inst_1(i__),Line(i__),lamobs(i__),
     +  lamlab(i__),l_F_F_Brg(i__),F_F_Brg(i__),e_F_F_Brg(i__),
     +  l_I_I_Brg(i__),I_I_Brg(i__),e_I_I_Brg(i__),Comm(i__)
        if(ar__1(33:39) .EQ. '') lamobs(i__) = rNULL__
        if(ar__1(52:59) .EQ. '') F_F_Brg(i__) = rNULL__
        if(ar__1(61:68) .EQ. '') e_F_F_Brg(i__) = rNULL__
        if(ar__1(75:82) .EQ. '') I_I_Brg(i__) = rNULL__
        if(ar__1(84:91) .EQ. '') e_I_I_Brg(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Gal_1(i__),SMP_1(i__),Inst_1(i__),Line(i__),lamobs(i__),
     +  lamlab(i__),l_F_F_Brg(i__),F_F_Brg(i__),e_F_F_Brg(i__),
     +  l_I_I_Brg(i__),I_I_Brg(i__),e_I_I_Brg(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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