FORTRAN Generation
(/./ftp/cats/J/A_A/589/A139)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/589/A139 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-Jun-08
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/589/A139      Stellar populations of NGC 3311            (Barbosa+, 2016)
*================================================================================
*The Hydra I cluster core. I. Stellar populations in the cD galaxy NGC 3311.
*    Barbosa C.E., Arnaboldi M., Coccato L., Hilker M., Mendes de Oliveira C.,
*    Richtler T.
*    <Astron. Astrophys. 589, A139 (2016)>
*    =2016A&A...589A.139B        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Equivalent width of Lick indices for the Hydra I
                                 cluster core and stellar population parameters

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

      character*8   ID         (nr__) ! Identification of the slit (1)
      real*4        Rad        (nr__) ! (kpc) Projected distance to the center of
*                                     NGC 3311 (2)
      real*4        PA         (nr__) ! (deg) [] Position angle in relation to NGC 3311
      real*4        S_N        (nr__) ! (10/nm) Signal-to-Noise ratio
      real*4        Hbeta      (nr__) ! (0.1nm) ?=- H{beta} equivalent width
      real*4        e_Hbeta    (nr__) ! (0.1nm) ?=- Error in the Hbeta equivalent width
      real*4        Fe5015     (nr__) ! (0.1nm) ?=- Fe 5015{AA} equivalent width
      real*4        e_Fe5015   (nr__) ! (0.1nm) ?=- Error in the Fe5015 equivalent width
      real*4        Mgb        (nr__) ! (0.1nm) ?=- Mgb equivalent width
      real*4        e_Mgb      (nr__) ! (0.1nm) ?=- Error in the Mgb equivalent width
      real*4        Fe5270     (nr__) ! (0.1nm) ?=- Fe 5270{AA} equivalent width
      real*4        e_Fe5270   (nr__) ! (0.1nm) ?=- Error in the Fe5270 equivalent width
      real*4        Fe5335     (nr__) ! (0.1nm) ?=- Fe 5335{AA} equivalent width
      real*4        e_Fe5335   (nr__) ! (0.1nm) ?=- Error in the Fe5335 equivalent width
      real*4        Fe5406     (nr__) ! (0.1nm) ?=- Fe 5406{AA} equivalent width
      real*4        e_Fe5406   (nr__) ! (0.1nm) ?=- Error in the Fe5406 equivalent width
      real*4        Fe5709     (nr__) ! (0.1nm) ?=- Fe 5709{AA} equivalent width
      real*4        e_Fe5709   (nr__) ! (0.1nm) ?=- Error in the Fe5709 equivalent width
      real*4        Age        (nr__) ! ([yr]) Luminosity weighted age
      real*4        e_Age      (nr__) ! ([yr]) Lower error in age
      real*4        E_Age_1    (nr__) ! ([yr]) Upper error in age
      real*4        v_Z_H_     (nr__) ! ([Sun]) Luminosity weighted total metallicity
      real*4        e__Z_H_    (nr__) ! ([Sun]) Lower error in total metallicity
      real*4        E__Z_H__1  (nr__) ! ([Sun]) Upper error in total metallicity
      real*4        v_alpha_Fe_(nr__) ! ([Sun]) Luminosity weighted alpha element abundance
      real*4        e__alpha_Fe_(nr__) ! ([Sun]) Lower error in alpha element abundance
      real*4        E__alpha_Fe__1(nr__) ! ([Sun]) Upper error in alpha element abundance
*Note (1): Slits names as cen1sNN, cen2sNN,inn1sNN, inn2sNN, out1sNN and out2sNN.
*Note (2): Radial distances to the center of NGC 3311 are calculated assuming
*          a distance of 50.7Mpc to the core of the Hydra I cluster.

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

C  Loading file 'table1.dat'	! Equivalent width of Lick indices for the Hydra I
*                                 cluster core and stellar population parameters

C  Format for file interpretation

    1 format(
     +  A8,1X,F4.1,1X,F6.1,1X,F5.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,119
        read(1,'(A152)')ar__
        read(ar__,1)
     +  ID(i__),Rad(i__),PA(i__),S_N(i__),Hbeta(i__),e_Hbeta(i__),
     +  Fe5015(i__),e_Fe5015(i__),Mgb(i__),e_Mgb(i__),Fe5270(i__),
     +  e_Fe5270(i__),Fe5335(i__),e_Fe5335(i__),Fe5406(i__),
     +  e_Fe5406(i__),Fe5709(i__),e_Fe5709(i__),Age(i__),e_Age(i__),
     +  E_Age_1(i__),v_Z_H_(i__),e__Z_H_(i__),E__Z_H__1(i__),
     +  v_alpha_Fe_(i__),e__alpha_Fe_(i__),E__alpha_Fe__1(i__)
        if (idig(ar__(28:32)).EQ.0) Hbeta(i__) =  rNULL__
        if (idig(ar__(34:38)).EQ.0) e_Hbeta(i__) =  rNULL__
        if (idig(ar__(40:44)).EQ.0) Fe5015(i__) =  rNULL__
        if (idig(ar__(46:49)).EQ.0) e_Fe5015(i__) =  rNULL__
        if (idig(ar__(51:55)).EQ.0) Mgb(i__) =  rNULL__
        if (idig(ar__(57:60)).EQ.0) e_Mgb(i__) =  rNULL__
        if (idig(ar__(62:66)).EQ.0) Fe5270(i__) =  rNULL__
        if (idig(ar__(68:71)).EQ.0) e_Fe5270(i__) =  rNULL__
        if (idig(ar__(73:77)).EQ.0) Fe5335(i__) =  rNULL__
        if (idig(ar__(79:82)).EQ.0) e_Fe5335(i__) =  rNULL__
        if (idig(ar__(84:88)).EQ.0) Fe5406(i__) =  rNULL__
        if (idig(ar__(90:93)).EQ.0) e_Fe5406(i__) =  rNULL__
        if (idig(ar__(95:99)).EQ.0) Fe5709(i__) =  rNULL__
        if (idig(ar__(101:104)).EQ.0) e_Fe5709(i__) =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),Rad(i__),PA(i__),S_N(i__),Hbeta(i__),e_Hbeta(i__),
     +  Fe5015(i__),e_Fe5015(i__),Mgb(i__),e_Mgb(i__),Fe5270(i__),
     +  e_Fe5270(i__),Fe5335(i__),e_Fe5335(i__),Fe5406(i__),
     +  e_Fe5406(i__),Fe5709(i__),e_Fe5709(i__),Age(i__),e_Age(i__),
     +  E_Age_1(i__),v_Z_H_(i__),e__Z_H_(i__),E__Z_H__1(i__),
     +  v_alpha_Fe_(i__),e__alpha_Fe_(i__),E__alpha_Fe__1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end