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

Conversion of standardized ReadMe file for file /./ftp/cats/II/267/ 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/267       The SDSS Photometric Catalog, Release 4   (Adelman-McCarthy+, 2006)
*================================================================================
*The SDSS Photometric Catalog, Release 4
*     Adelman-McCarthy J.K. et al.
*    <Astrophys. J. Suppl. Ser., 162, 38 (2006)>
*    =2006ApJS..162...38A
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'sdss4.sam'	! *Sample output of the SDSS-DR4 photometric catalog

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

C  Position composed of: RAdeg DEdeg
      integer*4     mode        ! [1,2] 1: primary, 2: secondary
*                                       3: family object (not here)
      integer*4     cl          ! Type (class) of the object (type) (1)
      character*19  SDSS        ! Recommended SDSS name, based on J2000 position
      character*1   m_SDSS      ! [*] The asterisk indicates that 2 different
*                                     SDSS objects share the same SDSS name
      character*1   R           ! [3] Release: blank for new, '3' included in DR3
      character*15  SDSS_ID     ! [0-9.] SDSS object identifier (2)
      character*18  objID       ! [0-9] SDSS unique object identifier (2)
      character*18  SpObjID     ! Pointer to the spectrum of object, or 0
      character*18  parentID    ! Pointer to parent (if object deblended)
      integer*4     rerun       ! [40,44] Re-run number (2)
      character*16  flags       ! [0-9A-F] Photo Object Attribute flags (3)
      character*8   Status      ! [0-9A-F] Hexadecimal  status (4)
      real*8        RAdeg       ! (deg) Right Ascension of the object (ICRS)
      real*8        DEdeg       ! (deg) Declination of the object (ICRS)
      real*4        e_RAdeg     ! (arcsec) Mean error on RAdeg
      real*4        e_DEdeg     ! (arcsec) Mean error on DEdeg
      real*4        ccPos       ! Correlation between e_RAdeg and e_DEdeg
      real*8        ObsDate     ! (yr) Observation date
      integer*4     Q           ! Quality of the observation:
*                                 1=bad 2=acceptable 3=good 4=missing 5=hole
      real*4        umag        ! (mag) Model magnitude in u filter (5)
      real*4        e_umag      ! (mag) Mean error on umag
      real*4        gmag        ! (mag) Model magnitude in g filter (5)
      real*4        e_gmag      ! (mag) Mean error on gmag
      real*4        rmag        ! (mag) Model magnitude in r filter (5)
      real*4        e_rmag      ! (mag) Mean error on rmag
      real*4        imag        ! (mag) Model magnitude in i filter (5)
      real*4        e_imag      ! (mag) Mean error on imag
      real*4        zmag        ! (mag) Model magnitude in z filter (5)
      real*4        e_zmag      ! (mag) Mean error on zmag
*--------------------------------------------------------------------------------
      character*16  uFlags      ! [0-9A-F] Detection flags, u band (flags_u)(3)
      integer*4     us          ! [0,1] 0=notStar, 1=Star in u band (probPSF_u)
      integer*4     uc          ! [0,9] Phototype in u band, 6=Star (type_u) (1)
      real*8        uDate       ! (yr) Date of observation in u band
      real*4        u_mag       ! (mag) ? Model magnitude in u filter (= umag) (5)
      real*4        e_u_mag     ! (mag) ? Mean error on u'mag (= e_umag)
      real*4        upmag       ! (mag) ? PSF magnitude in u band (psfMag_u)
      real*4        e_upmag     ! (mag) ? Mean error on upmag (psfMagErr_u)
      real*4        uPmag_1     ! (mag) ? Petrosian magnitude in u band (petroMag_u)
      real*4        e_uPmag_1   ! (mag) ? Mean error on  uPmag (petroMagErr_u)
      real*8        uPrad       ! (arcsec) ? Petrosian radius in u band (petroRad_u)
      real*8        e_uPrad     ! (arcsec) ? Mean error in uPrad (petroRadErr_u)
      real*4        uoRA        ! (arcsec) ? Offset of u position along RA (offsetRa_u)
      real*4        uoDE        ! (arcsec) ? Offset of u position along Dec (offsetDec_u)
      real*8        uisoA       ! (arcsec) ? Isophotal major diameter in u band (isoA_u)
      real*8        uisoB       ! (arcsec) ? Isophotal minor diameter in u band (isoB_u)
      real*4        uPA         ! (deg) ? Position angle of uisoA in u band (isoPhi_u)
*--------------------------------------------------------------------------------
      character*16  gFlags      ! [0-9A-F] Detection flags, g band (flags_g)(3)
      integer*4     gs          ! [0,1] 0=notStar, 1=Star in g band (probPSF_g)
      integer*4     gc          ! [0,9] Phototype in g band, 6=Star (type_g) (1)
      real*8        gDate       ! (yr) Date of observation in g band
      real*4        g_mag       ! (mag) ? Model magnitude in g filter (= gmag) (5)
      real*4        e_g_mag     ! (mag) ? Mean error on g'mag (= e_gmag)
      real*4        gpmag       ! (mag) ? PSF magnitude in g band (psfMag_g)
      real*4        e_gpmag     ! (mag) ? Mean error on gpmag (psfMagErr_g)
      real*4        gPmag_1     ! (mag) ? Petrosian magnitude in g band (petroMag_g)
      real*4        e_gPmag_1   ! (mag) ? Mean error on  gPmag (petroMagErr_g)
      real*8        gPrad       ! (arcsec) ? Petrosian radius in g band (petroMagErr_g)
      real*8        e_gPrad     ! (arcsec) ? Mean error in gPrad (petroMagErr_g)
      real*4        goRA        ! (arcsec) ? Offset of g position along RA (offsetRa_g)
      real*4        goDE        ! (arcsec) ? Offset of g position along Dec (offsetDec_g)
      real*8        gisoA       ! (arcsec) ? Isophotal major diameter in g band (isoA_g)
      real*8        gisoB       ! (arcsec) ? Isophotal minor diameter in g band (isoB_g)
      real*4        gPA         ! (deg) ? Position angle of gisoA in g band (isoPhi_g)
*--------------------------------------------------------------------------------
      character*16  rFlags      ! [0-9A-F]  Detection flags, r band (flags_r)(3)
      integer*4     rs          ! [0,1] 0=notStar, 1=Star in r band (probPSF_r)
      integer*4     rc          ! [0,9] Phototype in r band, 6=Star (type_r) (1)
      real*8        rDate       ! (yr) Date of observation in r band
      real*4        r_mag       ! (mag) ? Model magnitude in r filter (= rmag) (5)
      real*4        e_r_mag     ! (mag) ? Mean error on r'mag (= e_rmag)
      real*4        rpmag       ! (mag) ? PSF magnitude in r band (psfMag_r)
      real*4        e_rpmag     ! (mag) ? Mean error on rpmag (psfMagErr_r)
      real*4        rPmag_1     ! (mag) ? Petrosian magnitude in r band (petroMag_r)
      real*4        e_rPmag_1   ! (mag) ? Mean error on  rPmag (petroMagErr_r)
      real*8        rPrad       ! (arcsec) ? Petrosian radius in r band (petroMagErr_r)
      real*8        e_rPrad     ! (arcsec) ? Mean error in rPrad (petroMagErr_r)
      real*4        roRA        ! (arcsec) ? Offset of r position along RA (offsetRa_r)
      real*4        roDE        ! (arcsec) ? Offset of r position along Dec (offsetDec_r)
      real*8        risoA       ! (arcsec) ? Isophotal major diameter in r band (isoA_r)
      real*8        risoB       ! (arcsec) ? Isophotal minor diameter in r band (isoB_r)
      real*4        rPA         ! (deg) ? Position angle of risoA in g band (isoPhi_r)
*--------------------------------------------------------------------------------
      character*16  iFlags      ! [0-9A-F]  Detection flags, i band (flags_i)(3)
      integer*4     is          ! [0,1] 0=notStar, 1=Star in i band (probPSF_i)
      integer*4     ic          ! [0,9] Phototype in i band, 6=Star (type_i) (1)
      real*8        iDate       ! (yr) Date of observation in i band
      real*4        i_mag       ! (mag) ? Model magnitude in i filter (= imag) (5)
      real*4        e_i_mag     ! (mag) ? Mean error on i'mag (= e_imag)
      real*4        ipmag       ! (mag) ? PSF magnitude in i band (psfMag_i)
      real*4        e_ipmag     ! (mag) ? Mean error on ipmag (psfMagErr_i)
      real*4        iPmag_1     ! (mag) ? Petrosian magnitude in i band (petroMag_i)
      real*4        e_iPmag_1   ! (mag) ? Mean error on  iPmag (petroMagErr_i)
      real*8        iPrad       ! (arcsec) ? Petrosian radius in i band (petroMagErr_i)
      real*8        e_iPrad     ! (arcsec) ? Mean error in iPrad (petroMagErr_i)
      real*4        ioRA        ! (arcsec) ? Offset of i position along RA (offsetRa_i)
      real*4        ioDE        ! (arcsec) ? Offset of i position along Dec (offsetDec_i)
      real*8        iisoA       ! (arcsec) ? Isophotal major diameter in i band (isoA_i)
      real*8        iisoB       ! (arcsec) ? Isophotal minor diameter in i band (isoB_i)
      real*4        iPA         ! (deg) ? Position angle of iisoA in g band (isoPhi_i)
*--------------------------------------------------------------------------------
      character*16  zFlags      ! [0-9A-F]  Detection flags, z band (flags_z)(3)
      integer*4     zs          ! [0,1] 0=notStar, 1=Star in z band (probPSF_z)
      integer*4     zc          ! [0,9] Phototype in z band, 6=Star (type_z) (1)
      real*8        zDate       ! (yr) Date of observation in z band
      real*4        z_mag       ! (mag) ? Model magnitude in z filter (= zmag) (5)
      real*4        e_z_mag     ! (mag) ? Mean error on z'mag (= e_zmag)
      real*4        zpmag       ! (mag) ? PSF magnitude in z band (psfMag_z)
      real*4        e_zpmag     ! (mag) ? Mean error on zpmag (psfMagErr_z)
      real*4        zPmag_1     ! (mag) ? Petrosian magnitude in z band (petroMag_z)
      real*4        e_zPmag_1   ! (mag) ? Mean error on  zPmag (petroMagErr_z)
      real*8        zPrad       ! (arcsec) ? Petrosian radius in z band (petroMagErr_z)
      real*8        e_zPrad     ! (arcsec) ? Mean error in zPrad (petroMagErr_z)
      real*4        zoRA        ! (arcsec) ? Offset of z position along RA (offsetRa_z)
      real*4        zoDE        ! (arcsec) ? Offset of z position along Dec (offsetDec_z)
      real*8        zisoA       ! (arcsec) ? Isophotal major diameter in z band (isoA_z)
      real*8        zisoB       ! (arcsec) ? Isophotal minor diameter in z band (isoB_z)
      real*4        zPA         ! (deg) ? Position angle of zisoA in g band (isoPhi_z)
*Note (1): the type (class) is numbered as follows:
*      0 = Unknown: Object type is not known.
*      1 = Cosmic-ray track (not used).
*      2 = Defect (not used)
*      3 = Galaxy: An extended object composed of many stars and other matter.
*      4 = Ghost: Object created by reflected or refracted light. (not used)
*      5 = KnownObject: Object came from some other catalog (not yet used)
*      6 = Star: A a self-luminous gaseous celestial body.
*      7 = Trail: A satellite or asteroid or meteor trail. (not yet used)
*      8 = Sky: Blank sky spectrogram (no objects in this arcsecond area).
*      9 = NotAType
*Note (2): the SDSS identifier is composed of the skyVersion number,
*      run, rerun, camcol, field and object number. The SDSS-ID
*      is written run.cam.field.obj (4 numbers and 3 separator dots).
*      The objID is also combined from these numbers:
*      objID = skyVersion<<59 + rerun<<48 + run<<32
*            + camcol<<29 + field<<16 + obj
*Note (3): flags written in hexadecimal, full description available at
*      http://cas.sdss.org/astro/en/help/browser/enum.asp?n=PhotoFlags
* ------------------------------------------------------------------------------
* 0000000000000001 = CANONICAL_CENTER   Measurements used the center in r*
* 0000000000000002 = BRIGHT             Object detected in bright object-finding
* 0000000000000004 = EDGE               Object is too close to edge of frame
* 0000000000000008 = BLENDED            Multiple peaks detected
* 0000000000000010 = CHILD              Product of deblending a BLENDED object.
* 0000000000000020 = PEAKCENTER         Position of peak pixel rather than ML
* 0000000000000040 = NODEBLEND          No deblending attempted on BLENDED object
* 0000000000000080 = NOPROFILE          Radial profile could not be estimated
* 0000000000000100 = NOPETRO            No valid Petrosian radius found
* 0000000000000200 = MANYPETRO          More than one Petrosian radius found.
* 0000000000000400 = NOPETRO_BIG        Petrosian radius beyond radial profile.
* 0000000000000800 = DEBLEND_TOO_MANY_PEAKS More than 25 peaks in this object
* 0000000000001000 = COSMIC_RAY         Contains a pixel part of a cosmic ray.
* 0000000000002000 = MANYR50            More than one 50% light radius.
* 0000000000004000 = MANYR90            More than one 90% light radius.
* 0000000000008000 = BAD_RADIAL         Some points in radial profile S/N<0
* 0000000000010000 = INCOMPLETE_PROFILE Petrosian radius intersects frame edge
* 0000000000020000 = INTERP             Pixel(s) determined by interpolation.
* 0000000000040000 = SATURATED          Saturated pixel(s)
* 0000000000080000 = NOTCHECKED         Not checked for local peaks
* 0000000000100000 = SUBTRACTED         BRIGHT object, wings subtracted
* 0000000000200000 = NOSTOKES           No measured Stokes params
* 0000000000400000 = BADSKY             Sky level so bad that highest pixel <0.
* 0000000000800000 = PETROFAINT         Petrosian radius(es) rejected for low SB
* 0000000001000000 = TOO_LARGE          Object too large to measure its profile
* 0000000002000000 = DEBLENDED_AS_PSF   Deblender treated obj as PSF
* 0000000004000000 = DEBLEND_PRUNED     At least one child was removed
* 0000000008000000 = ELLIPFAINT         Object center fainter than the isophote
* 0000000010000000 = BINNED1            Object detected in 1x1 binned image
* 0000000020000000 = BINNED2            Object detected in 2x2 binned image
* 0000000040000000 = BINNED4            Object detected in 4x4 binned image
* 0000000080000000 = MOVED              Possibly moving object
* 0000000100000000 = DEBLENDED_AS_MOVING A MOVED object treated as moving.
* 0000000200000000 = NODEBLEND_MOVING    A MOVED object not treated as moving.
* 0000000400000000 = TOO_FEW_DETECTIONS Too few detections to deblend as moving.
* 0000000800000000 = BAD_MOVING_FIT     Moving fit too poor to be believable.
* 0000001000000000 = STATIONARY         Object consistent with being stationary.
* 0000002000000000 = PEAKS_TOO_CLOSE    Some peaks merged into a single peak.
* 0000004000000000 = MEDIAN_CENTER      Center given is of median-smoothed image.
* 0000008000000000 = LOCAL_EDGE         Center too close to an edge.
* 0000010000000000 = BAD_COUNTS_ERROR   Too few good pixels, error underestimated
* 0000020000000000 = BAD_MOVING_FIT_CHILD Moving child's velocity fit too poor
* 0000040000000000 = DEBLEND_UNASSIGNED_FLUX  After deblending flux not assigned
* 0000080000000000 = SATUR_CENTER        Center is close to saturated pixel(s)
* 0000100000000000 = INTERP_CENTER       Center is close to interpolated pixel(s)
* 0000200000000000 = DEBLENDED_AT_EDGE   Object close to edge was deblended
* 0000400000000000 = DEBLEND_NOPEAK      No detected peak within this child
* 0000800000000000 = PSF_FLUX_INTERP     PSF flux from interpolated pixels > 20%
* 0001000000000000 = TOO_FEW_GOOD_DETECTIONS Too few good detections
* 0002000000000000 = CENTER_OFF_AIMAGE   At least one peak's center lay off image
* 0004000000000000 = DEBLEND_DEGENERATE  Identical candidate children, 1 retained
* 0008000000000000 = BRIGHTEST_GALAXY_CHILD Brightest member classified as galaxy
* 0010000000000000 = CANONICAL_BAND       Band chosen as 'Canonical', generally r
* 0020000000000000 = AMOMENT_FAINT        Object too faint to measure moments
* 0040000000000000 = AMOMENT_SHIFT        Centroid shift too large
* 0080000000000000 = AMOMENT_MAXITER      Maximum number of iterations exceeded
* 0100000000000000 = MAYBE_CR             Object could actually be a cosmic ray
* 0200000000000000 = MAYBE_EGHOST         Object could actually be a ghost
* 0400000000000000 = NOTCHECKED_CENTER    Object center in a region not searched
* 0800000000000000 = OBJECT2_HAS_SATUR_DN Trails included in estimated flux.
* 1000000000000000 = OBJECT2_DEBLEND_PEEPHOLE Deblend modified by peephole
* 2000000000000000 = GROWN_MERGED         Growing led to a merger
* 4000000000000000 = HAS_CENTER           Object has a canonical center
* 8000000000000000 = RESERVED             Not used
* ------------------------------------------------------------------------------
*Note (4): status flags written in hexadecimal, full description available at
* http://cas.sdss.org/astro/en/help/browser/enum.asp?n=PhotoStatus
* ------------------------------------------------------------------------------
* 00000001 = SET         Object status has been set in reference to its own run
* 00000002 = GOOD        Object is good -- absence implies bad.
* 00000004 = DUPLICATE   Duplicate detections in an adjacent field
* 00000010 = OK_RUN      Object is usable
* 00000020 = RESOLVED    Object has been resolved against other runs.
* 00000040 = PSEGMENT    Object Belongs to a PRIMARY segment.
* 00000100 = FIRST_FIELD Object belongs to the first field in its segment.
* 00000200 = OK_SCANLINE Object lies within valid nu range for its scanline.
* 00000400 = OK_STRIPE   Object lies within valid eta range for its stripe.
* 00001000 = SECONDARY   This is a secondary survey object.
* 00002000 = PRIMARY     This is a primary survey object.
* 00004000 = TARGET      This is a spectroscopic target.
* ------------------------------------------------------------------------------
*Note (5):
*    The response curves of the SDSS filters (u g r i z) can be found on
*    the SDSS pages: http://www.sdss.org/dr4/instruments/imager/index.html
*    The central wavelength and FWHM are:
*    -----------------------------------------------
*     (nm)      u'      g'      r'      i'      z'
*    -----------------------------------------------
*    lambda   354.3    477.0   623.1   762.5   913.4
*    FWHM      56.7    138.7   137.3   152.6    95.0
*    -----------------------------------------------

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

C  Loading file 'sdss4.sam'	! *Sample output of the SDSS-DR4 photometric catalog

C  Format for file interpretation

    1 format(
     +  I1,1X,I1,1X,A19,A1,1X,A1,1X,A15,1X,A18,1X,A18,1X,A18,1X,I2,1X,
     +  A16,1X,A8,2X,F10.6,F10.6,1X,F5.3,1X,F5.3,1X,F5.2,1X,F9.4,1X,
     +  I1,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,A16,1X,I1,1X,I1,1X,F9.4,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F7.3,1X,F6.2,
     +  1X,F6.2,1X,F7.3,1X,F7.3,1X,F5.1,1X,A16,1X,I1,1X,I1,1X,F9.4,1X,
     +  F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F7.3,
     +  1X,F6.2,1X,F6.2,1X,F7.3,1X,F7.3,1X,F5.1,1X,A16,1X,I1,1X,I1,1X,
     +  F9.4,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F7.3,
     +  1X,F7.3,1X,F6.2,1X,F6.2,1X,F7.3,1X,F7.3,1X,F5.1,1X,A16,1X,I1,
     +  1X,I1,1X,F9.4,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,
     +  1X,F7.3,1X,F7.3,1X,F6.2,1X,F6.2,1X,F7.3,1X,F7.3,1X,F5.1,1X,
     +  A16,1X,I1,1X,I1,1X,F9.4,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F6.3,1X,F6.3,1X,F7.3,1X,F7.3,1X,F6.2,1X,F6.2,1X,F7.3,1X,F7.3,
     +  1X,F5.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'sdss4.sam')
      write(6,*) '....Loading file: sdss4.sam'
      do i__=1,100
        read(1,'(A871)')ar__
        read(ar__,1)
     +  mode,cl,SDSS,m_SDSS,R,SDSS_ID,objID,SpObjID,parentID,rerun,
     +  flags,Status,RAdeg,DEdeg,e_RAdeg,e_DEdeg,ccPos,ObsDate,Q,umag,
     +  e_umag,gmag,e_gmag,rmag,e_rmag,imag,e_imag,zmag,e_zmag,uFlags,
     +  us,uc,uDate,u_mag,e_u_mag,upmag,e_upmag,uPmag_1,e_uPmag_1,
     +  uPrad,e_uPrad,uoRA,uoDE,uisoA,uisoB,uPA,gFlags,gs,gc,gDate,
     +  g_mag,e_g_mag,gpmag,e_gpmag,gPmag_1,e_gPmag_1,gPrad,e_gPrad,
     +  goRA,goDE,gisoA,gisoB,gPA,rFlags,rs,rc,rDate,r_mag,e_r_mag,
     +  rpmag,e_rpmag,rPmag_1,e_rPmag_1,rPrad,e_rPrad,roRA,roDE,risoA,
     +  risoB,rPA,iFlags,is,ic,iDate,i_mag,e_i_mag,ipmag,e_ipmag,
     +  iPmag_1,e_iPmag_1,iPrad,e_iPrad,ioRA,ioDE,iisoA,iisoB,iPA,
     +  zFlags,zs,zc,zDate,z_mag,e_z_mag,zpmag,e_zpmag,zPmag_1,
     +  e_zPmag_1,zPrad,e_zPrad,zoRA,zoDE,zisoA,zisoB,zPA
        if(ar__(278:283) .EQ. '') u_mag = rNULL__
        if(ar__(285:290) .EQ. '') e_u_mag = rNULL__
        if(ar__(292:297) .EQ. '') upmag = rNULL__
        if(ar__(299:304) .EQ. '') e_upmag = rNULL__
        if(ar__(306:311) .EQ. '') uPmag_1 = rNULL__
        if(ar__(313:318) .EQ. '') e_uPmag_1 = rNULL__
        if(ar__(320:326) .EQ. '') uPrad = rNULL__
        if(ar__(328:334) .EQ. '') e_uPrad = rNULL__
        if(ar__(336:341) .EQ. '') uoRA = rNULL__
        if(ar__(343:348) .EQ. '') uoDE = rNULL__
        if(ar__(350:356) .EQ. '') uisoA = rNULL__
        if(ar__(358:364) .EQ. '') uisoB = rNULL__
        if(ar__(366:370) .EQ. '') uPA = rNULL__
        if(ar__(403:408) .EQ. '') g_mag = rNULL__
        if(ar__(410:415) .EQ. '') e_g_mag = rNULL__
        if(ar__(417:422) .EQ. '') gpmag = rNULL__
        if(ar__(424:429) .EQ. '') e_gpmag = rNULL__
        if(ar__(431:436) .EQ. '') gPmag_1 = rNULL__
        if(ar__(438:443) .EQ. '') e_gPmag_1 = rNULL__
        if(ar__(445:451) .EQ. '') gPrad = rNULL__
        if(ar__(453:459) .EQ. '') e_gPrad = rNULL__
        if(ar__(461:466) .EQ. '') goRA = rNULL__
        if(ar__(468:473) .EQ. '') goDE = rNULL__
        if(ar__(475:481) .EQ. '') gisoA = rNULL__
        if(ar__(483:489) .EQ. '') gisoB = rNULL__
        if(ar__(491:495) .EQ. '') gPA = rNULL__
        if(ar__(528:533) .EQ. '') r_mag = rNULL__
        if(ar__(535:540) .EQ. '') e_r_mag = rNULL__
        if(ar__(542:547) .EQ. '') rpmag = rNULL__
        if(ar__(549:554) .EQ. '') e_rpmag = rNULL__
        if(ar__(556:561) .EQ. '') rPmag_1 = rNULL__
        if(ar__(563:568) .EQ. '') e_rPmag_1 = rNULL__
        if(ar__(570:576) .EQ. '') rPrad = rNULL__
        if(ar__(578:584) .EQ. '') e_rPrad = rNULL__
        if(ar__(586:591) .EQ. '') roRA = rNULL__
        if(ar__(593:598) .EQ. '') roDE = rNULL__
        if(ar__(600:606) .EQ. '') risoA = rNULL__
        if(ar__(608:614) .EQ. '') risoB = rNULL__
        if(ar__(616:620) .EQ. '') rPA = rNULL__
        if(ar__(653:658) .EQ. '') i_mag = rNULL__
        if(ar__(660:665) .EQ. '') e_i_mag = rNULL__
        if(ar__(667:672) .EQ. '') ipmag = rNULL__
        if(ar__(674:679) .EQ. '') e_ipmag = rNULL__
        if(ar__(681:686) .EQ. '') iPmag_1 = rNULL__
        if(ar__(688:693) .EQ. '') e_iPmag_1 = rNULL__
        if(ar__(695:701) .EQ. '') iPrad = rNULL__
        if(ar__(703:709) .EQ. '') e_iPrad = rNULL__
        if(ar__(711:716) .EQ. '') ioRA = rNULL__
        if(ar__(718:723) .EQ. '') ioDE = rNULL__
        if(ar__(725:731) .EQ. '') iisoA = rNULL__
        if(ar__(733:739) .EQ. '') iisoB = rNULL__
        if(ar__(741:745) .EQ. '') iPA = rNULL__
        if(ar__(778:783) .EQ. '') z_mag = rNULL__
        if(ar__(785:790) .EQ. '') e_z_mag = rNULL__
        if(ar__(792:797) .EQ. '') zpmag = rNULL__
        if(ar__(799:804) .EQ. '') e_zpmag = rNULL__
        if(ar__(806:811) .EQ. '') zPmag_1 = rNULL__
        if(ar__(813:818) .EQ. '') e_zPmag_1 = rNULL__
        if(ar__(820:826) .EQ. '') zPrad = rNULL__
        if(ar__(828:834) .EQ. '') e_zPrad = rNULL__
        if(ar__(836:841) .EQ. '') zoRA = rNULL__
        if(ar__(843:848) .EQ. '') zoDE = rNULL__
        if(ar__(850:856) .EQ. '') zisoA = rNULL__
        if(ar__(858:864) .EQ. '') zisoB = rNULL__
        if(ar__(866:870) .EQ. '') zPA = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  mode,cl,SDSS,m_SDSS,R,SDSS_ID,objID,SpObjID,parentID,rerun,
     +  flags,Status,RAdeg,DEdeg,e_RAdeg,e_DEdeg,ccPos,ObsDate,Q,umag,
     +  e_umag,gmag,e_gmag,rmag,e_rmag,imag,e_imag,zmag,e_zmag,uFlags,
     +  us,uc,uDate,u_mag,e_u_mag,upmag,e_upmag,uPmag_1,e_uPmag_1,
     +  uPrad,e_uPrad,uoRA,uoDE,uisoA,uisoB,uPA,gFlags,gs,gc,gDate,
     +  g_mag,e_g_mag,gpmag,e_gpmag,gPmag_1,e_gPmag_1,gPrad,e_gPrad,
     +  goRA,goDE,gisoA,gisoB,gPA,rFlags,rs,rc,rDate,r_mag,e_r_mag,
     +  rpmag,e_rpmag,rPmag_1,e_rPmag_1,rPrad,e_rPrad,roRA,roDE,risoA,
     +  risoB,rPA,iFlags,is,ic,iDate,i_mag,e_i_mag,ipmag,e_ipmag,
     +  iPmag_1,e_iPmag_1,iPrad,e_iPrad,ioRA,ioDE,iisoA,iisoB,iPA,
     +  zFlags,zs,zc,zDate,z_mag,e_z_mag,zpmag,e_zpmag,zPmag_1,
     +  e_zPmag_1,zPrad,e_zPrad,zoRA,zoDE,zisoA,zisoB,zPA
c    .......End.of.Just test output...........
      end do
      close(1)

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