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

Conversion of standardized ReadMe file for file /./ftp/cats/II/341 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-16
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/341              VPHAS+ DR2 survey                            (Drew+, 2016)
*================================================================================
*The VST Photometric Halpha Survey of the Southern Galactic Plane and Bulge
*(VPHAS+). Second VPHAS data release (DR2).
*     Drew J.E., Gonzalez-Solares E., Greimel R., Irwin M.J., Kupcu Yoldas A.,
*     Lewis J., Barentsen G., Eisloffel J., Farnhill H.J., Martin W.E.,
*     Walsh J.R., Walton N.A., Mohr-Smith M., Raddi R., Sale S.E., Wright N.J.,
*     Groot P., Barlow M.J., Corradi R.L.M., Drake J.J., Fabregat J., Frew D.J.,
*     Gansicke B.T., Knigge C., Mampaso A., Morris R.A.H., Naylor T.,
*     Parker Q.A., Phillipps S., Ruhland C., Steeghs D., Unruh Y.C., Vink J.S.,
*     Wesson R., Zijlstra A.A.
*    <VPHAS+ survey, (2016)>
*    =2016yCat.2341....0D
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'vphasp.sam'	! VPHAS+ catalog, sample output
                                      (on a total of 319,126,837 sources)

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

C  J2000 position composed of: RAdeg DEdeg
      character*14  sourceID   (nr__) ! Unique VPHAS source identifier in the format
*                                   field-extension-number (sourceID)
      real*8        RAdeg      (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
*                                   (RAJ2000)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000)
*                                   (DEJ2000)
      character*27  VPHASDR2   (nr__) ! Position-based source name. Not a unique
*                                   identifier (VPHASDR2 JHHMMSS.s+DDMMSS.s)
*                                   (name)
      character*14  PrimaryID  (nr__) ! SourceID of the preferred catalogue entry
*                                   for this source (primaryID)
      integer*4     fPrimary   (nr__) ! [0/1] Flag =1 if the entry is the preferred
*                                   observation of a unique source
*                                   (Primary_source)
      integer*4     Nobs       (nr__) ! [1-6] Number of observations of this source
*                                   in the survey (nObs)
      integer*4     clean      (nr__) ! [0/1] =1 if cleang, cleanr2, cleanha, cleanr
*                                    and cleani are true (clean)
      real*4        u_g        (nr__) ! (mag) ? (u-g) colour index, formed by subtracting
*                                     columns u and g (u_g)
      real*4        g_r2       (nr__) ! (mag) ? (g-r2) colour index, formed by subtracting
*                                     columns g and r2 (g_r2)
      real*4        r_i        (nr__) ! (mag) ? (r-i) colour index, formed by subtracting
*                                     columns r and i (r_i)
      real*4        r_ha       (nr__) ! (mag) ? (r-Halpha) colour index, formed by
*                                    subtracting columns r and ha (r_ha)
      integer*4     cleanu     (nr__) ! [0/1] =1 if the u-band detection is
*                                    significant and the PSF fit good (clean_u)
      real*4        umag       (nr__) ! (mag) ? Default u-band VEGA magnitude obtained using
*                                   PSF fitting (u)
      real*4        umagAB     (nr__) ! (mag) ? Default u-band AB magnitude obtained using
*                                   PSF fitting (u_AB)
      real*4        e_umag     (nr__) ! (mag) ? Statistical uncertainty for u,
*                                   not including systematics (err_u)
      real*4        chiu       (nr__) ! ? Goodness of the PSF fit in u,
*                                    as determined by DAOPHOT ALLSTAR (chi_u)
      character*12  warningu   (nr__) ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_u) (1)
      real*4        umagap     (nr__) ! (mag) ? u-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_u)
      real*4        umagapAB   (nr__) ! (mag) ? u-band AB magnitude obtained using
*                                   aperture photometry (aperMag_u_AB)
      real*4        e_umagap   (nr__) ! (mag) ? Statistical uncertainty for umagap,
*                                   not including systematics (aperMagErr_u)
      real*4        snru       (nr__) ! ? Ratio of the signal in the aperture relative
*                                    to the background (snr_u)
      real*4        umaglim    (nr__) ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_u)
      real*4        PSFFWHMu   (nr__) ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_u)
      real*8        MJDu       (nr__) ! (d) ? Modified Julian Date at the start of the
*                                   u-band exposure (mjd_u)
      character*23  detIDu     (nr__) ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_u)
      integer*4     cleang     (nr__) ! [0/1] =1 if the g-band detection is
*                                   significant and the PSF fit good (clean_g)
      real*4        gmag       (nr__) ! (mag) ? Default g-band VEGA magnitude obtained using
*                                   PSF fitting (g)
      real*4        gmagAB     (nr__) ! (mag) ? Default g-band AB magnitude obtained using
*                                   PSF fitting (g_AB)
      real*4        e_gmag     (nr__) ! (mag) ? Statistical uncertainty for g,
*                                   not including systematics (err_g)
      real*4        chig       (nr__) ! ? Goodness of the PSF fit in g,
*                                   as determined by DAOPHOT ALLSTAR (chi_g)
      character*12  warningg   (nr__) ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_g) (1)
      real*4        gmagap     (nr__) ! (mag) ? g-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_g)
      real*4        gmagapAB   (nr__) ! (mag) ? g-band AB magnitude obtained using
*                                   aperture photometry (aperMag_g_AB)
      real*4        e_gmagap   (nr__) ! (mag) ? Statistical uncertainty for gmagap,
*                                   not including systematics (aperMagErr_g)
      real*4        snrg       (nr__) ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_g)
      real*4        gmaglim    (nr__) ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_g)
      real*4        PSFFWHMg   (nr__) ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_g)
      real*8        MJDg       (nr__) ! (d) ? Modified Julian Date at the start of the
*                                   g-band exposure (mjd_g)
      character*23  detIDg     (nr__) ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_g)
      integer*4     cleanr2    (nr__) ! [0/1] =1 if the r2-band detection is
*                                   significant and the PSF fit good (clean_r2)
      real*4        r2mag      (nr__) ! (mag) ? Second r-band VEGA magnitude obtained using
*                                   PSF fitting (r2)
      real*4        r2magAB    (nr__) ! (mag) ? Second r-band AB magnitude obtained using
*                                   PSF fitting (r2_AB)
      real*4        e_r2mag    (nr__) ! (mag) ? Statistical uncertainty for r2,
*                                   not including systematics (err_r2)
      real*4        chir2      (nr__) ! ? Goodness of the PSF fit in r2,
*                                   as determined by DAOPHOT ALLSTAR (chi_r2)
      character*12  warningr2  (nr__) ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_r2) (1)
      real*4        r2magap    (nr__) ! (mag) ? r2-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_r2)
      real*4        r2magapAB  (nr__) ! (mag) ? r2-band AB magnitude obtained using
*                                   aperture photometry (aperMag_r2_AB)
      real*4        e_r2magap  (nr__) ! (mag) ? Statistical uncertainty for r2magap,
*                                   not including systematics (aperMagErr_r2)
      real*4        snrr2      (nr__) ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_r2)
      real*4        r2maglim   (nr__) ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_r2)
      real*4        PSFFWHMr2  (nr__) ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_r2)
      real*8        MJDr2      (nr__) ! (d) ? Modified Julian Date at the start of the
*                                   r2-band exposure (mjd_r2)
      character*23  detIDr2    (nr__) ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_r2)
      integer*4     cleanHa    (nr__) ! [0/1] =1 if the H-alpha detection is
*                                   significant and the PSF fit good (clean_ha)
      real*4        Hamag      (nr__) ! (mag) ? Default H-alpha VEGA magnitude obtained
*                                   using PSF fitting (ha)
      real*4        e_Hamag    (nr__) ! (mag) ? Statistical uncertainty for ha,
*                                   not including systematics (err_ha)
      real*4        chiHa      (nr__) ! ? Goodness of the PSF fit in H-alpha,
*                                   as determined by DAOPHOT (chi_ha)
      character*9   warningHa  (nr__) ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_ha) (1)
      real*4        Hamagap    (nr__) ! (mag) ? H-alpha VEGA magnitude obtained using
*                                   aperture photometry (aperMag_ha)
      real*4        e_Hamagap  (nr__) ! (mag) ? Statistical uncertainty for Hamagap,
*                                   not including systematics (aperMagErr_ha)
      real*4        snrHa      (nr__) ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_ha)
      real*4        Hamaglim   (nr__) ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_ha)
      real*4        PSFFWHMHa  (nr__) ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_ha)
      real*8        MJDHa      (nr__) ! (d) ? Modified Julian Date at the start of the
*                                   H-alpha exposure (mjd_ha)
      character*23  detIDHa    (nr__) ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_ha)
      integer*4     cleanr     (nr__) ! [0/1] =1 if the r-band detection is
*                                   significant and the PSF fit good (clean_r)
      real*4        rmag       (nr__) ! (mag) ? Default r-band VEGA magnitude obtained using
*                                   PSF fitting (r)
      real*4        rmagAB     (nr__) ! (mag) ? Default r-band AB magnitude obtained using
*                                   PSF fitting (r_AB)
      real*4        e_rmag     (nr__) ! (mag) ? Statistical uncertainty for r,
*                                   not including systematics (err_r)
      real*4        chir       (nr__) ! ? Goodness of the PSF fit in r,
*                                   as determined by DAOPHOT ALLSTAR (chi_r)
      character*9   warningr   (nr__) ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_r) (1)
      real*4        rmagap     (nr__) ! (mag) ? r-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_r)
      real*4        rmagapAB   (nr__) ! (mag) ? r-band AB magnitude obtained using
*                                   aperture photometry (aperMag_r_AB)
      real*4        e_rmagap   (nr__) ! (mag) ? Statistical uncertainty for rmagap,
*                                   not including systematics (aperMagErr_r)
      real*4        snrr       (nr__) ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_r)
      real*4        rmaglim    (nr__) ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_r)
      real*4        PSFFWHMr   (nr__) ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_r)
      real*8        MJDr       (nr__) ! (d) ? Modified Julian Date at the start of the
*                                   r-band exposure (mjd_r)
      character*23  detIDr     (nr__) ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_r)
      integer*4     cleani     (nr__) ! [0/1] =1 if the i-band detection is
*                                   significant and the PSF fit good (clean_i)
      real*4        imag       (nr__) ! (mag) ? Default i-band VEGA magnitude obtained using
*                                   PSF fitting (i)
      real*4        imagAB     (nr__) ! (mag) ? Default i-band AB magnitude obtained using
*                                   PSF fitting (i_AB)
      real*4        e_imag     (nr__) ! (mag) ? Statistical uncertainty for i,
*                                   not including systematics (err_i)
      real*4        chii       (nr__) ! ? Goodness of the PSF fit in i,
*                                   as determined by DAOPHOT ALLSTAR (chi_i)
      character*9   warningi   (nr__) ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_i) (1)
      real*4        imagap     (nr__) ! (mag) ? i-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_i)
      real*4        imagapAB   (nr__) ! (mag) ? i-band AB magnitude obtained using
*                                   aperture photometry (aperMag_i_AB)
      real*4        e_imagap   (nr__) ! (mag) ? Statistical uncertainty for imagap,
*                                   not including systematics (aperMagErr_i)
      real*4        snri       (nr__) ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_i)
      real*4        imaglim    (nr__) ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_i)
      real*4        PSFFWHMi   (nr__) ! (arcsec) Typical PSF FWHM in the CCD frame (psffwhm_i)
      real*8        MJDi       (nr__) ! (d) Modified Julian Date at the start of the
*                                   i-band exposure (mjd_i)
      character*23  detIDi     (nr__) ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_i)
      character*5   Field      (nr__) ! Survey-specific identifier of the telescope
*                                   pointing (field)
      integer*4     Ext        (nr__) ! [1/32] OmegaCAM CCD extension used (ext)
      real*4        nbDist     (nr__) ! (arcsec) Distance to the nearest neighbour detected in
*Note (1): Possible values include
*   'No_error', 'Off_image', 'Too_faint', and 'Bad_fit'.
*   The magnitude values are left empty when an error occurred in catalogue
*    generation. This column is for information only.

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

C  Loading file 'vphasp.sam'	! VPHAS+ catalog, sample output
*                                      (on a total of 319,126,837 sources)

C  Format for file interpretation

    1 format(
     +  A14,1X,F9.5,1X,F9.5,1X,A27,1X,A14,1X,I1,1X,I1,1X,I1,1X,F5.2,
     +  1X,F5.2,1X,F5.2,1X,F5.2,1X,I1,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.1,
     +  1X,A12,1X,F5.2,1X,F5.2,1X,F4.2,1X,F6.1,1X,F4.1,1X,F4.2,1X,
     +  F8.2,1X,A23,1X,I1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.1,1X,A12,1X,
     +  F5.2,1X,F5.2,1X,F4.2,1X,F6.1,1X,F4.1,1X,F4.2,1X,F8.2,1X,A23,
     +  1X,I1,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.1,1X,A12,1X,F5.2,1X,F5.2,
     +  1X,F4.2,1X,F6.1,1X,F4.1,1X,F4.2,1X,F8.2,1X,A23,1X,I1,1X,F5.2,
     +  1X,F4.2,1X,F4.1,1X,A9,1X,F5.2,1X,F4.2,1X,F6.1,1X,F4.1,1X,F4.2,
     +  1X,F8.2,1X,A23,1X,I1,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.1,1X,A9,1X,
     +  F5.2,1X,F5.2,1X,F4.2,1X,F6.1,1X,F4.1,1X,F4.2,1X,F8.2,1X,A23,
     +  1X,I1,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.1,1X,A9,1X,F5.2,1X,F5.2,
     +  1X,F4.2,1X,F6.1,1X,F4.1,1X,F4.2,1X,F8.2,1X,A23,1X,A5,1X,I2,1X,
     +  F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'vphasp.sam')
      write(6,*) '....Loading file: vphasp.sam'
      do i__=1,1000
        read(1,'(A726)')ar__
        read(ar__,1)
     +  sourceID(i__),RAdeg(i__),DEdeg(i__),VPHASDR2(i__),
     +  PrimaryID(i__),fPrimary(i__),Nobs(i__),clean(i__),u_g(i__),
     +  g_r2(i__),r_i(i__),r_ha(i__),cleanu(i__),umag(i__),
     +  umagAB(i__),e_umag(i__),chiu(i__),warningu(i__),umagap(i__),
     +  umagapAB(i__),e_umagap(i__),snru(i__),umaglim(i__),
     +  PSFFWHMu(i__),MJDu(i__),detIDu(i__),cleang(i__),gmag(i__),
     +  gmagAB(i__),e_gmag(i__),chig(i__),warningg(i__),gmagap(i__),
     +  gmagapAB(i__),e_gmagap(i__),snrg(i__),gmaglim(i__),
     +  PSFFWHMg(i__),MJDg(i__),detIDg(i__),cleanr2(i__),r2mag(i__),
     +  r2magAB(i__),e_r2mag(i__),chir2(i__),warningr2(i__),
     +  r2magap(i__),r2magapAB(i__),e_r2magap(i__),snrr2(i__),
     +  r2maglim(i__),PSFFWHMr2(i__),MJDr2(i__),detIDr2(i__),
     +  cleanHa(i__),Hamag(i__),e_Hamag(i__),chiHa(i__),
     +  warningHa(i__),Hamagap(i__),e_Hamagap(i__),snrHa(i__),
     +  Hamaglim(i__),PSFFWHMHa(i__),MJDHa(i__),detIDHa(i__),
     +  cleanr(i__),rmag(i__),rmagAB(i__),e_rmag(i__),chir(i__),
     +  warningr(i__),rmagap(i__),rmagapAB(i__),e_rmagap(i__),
     +  snrr(i__),rmaglim(i__),PSFFWHMr(i__),MJDr(i__),detIDr(i__),
     +  cleani(i__),imag(i__),imagAB(i__),e_imag(i__),chii(i__),
     +  warningi(i__),imagap(i__),imagapAB(i__),e_imagap(i__),
     +  snri(i__),imaglim(i__),PSFFWHMi(i__),MJDi(i__),detIDi(i__),
     +  Field(i__),Ext(i__),nbDist(i__)
        if(ar__(85:89) .EQ. '') u_g(i__) = rNULL__
        if(ar__(91:95) .EQ. '') g_r2(i__) = rNULL__
        if(ar__(97:101) .EQ. '') r_i(i__) = rNULL__
        if(ar__(103:107) .EQ. '') r_ha(i__) = rNULL__
        if(ar__(111:115) .EQ. '') umag(i__) = rNULL__
        if(ar__(117:121) .EQ. '') umagAB(i__) = rNULL__
        if(ar__(123:126) .EQ. '') e_umag(i__) = rNULL__
        if(ar__(128:132) .EQ. '') chiu(i__) = rNULL__
        if(ar__(147:151) .EQ. '') umagap(i__) = rNULL__
        if(ar__(153:157) .EQ. '') umagapAB(i__) = rNULL__
        if(ar__(159:162) .EQ. '') e_umagap(i__) = rNULL__
        if(ar__(164:169) .EQ. '') snru(i__) = rNULL__
        if(ar__(171:174) .EQ. '') umaglim(i__) = rNULL__
        if(ar__(176:179) .EQ. '') PSFFWHMu(i__) = rNULL__
        if(ar__(181:188) .EQ. '') MJDu(i__) = rNULL__
        if(ar__(216:220) .EQ. '') gmag(i__) = rNULL__
        if(ar__(222:226) .EQ. '') gmagAB(i__) = rNULL__
        if(ar__(228:232) .EQ. '') e_gmag(i__) = rNULL__
        if(ar__(234:237) .EQ. '') chig(i__) = rNULL__
        if(ar__(252:256) .EQ. '') gmagap(i__) = rNULL__
        if(ar__(258:262) .EQ. '') gmagapAB(i__) = rNULL__
        if(ar__(264:267) .EQ. '') e_gmagap(i__) = rNULL__
        if(ar__(269:274) .EQ. '') snrg(i__) = rNULL__
        if(ar__(276:279) .EQ. '') gmaglim(i__) = rNULL__
        if(ar__(281:284) .EQ. '') PSFFWHMg(i__) = rNULL__
        if(ar__(286:293) .EQ. '') MJDg(i__) = rNULL__
        if(ar__(321:325) .EQ. '') r2mag(i__) = rNULL__
        if(ar__(327:331) .EQ. '') r2magAB(i__) = rNULL__
        if(ar__(333:336) .EQ. '') e_r2mag(i__) = rNULL__
        if(ar__(338:341) .EQ. '') chir2(i__) = rNULL__
        if(ar__(356:360) .EQ. '') r2magap(i__) = rNULL__
        if(ar__(362:366) .EQ. '') r2magapAB(i__) = rNULL__
        if(ar__(368:371) .EQ. '') e_r2magap(i__) = rNULL__
        if(ar__(373:378) .EQ. '') snrr2(i__) = rNULL__
        if(ar__(380:383) .EQ. '') r2maglim(i__) = rNULL__
        if(ar__(385:388) .EQ. '') PSFFWHMr2(i__) = rNULL__
        if(ar__(390:397) .EQ. '') MJDr2(i__) = rNULL__
        if(ar__(425:429) .EQ. '') Hamag(i__) = rNULL__
        if(ar__(431:434) .EQ. '') e_Hamag(i__) = rNULL__
        if(ar__(436:439) .EQ. '') chiHa(i__) = rNULL__
        if(ar__(451:455) .EQ. '') Hamagap(i__) = rNULL__
        if(ar__(457:460) .EQ. '') e_Hamagap(i__) = rNULL__
        if(ar__(462:467) .EQ. '') snrHa(i__) = rNULL__
        if(ar__(469:472) .EQ. '') Hamaglim(i__) = rNULL__
        if(ar__(474:477) .EQ. '') PSFFWHMHa(i__) = rNULL__
        if(ar__(479:486) .EQ. '') MJDHa(i__) = rNULL__
        if(ar__(514:518) .EQ. '') rmag(i__) = rNULL__
        if(ar__(520:524) .EQ. '') rmagAB(i__) = rNULL__
        if(ar__(526:529) .EQ. '') e_rmag(i__) = rNULL__
        if(ar__(531:534) .EQ. '') chir(i__) = rNULL__
        if(ar__(546:550) .EQ. '') rmagap(i__) = rNULL__
        if(ar__(552:556) .EQ. '') rmagapAB(i__) = rNULL__
        if(ar__(558:561) .EQ. '') e_rmagap(i__) = rNULL__
        if(ar__(563:568) .EQ. '') snrr(i__) = rNULL__
        if(ar__(570:573) .EQ. '') rmaglim(i__) = rNULL__
        if(ar__(575:578) .EQ. '') PSFFWHMr(i__) = rNULL__
        if(ar__(580:587) .EQ. '') MJDr(i__) = rNULL__
        if(ar__(615:619) .EQ. '') imag(i__) = rNULL__
        if(ar__(621:625) .EQ. '') imagAB(i__) = rNULL__
        if(ar__(627:630) .EQ. '') e_imag(i__) = rNULL__
        if(ar__(632:635) .EQ. '') chii(i__) = rNULL__
        if(ar__(647:651) .EQ. '') imagap(i__) = rNULL__
        if(ar__(653:657) .EQ. '') imagapAB(i__) = rNULL__
        if(ar__(659:662) .EQ. '') e_imagap(i__) = rNULL__
        if(ar__(664:669) .EQ. '') snri(i__) = rNULL__
        if(ar__(671:674) .EQ. '') imaglim(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  sourceID(i__),RAdeg(i__),DEdeg(i__),VPHASDR2(i__),
     +  PrimaryID(i__),fPrimary(i__),Nobs(i__),clean(i__),u_g(i__),
     +  g_r2(i__),r_i(i__),r_ha(i__),cleanu(i__),umag(i__),
     +  umagAB(i__),e_umag(i__),chiu(i__),warningu(i__),umagap(i__),
     +  umagapAB(i__),e_umagap(i__),snru(i__),umaglim(i__),
     +  PSFFWHMu(i__),MJDu(i__),detIDu(i__),cleang(i__),gmag(i__),
     +  gmagAB(i__),e_gmag(i__),chig(i__),warningg(i__),gmagap(i__),
     +  gmagapAB(i__),e_gmagap(i__),snrg(i__),gmaglim(i__),
     +  PSFFWHMg(i__),MJDg(i__),detIDg(i__),cleanr2(i__),r2mag(i__),
     +  r2magAB(i__),e_r2mag(i__),chir2(i__),warningr2(i__),
     +  r2magap(i__),r2magapAB(i__),e_r2magap(i__),snrr2(i__),
     +  r2maglim(i__),PSFFWHMr2(i__),MJDr2(i__),detIDr2(i__),
     +  cleanHa(i__),Hamag(i__),e_Hamag(i__),chiHa(i__),
     +  warningHa(i__),Hamagap(i__),e_Hamagap(i__),snrHa(i__),
     +  Hamaglim(i__),PSFFWHMHa(i__),MJDHa(i__),detIDHa(i__),
     +  cleanr(i__),rmag(i__),rmagAB(i__),e_rmag(i__),chir(i__),
     +  warningr(i__),rmagap(i__),rmagapAB(i__),e_rmagap(i__),
     +  snrr(i__),rmaglim(i__),PSFFWHMr(i__),MJDr(i__),detIDr(i__),
     +  cleani(i__),imag(i__),imagAB(i__),e_imag(i__),chii(i__),
     +  warningi(i__),imagap(i__),imagapAB(i__),e_imagap(i__),
     +  snri(i__),imaglim(i__),PSFFWHMi(i__),MJDi(i__),detIDi(i__),
     +  Field(i__),Ext(i__),nbDist(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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