FORTRAN Generation
(/./ftp/cats/J/MNRAS/402/803)

Conversion of standardized ReadMe file for file /./ftp/cats/J/MNRAS/402/803 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-Apr-19
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/MNRAS/402/803      M31 globular cluster system             (Peacock+, 2010)
*================================================================================
*The M31 globular cluster system: ugriz and K-band photometry and structural
*parameters.
*    Peacock M.B., MacCarone T.J., Knigge C., Kundu A., Waters C.Z., Zepf S.E.,
*    Zurek D.R.
*   <Mon. Not. R. Astron. Soc., 402, 803-818 (2010)>
*   =2010MNRAS.402..803P
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Photometry and structural parameters of M31 GCs

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  Name       (nr__) ! Name taken from the RBC, Galleti et al.,
*                                       2004, Cat. J/A+A/416/917 (1)
      character*11  OName      (nr__) ! Other names from RBV (J/A+A/416/917)
      real*8        RAdeg      (nr__) ! (deg) Right ascension in decimal
*                                       degrees (J2000) (2)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000) (2)
      integer*4     f          (nr__) ! Classification of the globular cluster (3)
      character*7   n_f        (nr__) ! Additional classification for galaxies
      integer*4     f1         (nr__) ! ?=- Classification from the RBC V3.5 (3)
      character*9   f2         (nr__) ! Classification from Calwell et al.,
*                                       2009AJ....137...94C
      real*4        Rg         (nr__) ! (arcsec) ?=99.00 Aperture size used to obtain
*                                       total g magnitude
      real*4        gmag       (nr__) ! (mag) ?=99.00 g magnitude
      real*4        u_g        (nr__) ! (mag) ?=99.00 u-g colour index
      real*4        g_r        (nr__) ! (mag) ?=99.00 g-r colour index
      real*4        r_i        (nr__) ! (mag) ?=99.00 r-i colour index
      real*4        i_z        (nr__) ! (mag) ?=99.00 i-z colour index
      real*4        e_gmag     (nr__) ! (mag) ?=99.00 rms uncertainty on g total
*                                       magnitude, includes the statistical,
*                                       calibration and systematic errors
      real*4        e_u_g      (nr__) ! (mag) ?=99.00 rms uncertainty on u-g
      real*4        e_g_r      (nr__) ! (mag) ?=99.00 rms uncertainty on g-r
      real*4        e_r_i      (nr__) ! (mag) ?=99.00 rms uncertainty on r-i
      real*4        e_i_z      (nr__) ! (mag) ?=99.00 rms uncertainty on i-z
      real*4        RKs        (nr__) ! (arcsec) ?=99.00 Aperture size used to obtain
*                                       total Ks magnitude
      real*4        Ksmag      (nr__) ! (mag) ?=99.00 Ks magnitude
      real*4        e_Ksmag    (nr__) ! (mag) ?=99.00 rms uncertainty on Ksmag
      real*4        c          (nr__) ! ?=99.00 Concentration parameter
*                                       (C=log(rt/r_0_), r_0_ = core radius)
      real*4        rh         (nr__) ! (pc) ?=99.00 Half-light radius (4)
      real*4        rt         (nr__) ! (pc) ?=99.0 Tidal radius (4)
      real*4        log_rho0   (nr__) ! ([solLum/pc3]) ?=99.00 Cluster core density
*                                       (in L_K_({sun})/pc^3^ units)
*Note (1): Sources identified as M31GC JHHMMSS+DDMMSS in Simbad.
*Note (2): Position of object in SDSS r-band image
*Note (3): Classification as follows:
*    1 = old globular cluster
*    2 = candidate cluster (21-candidate old cluster, 23-candidate young cluster)
*    3 = young cluster
*    4 = galaxy
*    5 = HII region
*    6 = stellar source
*Note (4): assuming the distance of M31 to be 780kpc: McConnachie et al.
*     (2005MNRAS.356..979M)

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

C  Loading file 'table1.dat'	! Photometry and structural parameters of M31 GCs

C  Format for file interpretation

    1 format(
     +  A18,1X,A11,1X,F8.5,1X,F8.5,1X,I1,A7,1X,I1,1X,A9,1X,F4.1,1X,
     +  F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,2X,F5.2,1X,F5.2,1X,F5.2,
     +  1X,F5.2,1X,F5.2,1X,F4.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
     +  F4.1,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,1595
        read(1,'(A175)')ar__
        read(ar__,1)
     +  Name(i__),OName(i__),RAdeg(i__),DEdeg(i__),f(i__),n_f(i__),
     +  f1(i__),f2(i__),Rg(i__),gmag(i__),u_g(i__),g_r(i__),r_i(i__),
     +  i_z(i__),e_gmag(i__),e_u_g(i__),e_g_r(i__),e_r_i(i__),
     +  e_i_z(i__),RKs(i__),Ksmag(i__),e_Ksmag(i__),c(i__),rh(i__),
     +  rt(i__),log_rho0(i__)
        if (f1(i__) .EQ. 45) f1(i__) =  iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),OName(i__),RAdeg(i__),DEdeg(i__),f(i__),n_f(i__),
     +  f1(i__),f2(i__),Rg(i__),gmag(i__),u_g(i__),g_r(i__),r_i(i__),
     +  i_z(i__),e_gmag(i__),e_u_g(i__),e_g_r(i__),e_r_i(i__),
     +  e_i_z(i__),RKs(i__),Ksmag(i__),e_Ksmag(i__),c(i__),rh(i__),
     +  rt(i__),log_rho0(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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