FORTRAN Generation
(/./ftp/cats/J/A_A/405/999)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/405/999 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-28
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/405/999      Unidentified infrared bands (UIR)    (Kahanpaeae+, 2003)
*================================================================================
*Unidentified infrared bands in the interstellar medium across the Galaxy.
*    Kahanpaeae J., Mattila K., Lehtinen K., Leinert C., Lemke D.
*   <Astron. Astrophys. 405, 999 (2003)>
*   =2003A&A...405..999K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Summary of observed positions (tables 1+2)

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*7   Name        ! Target name
      character*1   n_Name      ! [a-d] Individual Note (1)
      integer*4     RAh         ! (h) Right ascension (J2000)
      integer*4     RAm         ! (min) Right ascension (J2000)
      real*4        RAs         ! (s) Right ascension (J2000)
      character*1   DE_         ! Declination sign (J2000)
      integer*4     DEd         ! (deg) Declination (J2000)
      integer*4     DEm         ! (arcmin) Declination (J2000)
      real*4        DEs         ! (arcsec) Declination (J2000)
      real*4        GLON        ! (deg) Galactic longitude
      real*4        GLAT        ! (deg) Galactic latitude
      real*4        I6_2um      ! (10-10W/cm2/um/sr) ? UIR line intensity at 6.2um band
      real*4        e_I6_2um    ! (10-10W/cm2/um/sr) ? rms uncertainty on I6.2um
      real*4        I7_7um      ! (10-10W/cm2/um/sr) ? UIR line intensity at 7.7um band
      real*4        e_I7_7um    ! (10-10W/cm2/um/sr) ? rms uncertainty on I7.7um
      real*4        I8_6um      ! (10-10W/cm2/um/sr) ? UIR line intensity at 8.5um band
      real*4        e_I8_6um    ! (10-10W/cm2/um/sr) ? rms uncertainty on I8.6um
      real*4        I11_3um     ! (10-10W/cm2/um/sr) ? UIR line intensity at 11.3um band
      real*4        e_I11_3um   ! (10-10W/cm2/um/sr) ? rms uncertainty on I11.3um
      integer*4     TDT1        ! Target Dedicated Time (TDT) number
      integer*4     TDT2        ! ? ISO second TDT number
*Note (1): Notes:
*     a: Very high dark current.
*     b: Has 7.7{mu}m "ghost" band in the dark measurement.
*     c: Shows evidence of point source contamination in one of the four
*        pixels. This pixel was omitted from the average value calculation.
*     d: Observed 7 months later than other G+30 points.

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

C  Loading file 'table1.dat'	! Summary of observed positions (tables 1+2)

C  Format for file interpretation

    1 format(
     +  A7,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,F4.1,1X,F6.2,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,
     +  1X,F4.2,1X,I8,1X,I8)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,49
        read(1,'(A103)')ar__
        read(ar__,1)
     +  Name,n_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,I6_2um,
     +  e_I6_2um,I7_7um,e_I7_7um,I8_6um,e_I8_6um,I11_3um,e_I11_3um,
     +  TDT1,TDT2
        if(ar__(46:49) .EQ. '') I6_2um = rNULL__
        if(ar__(51:54) .EQ. '') e_I6_2um = rNULL__
        if(ar__(56:59) .EQ. '') I7_7um = rNULL__
        if(ar__(61:64) .EQ. '') e_I7_7um = rNULL__
        if(ar__(66:69) .EQ. '') I8_6um = rNULL__
        if(ar__(71:74) .EQ. '') e_I8_6um = rNULL__
        if(ar__(76:80) .EQ. '') I11_3um = rNULL__
        if(ar__(82:85) .EQ. '') e_I11_3um = rNULL__
        if(ar__(96:103) .EQ. '') TDT2 = iNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Name,n_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,I6_2um,
     +  e_I6_2um,I7_7um,e_I7_7um,I8_6um,e_I8_6um,I11_3um,e_I11_3um,
     +  TDT1,TDT2
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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