FORTRAN Generation
(/./ftp/cats/J/A_A/667/A76)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/667/A76 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-Sep-09
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/667/A76       64 Virgo galaxies radial surface brightness  (Junais+, 2022)
*================================================================================
*A Virgo Environmental Survey Tracing Ionised Gas Emission (VESTIGE). XIII.
*The role of ram-pressure stripping in transforming the diffuse and ultra-diffuse
*galaxies in the Virgo cluster .
*    Junais S., Boissier S., Boselli A., Ferrarese L., Cote P., Gwyn S.,
*    Roediger J., Lim S., Peng E.W., Cuillandre J.-C., Longobardi A., Fossati M.,
*    Hensler G., Koda J., Bautista J., Boquien M., Malek K., Amram P., Roehlly Y.
*    <Astron. Astrophys. 667, A76 (2022)>
*    =2022A&A...667A..76J        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'tabled1.dat'	! *Properties of the selected sample of galaxies

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     ID         (nr__) ! Source ID based on the position in the 
*                                       NGVS catalog
      character*27  NGVS       (nr__) ! Name  of the source in NGVS 
*                                       (NGVSJHH:MM:SS.ss+DD:MM:SS.s)
      character*8   VCC        (nr__) ! VCC name
      real*8        RAdeg      (nr__) ! (deg) Right ascension (J2000)
      real*4        DEdeg      (nr__) ! (deg) Declination (J2000)
      real*4        DM87       (nr__) ! (Mpc) Projected distance of the galaxy from M87
      real*4        gmag       (nr__) ! (mag) g-band magnitude
      real*4        M_         (nr__) ! (10+7Msun) Stellar mass of the galaxies 
*                                       (see Sect. 2.1)
      real*4        Reg        (nr__) ! (kpc) Effective radius of the source
      real*4        mu0g       (nr__) ! (mag/arcsec+2) g-band central surface brightness
      real*4        q          (nr__) ! Axis ratio
      real*4        PA         (nr__) ! (deg) [] Position angle
      real*4        E_B_V      (nr__) ! (mag) Foreground galactic extinction from 
*                                       Schlegel et al. (1998ApJ...500..525S)
      integer*4     Flag       (nr__) ! [1/2] Flag for the galaxy type (1)
*Note (1): Flag for the galaxy type as follows:
*    1 = UDGs from Lim et al. (2020ApJ...899...69L) 
*    2 = Diffuse galaxies

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

C  Declarations for 'appen.dat'	! Radial surface brightness profiles

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     ID_1       (nr__1) ! ID of the galaxy
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension of center (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination of center (J2000)
      real*4        Rarcsec    (nr__1) ! (arcsec) Radius in angular units
      real*4        Rkpc       (nr__1) ! (kpc) Radius in physical units
      character*1   l_SuBru    (nr__1) ! Limit flag on SuBru
      real*4        SuBru      (nr__1) ! (mag/arcsec2) Observed u band surface brightness (1)
      real*4        E_SuBru    (nr__1) ! (mag/arcsec2) Upper limit uncertainty in SuBru
      real*4        e_SuBru_1  (nr__1) ! (mag/arcsec2) Lower limit uncertainty in SuBru
      character*1   l_SuBrg    (nr__1) ! Limit flag on SuBrg (3)
      real*4        SuBrg      (nr__1) ! (mag/arcsec2) Observed g band surface brightness (1)
      real*4        E_SuBrg    (nr__1) ! (mag/arcsec2) Upper limit uncertainty in SuBrg
      real*4        e_SuBrg_1  (nr__1) ! (mag/arcsec2) Lower limit uncertainty in SuBrg
      character*1   l_SuBrr    (nr__1) ! Limit flag on SuBrr (3)
      real*4        SuBrr      (nr__1) ! (mag/arcsec2) Observed r band surface brightness (1)
      real*4        E_SuBrr    (nr__1) ! (mag/arcsec2) Upper limit uncertainty in SuBrr
      real*4        e_SuBrr_1  (nr__1) ! (mag/arcsec2) Lower limit uncertainty in SuBrr
      character*1   l_SuBri    (nr__1) ! Limit flag on SuBri (3)
      real*4        SuBri      (nr__1) ! (mag/arcsec2) Observed i band surface brightness (1)
      real*4        E_SuBri    (nr__1) ! (mag/arcsec2) Upper limit uncertainty in SuBri
      real*4        e_SuBri_1  (nr__1) ! (mag/arcsec2) Lower limit uncertainty in SuBri
      character*1   l_SuBrz    (nr__1) ! Limit flag on SuBrz (3)
      real*4        SuBrz      (nr__1) ! (mag/arcsec2) Observed z band surface brightness (1)
      real*4        E_SuBrz    (nr__1) ! (mag/arcsec2) Upper limit uncertainty in SuBrz
      real*4        e_SuBrz_1  (nr__1) ! (mag/arcsec2) Lower limit uncertainty in SuBrz
      character*1   l_SigmaHa  (nr__1) ! Limit flag on SigmaHa (3)
      real*4        SigmaHa    (nr__1) ! (mW/m2/arcsec2) Observed H{alpha} narrow-band surface
*                                          brightness (erg/cm2/s/arcsec2)
      real*4        E_SigmaHa  (nr__1) ! (mW/m2/arcsec2) Upper limit uncertainty in SigmaHa
*                                          (erg/cm2/s/arcsec2)
      real*4        e_SigmaHa_1(nr__1) ! (mW/m2/arcsec2) Lower limit uncertainty in SigmaHa
*                                          (erg/cm2/s/arcsec2)
      character*1   l_SuBrFUV  (nr__1) ! Limit flag on SuBrFUV (3)
      real*8        SuBrFUV    (nr__1) ! (mag/arcsec2) ?=-100 Observed FUV band surface
*                                          brightness (1) (4)
      real*8        E_SuBrFUV  (nr__1) ! (mag/arcsec2) ?=-100  Upper limit uncertainty in
*                                          SuBrFUV (4)
      real*8        e_SuBrFUV_1(nr__1) ! (mag/arcsec2) ?=-100 Lower limit uncertainty in
*                                          SuBrFUV (4)
      character*1   l_SuBrNUV  (nr__1) ! Limit flag on SuBrnuv
      real*8        SuBrNUV    (nr__1) ! (mag/arcsec2) ?=-100 Observed NUV band surface
*                                          brightness (1) (4)
      real*8        E_SuBrNUV  (nr__1) ! (mag/arcsec2) ?=-100  Upper limit uncertainty in
*                                          SuBrnuv (4)
      real*8        e_SuBrNUV_1(nr__1) ! (mag/arcsec2) ?=-100 Lower limit uncertainty in
*                                          SuBrnuv (4)
*Note (1): Corrected for Galactic extinction and inclination.
*          All the magnitudes are in AB units.
*Note (3): Three sigma upper limits in surface brightness are marked with > or <.
*Note (4): Unavailable data-points are marked as -100.
*          FUV and NUV data is not available for galaxy ID 892 and ID 935.
*           For ID 1164, FUV data is not available.

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

C  Declarations for 'list.dat'	! List of fits images

      integer*4 nr__2
      parameter (nr__2=507)	! Number of records
      character*96 ar__2  	! Full-size record

C  J2000 position composed of: RAdeg DEdeg
      integer*4     ID_2       (nr__2) ! ID of the galaxy
      real*8        RAdeg_2    (nr__2) ! (deg) Right Ascension of center (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Declination of center (J2000)
      real*4        scale      (nr__2) ! (arcsec/pix) Scale of the image
      integer*4     Nx         (nr__2) ! Number of pixels along X-axis
      integer*4     Ny         (nr__2) ! Number of pixels along Y-axis
      integer*4     size       (nr__2) ! (Kibyte) Size of FITS file
      character*19  FileName   (nr__2) ! Name of FITS file, in subdirectory fits
      character*29  Title      (nr__2) ! Title of the FITS file

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

C  Loading file 'tabled1.dat'	! *Properties of the selected sample of galaxies

C  Format for file interpretation

    1 format(
     +  I4,1X,A27,1X,A8,1X,F7.3,1X,F6.3,1X,F4.2,1X,F5.2,1X,F5.2,1X,
     +  F4.2,1X,F5.2,1X,F4.2,1X,F6.2,1X,F5.3,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tabled1.dat')
      write(6,*) '....Loading file: tabled1.dat'
      do i__=1,64
        read(1,'(A104)')ar__
        read(ar__,1)
     +  ID(i__),NGVS(i__),VCC(i__),RAdeg(i__),DEdeg(i__),DM87(i__),
     +  gmag(i__),M_(i__),Reg(i__),mu0g(i__),q(i__),PA(i__),
     +  E_B_V(i__),Flag(i__)
c    ..............Just test output...........
        write(6,1)
     +  ID(i__),NGVS(i__),VCC(i__),RAdeg(i__),DEdeg(i__),DM87(i__),
     +  gmag(i__),M_(i__),Reg(i__),mu0g(i__),q(i__),PA(i__),
     +  E_B_V(i__),Flag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'appen.dat'	! Radial surface brightness profiles

C  Format for file interpretation

    2 format(
     +  I4,1X,F9.5,1X,F8.5,1X,F6.2,1X,F5.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,
     +  1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,A1,
     +  F5.2,1X,F4.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,A1,E8.4,2X,
     +  E8.4,2X,E8.4,1X,A1,F7.2,1X,F7.2,1X,F7.2,1X,A1,F7.2,1X,F7.2,1X,
     +  F7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'appen.dat')
      write(6,*) '....Loading file: appen.dat'
      do i__=1,1270
        read(1,'(A201)')ar__1
        read(ar__1,2)
     +  ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),Rarcsec(i__),Rkpc(i__),
     +  l_SuBru(i__),SuBru(i__),E_SuBru(i__),e_SuBru_1(i__),
     +  l_SuBrg(i__),SuBrg(i__),E_SuBrg(i__),e_SuBrg_1(i__),
     +  l_SuBrr(i__),SuBrr(i__),E_SuBrr(i__),e_SuBrr_1(i__),
     +  l_SuBri(i__),SuBri(i__),E_SuBri(i__),e_SuBri_1(i__),
     +  l_SuBrz(i__),SuBrz(i__),E_SuBrz(i__),e_SuBrz_1(i__),
     +  l_SigmaHa(i__),SigmaHa(i__),E_SigmaHa(i__),e_SigmaHa_1(i__),
     +  l_SuBrFUV(i__),SuBrFUV(i__),E_SuBrFUV(i__),e_SuBrFUV_1(i__),
     +  l_SuBrNUV(i__),SuBrNUV(i__),E_SuBrNUV(i__),e_SuBrNUV_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),Rarcsec(i__),Rkpc(i__),
     +  l_SuBru(i__),SuBru(i__),E_SuBru(i__),e_SuBru_1(i__),
     +  l_SuBrg(i__),SuBrg(i__),E_SuBrg(i__),e_SuBrg_1(i__),
     +  l_SuBrr(i__),SuBrr(i__),E_SuBrr(i__),e_SuBrr_1(i__),
     +  l_SuBri(i__),SuBri(i__),E_SuBri(i__),e_SuBri_1(i__),
     +  l_SuBrz(i__),SuBrz(i__),E_SuBrz(i__),e_SuBrz_1(i__),
     +  l_SigmaHa(i__),SigmaHa(i__),E_SigmaHa(i__),e_SigmaHa_1(i__),
     +  l_SuBrFUV(i__),SuBrFUV(i__),E_SuBrFUV(i__),e_SuBrFUV_1(i__),
     +  l_SuBrNUV(i__),SuBrNUV(i__),E_SuBrNUV(i__),e_SuBrNUV_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'list.dat'	! List of fits images

C  Format for file interpretation

    3 format(I4,1X,F9.5,F9.5,1X,F5.3,1X,I4,1X,I4,1X,I6,1X,A19,1X,A29)

C  Effective file loading

      open(unit=1,status='old',file=
     +'list.dat')
      write(6,*) '....Loading file: list.dat'
      do i__=1,507
        read(1,'(A96)')ar__2
        read(ar__2,3)
     +  ID_2(i__),RAdeg_2(i__),DEdeg_2(i__),scale(i__),Nx(i__),
     +  Ny(i__),size(i__),FileName(i__),Title(i__)
c    ..............Just test output...........
        write(6,3)
     +  ID_2(i__),RAdeg_2(i__),DEdeg_2(i__),scale(i__),Nx(i__),
     +  Ny(i__),size(i__),FileName(i__),Title(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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