FORTRAN Generation
(/./ftp/cats/J/ApJ/705/1533)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/705/1533 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-29
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/705/1533        NGC 5128 globular clusters          (Chattopadhyay+, 2009)
*================================================================================
*Study of NGC 5128 globular clusters under multivariate statistical paradigm.
*    Chattopadhyay A.K., Chattopadhyay T., Davoust E., Mondal S., Sharina M.
*   <Astrophys. J., 705, 1533-1547 (2009)>
*   =2009ApJ...705.1533C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! List of parameters of the globular clusters of
                             NGC 5128 used in principal component analysis
                             (PCA) and cluster analysis (CA)

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

      character*9   Name       (nr__) ! Globular cluster name (G1)
      character*1   f_Name     (nr__) ! [a] discarded in CA (G2)
      integer*4     Group      (nr__) ! [1,3]? Group membership found by CA (G3)
      real*4        c          (nr__) ! Concentration (c~log(R_tid_/R_c_)) (G1)
      real*4        muV        (nr__) ! (mag/arcsec2) V surface brightness averaged over
*                                    half-light radius (<{mu}_v_>_h_) (G1)
      real*4        mu0        (nr__) ! (mag/arcsec2) Extinction-corrected central surface
*                                    brightness in the F606W bandpass (G1)
      real*4        W0         (nr__) ! Dimensionless central potential of the
*                                    best-fitting model (G1)
      real*4        logRc      (nr__) ! ([pc]) Log core radius (G1)
      real*4        logRtid    (nr__) ! ([pc]) Log tidal radius (G1)
      real*4        logrh      (nr__) ! ([pc]) Log half-light radius (G1)
      real*4        logM       (nr__) ! ([Msun]) Integrated model mass (logM_tot_) (G1)
      real*4        rho0       (nr__) ! (Msun/pc3) Central volume density (G1)
      real*4        sig_P      (nr__) ! (km/s) Predicted line of sight velocity dispersion
*                                    at the cluster center ({sigma}_p,0_) (G1)
      real*4        trh        (nr__) ! (10+9yr) Two-body relaxation time at the model
*                                    projected half-mass radius (G1)

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

C  Declarations for 'table2.dat'	! List of parameters as well as derived parameters
                             of the globular clusters of NGC 5128

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

      character*9   Name_1     (nr__1) ! Globular cluster name (G1)
      character*1   f_Name_1   (nr__1) ! [a] discarded in CA (G2)
      integer*4     Group_1    (nr__1) ! [1,3]? Group membership found by CA (G3)
      real*4        Rgc        (nr__1) ! (kpc) Galactocentric radius (G1)
      real*4        T1         (nr__1) ! (mag) Washington T1 magnitude (G1)
      real*4        v_Fe_H_    (nr__1) ! ([Sun]) Fe/H abundance (G1)
      real*4        v_C_T1_0   (nr__1) ! (mag) Extinction corrected (C-T1)_0_ color index (G1)
      integer*4     RV         (nr__1) ! (km/s) ? Weighted mean radial velocity from Woodley
*                                   et al. (2007, Cat. J/AJ/134/494)
      real*4        v_Z_H_     (nr__1) ! ([Sun]) ? Metallicity
      character*1   u__Z_H_    (nr__1) ! [:] Uncertainty flag on [Z/H] (1)
      real*4        e__Z_H_    (nr__1) ! ([Sun]) ? [Z/H] uncertainty
      real*4        Age        (nr__1) ! (Gyr) ? Age
      character*1   u_Age      (nr__1) ! [:] Uncertainty flag on Age (1)
      real*4        e_Age      (nr__1) ! (Gyr) ? Age uncertainty
      real*4        Vrot       (nr__1) ! (km/s) ? Rotational velocity
      real*4        sigma      (nr__1) ! (km/s) ? Velocity dispersion ({sigma}_V_)
*Note (1): ":" are used for the errors of [Z/H] and ages when the error of
*     [Z/H]>0.4dex and the error of Age>4Gyr.

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

C  Loading file 'table1.dat'	! List of parameters of the globular clusters of
*                             NGC 5128 used in principal component analysis
*                             (PCA) and cluster analysis (CA)

C  Format for file interpretation

    1 format(
     +  A9,A1,1X,I1,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F6.3,1X,F4.2,
     +  1X,F6.3,1X,F4.2,1X,E10.6,1X,F5.2,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,130
        read(1,'(A82)')ar__
        read(ar__,1)
     +  Name(i__),f_Name(i__),Group(i__),c(i__),muV(i__),mu0(i__),
     +  W0(i__),logRc(i__),logRtid(i__),logrh(i__),logM(i__),
     +  rho0(i__),sig_P(i__),trh(i__)
        if(ar__(12:12) .EQ. '') Group(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),f_Name(i__),Group(i__),c(i__),muV(i__),mu0(i__),
     +  W0(i__),logRc(i__),logRtid(i__),logrh(i__),logM(i__),
     +  rho0(i__),sig_P(i__),trh(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! List of parameters as well as derived parameters
*                             of the globular clusters of NGC 5128

C  Format for file interpretation

    2 format(
     +  A9,A1,1X,I1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,I3,1X,F5.2,A1,
     +  1X,F4.2,1X,F5.2,A1,1X,F4.2,1X,F5.2,1X,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,130
        read(1,'(A76)')ar__1
        read(ar__1,2)
     +  Name_1(i__),f_Name_1(i__),Group_1(i__),Rgc(i__),T1(i__),
     +  v_Fe_H_(i__),v_C_T1_0(i__),RV(i__),v_Z_H_(i__),u__Z_H_(i__),
     +  e__Z_H_(i__),Age(i__),u_Age(i__),e_Age(i__),Vrot(i__),
     +  sigma(i__)
        if(ar__1(12:12) .EQ. '') Group_1(i__) = iNULL__
        if(ar__1(37:39) .EQ. '') RV(i__) = iNULL__
        if(ar__1(41:45) .EQ. '') v_Z_H_(i__) = rNULL__
        if(ar__1(48:51) .EQ. '') e__Z_H_(i__) = rNULL__
        if(ar__1(53:57) .EQ. '') Age(i__) = rNULL__
        if(ar__1(60:63) .EQ. '') e_Age(i__) = rNULL__
        if(ar__1(65:69) .EQ. '') Vrot(i__) = rNULL__
        if(ar__1(71:76) .EQ. '') sigma(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1(i__),f_Name_1(i__),Group_1(i__),Rgc(i__),T1(i__),
     +  v_Fe_H_(i__),v_C_T1_0(i__),RV(i__),v_Z_H_(i__),u__Z_H_(i__),
     +  e__Z_H_(i__),Age(i__),u_Age(i__),e_Age(i__),Vrot(i__),
     +  sigma(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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