FORTRAN Generation
(/./ftp/cats/J/A_A/706/A261)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/706/A261 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-May-17
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/706/A261     EELG phys. prop. and ion. efficiency (Gimenez-Alcazar+, 2026)
*================================================================================
*J-PAS: First identification, physical properties and ionization efficiency of
*extreme emission line galaxies.
*    Gimenez-Alcazar A., Amorin R., Vilchez J.M., Hernan-Caballero A.,
*    Gonzalez-Otero M., Arroyo-Polonio A., Iglesias-Paramo J.,
*    Lumbreras-Calle A., Fernandez-Ontiveros J.A., Bonatto L.,
*    Gonzalez Delgado R.M., Kehrig C., Torralba A., Rahna P.T., Jimenez-Teja Y.,
*    Marquez I., Breda I., Alvarez-Candal A., Abramo R., Alcaniz J., Benitez N.,
*    Bonoli S., Carneiro S., Cenarro J., Cristobal-Hornillos D., Dupke R.,
*    Ederoclite A., Hernandez-Monteagudo C., Marin-Franch A.,
*    Mendes de Oliveira C., Moles M., Sodre L.Jr, Taylor K., Varela J.,
*    Vazquez Ramio.H.
*    <Astron. Astrophys. 706, A261 (2026)>
*    =2026A&A...706A.261G        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Derived SED Parameters

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     TileID     (nr__) ! J-PAS tile identifier
      integer*4     Number     (nr__) ! Source number within the tile
      real*8        RAdeg      (nr__) ! (deg) Right Ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination (J2000)
      real*8        z          (nr__) ! Spectroscopic or photometric redshift
      real*8        Reff       (nr__) ! (arcsec) Effective radius enclosing specified
*                                       fraction of flux
      real*8        logMstar   (nr__) ! ([Msun]) Stellar mass from SED fitting  (1)
      real*8        e_logMstar (nr__) ! ([Msun]) 1-sigma uncertainty from Bayesian
*                                       analysis on logMstar (1)
      real*8        logSFR10   (nr__) ! ([Msun/yr]) Star formation rate averaged over
*                                       last 10Myr (1)
      real*8        e_logSFR10 (nr__) ! ([Msun/yr]) 1-sigma uncertainty from Bayesian
*                                       analysis on logSFR10 (1)
      real*8        E_B_V      (nr__) ! (mag) Dust attenuation
      real*4        e_E_B_V    (nr__) ! (mag) 1-sigma uncertainty from Bayesian
*                                       analysis on E(B-V)
*Note (1): Stellar mass and SFR are derived using CIGALE with parametric SFHs,
*  Charlot & Bruzal (2019, see 2003MNRAS.344.1000B) stellar populations,
*  Chabrier (2003PASP..115..763C) IMF, and Dale et al. (2014ApJ...784...83D)
*  dust models.

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

C  Declarations for 'table2.dat'	! Emission line quantities

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

      integer*4     TileID_1   (nr__1) ! J-PAS tile identifier
      integer*4     Number_1   (nr__1) ! Source number within the tile
      real*8        z_1        (nr__1) ! Spectroscopic or photometric redshift
      real*8        xiHalpha   (nr__1) ! (10+7Hz/J) ? Ionizing photon production
*                                           efficiency from H-alpha (in Hz/erg
*                                           unit) (1)
      real*8        e_xiHalpha (nr__1) ! (10+7Hz/J) ? 1-sigma uncertainty on xi_Halpha
      real*8        xiHbeta    (nr__1) ! (10+7Hz/J) ? Ionizing photon production
*                                           efficiency from H-beta (in Hz/erg
*                                           unit) (1)
      real*8        e_xiHbeta  (nr__1) ! (10+7Hz/J) ? 1-sigma uncertainty on xi_Hbeta
      real*8        FHalpha    (nr__1) ! (10-20W/m2/nm) ? H-alpha line flux
*                                          (in 10^-19^erg/s/cm^2/{AA} unit)
      real*8        e_FHalpha  (nr__1) ! (10-20W/m2/nm) ? 1-sigma uncertainty on FHalpha
*                                          (in 10^-19^erg/s/cm^2/{AA} unit)
      real*8        FHbeta     (nr__1) ! (10-20W/m2/nm) ? H-beta line flux
*                                          (in 10^-19^erg/s/cm^2/{AA} unit)
      real*8        e_FHbeta   (nr__1) ! (10-20W/m2/nm) ? 1-sigma uncertainty on FHbetaa
*                                          (in 10^-19^erg/s/cm^2/{AA} unit)
      real*8        EWHalpha   (nr__1) ! (0.1nm) ? H-alpha equivalent width
*                                           (in {AA} unit)
      real*8        e_EWHalpha (nr__1) ! (0.1nm) ? 1-sigma uncertainty on EWHalpha
*                                           (in {AA} unit)
      real*8        EWHbeta    (nr__1) ! (0.1nm) ? H-beta equivalent width
*                                           (in {AA} unit)
      real*8        e_EWHbeta  (nr__1) ! (0.1nm) ? 1-sigma uncertainty on EWHbeta
*                                           (in {AA} unit)
      real*8        EWOIII5007 (nr__1) ! (0.1nm) ? Line [O III] 5007 A equivalent width
*                                           (in {AA} unit)
      real*8        e_EWOIII5007(nr__1) ! (0.1nm) ? 1-sigma uncertainty on EWOIII5007
*                                           (in {AA} unit)
*Note (1): xi values assume f_esc_ = 0, i.e., all Lyman continuum photons are
*  reprocessed into Balmer lines.

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

C  Loading file 'table1.dat'	! Derived SED Parameters

C  Format for file interpretation

    1 format(
     +  I4,1X,I5,1X,F18.14,1X,F18.15,1X,F11.9,1X,F10.8,2X,F19.16,1X,
     +  F20.18,1X,F21.18,1X,F19.17,1X,F19.17,1X,E22.19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,917
        read(1,'(A198)')ar__
        read(ar__,1)
     +  TileID(i__),Number(i__),RAdeg(i__),DEdeg(i__),z(i__),
     +  Reff(i__),logMstar(i__),e_logMstar(i__),logSFR10(i__),
     +  e_logSFR10(i__),E_B_V(i__),e_E_B_V(i__)
c    ..............Just test output...........
        write(6,1)
     +  TileID(i__),Number(i__),RAdeg(i__),DEdeg(i__),z(i__),
     +  Reff(i__),logMstar(i__),e_logMstar(i__),logSFR10(i__),
     +  e_logSFR10(i__),E_B_V(i__),e_E_B_V(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Emission line quantities

C  Format for file interpretation

    2 format(
     +  I4,1X,I5,1X,F11.9,1X,F18.15,1X,F20.18,1X,F18.15,1X,F20.18,1X,
     +  F21.16,1X,F19.16,1X,F24.18,1X,F20.16,1X,F21.16,1X,F20.16,1X,
     +  F23.17,1X,F21.16,1X,F23.16,1X,F19.15)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,917
        read(1,'(A323)')ar__1
        read(ar__1,2)
     +  TileID_1(i__),Number_1(i__),z_1(i__),xiHalpha(i__),
     +  e_xiHalpha(i__),xiHbeta(i__),e_xiHbeta(i__),FHalpha(i__),
     +  e_FHalpha(i__),FHbeta(i__),e_FHbeta(i__),EWHalpha(i__),
     +  e_EWHalpha(i__),EWHbeta(i__),e_EWHbeta(i__),EWOIII5007(i__),
     +  e_EWOIII5007(i__)
        if(ar__1(24:41) .EQ. '') xiHalpha(i__) = rNULL__
        if(ar__1(43:62) .EQ. '') e_xiHalpha(i__) = rNULL__
        if(ar__1(64:81) .EQ. '') xiHbeta(i__) = rNULL__
        if(ar__1(83:102) .EQ. '') e_xiHbeta(i__) = rNULL__
        if(ar__1(104:124) .EQ. '') FHalpha(i__) = rNULL__
        if(ar__1(126:144) .EQ. '') e_FHalpha(i__) = rNULL__
        if(ar__1(146:169) .EQ. '') FHbeta(i__) = rNULL__
        if(ar__1(171:190) .EQ. '') e_FHbeta(i__) = rNULL__
        if(ar__1(192:212) .EQ. '') EWHalpha(i__) = rNULL__
        if(ar__1(214:233) .EQ. '') e_EWHalpha(i__) = rNULL__
        if(ar__1(235:257) .EQ. '') EWHbeta(i__) = rNULL__
        if(ar__1(259:279) .EQ. '') e_EWHbeta(i__) = rNULL__
        if(ar__1(281:303) .EQ. '') EWOIII5007(i__) = rNULL__
        if(ar__1(305:323) .EQ. '') e_EWOIII5007(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  TileID_1(i__),Number_1(i__),z_1(i__),xiHalpha(i__),
     +  e_xiHalpha(i__),xiHbeta(i__),e_xiHbeta(i__),FHalpha(i__),
     +  e_FHalpha(i__),FHbeta(i__),e_FHbeta(i__),EWHalpha(i__),
     +  e_EWHalpha(i__),EWHbeta(i__),e_EWHbeta(i__),EWOIII5007(i__),
     +  e_EWOIII5007(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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