FORTRAN Generation
(/./ftp/cats/J/A_A/534/A102)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/534/A102 into FORTRAN code for reading data files line by line.

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/A+A/534/A102      AMIGA IX. Molecular gas properties       (Lisenfeld+, 2011)
*================================================================================
*The AMIGA sample of isolated galaxies.
*IX. Molecular gas properties.
*    Lisenfeld U., Espada D., Verdes-Montenegro L., Kuno N., Leon S., Sabater J.,
*    Sato N., Sulentic J., Verley S., Min Y.
*   <Astron. Astrophys. 534, A102 (2011)>
*   =2011A&A...534A.102L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! General data for the total CO sample

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

      integer*4     CIG         ! CIG-catalog galaxy number
      integer*4     Dist        ! (Mpc) Distance (based on H_0_=75km/s/Mpc)
      integer*4     Vel         ! (km/s) Recession velocity
      real*4        D25         ! (arcmin) Isophotal major diameter at the isophotal
*                                 level 25mag/arcsec^2^
      integer*4     i           ! (deg) Inclination from LEDA
      integer*4     TT          ! Morphological type, given in the RC3
*                                 numerical scale
      integer*4     Mi          ! Code for morphological signature of a possible
*                                 interaction (1)
      real*4        log_LB      ! ([Lsun]) Decimal logarithm of the blue luminosity in
*                                 units of bolometric solar luminosity
      character*1   l_log_LFIR  ! Limit flag on log(LFIR) (2)
      real*4        log_LFIR    ! ([Lsun]) ?=- Decimal logarithm of far-infrared luminosity
*                                 in units of bolometric solar luminosity
      real*4        log_LK      ! ([Lsun]) ?=- Decimal logarithm of K-band luminosity
*                                 from 2MASS in units of K-band solar luminosity
*Note (1): Code for interaction as follows:
*    0 = No signs of large distortions which cannot be clearly explained
*        by e.g. dust
*    1 = At least one of the following signs of possible interaction:
*        asymmetric, lopsided, distorted, confirmed pair in the Nasa
*        Extragalactic Database (NED) or in Verley et al. (2007,
*        Cat. J/A+A/472/121), integral sign shape, warp, tidal features
*        (tail, bridge, shell)
*    2 = Merger-like morphology or superposition of two galaxies
*Note (2): Upper limit code for FIR luminosity, LFIR
*     if <, LFIR is an upper limit, else LFIR is a detection

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

C  Declarations for 'table4.dat'	! Velocity integrated ^12^CO(1-0) line intensities,
                              mean velocities and line widths

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

      integer*4     CIG_1       ! CIG-catalog galaxy number
      character*1   u_CIG       ! [*] Code for data quality:
*                                 * means a marginal detections
      integer*4     oRA         ! (arcsec) RA offset from the center
      integer*4     oDE         ! (arcsec) Declination offset from the center
      real*4        rms         ! (mK) Root mean square noise for a velocity
*                                 resolution of 10.6km/s (IRAM), 13.1km/s (FCRAO)
      character*1   l_ICO       ! Limit flag on ICO (1)
      real*4        ICO         ! (K.km/s) Velocity integrated CO line temperature (Tmb)
      real*4        e_ICO       ! (K.km/s) ?=- rms uncertainty on ICO
      integer*4     VCO         ! (km/s) ?=- Mean velocity of the CO line
      integer*4     WCO         ! (km/s) ?=- Zero line-width of the CO spectrum
      integer*4     Tel         ! [1/2] Radio telescope code: 1=IRAM, 2=FCRAO
*Note (1): Upper limit code for ICO:
*     if <, ICO is an upper limit, else ICO is a detection

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

C  Declarations for 'table5.dat'	! Molecular gas mass

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

      integer*4     CIG_2       ! CIG-catalog galaxy number
      integer*4     Det         ! [0/2] Detection code (1)
      real*4        MH2c        ! ([Msun]) ?=0 Decimal logarithm of molecular gas mas
*                                  towards the central pointing (2)
      real*4        MH2m        ! ([Msun]) ?=0 Decimal logarithm of mapped molecular gas
*                                 mass (2)
      real*4        MH2e        ! ([Msun]) ?=0 Decimal logarithm of total molecular gas
*                                 mass, extrapolated from the central value (2)
      integer*4     Tel_1       ! [1/5] Radio-telescope code (3)
      integer*4     Ref         ! [1/6] Bibliographic code (4)
*Note (1): Detection code for the molecular mass as follows:
*   0 = Detection
*   1 = Upper limit
*   2 = Marginal detection
*Note (2): Zero indicates a blank value (no measurement taken).
*Note (3): Radio telescope code as follows:
*   1 = IRAM
*   2 = FCRAO
*   3 = SEST
*   4 = Onsala
*   5 = NRAO Kitt-Peak
*Note (4): Bibliographic code as follows:
*   1 = This work
*   2 = Young et al. (1995ApJS...98...219J)
*   3 = Braine & Combes (1993A&A...269....7B)
*   4 = Sauty et al. (2003A&A...411..318S)
*   5 = Elfhag et al. (1996, Cat. J/A+AS/115/439)
*   6 = Sage (1993A&A...272..123)

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

C  Loading file 'table1.dat'	! General data for the total CO sample

C  Format for file interpretation

    1 format(
     +  I4,1X,I3,1X,I5,1X,F5.2,1X,I2,1X,I2,1X,I1,1X,F5.2,1X,A1,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,273
        read(1,'(A47)')ar__
        read(ar__,1)
     +  CIG,Dist,Vel,D25,i,TT,Mi,log_LB,l_log_LFIR,log_LFIR,log_LK
        if (idig(ar__(37:41)).EQ.0) log_LFIR =  rNULL__
        if (idig(ar__(43:47)).EQ.0) log_LK =  rNULL__
c    ..............Just test output...........
        write(6,1)
     +  CIG,Dist,Vel,D25,i,TT,Mi,log_LB,l_log_LFIR,log_LFIR,log_LK
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Velocity integrated ^12^CO(1-0) line intensities,
*                              mean velocities and line widths

C  Format for file interpretation

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

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,285
        read(1,'(A43)')ar__1
        read(ar__1,2)
     +  CIG_1,u_CIG,oRA,oDE,rms,l_ICO,ICO,e_ICO,VCO,WCO,Tel
        if (idig(ar__1(29:32)).EQ.0) e_ICO =  rNULL__
        if (VCO .EQ. 45) VCO =  iNULL__
        if (WCO .EQ. 45) WCO =  iNULL__
c    ..............Just test output...........
        write(6,2)
     +  CIG_1,u_CIG,oRA,oDE,rms,l_ICO,ICO,e_ICO,VCO,WCO,Tel
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Molecular gas mass

C  Format for file interpretation

    3 format(I4,1X,I1,1X,F4.2,1X,F4.2,1X,F5.2,1X,I1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,273
        read(1,'(A26)')ar__2
        read(ar__2,3)CIG_2,Det,MH2c,MH2m,MH2e,Tel_1,Ref
c    ..............Just test output...........
        write(6,3)CIG_2,Det,MH2c,MH2m,MH2e,Tel_1,Ref
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