FORTRAN Generation
(/./ftp/cats/II/338)

Conversion of standardized ReadMe file for file /./ftp/cats/II/338 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-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. II/338          IRAS PSC/FSC Combined Catalogue               (Abrahamyan+ 2015)
*================================================================================
*The IRAS PSC/FSC Combined Catalogue
*     Abrahamyan H.V., Mickaelian A.M., Knyazyan A.V.
*    <Astron. & Computing 10, 99 (2015)>
*   =2015A&C....10...99A
*   =2015yCat.2338....0A
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! The IRAS PSC/FSC Combined Catalogue

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

C  J2000 position composed of: RAdeg DEdeg
      character*11  IRAS_PSC   (nr__) ! IRAS Point Source Catalog source name (II/125)
      integer*4     RApsc_h    (nr__) ! (h) ? RA hours (J2000)
      integer*4     RApsc_m    (nr__) ! (min) ? RA minutes (J2000)
      real*4        RApsc_s    (nr__) ! (sec) ? RA seconds (J2000)
      character*1   DEpsc__    (nr__) ! ? DEC sign (J2000)
      integer*4     DEpsc_d    (nr__) ! (deg) ? DEC degrees (J2000)
      integer*4     DEpsc_m    (nr__) ! (arcmin) ? DEC arcminutes (J2000)
      integer*4     DEpsc_s    (nr__) ! (arcsec) ? DEC arcseconds (J2000)
      real*8        RApsc_deg  (nr__) ! (deg) ? RA degrees (J2000)
      real*8        DEpsc_deg  (nr__) ! (deg) ? DEC degrees (J2000)
      integer*4     pscMaj     (nr__) ! (arcsec) ? Uncertainty ellipse semi-major axis
      integer*4     pscMin     (nr__) ! (arcsec) ? Uncertainty ellipse semi-minor axis
      integer*4     pscPA      (nr__) ! (deg) [0/180]? Uncertainty ellipse position angle
      real*4        psc12      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 12um
      real*4        psc25      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 25um
      real*4        psc60      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 60um
      real*4        psc100     (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 100um
      integer*4     q_psc12    (nr__) ! [0/3]? Flux density quality 12um (1)
      integer*4     q_psc25    (nr__) ! [0/3]? Flux density quality 25um (1)
      integer*4     q_psc60    (nr__) ! [0/3]? Flux density quality 60um (1)
      integer*4     q_psc100   (nr__) ! [0/3]? Flux density quality 100um (1)
      character*12  IRAS_FSC   (nr__) ! IRAS Faint Source Catalog source name (II/156)
      integer*4     RAfsc_h    (nr__) ! (h) ? RA hours (J2000)
      integer*4     RAfsc_m    (nr__) ! (min) ? RA minutes (J2000)
      real*4        RAfsc_s    (nr__) ! (sec) ? RA seconds (J2000)
      character*1   DEfsc__    (nr__) ! ? DEC sign (J2000)
      integer*4     DEfsc_d    (nr__) ! (deg) ? DEC degrees (J2000)
      integer*4     DEfsc_m    (nr__) ! (arcmin) ? DEC arcminutes (J2000)
      integer*4     DEfsc_s    (nr__) ! (arcsec) ? DEC arcseconds (J2000)
      real*8        RAfsc_deg  (nr__) ! (deg) ? RA degrees (J2000)
      real*8        DEfsc_deg  (nr__) ! (deg) ? DEC degrees (J2000)
      integer*4     fscMaj     (nr__) ! (arcsec) ? Uncertainty ellipse major axis
      integer*4     fscMin     (nr__) ! (arcsec) ? Uncertainty ellipse minor axis
      integer*4     fscPA      (nr__) ! (deg) [0/180]? Uncertainty ellipse position angle
      real*4        fsc12      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 12um
      real*4        fsc25      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 25um
      real*4        fsc60      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 60um
      real*4        fsc100     (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 100um
      integer*4     q_fsc12    (nr__) ! [0/3]? Flux density quality 12um (1)
      integer*4     q_fsc25    (nr__) ! [0/3]? Flux density quality 25um (1)
      integer*4     q_fsc60    (nr__) ! [0/3]? Flux density quality 60um (1)
      integer*4     q_fsc100   (nr__) ! [0/3]? Flux density quality 100um (1)
      real*8        RAir_deg   (nr__) ! (deg) PSC/FSC best RA degrees (J2000)
      real*8        DEir_deg   (nr__) ! (deg) PSC/FSC best DEC degrees (J2000)
      character*7   P_F        (nr__) ! Flag for combined coordinate:
*                                    "PSC", "FSC", or "PSC/FSC"
      real*4        e_RAir_deg (nr__) ! (deg) ? Positional error for RAJ2000
      real*4        e_DEir_deg (nr__) ! (deg) ? Positional error for DEJ2000
      integer*4     irMaj      (nr__) ! (arcsec) Uncertainty ellipse Major axis
      integer*4     irMin      (nr__) ! (arcsec) Uncertainty ellipse Minor axis
      integer*4     irPA       (nr__) ! (deg) [0/180] Uncertainty ellipse Position Angle
      real*4        d_ir       (nr__) ! (arcsec) [0/991]? Distance between PSC and FSC
      real*4        k_ir       (nr__) ! [0/3]? Relative distance K (in {sigma}'s)
      real*4        Fir12      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 12um
      real*4        Fir25      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 25um
      real*4        Fir60      (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 60um
      real*4        Fir100     (nr__) ! (Jy) ? Average non-color corrected flux density
*                                    at 100um
      character*3   f_ir       (nr__) ! Flag for flux ("PSC" or "FSC")
      character*7   cl         (nr__) ! Classification (type: "star" or "galaxy")
      character*14  Akari_IRC  (nr__) ! AKARI IRC source name (HHMMSSs+DDMMSS)(II/297)
      real*8        RAirc_deg  (nr__) ! (deg) ? AKARI-IRC RA degrees (J2000)
      real*8        DEirc_deg  (nr__) ! (deg) ? AKARI-IRC DEC degrees(J2000)
      real*4        ircMaj     (nr__) ! (arcsec) ? Major axis of position error ellipse
      real*4        ircMin     (nr__) ! (arcsec) ? Minor axis of position error ellipse
      real*4        ircPA      (nr__) ! (deg) [0/180]? Position angle of Major axis
      real*4        d_irc      (nr__) ! (arcsec) ? Distance between PSC/FSC and AKARI IRC
      real*4        k_irc      (nr__) ! [0/3]? Relative distance k (in {sigma}'s)
      real*4        S09        (nr__) ! (Jy) ? Flux density in AKARI/S9W filter
      real*4        S18        (nr__) ! (Jy) ? Flux density in AKARI/L18W filter
      real*4        e_S09      (nr__) ! (Jy) ? Flux error in S9W
      real*4        e_S18      (nr__) ! (Jy) ? Flux error in L18W
      integer*4     q_S09      (nr__) ! [0/3]? Flux quality flag for S9W (2)
      integer*4     q_S18      (nr__) ! [0/3]? Flux quality flag for L18W (2)
      real*4        r09        (nr__) ! (arcsec) ? Radius of source extent in S9W
      real*4        r18        (nr__) ! (arcsec) ? Radius of source extent in L18W
      character*14  Akari_FIS  (nr__) ! AKARI FIS source name (HHMMSSs+DDMMSS)(II/298)
      real*8        RAfis_deg  (nr__) ! (deg) ? AKARI-FIS RA degrees (J2000)
      real*8        DEfis_deg  (nr__) ! (deg) ? AKARI-FIS DEC degrees (J2000)
      integer*4     fisErr     (nr__) ! (arcsec) ? Error RA and Dec
      real*4        d_fis      (nr__) ! (arcsec) ? Distance between PSC/FSC and AKARI FIS
      real*4        k_fis      (nr__) ! [0/3]? Relative distance k (in {sigma}'s)
      real*4        S65        (nr__) ! (Jy) ? Flux density in N60
      real*4        S90        (nr__) ! (Jy) ? Flux density in WIDE-S
      real*4        S140       (nr__) ! (Jy) ? Flux density in WIDE-L
      real*4        S160       (nr__) ! (Jy) ? Flux density in N160
      real*4        e_S65      (nr__) ! (Jy) ? Uncertainty in N60
      real*4        e_S90      (nr__) ! (Jy) ? Uncertainty in WIDE-S
      real*4        e_S140     (nr__) ! (Jy) ? Uncertainty in WIDE-L
      real*4        e_S160     (nr__) ! (Jy) ? Uncertainty in N160
      integer*4     q_S65      (nr__) ! [0,3]? quality flag for N60 (2)
      integer*4     q_S90      (nr__) ! [0,3]? quality flag for WIDE-S (2)
      integer*4     q_S140     (nr__) ! [0,3]? quality flag for WIDE-L (2)
      integer*4     q_S160     (nr__) ! [0,3]? quality flag for N160 (2)
      character*19  WISE       (nr__) ! WISE All-Sky Release Catalog name (II/311)
      integer*4     RAwis_h    (nr__) ! (h) ? WISE RA hours (J2000)
      integer*4     RAwis_m    (nr__) ! (min) ? WISE RA minutes (J2000)
      real*4        RAwis_s    (nr__) ! (sec) ? WISE RA seconds (J2000)
      character*1   DEwis__    (nr__) ! ? WISE DEC sign (J2000)
      integer*4     DEwis_d    (nr__) ! (deg) ? WISE DEC degrees (J2000)
      integer*4     DEwis_m    (nr__) ! (arcmin) ? WISE DEC arcminutes (J2000)
      real*4        DEwis_s    (nr__) ! (arcsec) ? WISE DEC arcseconds (J2000)
      real*8        RAwis_deg  (nr__) ! (deg) ? WISE RA degrees (J2000)
      real*8        DEwis_deg  (nr__) ! (deg) ? WISE DEC degrees (J2000)
      real*4        wisMaj     (nr__) ! (arcsec) [0/2]? Semi-major axis of the error ellipse
      real*4        wisMin     (nr__) ! (arcsec) [0/2]? Semi-minor axis of the error ellipse
      integer*4     wisPA      (nr__) ! (deg) [0/180]? Position angle of the error ellipse
      real*4        d_wis      (nr__) ! (arcsec) ? Distance between IRAS/AKARI best and WISE
      real*4        Jmag       (nr__) ! (mag) ? 2MASS J magnitude
      real*8        Jflx       (nr__) ! (Jy) ? 2MASS J-band flux
      real*4        Hmag       (nr__) ! (mag) ? 2MASS H magnitude
      real*8        Hflx       (nr__) ! (Jy) ? 2MASS H-band flux
      real*4        Kmag       (nr__) ! (mag) ? 2MASS Ks magnitude
      real*8        Kflx       (nr__) ! (Jy) ? 2MASS Ks-band flux
      real*4        W1mag      (nr__) ! (mag) ? WISE W1 magnitude (3.35um)
      real*8        W1flx      (nr__) ! (Jy) ? WISE W1-band flux (3.35um)
      real*4        W2mag      (nr__) ! (mag) ? WISE W2 magnitude (4.6um)
      real*8        W2flx      (nr__) ! (Jy) ? WISE W2-band flux (4.6um)
      real*4        W3mag      (nr__) ! (mag) ? WISE W3 magnitude (11.6um)
      real*8        W3flx      (nr__) ! (Jy) ? WISE W3-band flux (11.6um)
      real*4        W4mag      (nr__) ! (mag) ? WISE W4 magnitude (22.1um)
      real*8        W4flx      (nr__) ! (Jy) ? WISE W4-band flux (22.1um)
      real*4        e_Jmag     (nr__) ! (mag) ? Mean error on J magnitude
      real*4        e_Hmag     (nr__) ! (mag) ? Mean error on H magnitude
      real*4        e_Kmag     (nr__) ! (mag) ? Mean error on Ks magnitude
      real*4        e_W1mag    (nr__) ! (mag) ? Mean error on W1 magnitude (3.35um)
      real*4        e_W2mag    (nr__) ! (mag) ? Mean error on W2 magnitude (4.6um)
      real*4        e_W3mag    (nr__) ! (mag) ? Mean error on W3 magnitude (11.6um)
      real*4        e_W4mag    (nr__) ! (mag) ? Mean error on W4 magnitude (22.1um)
      real*4        d2M        (nr__) ! (arcsec) [0/3]? Distance between WISE and 2MASS source
      integer*4     cat        (nr__) ! [1/3]? Identfication Category: 1=Best,
*                                         2=Reliable, 3=Probable
      real*8        RAdeg      (nr__) ! (deg) Best coordinate RA degrees (J2000)
      real*8        DEdeg      (nr__) ! (deg) Best coordinate DEC degrees (J2000)
      real*8        Glon       (nr__) ! (deg) Best galactic longitude
      real*8        Glat       (nr__) ! (deg) Best galactic latitude
      character*9   flag       (nr__) ! Best coordinate flag (source catalogue)
*Note (1): IRAS qualities are 1=upper limit, 2=moderate, 3=high
*Note (2): AKARI qualities are:
*     0 = not observed
*     1 = the source is not confirmed
*     2 = source is confirmed but the flux is not reliable
*     3 = high quality (source confirmed and flux is reliable)

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

C  Loading file 'catalog.dat'	! The IRAS PSC/FSC Combined Catalogue

C  Format for file interpretation

    1 format(
     +  A11,2X,I2,1X,I2,1X,E4.1,3X,A1,I2,1X,I2,1X,I2,3X,F10.6,2X,
     +  F10.6,4X,I3,5X,I2,6X,I3,5X,E10.3,3X,E10.3,3X,E10.3,3X,E10.3,
     +  6X,I1,10X,I1,10X,I1,10X,I1,8X,A12,2X,I2,1X,I2,1X,E4.1,3X,A1,
     +  I2,1X,I2,1X,I2,3X,F10.6,3X,F10.6,6X,I2,6X,I2,7X,I3,5X,E9.4,3X,
     +  E9.4,3X,E10.4,3X,E10.4,6X,I1,9X,I1,9X,I1,9X,I1,9X,F10.6,5X,
     +  F10.6,7X,A7,4X,E8.3,3X,E8.3,5X,I3,6X,I2,7X,I3,5X,E6.2,7X,E4.2,
     +  8X,E10.3,3X,E10.3,3X,E10.3,3X,E10.3,3X,A3,4X,A7,2X,A14,3X,
     +  F9.5,3X,F9.5,4X,F4.2,5X,F4.2,4X,F6.2,6X,F6.2,11X,F4.2,10X,
     +  E9.4,3X,E9.4,3X,E11.3,3X,E10.3,5X,I1,7X,I1,6X,F5.2,4X,F5.2,3X,
     +  A14,3X,F8.4,3X,F8.4,4X,I1,8X,F6.2,12X,F4.2,9X,E11.3,3X,E11.3,
     +  3X,E11.3,3X,E11.3,3X,E10.3,3X,E10.3,3X,E10.3,3X,E10.3,6X,I1,
     +  7X,I1,7X,I1,8X,I1,6X,A19,3X,I2,1X,I2,1X,F6.3,3X,A1,I2,1X,I2,
     +  1X,F5.2,3X,F10.6,3X,F10.6,3X,F5.3,3X,F5.3,3X,I3,4X,F5.2,3X,
     +  F6.3,3X,F11.6,3X,F6.3,3X,F10.6,3X,F6.3,3X,F10.6,3X,F6.3,3X,
     +  F10.6,3X,F6.3,3X,F10.6,3X,F6.3,3X,F9.6,3X,F6.3,3X,F8.6,4X,
     +  F5.3,6X,F5.3,6X,F5.3,6X,F5.3,8X,F5.3,7X,F5.3,7X,F5.3,5X,F5.3,
     +  6X,I1,7X,F10.6,3X,F10.6,3X,F8.4,3X,F8.4,3X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,345162
        read(1,'(A1385)')ar__
        read(ar__,1)
     +  IRAS_PSC(i__),RApsc_h(i__),RApsc_m(i__),RApsc_s(i__),
     +  DEpsc__(i__),DEpsc_d(i__),DEpsc_m(i__),DEpsc_s(i__),
     +  RApsc_deg(i__),DEpsc_deg(i__),pscMaj(i__),pscMin(i__),
     +  pscPA(i__),psc12(i__),psc25(i__),psc60(i__),psc100(i__),
     +  q_psc12(i__),q_psc25(i__),q_psc60(i__),q_psc100(i__),
     +  IRAS_FSC(i__),RAfsc_h(i__),RAfsc_m(i__),RAfsc_s(i__),
     +  DEfsc__(i__),DEfsc_d(i__),DEfsc_m(i__),DEfsc_s(i__),
     +  RAfsc_deg(i__),DEfsc_deg(i__),fscMaj(i__),fscMin(i__),
     +  fscPA(i__),fsc12(i__),fsc25(i__),fsc60(i__),fsc100(i__),
     +  q_fsc12(i__),q_fsc25(i__),q_fsc60(i__),q_fsc100(i__),
     +  RAir_deg(i__),DEir_deg(i__),P_F(i__),e_RAir_deg(i__),
     +  e_DEir_deg(i__),irMaj(i__),irMin(i__),irPA(i__),d_ir(i__),
     +  k_ir(i__),Fir12(i__),Fir25(i__),Fir60(i__),Fir100(i__),
     +  f_ir(i__),cl(i__),Akari_IRC(i__),RAirc_deg(i__),
     +  DEirc_deg(i__),ircMaj(i__),ircMin(i__),ircPA(i__),d_irc(i__),
     +  k_irc(i__),S09(i__),S18(i__),e_S09(i__),e_S18(i__),q_S09(i__),
     +  q_S18(i__),r09(i__),r18(i__),Akari_FIS(i__),RAfis_deg(i__),
     +  DEfis_deg(i__),fisErr(i__),d_fis(i__),k_fis(i__),S65(i__),
     +  S90(i__),S140(i__),S160(i__),e_S65(i__),e_S90(i__),
     +  e_S140(i__),e_S160(i__),q_S65(i__),q_S90(i__),q_S140(i__),
     +  q_S160(i__),WISE(i__),RAwis_h(i__),RAwis_m(i__),RAwis_s(i__),
     +  DEwis__(i__),DEwis_d(i__),DEwis_m(i__),DEwis_s(i__),
     +  RAwis_deg(i__),DEwis_deg(i__),wisMaj(i__),wisMin(i__),
     +  wisPA(i__),d_wis(i__),Jmag(i__),Jflx(i__),Hmag(i__),Hflx(i__),
     +  Kmag(i__),Kflx(i__),W1mag(i__),W1flx(i__),W2mag(i__),
     +  W2flx(i__),W3mag(i__),W3flx(i__),W4mag(i__),W4flx(i__),
     +  e_Jmag(i__),e_Hmag(i__),e_Kmag(i__),e_W1mag(i__),e_W2mag(i__),
     +  e_W3mag(i__),e_W4mag(i__),d2M(i__),cat(i__),RAdeg(i__),
     +  DEdeg(i__),Glon(i__),Glat(i__),flag(i__)
        if(ar__(14:15) .EQ. '') RApsc_h(i__) = iNULL__
        if(ar__(17:18) .EQ. '') RApsc_m(i__) = iNULL__
        if(ar__(20:23) .EQ. '') RApsc_s(i__) = rNULL__
        if(ar__(28:29) .EQ. '') DEpsc_d(i__) = iNULL__
        if(ar__(31:32) .EQ. '') DEpsc_m(i__) = iNULL__
        if(ar__(34:35) .EQ. '') DEpsc_s(i__) = iNULL__
        if(ar__(39:48) .EQ. '') RApsc_deg(i__) = rNULL__
        if(ar__(51:60) .EQ. '') DEpsc_deg(i__) = rNULL__
        if(ar__(65:67) .EQ. '') pscMaj(i__) = iNULL__
        if(ar__(73:74) .EQ. '') pscMin(i__) = iNULL__
        if(ar__(81:83) .EQ. '') pscPA(i__) = iNULL__
        if(ar__(89:98) .EQ. '') psc12(i__) = rNULL__
        if(ar__(102:111) .EQ. '') psc25(i__) = rNULL__
        if(ar__(115:124) .EQ. '') psc60(i__) = rNULL__
        if(ar__(128:137) .EQ. '') psc100(i__) = rNULL__
        if(ar__(144:144) .EQ. '') q_psc12(i__) = iNULL__
        if(ar__(155:155) .EQ. '') q_psc25(i__) = iNULL__
        if(ar__(166:166) .EQ. '') q_psc60(i__) = iNULL__
        if(ar__(177:177) .EQ. '') q_psc100(i__) = iNULL__
        if(ar__(200:201) .EQ. '') RAfsc_h(i__) = iNULL__
        if(ar__(203:204) .EQ. '') RAfsc_m(i__) = iNULL__
        if(ar__(206:209) .EQ. '') RAfsc_s(i__) = rNULL__
        if(ar__(214:215) .EQ. '') DEfsc_d(i__) = iNULL__
        if(ar__(217:218) .EQ. '') DEfsc_m(i__) = iNULL__
        if(ar__(220:221) .EQ. '') DEfsc_s(i__) = iNULL__
        if(ar__(225:234) .EQ. '') RAfsc_deg(i__) = rNULL__
        if(ar__(238:247) .EQ. '') DEfsc_deg(i__) = rNULL__
        if(ar__(254:255) .EQ. '') fscMaj(i__) = iNULL__
        if(ar__(262:263) .EQ. '') fscMin(i__) = iNULL__
        if(ar__(271:273) .EQ. '') fscPA(i__) = iNULL__
        if(ar__(279:287) .EQ. '') fsc12(i__) = rNULL__
        if(ar__(291:299) .EQ. '') fsc25(i__) = rNULL__
        if(ar__(303:312) .EQ. '') fsc60(i__) = rNULL__
        if(ar__(316:325) .EQ. '') fsc100(i__) = rNULL__
        if(ar__(332:332) .EQ. '') q_fsc12(i__) = iNULL__
        if(ar__(342:342) .EQ. '') q_fsc25(i__) = iNULL__
        if(ar__(352:352) .EQ. '') q_fsc60(i__) = iNULL__
        if(ar__(362:362) .EQ. '') q_fsc100(i__) = iNULL__
        if(ar__(415:422) .EQ. '') e_RAir_deg(i__) = rNULL__
        if(ar__(426:433) .EQ. '') e_DEir_deg(i__) = rNULL__
        if(ar__(465:470) .EQ. '') d_ir(i__) = rNULL__
        if(ar__(478:481) .EQ. '') k_ir(i__) = rNULL__
        if(ar__(490:499) .EQ. '') Fir12(i__) = rNULL__
        if(ar__(503:512) .EQ. '') Fir25(i__) = rNULL__
        if(ar__(516:525) .EQ. '') Fir60(i__) = rNULL__
        if(ar__(529:538) .EQ. '') Fir100(i__) = rNULL__
        if(ar__(575:583) .EQ. '') RAirc_deg(i__) = rNULL__
        if(ar__(587:595) .EQ. '') DEirc_deg(i__) = rNULL__
        if(ar__(600:603) .EQ. '') ircMaj(i__) = rNULL__
        if(ar__(609:612) .EQ. '') ircMin(i__) = rNULL__
        if(ar__(617:622) .EQ. '') ircPA(i__) = rNULL__
        if(ar__(629:634) .EQ. '') d_irc(i__) = rNULL__
        if(ar__(646:649) .EQ. '') k_irc(i__) = rNULL__
        if(ar__(660:668) .EQ. '') S09(i__) = rNULL__
        if(ar__(672:680) .EQ. '') S18(i__) = rNULL__
        if(ar__(684:694) .EQ. '') e_S09(i__) = rNULL__
        if(ar__(698:707) .EQ. '') e_S18(i__) = rNULL__
        if(ar__(713:713) .EQ. '') q_S09(i__) = iNULL__
        if(ar__(721:721) .EQ. '') q_S18(i__) = iNULL__
        if(ar__(728:732) .EQ. '') r09(i__) = rNULL__
        if(ar__(737:741) .EQ. '') r18(i__) = rNULL__
        if(ar__(762:769) .EQ. '') RAfis_deg(i__) = rNULL__
        if(ar__(773:780) .EQ. '') DEfis_deg(i__) = rNULL__
        if(ar__(785:785) .EQ. '') fisErr(i__) = iNULL__
        if(ar__(794:799) .EQ. '') d_fis(i__) = rNULL__
        if(ar__(812:815) .EQ. '') k_fis(i__) = rNULL__
        if(ar__(825:835) .EQ. '') S65(i__) = rNULL__
        if(ar__(839:849) .EQ. '') S90(i__) = rNULL__
        if(ar__(853:863) .EQ. '') S140(i__) = rNULL__
        if(ar__(867:877) .EQ. '') S160(i__) = rNULL__
        if(ar__(881:890) .EQ. '') e_S65(i__) = rNULL__
        if(ar__(894:903) .EQ. '') e_S90(i__) = rNULL__
        if(ar__(907:916) .EQ. '') e_S140(i__) = rNULL__
        if(ar__(920:929) .EQ. '') e_S160(i__) = rNULL__
        if(ar__(936:936) .EQ. '') q_S65(i__) = iNULL__
        if(ar__(944:944) .EQ. '') q_S90(i__) = iNULL__
        if(ar__(952:952) .EQ. '') q_S140(i__) = iNULL__
        if(ar__(961:961) .EQ. '') q_S160(i__) = iNULL__
        if(ar__(990:991) .EQ. '') RAwis_h(i__) = iNULL__
        if(ar__(993:994) .EQ. '') RAwis_m(i__) = iNULL__
        if(ar__(996:1001) .EQ. '') RAwis_s(i__) = rNULL__
        if(ar__(1006:1007) .EQ. '') DEwis_d(i__) = iNULL__
        if(ar__(1009:1010) .EQ. '') DEwis_m(i__) = iNULL__
        if(ar__(1012:1016) .EQ. '') DEwis_s(i__) = rNULL__
        if(ar__(1020:1029) .EQ. '') RAwis_deg(i__) = rNULL__
        if(ar__(1033:1042) .EQ. '') DEwis_deg(i__) = rNULL__
        if(ar__(1046:1050) .EQ. '') wisMaj(i__) = rNULL__
        if(ar__(1054:1058) .EQ. '') wisMin(i__) = rNULL__
        if(ar__(1062:1064) .EQ. '') wisPA(i__) = iNULL__
        if(ar__(1069:1073) .EQ. '') d_wis(i__) = rNULL__
        if(ar__(1077:1082) .EQ. '') Jmag(i__) = rNULL__
        if(ar__(1086:1096) .EQ. '') Jflx(i__) = rNULL__
        if(ar__(1100:1105) .EQ. '') Hmag(i__) = rNULL__
        if(ar__(1109:1118) .EQ. '') Hflx(i__) = rNULL__
        if(ar__(1122:1127) .EQ. '') Kmag(i__) = rNULL__
        if(ar__(1131:1140) .EQ. '') Kflx(i__) = rNULL__
        if(ar__(1144:1149) .EQ. '') W1mag(i__) = rNULL__
        if(ar__(1153:1162) .EQ. '') W1flx(i__) = rNULL__
        if(ar__(1166:1171) .EQ. '') W2mag(i__) = rNULL__
        if(ar__(1175:1184) .EQ. '') W2flx(i__) = rNULL__
        if(ar__(1188:1193) .EQ. '') W3mag(i__) = rNULL__
        if(ar__(1197:1205) .EQ. '') W3flx(i__) = rNULL__
        if(ar__(1209:1214) .EQ. '') W4mag(i__) = rNULL__
        if(ar__(1218:1225) .EQ. '') W4flx(i__) = rNULL__
        if(ar__(1230:1234) .EQ. '') e_Jmag(i__) = rNULL__
        if(ar__(1241:1245) .EQ. '') e_Hmag(i__) = rNULL__
        if(ar__(1252:1256) .EQ. '') e_Kmag(i__) = rNULL__
        if(ar__(1263:1267) .EQ. '') e_W1mag(i__) = rNULL__
        if(ar__(1276:1280) .EQ. '') e_W2mag(i__) = rNULL__
        if(ar__(1288:1292) .EQ. '') e_W3mag(i__) = rNULL__
        if(ar__(1300:1304) .EQ. '') e_W4mag(i__) = rNULL__
        if(ar__(1310:1314) .EQ. '') d2M(i__) = rNULL__
        if(ar__(1321:1321) .EQ. '') cat(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  IRAS_PSC(i__),RApsc_h(i__),RApsc_m(i__),RApsc_s(i__),
     +  DEpsc__(i__),DEpsc_d(i__),DEpsc_m(i__),DEpsc_s(i__),
     +  RApsc_deg(i__),DEpsc_deg(i__),pscMaj(i__),pscMin(i__),
     +  pscPA(i__),psc12(i__),psc25(i__),psc60(i__),psc100(i__),
     +  q_psc12(i__),q_psc25(i__),q_psc60(i__),q_psc100(i__),
     +  IRAS_FSC(i__),RAfsc_h(i__),RAfsc_m(i__),RAfsc_s(i__),
     +  DEfsc__(i__),DEfsc_d(i__),DEfsc_m(i__),DEfsc_s(i__),
     +  RAfsc_deg(i__),DEfsc_deg(i__),fscMaj(i__),fscMin(i__),
     +  fscPA(i__),fsc12(i__),fsc25(i__),fsc60(i__),fsc100(i__),
     +  q_fsc12(i__),q_fsc25(i__),q_fsc60(i__),q_fsc100(i__),
     +  RAir_deg(i__),DEir_deg(i__),P_F(i__),e_RAir_deg(i__),
     +  e_DEir_deg(i__),irMaj(i__),irMin(i__),irPA(i__),d_ir(i__),
     +  k_ir(i__),Fir12(i__),Fir25(i__),Fir60(i__),Fir100(i__),
     +  f_ir(i__),cl(i__),Akari_IRC(i__),RAirc_deg(i__),
     +  DEirc_deg(i__),ircMaj(i__),ircMin(i__),ircPA(i__),d_irc(i__),
     +  k_irc(i__),S09(i__),S18(i__),e_S09(i__),e_S18(i__),q_S09(i__),
     +  q_S18(i__),r09(i__),r18(i__),Akari_FIS(i__),RAfis_deg(i__),
     +  DEfis_deg(i__),fisErr(i__),d_fis(i__),k_fis(i__),S65(i__),
     +  S90(i__),S140(i__),S160(i__),e_S65(i__),e_S90(i__),
     +  e_S140(i__),e_S160(i__),q_S65(i__),q_S90(i__),q_S140(i__),
     +  q_S160(i__),WISE(i__),RAwis_h(i__),RAwis_m(i__),RAwis_s(i__),
     +  DEwis__(i__),DEwis_d(i__),DEwis_m(i__),DEwis_s(i__),
     +  RAwis_deg(i__),DEwis_deg(i__),wisMaj(i__),wisMin(i__),
     +  wisPA(i__),d_wis(i__),Jmag(i__),Jflx(i__),Hmag(i__),Hflx(i__),
     +  Kmag(i__),Kflx(i__),W1mag(i__),W1flx(i__),W2mag(i__),
     +  W2flx(i__),W3mag(i__),W3flx(i__),W4mag(i__),W4flx(i__),
     +  e_Jmag(i__),e_Hmag(i__),e_Kmag(i__),e_W1mag(i__),e_W2mag(i__),
     +  e_W3mag(i__),e_W4mag(i__),d2M(i__),cat(i__),RAdeg(i__),
     +  DEdeg(i__),Glon(i__),Glat(i__),flag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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