Conversion of standardized ReadMe file for
file /./ftp/cats/II/351 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-19
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/351 VISTA Magellanic Survey (VMC) catalog (Cioni+, 2011)
*================================================================================
*The VISTA Magellanic Survey (VMC) survey.
* Cioni M.-R.L., Clementini G., Girardi L., Guandalini R., Gullieuszik M.,
* Miszalski B., Moretti M.-I., Ripepi V., Rubele S., Bagheri G., Bekki K.,
* Cross N., de Blok W.J.G., de Grijs R., Emerson J.P., Evans C.J., Gibson B.,
* Gonzales-Solares E., Groenewegen M.A.T., Irwin M., Ivanov V.D., Lewis J.,
* Marconi M., Marquette J.-B., Mastropietro C., Moore B., Napiwotzki R.,
* Naylor T., Oliveira J.M., Read M., Sutorius E., van Loon J.Th.,
* Wilkinson M.I., Wood P.R.
* <Astron. Astrophys. 527, A116 (2011) >
* =2011A&A...527A.116C
* =2017yCat.2351....0C
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'vmc_dr4.sam' ! VISTA Magellanic Survey: band-merged master
source catalog (YJKs) DR4; sample output
(on a total of 13,857,646 sources)
integer*4 nr__
parameter (nr__=1000) ! Number of records
character*796 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*8 SrcID (nr__) ! [558345748481/558359606126] Unique ID of this
* merged detection as assigned by merge algorithm
* (SOURCEID)
real*8 RAdeg (nr__) ! (deg) Celestial Right Ascension (J2000) (RA2000)
real*8 DEdeg (nr__) ! (deg) Celestial Declination (J2000) (DEC2000)
character*24 Name (nr__) ! IAU Name (not unique; VMC JHHMMSS.ss+DDMMSS.ss)
* (IAUNAME) (1)
integer*4 CuID (nr__) ! [530311] Unique ID of curation event giving rise
* to this record (always =530311) (CUEVENTID)
integer*8 FrID (nr__) ! [558345748481/558345748499] Unique ID of the set
* of frames that this merged source comes from
* (FRAMESETID)
real*8 GLON (nr__) ! (deg) Galactic longitude (L)
real*8 GLAT (nr__) ! (deg) Galactic latitude (B)
integer*8 priOrSec (nr__) ! Seam code for a unique (=0) or duplicated (!=0)
* source (PRIORSEC) (2)
real*8 Y_Jpnt (nr__) ! (mag) [-10/10]? Point source colour index Y-J
* (using APERMAG3) (YMJPNT)
real*8 e_Y_Jpnt (nr__) ! (mag) [ /1016]? Error on Y-Jpnt (YMJPNTERR)
real*8 J_Kspnt (nr__) ! (mag) [-7/12]? Point source colour index J-Ks
* (using APERMAG3) (JMKSPNT)
real*8 e_J_Kspnt (nr__) ! (mag) [/1482]? Error on J-Kspnt (JMKSPNTERR)
real*8 Y_Jext (nr__) ! (mag) [-10/10]? Extended source colour index Y-J
* (using APERMAGNOAPERCORR3) (YMJEXT)
real*8 e_Y_Jext (nr__) ! (mag) [/1016]? Error on Y-Jext (YMJEXTERR)
real*8 J_Ksext (nr__) ! (mag) [-7/12]? Extended source colour index J-Ks
* (using APERMAGNOAPERCORR3) (JMKSEXT)
real*8 e_J_Ksext (nr__) ! (mag) [/1482]? Error on J-Ksext (JMKSEXTERR)
real*8 Clst (nr__) ! [-128/359] Merged N(0,1) stellarness-of-profile
* statistic (MERGEDCLASSSTAT)
integer*4 Cl (nr__) ! [-3/1]? Class flag (MERGEDCLASS) (3)
real*8 pStar (nr__) ! Probability that the source is a star (PSTAR)
real*8 pGal (nr__) ! Probability that the source is a galaxy (PGALAXY)
real*8 pNoise (nr__) ! Probability that the source is noise (PNOISE)
real*8 KsMJD (nr__) ! (d) ? Modified Julian Day in Ks band (KSMJD)
real*8 Kspmag (nr__) ! (mag) [3/34]? Extended source VISTA Ks mag (Petrosian)
* (KSPETROMAG)
real*8 e_Kspmag (nr__) ! (mag) [/88013]? Error in Kspmag (KSPETROMAGERR)
real*8 Ksap3 (nr__) ! (mag) [9/29]? Default point source Ks aperture
* corrected mag (2.0 arcsec diameter) (KSAPERMAG3)
real*8 e_Ksap3 (nr__) ! (mag) [/443]? Error in Ksap3 (KSAPERMAG3ERR)
real*8 Ksap4 (nr__) ! (mag) [9/34]? Point source Ks aperture corrected mag
* (2.8 arcsec aperture diameter) (KSAPERMAG4)
real*8 e_Ksap4 (nr__) ! (mag) [/76448]? Error in Ksap4 (KSAPERMAG4ERR)
real*8 Ksap6 (nr__) ! (mag) [7/36]? Point source Ks aperture corrected mag
* (5.7 arcsec aperture diameter) (KSAPERMAG6)
real*8 e_Ksap6 (nr__) ! (mag) [/4171555]? Error in Ksap6 (KSAPERMAG6ERR)
real*8 Ksapc3 (nr__) ! (mag) [10/29]? Default extended source Ks aperture
* mag (2.0 arcsec diameter) (KSAPERMAGNOAPERCORR3)
real*8 Ksapc4 (nr__) ! (mag) [9/34]? Extended source Ks aperture mag
* (2.8" aperture diameter) (KSAPERMAGNOAPERCORR4)
real*8 Ksapc6 (nr__) ! (mag) [7/37]? Extended source Ks aperture mag
* (5.7" aperture diameter) (KSAPERMAGNOAPERCORR6)
real*8 Ksrms (nr__) ! (pix) [0.7/53]? RMS of axes of ellipse fit in Ks
* (KSGAUSIG)
real*4 Ksell (nr__) ! ? 1-b/a, where a/b=semi-major/minor axes in Ks
* (KSELL)
real*4 KsPA (nr__) ! (deg) ? Ellipse fit celestial orientation in Ks (KSPA)
integer*4 Kserrb (nr__) ! [/27]? Processing warning/error bitwise flags
* in Ks (KSERRBITS)
real*8 Ksavg (nr__) ! [0/112]? Average confidence in 2" diam default
* aperture 3 Ks (KSAVERAGECONF)
integer*4 KsCl (nr__) ! [-2/1]? Discrete image classification flag in
* Ks (KSCLASS) (3)
real*8 KsClst (nr__) ! [-92/298]? N(0,1) stellarness-of-profile
* statistic in Ks (KSCLASSSTAT)
integer*4 Ksperrb (nr__) ! [16/]?=0 Additional Wide Field Astronomy Unit
* (WFAU) post-processing error bits in Ks
* (<256=minor issues) (KSPPERRBITS) (4)
integer*4 Ksseq (nr__) ! ? The running number of the Ks detection
* (KSSEQNUM)
real*8 JMJD (nr__) ! (d) ? Modified Julian Day in J band (JMJD)
real*8 Jpmag (nr__) ! (mag) [6/38]? Extended source VISTA J mag (Petrosian)
* (JPETROMAG)
real*8 e_Jpmag (nr__) ! (mag) [/1108416]? Error in Jpmag (JPETROMAGERR)
real*8 Jap3 (nr__) ! (mag) [10/32]? Default point source J aperture
* corrected mag (2.0 arcsec diam) (JAPERMAG3)
real*8 e_Jap3 (nr__) ! (mag) [/1482]? Error in Jap3 (JAPERMAG3ERR)
real*8 Jap4 (nr__) ! (mag) [10/33]? Point source J aperture corrected mag
* (2.8 arcsec aperture diam) (JAPERMAG4)
real*8 e_Jap4 (nr__) ! (mag) [/11944]? Error in Jap4 (JAPERMAG4ERR)
real*8 Jap6 (nr__) ! (mag) [8/36]? Point source J aperture corrected mag
* (5.7 arcsec aperture diam) (JAPERMAG6)
real*8 e_Jap6 (nr__) ! (mag) [/242845]? Error in Jap6 (JAPERMAG6ERR)
real*8 Japc3 (nr__) ! (mag) [10/32]? Default extended source J aperture
* mag (2.0 arcsec diam) (JAPERMAGNOAPERCORR3)
real*8 Japc4 (nr__) ! (mag) [10/33]? Extended source J aperture mag
* (2.8" aperture diameter) (JAPERMAGNOAPERCORR4)
real*8 Japc6 (nr__) ! (mag) [8/36]? Extended source J aperture mag
* (5.7" aperture diameter) (JAPERMAGNOAPERCORR6)
real*8 Jrms (nr__) ! (pix) [0.7/52]? RMS of axes of ellipse fit in J
* (JGAUSIG)
real*4 Jell (nr__) ! ? 1-b/a, where a/b=semi-major/minor axes in J
* (JELL)
real*4 JPA (nr__) ! (deg) ? Ellipse fit celestial orientation in J (JPA)
integer*4 Jerrb (nr__) ! [0/24]? Processing warning/error bitwise flags
* in J (JERRBITS)
real*8 Javg (nr__) ! [0/112]? Average confidence in 2 arcsec diam
* default aperture 3 J (JAVERAGECONF)
integer*4 JCl (nr__) ! [-2/1]? Discrete image classification flag in J
* (JCLASS) (3)
real*8 JClst (nr__) ! [-76/333]? N(0,1) stellarness-of-profile
* statistic in J (JCLASSSTAT)
integer*4 Jperrb (nr__) ! [16/]?=0 Additional WFAU post-processing error
* bits in J (<256=minor issues) (JPPERRBITS) (4)
integer*4 Jseq (nr__) ! ? The running number of the J detection (JSEQNUM)
real*8 Jxi (nr__) ! (arcsec) ? Offset of J detection from master Ks position
* (+east/-west) (JXI)
real*8 Jeta (nr__) ! (arcsec) ? Offset of J detection from master Ks position
* (+north/-south) (JETA)
real*8 YMJD (nr__) ! (d) ? Modified Julian Day in Y band (YMJD)
real*8 Ypmag (nr__) ! (mag) [6/38]? Extended source VISTA Y mag (Petrosian)
* (YPETROMAG)
real*8 e_Ypmag (nr__) ! (mag) [/296483]? Error in extended source Y mag
* (Petrosian) (YPETROMAGERR)
real*8 Yap3 (nr__) ! (mag) [10/31]? Default point source Y aperture
* corrected mag (2.0 arcsec diam) (YAPERMAG3)
real*8 e_Yap3 (nr__) ! (mag) [/524]? Error in Yap3 (YAPERMAG3ERR)
real*8 Yap4 (nr__) ! (mag) [9/33]? Point source Y aperture corrected mag
* (2.8 arcsec aperture diam) (YAPERMAG4)
real*8 e_Yap4 (nr__) ! (mag) [/5766]? Error in Yap4 (YAPERMAG4ERR)
real*8 Yap6 (nr__) ! (mag) [8/36]? Point source Y aperture corrected mag
* (5.7 arcsec aperture diam) (YAPERMAG6)
real*8 e_Yap6 (nr__) ! (mag) [/154823]? Error in Yap6 (YAPERMAG6ERR)
real*8 Yapc3 (nr__) ! (mag) [10/31]? Default extended source Y aperture
* mag (2.0 arcsec diam) (YAPERMAGNOAPERCORR3)
real*8 Yapc4 (nr__) ! (mag) [10/33]? Extended source Y aperture mag
* (2.8" aperture diameter) (YAPERMAGNOAPERCORR4)
real*8 Yapc6 (nr__) ! (mag) [9/36]? Extended source Y aperture mag
* (5.7" aperture diameter) (YAPERMAGNOAPERCORR6)
real*8 Yrms (nr__) ! (pix) [0.7/53]? RMS of axes of ellipse fit in Y
* (YGAUSIG)
real*4 Yell (nr__) ! ? 1-b/a, where a/b=semi-major/minor axes in Y
* (YELL)
real*4 YPA (nr__) ! (deg) ? Ellipse fit celestial orientation in Y (YPA)
integer*4 Yerrb (nr__) ! [0/24]? Processing warning/error bitwise flags
* in Y (YERRBITS)
real*8 Yavg (nr__) ! [0/112]? Average confidence in 2 arcsec diam
* default aperture 3 Y (YAVERAGECONF)
integer*4 YCl (nr__) ! [-2/1]? Discrete image classification flag in Y
* (YCLASS) (3)
real*8 YClst (nr__) ! [-100/263]? N(0,1) stellarness-of-profile
* statistic in Y (YCLASSSTAT)
integer*4 Yperrb (nr__) ! [16/]?=0 Additional WFAU post-processing error
* bits in Y (<256=minor issues) (YPPERRBITS) (4)
integer*4 Yseq (nr__) ! ? The running number of the Y detection (YSEQNUM)
real*8 Yxi (nr__) ! (arcsec) ? Offset of Y detection from master Ks position
* (+east/-west) (YXI)
real*8 Yeta (nr__) ! (arcsec) ? Offset of Y detection from master Ks position
* (+north/-south) (YETA)
integer*4 Var (nr__) ! [0/1]? Classification of variability across all
* bands (VARFLAG)
integer*4 Prim (nr__) ! [0/1]? 1=Primary source; 0=secondary source
* (PRIMARY_SOURCE)
*Note (1): The IAUNAME of sources in the PSF catalogues may not be
* unique. At this stage, sources in the overlap of tiles will appear
* with the same IAUNAME. Furthermore, the IAUNAME is rounded to two
* decimal points in arcsec, hence, it may be possible that two
* sufficiently close extractions result in two sources with the
* same IAUNAME.
*Note (2): Sources are unique within each tile. Where PRIORSEC>0
* signifies that a source is located in a region of overlap with
* an adjacent tile that is not yet part of the current release.
*Note (3): Class flag as follows:
* 1 = gal (6605196 occurrences in merged bands)
* 0 = noise (36890 occurrences in merged bands)
* -1 = star (6473980 occurrences in merged bands)
* -2 = probStar (702149 occurrences in merged bands)
* -3 = probGal (39431 occurrences in merged bands)
* -9 = saturated (0 occurrence)
*Note (4): Quality error bit flags assigned during post processing
* are listed at http://horus.roe.ac.uk/vsa/ppErrBits.html
* These flags refer to quality issues of varying severity. For each
* pass-band nine quality issues are implemented as follows, where the
* corresponding value of the ppErrBit is given in parenthesis.
* Source is deblended (16), has bad pixel(s) in default aperture (64),
* has low confidence in default aperture (128), lies within detector #16
* region of a tile (4096), is close to saturation (65536), has
* photometric calibration probably subject to systematic errors
* (131072), lies within a dither offset of the stacked frame boundary
* (4194304), lies within the underexposed strip of a tile (8388606), and
* lies within an underexposed region of a tile due to missing detector
* (16777216). To select only sources without quality issues the user can
* filter on ppErrBits=0, but note that all sources in this catalog
* have at least ppErrBits=16 due to the dense stellar field, and to
* include only sources with minor quality issues use ppErrBits<256.
C=============================================================================
C Loading file 'vmc_dr4.sam' ! VISTA Magellanic Survey: band-merged master
* source catalog (YJKs) DR4; sample output
* (on a total of 13,857,646 sources)
C Format for file interpretation
1 format(
+ I12,1X,F9.6,1X,F10.6,1X,A24,1X,I6,1X,I12,1X,F10.6,1X,F10.6,1X,
+ I12,1X,F7.4,1X,F9.4,1X,F7.4,1X,F9.4,1X,F7.4,1X,F9.4,1X,F7.4,
+ 1X,F9.4,1X,F9.4,1X,I2,1X,F8.6,1X,F8.6,1X,F8.6,1X,F14.8,1X,
+ F7.4,1X,F10.4,1X,F7.4,1X,F8.4,1X,F7.4,1X,F10.4,1X,F7.4,1X,
+ F12.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F6.4,1X,F5.1,1X,I2,
+ 1X,F7.3,1X,I2,1X,F7.3,1X,I8,1X,I7,1X,F14.8,1X,F7.4,1X,F12.4,
+ 1X,F7.4,1X,F9.4,1X,F7.4,1X,F10.4,1X,F7.4,1X,F11.4,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F6.4,1X,F5.1,1X,I2,1X,F7.3,1X,I2,1X,
+ F7.3,1X,I8,1X,I7,1X,F9.6,1X,F9.6,1X,F14.8,1X,F7.4,1X,F11.4,1X,
+ F7.4,1X,F8.4,1X,F7.4,1X,F9.4,1X,F7.4,1X,F11.4,1X,F7.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F6.4,1X,F5.1,1X,I2,1X,F7.3,1X,I2,1X,F7.3,
+ 1X,I8,1X,I7,1X,F9.6,1X,F9.6,2X,I1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'vmc_dr4.sam')
write(6,*) '....Loading file: vmc_dr4.sam'
do i__=1,1000
read(1,'(A796)')ar__
read(ar__,1)
+ SrcID(i__),RAdeg(i__),DEdeg(i__),Name(i__),CuID(i__),
+ FrID(i__),GLON(i__),GLAT(i__),priOrSec(i__),Y_Jpnt(i__),
+ e_Y_Jpnt(i__),J_Kspnt(i__),e_J_Kspnt(i__),Y_Jext(i__),
+ e_Y_Jext(i__),J_Ksext(i__),e_J_Ksext(i__),Clst(i__),Cl(i__),
+ pStar(i__),pGal(i__),pNoise(i__),KsMJD(i__),Kspmag(i__),
+ e_Kspmag(i__),Ksap3(i__),e_Ksap3(i__),Ksap4(i__),e_Ksap4(i__),
+ Ksap6(i__),e_Ksap6(i__),Ksapc3(i__),Ksapc4(i__),Ksapc6(i__),
+ Ksrms(i__),Ksell(i__),KsPA(i__),Kserrb(i__),Ksavg(i__),
+ KsCl(i__),KsClst(i__),Ksperrb(i__),Ksseq(i__),JMJD(i__),
+ Jpmag(i__),e_Jpmag(i__),Jap3(i__),e_Jap3(i__),Jap4(i__),
+ e_Jap4(i__),Jap6(i__),e_Jap6(i__),Japc3(i__),Japc4(i__),
+ Japc6(i__),Jrms(i__),Jell(i__),JPA(i__),Jerrb(i__),Javg(i__),
+ JCl(i__),JClst(i__),Jperrb(i__),Jseq(i__),Jxi(i__),Jeta(i__),
+ YMJD(i__),Ypmag(i__),e_Ypmag(i__),Yap3(i__),e_Yap3(i__),
+ Yap4(i__),e_Yap4(i__),Yap6(i__),e_Yap6(i__),Yapc3(i__),
+ Yapc4(i__),Yapc6(i__),Yrms(i__),Yell(i__),YPA(i__),Yerrb(i__),
+ Yavg(i__),YCl(i__),YClst(i__),Yperrb(i__),Yseq(i__),Yxi(i__),
+ Yeta(i__),Var(i__),Prim(i__)
if(ar__(115:121) .EQ. '') Y_Jpnt(i__) = rNULL__
if(ar__(123:131) .EQ. '') e_Y_Jpnt(i__) = rNULL__
if(ar__(133:139) .EQ. '') J_Kspnt(i__) = rNULL__
if(ar__(141:149) .EQ. '') e_J_Kspnt(i__) = rNULL__
if(ar__(151:157) .EQ. '') Y_Jext(i__) = rNULL__
if(ar__(159:167) .EQ. '') e_Y_Jext(i__) = rNULL__
if(ar__(169:175) .EQ. '') J_Ksext(i__) = rNULL__
if(ar__(177:185) .EQ. '') e_J_Ksext(i__) = rNULL__
if(ar__(197:198) .EQ. '') Cl(i__) = iNULL__
if(ar__(227:240) .EQ. '') KsMJD(i__) = rNULL__
if(ar__(242:248) .EQ. '') Kspmag(i__) = rNULL__
if(ar__(250:259) .EQ. '') e_Kspmag(i__) = rNULL__
if(ar__(261:267) .EQ. '') Ksap3(i__) = rNULL__
if(ar__(269:276) .EQ. '') e_Ksap3(i__) = rNULL__
if(ar__(278:284) .EQ. '') Ksap4(i__) = rNULL__
if(ar__(286:295) .EQ. '') e_Ksap4(i__) = rNULL__
if(ar__(297:303) .EQ. '') Ksap6(i__) = rNULL__
if(ar__(305:316) .EQ. '') e_Ksap6(i__) = rNULL__
if(ar__(318:324) .EQ. '') Ksapc3(i__) = rNULL__
if(ar__(326:332) .EQ. '') Ksapc4(i__) = rNULL__
if(ar__(334:340) .EQ. '') Ksapc6(i__) = rNULL__
if(ar__(342:348) .EQ. '') Ksrms(i__) = rNULL__
if(ar__(350:355) .EQ. '') Ksell(i__) = rNULL__
if(ar__(357:361) .EQ. '') KsPA(i__) = rNULL__
if(ar__(363:364) .EQ. '') Kserrb(i__) = iNULL__
if(ar__(366:372) .EQ. '') Ksavg(i__) = rNULL__
if(ar__(374:375) .EQ. '') KsCl(i__) = iNULL__
if(ar__(377:383) .EQ. '') KsClst(i__) = rNULL__
if(ar__(394:400) .EQ. '') Ksseq(i__) = iNULL__
if(ar__(402:415) .EQ. '') JMJD(i__) = rNULL__
if(ar__(417:423) .EQ. '') Jpmag(i__) = rNULL__
if(ar__(425:436) .EQ. '') e_Jpmag(i__) = rNULL__
if(ar__(438:444) .EQ. '') Jap3(i__) = rNULL__
if(ar__(446:454) .EQ. '') e_Jap3(i__) = rNULL__
if(ar__(456:462) .EQ. '') Jap4(i__) = rNULL__
if(ar__(464:473) .EQ. '') e_Jap4(i__) = rNULL__
if(ar__(475:481) .EQ. '') Jap6(i__) = rNULL__
if(ar__(483:493) .EQ. '') e_Jap6(i__) = rNULL__
if(ar__(495:501) .EQ. '') Japc3(i__) = rNULL__
if(ar__(503:509) .EQ. '') Japc4(i__) = rNULL__
if(ar__(511:517) .EQ. '') Japc6(i__) = rNULL__
if(ar__(519:525) .EQ. '') Jrms(i__) = rNULL__
if(ar__(527:532) .EQ. '') Jell(i__) = rNULL__
if(ar__(534:538) .EQ. '') JPA(i__) = rNULL__
if(ar__(540:541) .EQ. '') Jerrb(i__) = iNULL__
if(ar__(543:549) .EQ. '') Javg(i__) = rNULL__
if(ar__(551:552) .EQ. '') JCl(i__) = iNULL__
if(ar__(554:560) .EQ. '') JClst(i__) = rNULL__
if(ar__(571:577) .EQ. '') Jseq(i__) = iNULL__
if(ar__(579:587) .EQ. '') Jxi(i__) = rNULL__
if(ar__(589:597) .EQ. '') Jeta(i__) = rNULL__
if(ar__(599:612) .EQ. '') YMJD(i__) = rNULL__
if(ar__(614:620) .EQ. '') Ypmag(i__) = rNULL__
if(ar__(622:632) .EQ. '') e_Ypmag(i__) = rNULL__
if(ar__(634:640) .EQ. '') Yap3(i__) = rNULL__
if(ar__(642:649) .EQ. '') e_Yap3(i__) = rNULL__
if(ar__(651:657) .EQ. '') Yap4(i__) = rNULL__
if(ar__(659:667) .EQ. '') e_Yap4(i__) = rNULL__
if(ar__(669:675) .EQ. '') Yap6(i__) = rNULL__
if(ar__(677:687) .EQ. '') e_Yap6(i__) = rNULL__
if(ar__(689:695) .EQ. '') Yapc3(i__) = rNULL__
if(ar__(697:703) .EQ. '') Yapc4(i__) = rNULL__
if(ar__(705:711) .EQ. '') Yapc6(i__) = rNULL__
if(ar__(713:719) .EQ. '') Yrms(i__) = rNULL__
if(ar__(721:726) .EQ. '') Yell(i__) = rNULL__
if(ar__(728:732) .EQ. '') YPA(i__) = rNULL__
if(ar__(734:735) .EQ. '') Yerrb(i__) = iNULL__
if(ar__(737:743) .EQ. '') Yavg(i__) = rNULL__
if(ar__(745:746) .EQ. '') YCl(i__) = iNULL__
if(ar__(748:754) .EQ. '') YClst(i__) = rNULL__
if(ar__(765:771) .EQ. '') Yseq(i__) = iNULL__
if(ar__(773:781) .EQ. '') Yxi(i__) = rNULL__
if(ar__(783:791) .EQ. '') Yeta(i__) = rNULL__
if(ar__(794:794) .EQ. '') Var(i__) = iNULL__
if(ar__(796:796) .EQ. '') Prim(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ SrcID(i__),RAdeg(i__),DEdeg(i__),Name(i__),CuID(i__),
+ FrID(i__),GLON(i__),GLAT(i__),priOrSec(i__),Y_Jpnt(i__),
+ e_Y_Jpnt(i__),J_Kspnt(i__),e_J_Kspnt(i__),Y_Jext(i__),
+ e_Y_Jext(i__),J_Ksext(i__),e_J_Ksext(i__),Clst(i__),Cl(i__),
+ pStar(i__),pGal(i__),pNoise(i__),KsMJD(i__),Kspmag(i__),
+ e_Kspmag(i__),Ksap3(i__),e_Ksap3(i__),Ksap4(i__),e_Ksap4(i__),
+ Ksap6(i__),e_Ksap6(i__),Ksapc3(i__),Ksapc4(i__),Ksapc6(i__),
+ Ksrms(i__),Ksell(i__),KsPA(i__),Kserrb(i__),Ksavg(i__),
+ KsCl(i__),KsClst(i__),Ksperrb(i__),Ksseq(i__),JMJD(i__),
+ Jpmag(i__),e_Jpmag(i__),Jap3(i__),e_Jap3(i__),Jap4(i__),
+ e_Jap4(i__),Jap6(i__),e_Jap6(i__),Japc3(i__),Japc4(i__),
+ Japc6(i__),Jrms(i__),Jell(i__),JPA(i__),Jerrb(i__),Javg(i__),
+ JCl(i__),JClst(i__),Jperrb(i__),Jseq(i__),Jxi(i__),Jeta(i__),
+ YMJD(i__),Ypmag(i__),e_Ypmag(i__),Yap3(i__),e_Yap3(i__),
+ Yap4(i__),e_Yap4(i__),Yap6(i__),e_Yap6(i__),Yapc3(i__),
+ Yapc4(i__),Yapc6(i__),Yrms(i__),Yell(i__),YPA(i__),Yerrb(i__),
+ Yavg(i__),YCl(i__),YClst(i__),Yperrb(i__),Yseq(i__),Yxi(i__),
+ Yeta(i__),Var(i__),Prim(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end