CDS.headStuff2 FORTRAN Generation
()

Conversion of standardized ReadMe file for catalog 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-25
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. II/297           AKARI/IRC mid-IR all-sky Survey    (ISAS/JAXA, 2010)
*================================================================================
*The AKARI/IRC Mid-Infrared All-Sky Survey (Version 1)
*    Ishihara D., Onaka T., Kataza H., Salama A., Alfageme C., Cassatella A.,
*    Cox N., Garcia-Lario P., Stephenson C., Cohen M., Fujishiro N., 
*    Fujiwara H., Hasegawa S., Ita Y., Kim W., Matsuhara H., Murakami H.,
*    Muller T.G., Nakagawa T., Ohyama Y., Oyabu S., Pyo J., Sakon I.,
*    Shibai H., Takita S., Tanab T., Uemizu K., Ueno M., Usui F., Wada T.,
*    Watarai H., Yamamura I., Yamauchi C. 
*   <Astron. Astrophys. 514, A1 (2010)>
*   =2010A&A...514A...1I
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'irc.dat'	! AKARI/IRC All-Sky Survey Point Source Catalogue
                                  (Version 1.0)

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     objID      (nr__) ! [200000001/200870973] AKARI source ID number.
      character*14  objName    (nr__) ! AKARI source name (HHMMSSs+DDMMSS) (2)
      real*8        RAdeg      (nr__) ! (deg) Right Ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination (J2000)
      real*4        errMaj     (nr__) ! (arcsec) Major axis of position error ellipse
      real*4        errMin     (nr__) ! (arcsec) Minor axis of position error ellipse
      real*4        errPA      (nr__) ! (deg) Position angle of Major axis
      real*4        S09        (nr__) ! (Jy) ?=-999.9 Flux density in AKARI/S9W filter (1)
      real*4        S18        (nr__) ! (Jy) ?=-999.9 Flux density in AKARI/L18W filter (1)
      real*4        e_S09      (nr__) ! (Jy) ?=-1000 Flux error in S9W
      real*4        e_S18      (nr__) ! (Jy) ?=-1000 Flux error in L18W
      integer*4     q_S09      (nr__) ! [0,3] Flux quality flag for S9W (3)
      integer*4     q_S18      (nr__) ! [0,3] Flux quality flag for L18W (3)
      character*1   f09        (nr__) ! [0-9a-f] Bit flags for S9W (4)
      character*1   f18        (nr__) ! [0-9a-f] Bit flags for L18W (4)
      integer*4     Ns09       (nr__) ! Number of scans with source detection in S9W
      integer*4     Ns18       (nr__) ! Number of scans with source detection in L18W
      integer*4     Np09       (nr__) ! Number of scans with possible detection in S09W
      integer*4     Np18       (nr__) ! Number of scans with possible detection in L18W
      integer*4     M09        (nr__) ! [0,1]?=-1 1 is month confirmed and 0 is not.
*                                    (inverted value of lower bit of f09)
      integer*4     M18        (nr__) ! [0,1]?=-1 1 is month confirmed and 0 is not.
*                                    (inverted value of lower bit of f18)
      integer*4     Nd09       (nr__) ! ?=-1 Number of sources in 45" radius in S9W
      integer*4     Nd18       (nr__) ! ?=-1 Number of sources in 45" radius in L18W
      integer*4     X09        (nr__) ! [0,1]?=-1 Extended source flag (5)
      integer*4     X18        (nr__) ! [0,1]?=-1 Extended source flag (5)
      real*8        r09        (nr__) ! (arcsec) ?=-999.9 Radius of source extent in S9W (6)
      real*8        r18        (nr__) ! (arcsec) ?=-999.9 Radius of source extent in L18W (6)
      integer*4     Ndet       (nr__) ! Number of events used for position computation
      integer*4     N09        (nr__) ! Number of events used for S09 measurement
      integer*4     N18        (nr__) ! Number of events used for S18 measurement
*Note (1): the filter characteristics are:
*     -------------------------------------------------
*           Filter:    S9W      L18W
*     -------------------------------------------------
*       Center(um):     9         18
*        Width(um):    4.10       9.97
*        Range(um):  6.7-11.6  13.9-25.6
*         Pixel("):  9.4x9.4   10.4x9.4
*       Limit(mJy):     50       120     [at 5{sigma}]
*     -------------------------------------------------
*Note (2):
*    The sources should be referred in the literatures by their full name
*    (AKARI-IRC-V1) followed by the letter 'J' and the objName, e.g.
*    AKARI-IRC-V1 J0123498-025805
*Note (3): Four-level flux quality indicator:
*    3 = high quality (source confirmed and flux is reliable)
*    2 = source is confirmed but the flux is not reliable (see the flags)
*    1 = the source is not confirmed
*    0 = not observed
*Note (4): Bit flags of data quality in hexadecimal:
*      1 = not month confirmed.
*         This means that the period between the first detection and
*                the last detection is shorter than a month.
*      2 = saturated (not used this version);
*      4 = use events affected by the South Atlantic Anomaly (SAA)
*          (not used in this version),
*      8 = use edge events
*      Combined values are represented by the sum, e.g. 9 = 8 + 1 =
*      not month confirmed AND used edge events.
*Note (5): the value 1 means that the source is possibly more extended than
*     the point spread function (>15.6arcsec)
*Note (6): average of major and minor axes of source extent = (a+b)/2

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

C  Loading file 'irc.dat'	! AKARI/IRC All-Sky Survey Point Source Catalogue
*                                  (Version 1.0)

C  Format for file interpretation

    1 format(
     +  1X,I9,1X,A14,1X,F9.5,1X,F9.5,3X,F5.2,3X,F5.2,2X,F6.2,E11.4,
     +  E11.4,E10.3,E10.3,1X,I1,1X,I1,4X,A1,4X,A1,1X,I4,1X,I4,1X,I4,
     +  1X,I4,1X,I2,1X,I2,1X,I3,1X,I3,1X,I2,1X,I2,1X,F7.2,1X,F7.2,1X,
     +  I4,1X,I4,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'irc.dat')
      write(6,*) '....Loading file: irc.dat'
      do i__=1,870973
        read(1,'(A196)')ar__
        read(ar__,1)
     +  objID(i__),objName(i__),RAdeg(i__),DEdeg(i__),errMaj(i__),
     +  errMin(i__),errPA(i__),S09(i__),S18(i__),e_S09(i__),
     +  e_S18(i__),q_S09(i__),q_S18(i__),f09(i__),f18(i__),Ns09(i__),
     +  Ns18(i__),Np09(i__),Np18(i__),M09(i__),M18(i__),Nd09(i__),
     +  Nd18(i__),X09(i__),X18(i__),r09(i__),r18(i__),Ndet(i__),
     +  N09(i__),N18(i__)
c    ..............Just test output...........
        write(6,1)
     +  objID(i__),objName(i__),RAdeg(i__),DEdeg(i__),errMaj(i__),
     +  errMin(i__),errPA(i__),S09(i__),S18(i__),e_S09(i__),
     +  e_S18(i__),q_S09(i__),q_S18(i__),f09(i__),f18(i__),Ns09(i__),
     +  Ns18(i__),Np09(i__),Np18(i__),M09(i__),M18(i__),Nd09(i__),
     +  Nd18(i__),X09(i__),X18(i__),r09(i__),r18(i__),Ndet(i__),
     +  N09(i__),N18(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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