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