FORTRAN Generation
(/./ftp/cats/J/ApJS/265/40)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/265/40 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 2026-Aug-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. J/ApJS/265/40    MUDF III. HST grism spectroscopy & imaging    (Revalski+, 2023)
*================================================================================
*The MUSE Ultra Deep Field (MUDF).
*III. Hubble Space Telescope WFC3 grism spectroscopy and imaging.
*    Revalski M., Rafelski M., Fumagalli M., Fossati M., Pirzkal N.,
*    Sunnquist B., Prichard L.J., Henry A., Bagley M., Dutta R., Papini G.,
*    Battaia F.A., D'Odorico V., Dayal P., Estrada-Carpenter V.,
*    Lofthouse E.K., Lusso E., Morris S.L., Nedkova K.V., Papovich C.,
*    Peroux C.
*   <Astrophys. J. Suppl. Ser., 265, 40 (2023)>
*   =2023ApJS..265...40R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! *Multi-filter source Extractor catalog

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     Seq        (nr__) ! [1/21694] Object ID number (NUMBER)
      real*8        Xpos       (nr__) ! (pix) [1463/6723] Source position along X (X_IMAGE)
      real*8        Ypos       (nr__) ! (pix) [1137/6694] Source position along Y (Y_IMAGE)
      integer*4     b_Xpos     (nr__) ! (pix) [1462/6721] Minimum x-coordinate among
*                                    detected pixels (XMIN_IMAGE)
      integer*4     b_Ypos     (nr__) ! (pix) [1131/6690] Minimum y-coordinate among
*                                    detected pixels (YMIN_IMAGE)
      integer*4     B_Xpos_1   (nr__) ! (pix) [1466/6725] Maximum x-coordinate among
*                                    detected pixels (XMAX_IMAGE)
      integer*4     B_Ypos_1   (nr__) ! (pix) [1144/6697] Maximum y-coordinate among
*                                    detected pixels (YMAX_IMAGE)
      real*8        RAdeg      (nr__) ! (deg) [325.54/325.67] Right ascension of
*                                    barycenter (J2000) (ALPHA_J2000)
      real*8        DEdeg      (nr__) ! (deg) [-44.4/-44.27] Declination of barycenter
*                                    (J2000) (DELTA_J2000)
      real*4        Aimg       (nr__) ! (pix) [0.8/29.6] Profile RMS along major axis
*                                    (A_IMAGE)
      real*4        Bimg       (nr__) ! (pix) [0.4/18.8] Profile RMS along minor axis
*                                    (B_IMAGE)
      real*4        PA         (nr__) ! (deg) [-89.8/90] Position angle (THETA_IMAGE)
      real*4        Elong      (nr__) ! [1/5.2] Elongation parameter (Aimg/Bimg)
*                                    (ELONGATION)
      real*4        Ell        (nr__) ! [0.004/0.81] Ellipticity (ELLIPTICITY)
      real*4        RKF140W    (nr__) ! [0/11.71] HST WFC3/IR F140W Kron radius in
*                                    units of A or B (KRON_RADIUS_F140W)
      real*8        RF140W     (nr__) ! (pix) [-74.9/33.5] HST WFC3/IR F140W effective
*                                    radius (FLUX_RADIUS_F140W)
      real*4        FWHMF140W  (nr__) ! [2.58/92.3] HST WFC3/IR F140W full width at
*                                    half maximum (FWHM_IMAGE_F140W)
      real*4        S_GF140W   (nr__) ! [0/1] HST WFC3/IR F140W S/G classifier
*                                    (0=galaxy; 1=star) (CLASS_STAR_F140W)
      integer*4     AreaF140W  (nr__) ! (arcmin2) [6/23405] HST WFC3/IR F140W isophotal area
*                                    (ISOAREA_IMAGE_F140W)
      real*8        F140Wiso   (nr__) ! (mag) [15/29.1] HST WFC3/IR F140W isophotal
*                                    magnitude (MAG_ISO_F140W)
      real*4        e_F140Wiso (nr__) ! (mag) [0.0001/0.7] F140Wiso uncertainty
*                                    (MAGERR_ISO_F140W)
      real*8        FF140Wi    (nr__) ! (ct) [0.09/34719] HST WFC3/IR F140W isophotal
*                                    flux (FLUX_ISO_F140W)
      real*8        e_FF140Wi  (nr__) ! (ct) [0.03/4.4] FF140Wi uncertainty
*                                    (FLUXERR_ISO_F140W)
      real*8        F140Wmag   (nr__) ! (mag) [15/29.8]?=99 Kron-like elliptical aperture
*                                    HST WFC3/IR F140W magnitude (MAG_AUTO_F140W)
      real*8        e_F140Wmag (nr__) ! (mag) [0.0001/2.9e+10] F140Wmag uncertainty
*                                    (MAGERR_AUTO_F140W)
      real*8        FF140W     (nr__) ! (ct) [-2.3/33417] HST WFC3/IR F140W flux within a
*                                    Kron-like elliptical aperture
*                                    (FLUX_AUTO_F140W)
      real*8        e_FF140W   (nr__) ! (ct) [0.07/6.1e+11] FF140W uncertainty
*                                    (FLUXERR_AUTO_F140W)
      real*4        BkgF140W   (nr__) ! (ct) [-0.004/0.031] HST WFC3/IR F140W background
*                                    at centroid position (BACKGROUND_F140W)
      integer*4     f_F140W    (nr__) ! [0/19] Sextractor F140W flags (FLAGS_F140W)
      real*4        RKF125W    (nr__) ! [0/9.1] HST WFC3/IR F125W Kron radius in
*                                    units of A or B (KRON_RADIUS_F125W)
      real*8        RF125W     (nr__) ! (pix) [-255.2/548.5] HST WFC3/IR F125W effective
*                                    radius (FLUX_RADIUS_F125W)
      real*4        FWHMF125W  (nr__) ! [0/92.3] HST WFC3/IR F125W full width at
*                                    half maximum (FWHM_IMAGE_F125W)
      real*4        S_GF125W   (nr__) ! [0/1] HST WFC3/IR F125W S/G classifier
*                                    (0=galaxy; 1=star) (CLASS_STAR_F125W)
      integer*4     AreaF125W  (nr__) ! (arcmin2) [0/23377] HST WFC3/IR F125W isophotal area
*                                    (ISOAREA_IMAGE_F125W)
      real*8        F125Wiso   (nr__) ! (mag) [15.5/46.8]?=99 HST WFC3/IR F125W isophotal
*                                    magnitude (MAG_ISO_F125W)
      real*4        e_F125Wiso (nr__) ! (mag) [0.0002/1.6e+19]?=99 F125Wiso uncertainty
*                                    (MAGERR_ISO_F125W)
      real*4        FF125Wi    (nr__) ! (ct) [-0.3/18653] HST WFC3/IR F125W isophotal
*                                    flux (FLUX_ISO_F125W)
      real*8        e_FF125Wi  (nr__) ! (ct) [0.06/2.3e+12] FF125Wi uncertainty
*                                    (FLUXERR_ISO_F125W)
      real*8        F125Wmag   (nr__) ! (mag) [15.6/44.6]?=99 Kron-like elliptical
*                                    aperture HST WFC3/IR F125W magnitude
*                                    (MAG_AUTO_F125W)
      real*4        e_F125Wmag (nr__) ! (mag) [0.0002/8.1e+18]?=99 F125Wmag uncertainty
*                                    (MAGERR_AUTO_F250W)
      real*4        FF125W     (nr__) ! (ct) [-1.43/17833] HST WFC3/IR F125W flux within
*                                    a Kron-like elliptical aperture
*                                    (FLUX_AUTO_F125W)
      real*8        e_FF125W   (nr__) ! (ct) [0.14/2.5e+12] FF125W uncertainty
*                                    (FLUXERR_AUTO_F125W)
      real*4        BkgF125W   (nr__) ! (ct) [-0.006/0.03] HST WFC3/IR F125W background
*                                    at centroid position (BACKGROUND_F125W)
      integer*4     f_F125W    (nr__) ! [0/19] Sextractor HST WFC3/IR F125W flags
*                                    (FLAGS_F125W)
      real*4        RKF702W    (nr__) ! [0/9.1] HST WFPC2 F702W Kron radius in units
*                                    of A or B (KRON_RADIUS_F702W)
      real*8        RF702W     (nr__) ! (pix) [-8828/4364] HST WFPC2 F702W effective
*                                    radius (FLUX_RADIUS_F702W)
      real*4        FWHMF702W  (nr__) ! [0/50.4] HST WFPC2 F702W full width at half
*                                    maximum (FWHM_IMAGE_F702W)
      real*4        S_GF702W   (nr__) ! [0/1] HST WFPC2 F702W S/G classifier
*                                    (0=galaxy; 1=star) (CLASS_STAR_F702W)
      integer*4     AreaF702W  (nr__) ! (arcmin2) [0/10342] HST WFPC2 F702W isophotal area
*                                    (ISOAREA_IMAGE_F702W)
      real*8        F702Wiso   (nr__) ! (mag) [18.2/51.6]?=99 HST WFPC2 F702W isophotal
*                                    magnitude (MAG_ISO_F702W)
      real*8        e_F702Wiso (nr__) ! (mag) [0.002/5.6e+11]?=99 F702Wiso uncertainty
*                                    (MAGERR_ISO_F702W)
      real*4        FF702Wi    (nr__) ! (ct) [-0.021/64] HST WFPC2 F702W isophotal flux
*                                    (FLUX_ISO_F702W)
      real*4        e_FF702Wi  (nr__) ! (ct) [0/2.4e+11] FF702Wi uncertainty
*                                    (FLUXERR_ISO_F702W)
      real*8        F702Wmag   (nr__) ! (mag) [18.2/50.3]?=99 Kron-like elliptical
*                                    aperture HST WFPC2 F702W magnitude
*                                    (MAG_AUTO_F702W)
      real*8        e_F702Wmag (nr__) ! (mag) [0.0015/9.1e+11]?=99 F702Wmag uncertainty
*                                    (MAGERR_AUTO_F702W)
      real*4        FF702W     (nr__) ! (ct) [-0.13/62.5] HST WFPC2 F702W flux within a
*                                    Kron-like elliptical aperture
*                                    (FLUX_AUTO_F702W)
      real*4        e_FF702W   (nr__) ! (ct) [0/2.1e+11] FF702W uncertainty
*                                    (FLUXERR_AUTO_F702W)
      real*4        BkgF702W   (nr__) ! (ct) [-3.8e-5/0.0009] HST WFPC2 F702W background
*                                    at centroid position (BACKGROUND_F702W)
      integer*4     f_F702W    (nr__) ! [0/19] Sextractor HST WFPC2 F702W flags
*                                    (FLAGS_F702W)
      real*4        RKF450W    (nr__) ! [0/9.1] HST WFPC2 F450W Kron radius in units
*                                    of A or B (KRON_RADIUS_F450W)
      real*8        RF450W     (nr__) ! (pix) [-1155/22403] HST WFPC2 F450W effective
*                                    radius (FLUX_RADIUS_F450W)
      real*4        FWHMF450W  (nr__) ! [0/52.4] HST WFPC2 F450W full width at half
*                                    maximum (FWHM_IMAGE_F450W)
      real*4        S_GF450W   (nr__) ! [0/1] HST WFPC2 F450W S/G classifier
*                                    (0=galaxy; 1=star) (CLASS_STAR_F450W)
      integer*4     AreaF450W  (nr__) ! (arcmin2) [0/9435] HST WFPC2 F450W isophotal area
*                                    (ISOAREA_IMAGE_F450W)
      real*8        F450Wiso   (nr__) ! (mag) [18.5/49.4]?=99 HST WFPC2 F450W isophotal
*                                    magnitude (MAG_ISO_F450W)
      real*8        e_F450Wiso (nr__) ! (mag) [0.003/7e+11]?=99 F450Wiso uncertainty
*                                    (MAGERR_ISO_F450W)
      real*4        FF450Wi    (nr__) ! (ct) [-0.8/22.4] HST WFPC2 F450W isophotal flux
*                                    (FLUX_ISO_F450W)
      real*4        e_FF450Wi  (nr__) ! (ct) [0/1.6e+11] FF450Wi uncertainty
*                                    (FLUXERR_ISO_F450W)
      real*8        F450Wmag   (nr__) ! (mag) [18.3/48.4]?=99 Kron-like elliptical
*                                    aperture HST WFPC2 F450W magnitude
*                                    (MAG_AUTO_F450W)
      real*8        e_F450Wmag (nr__) ! (mag) [0.002/1.2e+12]?=99 F450Wmag uncertainty
*                                    (MAGERR_AUTO_F450W)
      real*4        FF450W     (nr__) ! (ct) [-0.8/26.7] HST WFPC2 F450W flux within a
*                                    Kron-like elliptical aperture
*                                    (FLUX_AUTO_F450W)
      real*4        e_FF450W   (nr__) ! (ct) [0/1.3e+11] FF450Wmag uncertainty
*                                    (FLUXERR_AUTO_F450W)
      real*4        BkgF450W   (nr__) ! (ct) [-4.2e-5/0.0004] HST WFPC2 F450W background
*                                    at centroid position (BACKGROUND_F450W)
      integer*4     f_F450W    (nr__) ! [0/19] Sextractor HST WFPC2 F450W flags
*                                    (FLAGS_F450W)
      real*4        RKF336W    (nr__) ! [0/9.1] HST WFC3/UVIS F336W Kron radius in
*                                    units of A or B (KRON_RADIUS_F336W)
      real*8        RF336W     (nr__) ! (pix) [-3506/22075] HST WFC3/UVIS F336W effective
*                                    radius (FLUX_RADIUS_F336W)
      real*4        FWHMF336W  (nr__) ! [0/128.2] HST WFC3/UVIS F336W full width at
*                                    half maximum (FWHM_IMAGE_F336W)
      real*4        S_GF336W   (nr__) ! [0/1] HST WFC3/UVIS F336W S/G classifier
*                                    (0=galaxy; 1=star) (CLASS_STAR_F336W)
      integer*4     AreaF336W  (nr__) ! (arcmin2) [0/9160] HST WFC3/UVIS F336W isophotal area
*                                    (ISOAREA_IMAGE_F336W)
      real*8        F336Wiso   (nr__) ! (mag) [19.97/50.1]?=99 HST WFC3/UVIS F336W
*                                    isophotal magnitude (MAG_ISO_F336W)
      real*8        e_F336Wiso (nr__) ! (mag) [0.0009/1311]?=99 F336Wiso uncertainty
*                                    (MAGERR_ISO_F336W)
      real*4        FF336Wi    (nr__) ! (ct) [-0.11/76.6] HST WFC3/UVIS F336W isophotal
*                                    flux (FLUX_ISO_F336W)
      real*4        e_FF336Wi  (nr__) ! (ct) [0/0.09] FF336Wi uncertainty
*                                    (FLUXERR_ISO_F336W)
      real*8        F336Wmag   (nr__) ! (mag) [20/48.5]?=99 Kron-like elliptical aperture
*                                    HST WFC3/UVIS F336W magnitude
*                                    (MAG_AUTO_F336W)
      real*8        e_F336Wmag (nr__) ! (mag) [0.0009/428]?=99 F336Wmag uncertainty
*                                    (MAGERR_AUTO_F336W)
      real*4        FF336W     (nr__) ! (ct) [-0.16/74.3] HST WFC3/UVIS F336W flux within
*                                    a Kron-like elliptical aperture
*                                    (FLUX_AUTO_F336W)
      real*4        e_FF336W   (nr__) ! (ct) [0/0.1] FF336Wmag uncertainty
*                                    (FLUXERR_AUTO_F336W)
      real*4        BkgF336W   (nr__) ! (ct) [-0.0004/0.0008] HST WFC3/UVIS F336W
*                                    background at centroid position
*                                    (BACKGROUND_F336W)
      integer*4     f_F336W    (nr__) ! [0/19] Sextractor HST WFC3/UVIS F336W flags
*                                    (FLAGS_F336W)
      real*8        zsp        (nr__) ! [0.06/5.4]?=-99 Derived spectroscopic
*                                    redshift (SPEC_Z)
      real*4        e_zsp      (nr__) ! [1e-5/0.11]?=-99 zsp uncertainty (SPEC_Z_ERR)
      integer*4     q_zsp      (nr__) ! [0/4]? Redshift quality flag (4=best)
*                                    (SPEC_Z_FLAG) (1)
      character*2   f_zsp      (nr__) ! sp: Indicates 1D FITS spectrum available
*                                    in subdirectory "sp"; flag added by CDS
*Note (1): Redshift quality flag as follows:
*  4 = 2+ lines with S/N>=5 (366 occurrences)
*  3 = 2+ lines with S/N>=3 (22 occurrences)
*  2 = 1 line with S/N>=5 (31 occurrences)
*  1 = spectral coverage, but no redshift (1088 occurrences)
*  0 = no spectral coverage for the source (1868 occurrences)
*
*          See Section 5.3.

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

C  Loading file 'catalog.dat'	! *Multi-filter source Extractor catalog

C  Format for file interpretation

    1 format(
     +  I5,1X,F8.3,1X,F8.3,1X,I4,1X,I4,1X,I4,1X,I4,1X,F11.7,1X,F11.7,
     +  1X,F6.3,1X,F6.3,1X,F5.1,1X,F5.3,1X,F5.3,1X,F5.2,1X,F7.3,1X,
     +  F5.2,1X,F4.2,1X,I5,1X,F9.6,1X,F6.4,1X,F14.8,1X,F10.8,1X,F7.4,
     +  1X,F16.4,1X,F16.9,1X,F21.8,1X,E13.10,1X,I2,1X,F4.2,1X,F8.3,1X,
     +  F5.2,1X,F4.2,1X,I5,1X,F7.4,1X,E22.17,1X,E13.10,1X,F22.8,1X,
     +  F7.4,1X,E22.17,1X,E13.10,1X,F21.7,1X,E13.10,1X,I2,1X,F4.2,1X,
     +  F9.3,1X,F5.2,1X,F4.2,1X,I5,1X,F7.4,1X,F18.4,1X,E13.10,1X,
     +  E14.10,1X,F7.4,1X,F18.4,1X,E13.10,1X,E14.10,1X,E13.10,1X,I2,
     +  1X,F4.2,1X,F9.3,1X,F5.2,1X,F4.2,1X,I4,1X,F7.4,1X,F18.4,1X,
     +  E13.10,1X,E14.9,1X,F7.4,1X,F18.4,1X,E13.10,1X,E14.9,1X,E13.10,
     +  1X,I2,1X,F4.2,1X,F9.3,1X,F6.2,1X,F4.2,1X,I4,1X,F7.4,1X,F9.4,
     +  1X,E13.10,1X,E12.10,1X,F7.4,1X,F8.4,1X,E13.10,1X,E12.10,1X,
     +  E13.10,1X,I2,1X,F9.5,1X,E7.5,1X,I1,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,3375
        read(1,'(A919)')ar__
        read(ar__,1)
     +  Seq(i__),Xpos(i__),Ypos(i__),b_Xpos(i__),b_Ypos(i__),
     +  B_Xpos_1(i__),B_Ypos_1(i__),RAdeg(i__),DEdeg(i__),Aimg(i__),
     +  Bimg(i__),PA(i__),Elong(i__),Ell(i__),RKF140W(i__),
     +  RF140W(i__),FWHMF140W(i__),S_GF140W(i__),AreaF140W(i__),
     +  F140Wiso(i__),e_F140Wiso(i__),FF140Wi(i__),e_FF140Wi(i__),
     +  F140Wmag(i__),e_F140Wmag(i__),FF140W(i__),e_FF140W(i__),
     +  BkgF140W(i__),f_F140W(i__),RKF125W(i__),RF125W(i__),
     +  FWHMF125W(i__),S_GF125W(i__),AreaF125W(i__),F125Wiso(i__),
     +  e_F125Wiso(i__),FF125Wi(i__),e_FF125Wi(i__),F125Wmag(i__),
     +  e_F125Wmag(i__),FF125W(i__),e_FF125W(i__),BkgF125W(i__),
     +  f_F125W(i__),RKF702W(i__),RF702W(i__),FWHMF702W(i__),
     +  S_GF702W(i__),AreaF702W(i__),F702Wiso(i__),e_F702Wiso(i__),
     +  FF702Wi(i__),e_FF702Wi(i__),F702Wmag(i__),e_F702Wmag(i__),
     +  FF702W(i__),e_FF702W(i__),BkgF702W(i__),f_F702W(i__),
     +  RKF450W(i__),RF450W(i__),FWHMF450W(i__),S_GF450W(i__),
     +  AreaF450W(i__),F450Wiso(i__),e_F450Wiso(i__),FF450Wi(i__),
     +  e_FF450Wi(i__),F450Wmag(i__),e_F450Wmag(i__),FF450W(i__),
     +  e_FF450W(i__),BkgF450W(i__),f_F450W(i__),RKF336W(i__),
     +  RF336W(i__),FWHMF336W(i__),S_GF336W(i__),AreaF336W(i__),
     +  F336Wiso(i__),e_F336Wiso(i__),FF336Wi(i__),e_FF336Wi(i__),
     +  F336Wmag(i__),e_F336Wmag(i__),FF336W(i__),e_FF336W(i__),
     +  BkgF336W(i__),f_F336W(i__),zsp(i__),e_zsp(i__),q_zsp(i__),
     +  f_zsp(i__)
        if(ar__(916:916) .EQ. '') q_zsp(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Seq(i__),Xpos(i__),Ypos(i__),b_Xpos(i__),b_Ypos(i__),
     +  B_Xpos_1(i__),B_Ypos_1(i__),RAdeg(i__),DEdeg(i__),Aimg(i__),
     +  Bimg(i__),PA(i__),Elong(i__),Ell(i__),RKF140W(i__),
     +  RF140W(i__),FWHMF140W(i__),S_GF140W(i__),AreaF140W(i__),
     +  F140Wiso(i__),e_F140Wiso(i__),FF140Wi(i__),e_FF140Wi(i__),
     +  F140Wmag(i__),e_F140Wmag(i__),FF140W(i__),e_FF140W(i__),
     +  BkgF140W(i__),f_F140W(i__),RKF125W(i__),RF125W(i__),
     +  FWHMF125W(i__),S_GF125W(i__),AreaF125W(i__),F125Wiso(i__),
     +  e_F125Wiso(i__),FF125Wi(i__),e_FF125Wi(i__),F125Wmag(i__),
     +  e_F125Wmag(i__),FF125W(i__),e_FF125W(i__),BkgF125W(i__),
     +  f_F125W(i__),RKF702W(i__),RF702W(i__),FWHMF702W(i__),
     +  S_GF702W(i__),AreaF702W(i__),F702Wiso(i__),e_F702Wiso(i__),
     +  FF702Wi(i__),e_FF702Wi(i__),F702Wmag(i__),e_F702Wmag(i__),
     +  FF702W(i__),e_FF702W(i__),BkgF702W(i__),f_F702W(i__),
     +  RKF450W(i__),RF450W(i__),FWHMF450W(i__),S_GF450W(i__),
     +  AreaF450W(i__),F450Wiso(i__),e_F450Wiso(i__),FF450Wi(i__),
     +  e_FF450Wi(i__),F450Wmag(i__),e_F450Wmag(i__),FF450W(i__),
     +  e_FF450W(i__),BkgF450W(i__),f_F450W(i__),RKF336W(i__),
     +  RF336W(i__),FWHMF336W(i__),S_GF336W(i__),AreaF336W(i__),
     +  F336Wiso(i__),e_F336Wiso(i__),FF336Wi(i__),e_FF336Wi(i__),
     +  F336Wmag(i__),e_F336Wmag(i__),FF336W(i__),e_FF336W(i__),
     +  BkgF336W(i__),f_F336W(i__),zsp(i__),e_zsp(i__),q_zsp(i__),
     +  f_zsp(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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