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

Conversion of standardized ReadMe file for file /./ftp/cats/II/341 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. 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    ! Unique VPHAS source identifier in the format
*                                   field-extension-number (sourceID)
      real*8        RAdeg       ! (deg) Right Ascension in decimal degrees (J2000)
*                                   (RAJ2000)
      real*8        DEdeg       ! (deg) Declination in decimal degrees (J2000)
*                                   (DEJ2000)
      character*27  VPHASDR2    ! Position-based source name. Not a unique
*                                   identifier (VPHASDR2 JHHMMSS.s+DDMMSS.s)
*                                   (name)
      character*14  PrimaryID   ! SourceID of the preferred catalogue entry
*                                   for this source (primaryID)
      integer*4     fPrimary    ! [0/1] Flag =1 if the entry is the preferred
*                                   observation of a unique source
*                                   (Primary_source)
      integer*4     Nobs        ! [1-6] Number of observations of this source
*                                   in the survey (nObs)
      integer*4     clean       ! [0/1] =1 if cleang, cleanr2, cleanha, cleanr
*                                    and cleani are true (clean)
      real*4        u_g         ! (mag) ? (u-g) colour index, formed by subtracting
*                                     columns u and g (u_g)
      real*4        g_r2        ! (mag) ? (g-r2) colour index, formed by subtracting
*                                     columns g and r2 (g_r2)
      real*4        r_i         ! (mag) ? (r-i) colour index, formed by subtracting
*                                     columns r and i (r_i)
      real*4        r_ha        ! (mag) ? (r-Halpha) colour index, formed by
*                                    subtracting columns r and ha (r_ha)
      integer*4     cleanu      ! [0/1] =1 if the u-band detection is
*                                    significant and the PSF fit good (clean_u)
      real*4        umag        ! (mag) ? Default u-band VEGA magnitude obtained using
*                                   PSF fitting (u)
      real*4        umagAB      ! (mag) ? Default u-band AB magnitude obtained using
*                                   PSF fitting (u_AB)
      real*4        e_umag      ! (mag) ? Statistical uncertainty for u,
*                                   not including systematics (err_u)
      real*4        chiu        ! ? Goodness of the PSF fit in u,
*                                    as determined by DAOPHOT ALLSTAR (chi_u)
      character*12  warningu    ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_u) (1)
      real*4        umagap      ! (mag) ? u-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_u)
      real*4        umagapAB    ! (mag) ? u-band AB magnitude obtained using
*                                   aperture photometry (aperMag_u_AB)
      real*4        e_umagap    ! (mag) ? Statistical uncertainty for umagap,
*                                   not including systematics (aperMagErr_u)
      real*4        snru        ! ? Ratio of the signal in the aperture relative
*                                    to the background (snr_u)
      real*4        umaglim     ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_u)
      real*4        PSFFWHMu    ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_u)
      real*8        MJDu        ! (d) ? Modified Julian Date at the start of the
*                                   u-band exposure (mjd_u)
      character*23  detIDu      ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_u)
      integer*4     cleang      ! [0/1] =1 if the g-band detection is
*                                   significant and the PSF fit good (clean_g)
      real*4        gmag        ! (mag) ? Default g-band VEGA magnitude obtained using
*                                   PSF fitting (g)
      real*4        gmagAB      ! (mag) ? Default g-band AB magnitude obtained using
*                                   PSF fitting (g_AB)
      real*4        e_gmag      ! (mag) ? Statistical uncertainty for g,
*                                   not including systematics (err_g)
      real*4        chig        ! ? Goodness of the PSF fit in g,
*                                   as determined by DAOPHOT ALLSTAR (chi_g)
      character*12  warningg    ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_g) (1)
      real*4        gmagap      ! (mag) ? g-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_g)
      real*4        gmagapAB    ! (mag) ? g-band AB magnitude obtained using
*                                   aperture photometry (aperMag_g_AB)
      real*4        e_gmagap    ! (mag) ? Statistical uncertainty for gmagap,
*                                   not including systematics (aperMagErr_g)
      real*4        snrg        ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_g)
      real*4        gmaglim     ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_g)
      real*4        PSFFWHMg    ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_g)
      real*8        MJDg        ! (d) ? Modified Julian Date at the start of the
*                                   g-band exposure (mjd_g)
      character*23  detIDg      ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_g)
      integer*4     cleanr2     ! [0/1] =1 if the r2-band detection is
*                                   significant and the PSF fit good (clean_r2)
      real*4        r2mag       ! (mag) ? Second r-band VEGA magnitude obtained using
*                                   PSF fitting (r2)
      real*4        r2magAB     ! (mag) ? Second r-band AB magnitude obtained using
*                                   PSF fitting (r2_AB)
      real*4        e_r2mag     ! (mag) ? Statistical uncertainty for r2,
*                                   not including systematics (err_r2)
      real*4        chir2       ! ? Goodness of the PSF fit in r2,
*                                   as determined by DAOPHOT ALLSTAR (chi_r2)
      character*12  warningr2   ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_r2) (1)
      real*4        r2magap     ! (mag) ? r2-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_r2)
      real*4        r2magapAB   ! (mag) ? r2-band AB magnitude obtained using
*                                   aperture photometry (aperMag_r2_AB)
      real*4        e_r2magap   ! (mag) ? Statistical uncertainty for r2magap,
*                                   not including systematics (aperMagErr_r2)
      real*4        snrr2       ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_r2)
      real*4        r2maglim    ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_r2)
      real*4        PSFFWHMr2   ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_r2)
      real*8        MJDr2       ! (d) ? Modified Julian Date at the start of the
*                                   r2-band exposure (mjd_r2)
      character*23  detIDr2     ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_r2)
      integer*4     cleanHa     ! [0/1] =1 if the H-alpha detection is
*                                   significant and the PSF fit good (clean_ha)
      real*4        Hamag       ! (mag) ? Default H-alpha VEGA magnitude obtained
*                                   using PSF fitting (ha)
      real*4        e_Hamag     ! (mag) ? Statistical uncertainty for ha,
*                                   not including systematics (err_ha)
      real*4        chiHa       ! ? Goodness of the PSF fit in H-alpha,
*                                   as determined by DAOPHOT (chi_ha)
      character*9   warningHa   ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_ha) (1)
      real*4        Hamagap     ! (mag) ? H-alpha VEGA magnitude obtained using
*                                   aperture photometry (aperMag_ha)
      real*4        e_Hamagap   ! (mag) ? Statistical uncertainty for Hamagap,
*                                   not including systematics (aperMagErr_ha)
      real*4        snrHa       ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_ha)
      real*4        Hamaglim    ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_ha)
      real*4        PSFFWHMHa   ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_ha)
      real*8        MJDHa       ! (d) ? Modified Julian Date at the start of the
*                                   H-alpha exposure (mjd_ha)
      character*23  detIDHa     ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_ha)
      integer*4     cleanr      ! [0/1] =1 if the r-band detection is
*                                   significant and the PSF fit good (clean_r)
      real*4        rmag        ! (mag) ? Default r-band VEGA magnitude obtained using
*                                   PSF fitting (r)
      real*4        rmagAB      ! (mag) ? Default r-band AB magnitude obtained using
*                                   PSF fitting (r_AB)
      real*4        e_rmag      ! (mag) ? Statistical uncertainty for r,
*                                   not including systematics (err_r)
      real*4        chir        ! ? Goodness of the PSF fit in r,
*                                   as determined by DAOPHOT ALLSTAR (chi_r)
      character*9   warningr    ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_r) (1)
      real*4        rmagap      ! (mag) ? r-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_r)
      real*4        rmagapAB    ! (mag) ? r-band AB magnitude obtained using
*                                   aperture photometry (aperMag_r_AB)
      real*4        e_rmagap    ! (mag) ? Statistical uncertainty for rmagap,
*                                   not including systematics (aperMagErr_r)
      real*4        snrr        ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_r)
      real*4        rmaglim     ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_r)
      real*4        PSFFWHMr    ! (arcsec) ? Typical PSF FWHM in the CCD frame
*                                    (psffwhm_r)
      real*8        MJDr        ! (d) ? Modified Julian Date at the start of the
*                                   r-band exposure (mjd_r)
      character*23  detIDr      ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_r)
      integer*4     cleani      ! [0/1] =1 if the i-band detection is
*                                   significant and the PSF fit good (clean_i)
      real*4        imag        ! (mag) ? Default i-band VEGA magnitude obtained using
*                                   PSF fitting (i)
      real*4        imagAB      ! (mag) ? Default i-band AB magnitude obtained using
*                                   PSF fitting (i_AB)
      real*4        e_imag      ! (mag) ? Statistical uncertainty for i,
*                                   not including systematics (err_i)
      real*4        chii        ! ? Goodness of the PSF fit in i,
*                                   as determined by DAOPHOT ALLSTAR (chi_i)
      character*9   warningi    ! Flags problems encountered while fitting the
*                                   PSF, if any (warning_i) (1)
      real*4        imagap      ! (mag) ? i-band VEGA magnitude obtained using
*                                   aperture photometry (aperMag_i)
      real*4        imagapAB    ! (mag) ? i-band AB magnitude obtained using
*                                   aperture photometry (aperMag_i_AB)
      real*4        e_imagap    ! (mag) ? Statistical uncertainty for imagap,
*                                   not including systematics (aperMagErr_i)
      real*4        snri        ! ? Ratio of the signal in the aperture relative
*                                   to the background (snr_i)
      real*4        imaglim     ! (mag) ? The Vega magnitude a SNR=3 source would have
*                                   at this position (magLim_i)
      real*4        PSFFWHMi    ! (arcsec) Typical PSF FWHM in the CCD frame (psffwhm_i)
      real*8        MJDi        ! (d) Modified Julian Date at the start of the
*                                   i-band exposure (mjd_i)
      character*23  detIDi      ! Detection identifier:
*                                   night-#exposure-#extension-#number
*                                   (detectionID_i)
      character*5   Field       ! Survey-specific identifier of the telescope
*                                   pointing (field)
      integer*4     Ext         ! [1/32] OmegaCAM CCD extension used (ext)
      real*4        nbDist      ! (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,RAdeg,DEdeg,VPHASDR2,PrimaryID,fPrimary,Nobs,clean,
     +  u_g,g_r2,r_i,r_ha,cleanu,umag,umagAB,e_umag,chiu,warningu,
     +  umagap,umagapAB,e_umagap,snru,umaglim,PSFFWHMu,MJDu,detIDu,
     +  cleang,gmag,gmagAB,e_gmag,chig,warningg,gmagap,gmagapAB,
     +  e_gmagap,snrg,gmaglim,PSFFWHMg,MJDg,detIDg,cleanr2,r2mag,
     +  r2magAB,e_r2mag,chir2,warningr2,r2magap,r2magapAB,e_r2magap,
     +  snrr2,r2maglim,PSFFWHMr2,MJDr2,detIDr2,cleanHa,Hamag,e_Hamag,
     +  chiHa,warningHa,Hamagap,e_Hamagap,snrHa,Hamaglim,PSFFWHMHa,
     +  MJDHa,detIDHa,cleanr,rmag,rmagAB,e_rmag,chir,warningr,rmagap,
     +  rmagapAB,e_rmagap,snrr,rmaglim,PSFFWHMr,MJDr,detIDr,cleani,
     +  imag,imagAB,e_imag,chii,warningi,imagap,imagapAB,e_imagap,
     +  snri,imaglim,PSFFWHMi,MJDi,detIDi,Field,Ext,nbDist
        if(ar__(85:89) .EQ. '') u_g = rNULL__
        if(ar__(91:95) .EQ. '') g_r2 = rNULL__
        if(ar__(97:101) .EQ. '') r_i = rNULL__
        if(ar__(103:107) .EQ. '') r_ha = rNULL__
        if(ar__(111:115) .EQ. '') umag = rNULL__
        if(ar__(117:121) .EQ. '') umagAB = rNULL__
        if(ar__(123:126) .EQ. '') e_umag = rNULL__
        if(ar__(128:132) .EQ. '') chiu = rNULL__
        if(ar__(147:151) .EQ. '') umagap = rNULL__
        if(ar__(153:157) .EQ. '') umagapAB = rNULL__
        if(ar__(159:162) .EQ. '') e_umagap = rNULL__
        if(ar__(164:169) .EQ. '') snru = rNULL__
        if(ar__(171:174) .EQ. '') umaglim = rNULL__
        if(ar__(176:179) .EQ. '') PSFFWHMu = rNULL__
        if(ar__(181:188) .EQ. '') MJDu = rNULL__
        if(ar__(216:220) .EQ. '') gmag = rNULL__
        if(ar__(222:226) .EQ. '') gmagAB = rNULL__
        if(ar__(228:232) .EQ. '') e_gmag = rNULL__
        if(ar__(234:237) .EQ. '') chig = rNULL__
        if(ar__(252:256) .EQ. '') gmagap = rNULL__
        if(ar__(258:262) .EQ. '') gmagapAB = rNULL__
        if(ar__(264:267) .EQ. '') e_gmagap = rNULL__
        if(ar__(269:274) .EQ. '') snrg = rNULL__
        if(ar__(276:279) .EQ. '') gmaglim = rNULL__
        if(ar__(281:284) .EQ. '') PSFFWHMg = rNULL__
        if(ar__(286:293) .EQ. '') MJDg = rNULL__
        if(ar__(321:325) .EQ. '') r2mag = rNULL__
        if(ar__(327:331) .EQ. '') r2magAB = rNULL__
        if(ar__(333:336) .EQ. '') e_r2mag = rNULL__
        if(ar__(338:341) .EQ. '') chir2 = rNULL__
        if(ar__(356:360) .EQ. '') r2magap = rNULL__
        if(ar__(362:366) .EQ. '') r2magapAB = rNULL__
        if(ar__(368:371) .EQ. '') e_r2magap = rNULL__
        if(ar__(373:378) .EQ. '') snrr2 = rNULL__
        if(ar__(380:383) .EQ. '') r2maglim = rNULL__
        if(ar__(385:388) .EQ. '') PSFFWHMr2 = rNULL__
        if(ar__(390:397) .EQ. '') MJDr2 = rNULL__
        if(ar__(425:429) .EQ. '') Hamag = rNULL__
        if(ar__(431:434) .EQ. '') e_Hamag = rNULL__
        if(ar__(436:439) .EQ. '') chiHa = rNULL__
        if(ar__(451:455) .EQ. '') Hamagap = rNULL__
        if(ar__(457:460) .EQ. '') e_Hamagap = rNULL__
        if(ar__(462:467) .EQ. '') snrHa = rNULL__
        if(ar__(469:472) .EQ. '') Hamaglim = rNULL__
        if(ar__(474:477) .EQ. '') PSFFWHMHa = rNULL__
        if(ar__(479:486) .EQ. '') MJDHa = rNULL__
        if(ar__(514:518) .EQ. '') rmag = rNULL__
        if(ar__(520:524) .EQ. '') rmagAB = rNULL__
        if(ar__(526:529) .EQ. '') e_rmag = rNULL__
        if(ar__(531:534) .EQ. '') chir = rNULL__
        if(ar__(546:550) .EQ. '') rmagap = rNULL__
        if(ar__(552:556) .EQ. '') rmagapAB = rNULL__
        if(ar__(558:561) .EQ. '') e_rmagap = rNULL__
        if(ar__(563:568) .EQ. '') snrr = rNULL__
        if(ar__(570:573) .EQ. '') rmaglim = rNULL__
        if(ar__(575:578) .EQ. '') PSFFWHMr = rNULL__
        if(ar__(580:587) .EQ. '') MJDr = rNULL__
        if(ar__(615:619) .EQ. '') imag = rNULL__
        if(ar__(621:625) .EQ. '') imagAB = rNULL__
        if(ar__(627:630) .EQ. '') e_imag = rNULL__
        if(ar__(632:635) .EQ. '') chii = rNULL__
        if(ar__(647:651) .EQ. '') imagap = rNULL__
        if(ar__(653:657) .EQ. '') imagapAB = rNULL__
        if(ar__(659:662) .EQ. '') e_imagap = rNULL__
        if(ar__(664:669) .EQ. '') snri = rNULL__
        if(ar__(671:674) .EQ. '') imaglim = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  sourceID,RAdeg,DEdeg,VPHASDR2,PrimaryID,fPrimary,Nobs,clean,
     +  u_g,g_r2,r_i,r_ha,cleanu,umag,umagAB,e_umag,chiu,warningu,
     +  umagap,umagapAB,e_umagap,snru,umaglim,PSFFWHMu,MJDu,detIDu,
     +  cleang,gmag,gmagAB,e_gmag,chig,warningg,gmagap,gmagapAB,
     +  e_gmagap,snrg,gmaglim,PSFFWHMg,MJDg,detIDg,cleanr2,r2mag,
     +  r2magAB,e_r2mag,chir2,warningr2,r2magap,r2magapAB,e_r2magap,
     +  snrr2,r2maglim,PSFFWHMr2,MJDr2,detIDr2,cleanHa,Hamag,e_Hamag,
     +  chiHa,warningHa,Hamagap,e_Hamagap,snrHa,Hamaglim,PSFFWHMHa,
     +  MJDHa,detIDHa,cleanr,rmag,rmagAB,e_rmag,chir,warningr,rmagap,
     +  rmagapAB,e_rmagap,snrr,rmaglim,PSFFWHMr,MJDr,detIDr,cleani,
     +  imag,imagAB,e_imag,chii,warningi,imagap,imagapAB,e_imagap,
     +  snri,imaglim,PSFFWHMi,MJDi,detIDi,Field,Ext,nbDist
c    .......End.of.Just test output...........
      end do
      close(1)

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