Conversion of standardized ReadMe file for
file /./ftp/cats/II/367 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 2026-Oct-11
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/367 The VISTA Hemisphere Survey (VHS) catalog DR5 (McMahon+, 2020)
*================================================================================
*VISTA Hemisphere Survey band-merged multi-waveband catalogues (VHS) DR5
* McMahon, R. G., Banerji, M., Gonzalez, E., Koposov, S. E., Bejar, V. J.,
* Lodieu, N., Rebolo, R., The VHS Collaboration
* <The Messenger. 154, 35-37 (2013)>
* =2021yCat.2367....0M
* =2013Msngr.154...35M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'vhs_dr5.sam' ! VISTA Hemisphere Survey band-merged
multi-waveband catalogues (VHS) DR5 (version 3);
sample output (on a total of 1,374,207,485 sources)
integer*4 nr__
parameter (nr__=1000) ! Number of records
character*1110 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*8 SrcID ! Unique ID of this merged detection as
* assigned by merge algorithm
* unless PrioOrSec!=0 (SOURCEID)
real*8 RAdeg ! (deg) Celestial Right Ascension (J2000)
* (RA2000)
real*8 DEdeg ! (deg) [-90/1.5] Celestial Declination (J2000)
* (DEC2000)
character*23 Name ! Source name in IAU convention
* (VHS JHHMMSS.ss+DDMMSS.s) (SOURCENAME)
integer*4 cuEventID ! [587931] UID of curation event giving
* rise to this record (always "587931")
integer*8 FrameID ! UID of the set of frames that this
* merged source comes from (FRAMESETID)
real*8 GLON ! (deg) Galactic longitude (L)
real*8 GLAT ! (deg) Galactic latitude (B)
integer*8 PriOrSec ! Seam code for a unique (=0) or
* duplicated (!=0) source (5)
real*8 Y_Jpnt ! (mag) [-12.1/10.3]? Point source color Y-J
* (using aperMag3) (YMJPNT)
real*8 e_Y_Jpnt ! (mag) [4E-4/5999]? Error on point source
* color Y-J (YMJPNTERR)
real*8 J_Hpnt ! (mag) [-11.3/9.7]? Point source color J-H
* (using aperMag3) (JMHPNT)
real*8 e_J_Hpnt ! (mag) [5E-4/744]? Error on point source
* color J-H (JMHPNTERR)
real*8 H_Kspnt ! (mag) [-10.7/11.3]? Point source color H-Ks
* (using aperMag3) (HMKSPNT)
real*8 e_H_Kspnt ! (mag) [5E-4/627]? Error on point source
* color H-Ks (HMKSPNTERR)
real*8 J_Kspnt ! (mag) [-12/21.4]? Point source color J-Ks
* (using aperMag3) (JMKSPNT)
real*8 e_J_Kspnt ! (mag) [5E-4/13621]? Error on point source
* color J-Ks (JMKSPNTERR)
real*8 Y_Jext ! (mag) [-12/10.3]? Extended source color Y-J
* (using aperMagNoAperCorr3) (YMJEXT)
real*8 e_Y_Jext ! (mag) [4E-4/5999]? Error on extended source
* color Y-J (YMJEXTERR)
real*8 J_Hext ! (mag) [-11.2/9.8]? Extended source color J-H
* (using aperMagNoAperCorr3) (JMHEXT)
real*8 e_J_Hext ! (mag) [5E-4/744]? Error on extended source
* color J-H (JMHEXTERR)
real*8 H_Ksext ! (mag) [-10.6/11.2]? Extended source color
* H-Ks (using aperMagNoAperCorr3)
* (HMKSEXT)
real*8 e_H_Ksext ! (mag) [5E-4/627]? Error on extended source
* color H-Ks (HMKSEXTERR)
real*8 J_Ksext ! (mag) [-11.8/22.1]? Extended source color
* J-Ks (using aperMagNoAperCorr3)
* (JMKSEXT)
real*8 e_J_Ksext ! (mag) [5E-4/13621]? Error on extended source
* color J-Ks (JMKSEXTERR)
real*8 Mstat ! [-373/538] Merged N(0,1)
* stellarness-of-profile statistic
* (MERGEDCLASSSTAT) (1)
integer*4 Mclass ! [-3/1] Class flag
* (MERGEDCLASS) (2)
real*8 pStar ! [0/1] Probability that the source is a
* star
real*8 pGal ! [0/1] Probability that the source is a
* galaxy (PGALAXY)
real*8 pNoise ! [0/1] Probability that the source is
* noise
real*8 E_B_V ! (mag) [4E-4/41] The galactic dust extinction
* value measured from the Schlegel
* (1998ApJ...500..525S) maps (EBV)
real*8 AY ! (mag) [5E-4/53] The galactic extinction
* correction in the Y band
real*8 AJ ! (mag) [4E-4/38] The galactic extinction
* correction in the J band
real*8 AH ! (mag) [2E-4/24] The galactic extinction
* correction in the H band
real*8 AKs ! (mag) [1E-4/16] The galactic extinction
* correction in the Ks band
real*8 YMJD ! (d) [55160.2/57811.2]? Modified Julian Day
* in Y band
real*8 Ypmag ! (mag) [3.6/37]? Extended source Y (Vega) mag
* (Petrosian) (YPETROMAG) (3)
real*8 e_Ypmag ! (mag) [0/100]? Error in extended source Y mag
* (Petrosian) (YPETROMAGERR) (3)
real*8 Yap3 ! (mag) [4.2/34]? Default point source Y
* aperture corrected Vega mag (2.0" diam)
* (YAPERMAG3) (3)
real*8 e_Yap3 ! (mag) [1E-4/100]? Error in default
* point/extended source Y mag
* (YAPERMAG3ERR) (3)
real*8 Yap4 ! (mag) [4.1/34]? Point source Y aperture
* corrected Vega mag (2.8" aperture diam)
* (YAPERMAG4) (3)
real*8 e_Yap4 ! (mag) [1E-4/100]? Error in point/extended
* source Y mag (2.8 arcsec aperture diam)
* (YAPERMAG4ERR) (3)
real*8 Yap6 ! (mag) [3.1/38]? Point source Y aperture
* corrected Vega mag (5.7" aperture diam)
* (YAPERMAG6) (3)
real*8 e_Yap6 ! (mag) [1E-4/100]? Error in point/extended
* source Y mag (5.7 arcsec aperture diam)
* (YAPERMAG6ERR) (3)
real*8 Yapc3 ! [4.6/35]? Default extended source Y
* aperture Vega mag (2.0 arcsec diam)
* (YAPERMAGNOAPERCORR3)
real*8 Yapc4 ! (mag) [4.4/35]? Extended source Y aperture
* Vega mag (2.8 arcsec aperture diameter)
* (YAPERMAGNOAPERCORR4)
real*8 Yapc6 ! (mag) [3.2/39]? Extended source Y aperture
* Vega mag (5.7 arcsec aperture diameter)
* (YAPERMAGNOAPERCORR6)
real*8 Yamaj ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in Y band
* (YHLCORSMJRADAS)
real*8 Yrms ! (pix) [0.7/530]? RMS of axes of ellipse fit
* in Y (YGAUSIG)
real*4 Yell ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in Y
real*4 YPA ! (deg) [0/180]? Ellipse fit celestial
* orientation in Y
integer*4 Yerrbits ! [-1/28] Processing warning/error
* bitwise flags in Y
real*8 Yavgap3 ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 Y (YAVERAGECONF)
integer*4 Yclass ! [-2/1] Discrete image classification
* flag in Y
real*8 YclasSt ! [-156/444]? N(0,1)
* stellarness-of-profile statistic in Y
* (YCLASSSTAT)
integer*4 Yperrbits ! [0/12587216] Additional WFAU
* post-processing error bits in Y (4)
integer*4 Yseq ! [0/1413053]? The running number of the
* Y detection (YSEQNUM)
real*8 JMJD ! (d) [55139/57844.2]? Modified Julian Day in
* J band
real*8 Jpmag ! (mag) [-24/38]? Extended source J (Vega) mag
* (Petrosian) (JPETROMAG) (3)
real*8 e_Jpmag ! (mag) [0/100]? Error in extended source J mag
* (Petrosian) (JPETROMAGERR) (3)
real*8 Jap3 ! (mag) [-21.1/34.1]? Default point source J
* aperture corrected Vega mag (2.0" diam)
* (JAPERMAG3) (3)
real*8 e_Jap3 ! (mag) [0/100]? Error in default
* point/extended source J mag
* (JAPERMAG3ERR) (3)
real*8 Jap4 ! (mag) [-21.2/39]? Point source J aperture
* corrected Vega mag (2.8" aperture diam)
* (JAPERMAG4) (3)
real*8 e_Jap4 ! (mag) [0/100]? Error in point/extended source
* J mag (2.8 arcsec aperture diam)
* (JAPERMAG4ERR) (3)
real*8 Jap6 ! (mag) [-21.7/39]? Point source J aperture
* corrected Vega mag (5.7" aperture diam)
* (JAPERMAG6) (3)
real*8 e_Jap6 ! (mag) [0/100]? Error in point/extended source
* J mag (5.7 arcsec aperture diam)
* (JAPERMAG6ERR) (3)
real*8 Japc3 ! (mag) [-20.5/35]? Default extended source J
* aperture Vega mag (2.0 arcsec diam)
* (JAPERMAGNOAPERCORR3)
real*8 Japc4 ! (mag) [-21/39]? Extended source J aperture
* Vega mag (2.8 arcsec aperture diameter)
* (JAPERMAGNOAPERCORR4)
real*8 Japc6 ! (mag) [-22/39]? Extended source J aperture
* Vega mag (5.7 arcsec aperture diameter)
* (JAPERMAGNOAPERCORR6)
real*8 Jamaj ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in J band
* (JHLCORSMJRADAS)
real*8 Jrms ! (pix) [0.7/584]? RMS of axes of ellipse fit
* in J (JGAUSIG)
real*4 Jell ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in J
real*4 JPA ! (deg) [0/180]? Ellipse fit celestial
* orientation in J
integer*4 Jerrbits ! [-1/28]? Processing warning/error
* bitwise flags in J
real*8 Javgap3 ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 J (JAVERAGECONF)
integer*4 Jclass ! [-2/1]? Discrete image classification
* flag in J
real*8 JclasStat ! [-373/458]? N(0,1)
* stellarness-of-profile statistic in J
* (JCLASSSTAT)
integer*4 Jperrbits ! [0/12587216]? Additional WFAU
* post-processing error bits in J (4)
integer*4 Jseq ! [1/1661208]? The running number of the
* J detection (JSEQNUM)
real*8 Jxi ! (arcsec) [-2/2]? Offset of J detection from
* master position (+east/-west)
real*8 Jeta ! (arcsec) [-2/2]? Offset of J detection from
* master position (+north/-south)
real*8 HMJD ! (d) [55139/56074.2]? Modified Julian Day in
* H band
real*8 Hpmag ! (mag) [-24/36]? Extended source H (Vega) mag
* (Petrosian) (HPETROMAG) (3)
real*8 e_Hpmag ! (mag) [0/100]? Error in extended source H mag
* (Petrosian) (HPETROMAGERR) (3)
real*8 Hap3 ! (mag) [-21.5/29]? Default point source H
* aperture corrected Vega mag (2.0" diam)
* (HAPERMAG3) (3)
real*8 e_Hap3 ! (mag) [0/100]? Error in default
* point/extended source H mag
* (HAPERMAG3ERR) (3)
real*8 Hap4 ! (mag) [-22/30.1]? Point source H aperture
* corrected Vega mag (2.8" aperture diam)
* (HAPERMAG4) (3)
real*8 e_Hap4 ! (mag) [0/100]? Error in point/extended source
* H mag (2.8 arcsec aperture diam)
* (HAPERMAG4ERR) (3)
real*8 Hap6 ! (mag) [-22.2/35]? Point source H aperture
* corrected Vega mag (5.7" aperture diam)
* (HAPERMAG6) (3)
real*8 e_Hap6 ! (mag) [0/100]? Error in point/extended source
* H mag (5.7 arcsec aperture diam)
* (HAPERMAG6ERR) (3)
real*8 Hapc3 ! (mag) [-21/29.1]? Default extended source H
* aperture Vega mag (2.0 arcsec diam)
* (HAPERMAGNOAPERCORR3)
real*8 Hapc4 ! (mag) [-21.4/30.3]? Extended source H
* aperture Vega mag
* (2.8 arcsec aperture diameter)
* (HAPERMAGNOAPERCORR4)
real*8 Hapc6 ! (mag) [-22.2/35]? Extended source H aperture
* Vega mag (5.7 arcsec aperture diameter)
* (HAPERMAGNOAPERCORR6)
real*8 Hamaj ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in H band
* (HHLCORSMJRADAS)
real*8 Hrms ! (pix) [0.7/584]? RMS of axes of ellipse fit
* in H (HGAUSIG)
real*4 Hell ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in H
real*4 HPA ! (deg) [0/180]? Ellipse fit celestial
* orientation in H
integer*4 Herrbits ! [-1/28]? Processing warning/error
* bitwise flags in H
real*8 Havgap3 ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 H (HAVERAGECONF)
integer*4 Hclass ! [-2/1]? Discrete image classification
* flag in H
real*8 HclasSt ! [-96/383]? N(0,1)
* stellarness-of-profile statistic in H
* (HCLASSSTAT)
integer*4 Hperrbits ! [0/12587216]? Additional WFAU
* post-processing error bits in H (4)
integer*4 Hseq ! [1/125514]? The running number of the H
* detection (HSEQNUM)
real*8 Hxi ! (arcsec) [-2/2]? Offset of H detection from
* master position (+east/-west)
real*8 Heta ! (arcsec) [-2/2]? Offset of H detection from
* master position (+north/-south)
real*8 KsMJD ! (d) [55139/57844.2]? Modified Julian Day in
* Ks band
real*8 Kspmag ! (mag) [-20.3/34]? Extended source Ks (Vega) mag
* (Petrosian) (KSPETROMAG) (3)
real*8 e_Kspmag ! (mag) [0/100]? Error in extended source Ks
* mag (Petrosian) (KSPETROMAGERR) (3)
real*8 Ksap3 ! (mag) [-33.5/31]? Default point source Ks
* aperture corrected Vega mag (2.0" diam)
* (KSAPERMAG3) (3)
real*8 e_Ksap3 ! (mag) [0/96]? Error in default point/extended
* source Ks mag (KSAPERMAG3ERR) (3)
real*8 Ksap4 ! (mag) [-33.4/36]? Point source Ks aperture
* corrected Vega mag (2.8" aperture diam)
* (KSAPERMAG4) (3)
real*8 e_Ksap4 ! (mag) [0/100]? Error in point/extended source
* Ks mag (2.8 arcsec aperture diam)
* (KSAPERMAG4ERR) (3)
real*8 Ksap6 ! (mag) [-33.3/37]? Point source Ks aperture
* corrected Vega mag (5.7" aperture diam)
* (KSAPERMAG6) (3)
real*8 e_Ksap6 ! (mag) [0/100]? Error in point/extended source
* Ks mag (5.7 arcsec aperture diam)
* (KSAPERMAG6ERR) (3)
real*8 Ksapc3 ! (mag) [-33.4/32]? Default extended source Ks
* aperture Vega mag (2.0 arcsec diam)
* (KSAPERMAGNOAPERCORR3)
real*8 Ksapc4 ! (mag) [-33.3/36]? Extended source Ks aperture
* Vega mag (2.8 arcsec aperture diameter)
* (KSAPERMAGNOAPERCORR4)
real*8 Ksapc6 ! (mag) [-33.3/37.1]? Extended source Ks
* aperture Vega mag
* (5.7 arcsec aperture diameter)
* (KSAPERMAGNOAPERCORR6)
real*8 Ksamaj ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in Ks band
* (KSHLCORSMJRADAS)
real*8 Ksrms ! (pix) [0.7/549]? RMS of axes of ellipse fit
* in Ks (KSGAUSIG)
real*4 Ksell ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in Ks
real*4 KsPA ! (deg) [0/180]? Ellipse fit celestial
* orientation in Ks
integer*4 Kserrbits ! [-1/28]? Processing warning/error
* bitwise flags in Ks
real*8 Ksavgap3 ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 Ks
* (KSAVERAGECONF)
integer*4 Ksclass ! [-2/1]? Discrete image classification
* flag in Ks
real*8 KsclasSt ! [-108/385]? N(0,1)
* stellarness-of-profile statistic in Ks
* (KSCLASSSTAT)
integer*4 Ksperrbits ! [0/12587216]? Additional WFAU
* post-processing error bits in Ks (4)
integer*4 Ksseq ! [1/1455550]? The running number of the
* Ks detection (KSSEQNUM)
real*8 Ksxi ! (arcsec) [-2/2]? Offset of Ks detection from
* master position (+east/-west)
real*8 Kseta ! (arcsec) [-2/2]? Offset of Ks detection from
* master position (+north/-south)
real*4 prim ! [0/1]? Flag indicates primary sources
* (=1), otherwise (=0) (PRIMARY_SOURCE)
*Note (1): This profile classification statistic is a continuously distributed,
* Gaussian N(0,1) (i.e. zero mean, unit variance) statistic formed from
* the available individual classification statistics by averaging them
* and multiplying by N1/2 such that cuts on mergedClassStat result in
* completeness being independent of number of frames an object appears
* on, but with reliability improving with the number of frames.
*Note (2): Class flag from available measurements as follows:
* 1 = Galaxy (522,436,266 occurrences);
* 0 = Noise (25,899,664 occurrences);
* -1 = Star (631,842,962 occurrences);
* -2 = Probable star (180,291,372 occurrences);
* -3 = Probable galaxy (13,737,221 occurrences);
*Note (3): Magnitudes with errors >100mag have been set to null value by CDS;
* their errors have been set to 100mag.
*Note (4): Post-processing error quality bit flags assigned to detections in the
* archive curation procedure for survey data. From least to most
* significant byte in the 4-byte integer attribute byte 0 (bits 0 to 7)
* corresponds to information on generally innocuous conditions that are
* nonetheless potentially significant as regards the integrity of that
* detection; byte 1 (bits 8 to 15) corresponds to warnings; byte 2 (bits
* -3 to 23) corresponds to important warnings; and finally byte 3 (bits
* 0 to 31) corresponds to severe warnings.
* See http://horus.roe.ac.uk/vsa/ppErrBits.html
*Note (5): If a source is not found to be duplicated in overlap regions, then
* priOrSec=0; if a source is duplicated, then priOrSec will be set to
* the frameSetID of the source that should be considered the best one to
* use out of the set of duplicates.
C=============================================================================
C Loading file 'vhs_dr5.sam' ! VISTA Hemisphere Survey band-merged
* multi-waveband catalogues (VHS) DR5 (version 3);
* sample output (on a total of 1,374,207,485 sources)
C Format for file interpretation
1 format(
+ I12,1X,F10.6,1X,F10.6,1X,A23,1X,I6,1X,I12,1X,F10.6,1X,F10.6,
+ 1X,I12,1X,F8.4,1X,F9.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F10.4,1X,F8.4,1X,F9.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F10.4,1X,F9.4,1X,I2,1X,F8.6,1X,F8.6,1X,F8.6,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F12.6,1X,F7.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F8.4,1X,F6.4,1X,F5.1,1X,I2,1X,F7.3,1X,
+ I2,1X,F8.3,1X,I8,1X,I7,1X,F12.6,1X,F8.4,1X,F7.4,1X,F8.4,1X,
+ F7.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F8.4,1X,F8.4,
+ 1X,F7.4,1X,F8.4,1X,F6.4,1X,F5.1,1X,I2,1X,F7.3,1X,I2,1X,F8.3,
+ 1X,I8,1X,I7,1X,F9.6,1X,F9.6,1X,F12.6,1X,F8.4,1X,F7.4,1X,F8.4,
+ 1X,F7.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F8.4,1X,
+ F8.4,1X,F7.4,1X,F8.4,1X,F6.4,1X,F5.1,1X,I2,1X,F7.3,1X,I2,1X,
+ F7.3,1X,I8,1X,I6,1X,F9.6,1X,F9.6,1X,F12.6,1X,F8.4,1X,F7.4,1X,
+ F8.4,1X,F7.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F8.4,
+ 1X,F8.4,1X,F7.4,1X,F8.4,1X,F6.4,1X,F5.1,1X,I2,1X,F7.3,1X,I2,
+ 1X,F8.3,1X,I8,1X,I7,1X,F9.6,1X,F9.6,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'vhs_dr5.sam')
write(6,*) '....Loading file: vhs_dr5.sam'
do i__=1,1000
read(1,'(A1110)')ar__
read(ar__,1)
+ SrcID,RAdeg,DEdeg,Name,cuEventID,FrameID,GLON,GLAT,PriOrSec,
+ Y_Jpnt,e_Y_Jpnt,J_Hpnt,e_J_Hpnt,H_Kspnt,e_H_Kspnt,J_Kspnt,
+ e_J_Kspnt,Y_Jext,e_Y_Jext,J_Hext,e_J_Hext,H_Ksext,e_H_Ksext,
+ J_Ksext,e_J_Ksext,Mstat,Mclass,pStar,pGal,pNoise,E_B_V,AY,AJ,
+ AH,AKs,YMJD,Ypmag,e_Ypmag,Yap3,e_Yap3,Yap4,e_Yap4,Yap6,e_Yap6,
+ Yapc3,Yapc4,Yapc6,Yamaj,Yrms,Yell,YPA,Yerrbits,Yavgap3,Yclass,
+ YclasSt,Yperrbits,Yseq,JMJD,Jpmag,e_Jpmag,Jap3,e_Jap3,Jap4,
+ e_Jap4,Jap6,e_Jap6,Japc3,Japc4,Japc6,Jamaj,Jrms,Jell,JPA,
+ Jerrbits,Javgap3,Jclass,JclasStat,Jperrbits,Jseq,Jxi,Jeta,
+ HMJD,Hpmag,e_Hpmag,Hap3,e_Hap3,Hap4,e_Hap4,Hap6,e_Hap6,Hapc3,
+ Hapc4,Hapc6,Hamaj,Hrms,Hell,HPA,Herrbits,Havgap3,Hclass,
+ HclasSt,Hperrbits,Hseq,Hxi,Heta,KsMJD,Kspmag,e_Kspmag,Ksap3,
+ e_Ksap3,Ksap4,e_Ksap4,Ksap6,e_Ksap6,Ksapc3,Ksapc4,Ksapc6,
+ Ksamaj,Ksrms,Ksell,KsPA,Kserrbits,Ksavgap3,Ksclass,KsclasSt,
+ Ksperrbits,Ksseq,Ksxi,Kseta,prim
if(ar__(115:122) .EQ. '') Y_Jpnt = rNULL__
if(ar__(124:132) .EQ. '') e_Y_Jpnt = rNULL__
if(ar__(134:141) .EQ. '') J_Hpnt = rNULL__
if(ar__(143:150) .EQ. '') e_J_Hpnt = rNULL__
if(ar__(152:159) .EQ. '') H_Kspnt = rNULL__
if(ar__(161:168) .EQ. '') e_H_Kspnt = rNULL__
if(ar__(170:177) .EQ. '') J_Kspnt = rNULL__
if(ar__(179:188) .EQ. '') e_J_Kspnt = rNULL__
if(ar__(190:197) .EQ. '') Y_Jext = rNULL__
if(ar__(199:207) .EQ. '') e_Y_Jext = rNULL__
if(ar__(209:216) .EQ. '') J_Hext = rNULL__
if(ar__(218:225) .EQ. '') e_J_Hext = rNULL__
if(ar__(227:234) .EQ. '') H_Ksext = rNULL__
if(ar__(236:243) .EQ. '') e_H_Ksext = rNULL__
if(ar__(245:252) .EQ. '') J_Ksext = rNULL__
if(ar__(254:263) .EQ. '') e_J_Ksext = rNULL__
if(ar__(345:356) .EQ. '') YMJD = rNULL__
if(ar__(358:364) .EQ. '') Ypmag = rNULL__
if(ar__(366:372) .EQ. '') e_Ypmag = rNULL__
if(ar__(374:380) .EQ. '') Yap3 = rNULL__
if(ar__(382:388) .EQ. '') e_Yap3 = rNULL__
if(ar__(390:396) .EQ. '') Yap4 = rNULL__
if(ar__(398:404) .EQ. '') e_Yap4 = rNULL__
if(ar__(406:412) .EQ. '') Yap6 = rNULL__
if(ar__(414:420) .EQ. '') e_Yap6 = rNULL__
if(ar__(422:428) .EQ. '') Yapc3 = rNULL__
if(ar__(430:436) .EQ. '') Yapc4 = rNULL__
if(ar__(438:444) .EQ. '') Yapc6 = rNULL__
if(ar__(446:452) .EQ. '') Yamaj = rNULL__
if(ar__(454:461) .EQ. '') Yrms = rNULL__
if(ar__(463:468) .EQ. '') Yell = rNULL__
if(ar__(470:474) .EQ. '') YPA = rNULL__
if(ar__(479:485) .EQ. '') Yavgap3 = rNULL__
if(ar__(490:497) .EQ. '') YclasSt = rNULL__
if(ar__(508:514) .EQ. '') Yseq = iNULL__
if(ar__(516:527) .EQ. '') JMJD = rNULL__
if(ar__(529:536) .EQ. '') Jpmag = rNULL__
if(ar__(538:544) .EQ. '') e_Jpmag = rNULL__
if(ar__(546:553) .EQ. '') Jap3 = rNULL__
if(ar__(555:561) .EQ. '') e_Jap3 = rNULL__
if(ar__(563:570) .EQ. '') Jap4 = rNULL__
if(ar__(572:578) .EQ. '') e_Jap4 = rNULL__
if(ar__(580:587) .EQ. '') Jap6 = rNULL__
if(ar__(589:595) .EQ. '') e_Jap6 = rNULL__
if(ar__(597:604) .EQ. '') Japc3 = rNULL__
if(ar__(606:613) .EQ. '') Japc4 = rNULL__
if(ar__(615:622) .EQ. '') Japc6 = rNULL__
if(ar__(624:630) .EQ. '') Jamaj = rNULL__
if(ar__(632:639) .EQ. '') Jrms = rNULL__
if(ar__(641:646) .EQ. '') Jell = rNULL__
if(ar__(648:652) .EQ. '') JPA = rNULL__
if(ar__(654:655) .EQ. '') Jerrbits = iNULL__
if(ar__(657:663) .EQ. '') Javgap3 = rNULL__
if(ar__(665:666) .EQ. '') Jclass = iNULL__
if(ar__(668:675) .EQ. '') JclasStat = rNULL__
if(ar__(677:684) .EQ. '') Jperrbits = iNULL__
if(ar__(686:692) .EQ. '') Jseq = iNULL__
if(ar__(694:702) .EQ. '') Jxi = rNULL__
if(ar__(704:712) .EQ. '') Jeta = rNULL__
if(ar__(714:725) .EQ. '') HMJD = rNULL__
if(ar__(727:734) .EQ. '') Hpmag = rNULL__
if(ar__(736:742) .EQ. '') e_Hpmag = rNULL__
if(ar__(744:751) .EQ. '') Hap3 = rNULL__
if(ar__(753:759) .EQ. '') e_Hap3 = rNULL__
if(ar__(761:768) .EQ. '') Hap4 = rNULL__
if(ar__(770:776) .EQ. '') e_Hap4 = rNULL__
if(ar__(778:785) .EQ. '') Hap6 = rNULL__
if(ar__(787:793) .EQ. '') e_Hap6 = rNULL__
if(ar__(795:802) .EQ. '') Hapc3 = rNULL__
if(ar__(804:811) .EQ. '') Hapc4 = rNULL__
if(ar__(813:820) .EQ. '') Hapc6 = rNULL__
if(ar__(822:828) .EQ. '') Hamaj = rNULL__
if(ar__(830:837) .EQ. '') Hrms = rNULL__
if(ar__(839:844) .EQ. '') Hell = rNULL__
if(ar__(846:850) .EQ. '') HPA = rNULL__
if(ar__(852:853) .EQ. '') Herrbits = iNULL__
if(ar__(855:861) .EQ. '') Havgap3 = rNULL__
if(ar__(863:864) .EQ. '') Hclass = iNULL__
if(ar__(866:872) .EQ. '') HclasSt = rNULL__
if(ar__(874:881) .EQ. '') Hperrbits = iNULL__
if(ar__(883:888) .EQ. '') Hseq = iNULL__
if(ar__(890:898) .EQ. '') Hxi = rNULL__
if(ar__(900:908) .EQ. '') Heta = rNULL__
if(ar__(910:921) .EQ. '') KsMJD = rNULL__
if(ar__(923:930) .EQ. '') Kspmag = rNULL__
if(ar__(932:938) .EQ. '') e_Kspmag = rNULL__
if(ar__(940:947) .EQ. '') Ksap3 = rNULL__
if(ar__(949:955) .EQ. '') e_Ksap3 = rNULL__
if(ar__(957:964) .EQ. '') Ksap4 = rNULL__
if(ar__(966:972) .EQ. '') e_Ksap4 = rNULL__
if(ar__(974:981) .EQ. '') Ksap6 = rNULL__
if(ar__(983:989) .EQ. '') e_Ksap6 = rNULL__
if(ar__(991:998) .EQ. '') Ksapc3 = rNULL__
if(ar__(1000:1007) .EQ. '') Ksapc4 = rNULL__
if(ar__(1009:1016) .EQ. '') Ksapc6 = rNULL__
if(ar__(1018:1024) .EQ. '') Ksamaj = rNULL__
if(ar__(1026:1033) .EQ. '') Ksrms = rNULL__
if(ar__(1035:1040) .EQ. '') Ksell = rNULL__
if(ar__(1042:1046) .EQ. '') KsPA = rNULL__
if(ar__(1048:1049) .EQ. '') Kserrbits = iNULL__
if(ar__(1051:1057) .EQ. '') Ksavgap3 = rNULL__
if(ar__(1059:1060) .EQ. '') Ksclass = iNULL__
if(ar__(1062:1069) .EQ. '') KsclasSt = rNULL__
if(ar__(1071:1078) .EQ. '') Ksperrbits = iNULL__
if(ar__(1080:1086) .EQ. '') Ksseq = iNULL__
if(ar__(1088:1096) .EQ. '') Ksxi = rNULL__
if(ar__(1098:1106) .EQ. '') Kseta = rNULL__
if(ar__(1108:1110) .EQ. '') prim = rNULL__
c ..............Just test output...........
write(6,1)
+ SrcID,RAdeg,DEdeg,Name,cuEventID,FrameID,GLON,GLAT,PriOrSec,
+ Y_Jpnt,e_Y_Jpnt,J_Hpnt,e_J_Hpnt,H_Kspnt,e_H_Kspnt,J_Kspnt,
+ e_J_Kspnt,Y_Jext,e_Y_Jext,J_Hext,e_J_Hext,H_Ksext,e_H_Ksext,
+ J_Ksext,e_J_Ksext,Mstat,Mclass,pStar,pGal,pNoise,E_B_V,AY,AJ,
+ AH,AKs,YMJD,Ypmag,e_Ypmag,Yap3,e_Yap3,Yap4,e_Yap4,Yap6,e_Yap6,
+ Yapc3,Yapc4,Yapc6,Yamaj,Yrms,Yell,YPA,Yerrbits,Yavgap3,Yclass,
+ YclasSt,Yperrbits,Yseq,JMJD,Jpmag,e_Jpmag,Jap3,e_Jap3,Jap4,
+ e_Jap4,Jap6,e_Jap6,Japc3,Japc4,Japc6,Jamaj,Jrms,Jell,JPA,
+ Jerrbits,Javgap3,Jclass,JclasStat,Jperrbits,Jseq,Jxi,Jeta,
+ HMJD,Hpmag,e_Hpmag,Hap3,e_Hap3,Hap4,e_Hap4,Hap6,e_Hap6,Hapc3,
+ Hapc4,Hapc6,Hamaj,Hrms,Hell,HPA,Herrbits,Havgap3,Hclass,
+ HclasSt,Hperrbits,Hseq,Hxi,Heta,KsMJD,Kspmag,e_Kspmag,Ksap3,
+ e_Ksap3,Ksap4,e_Ksap4,Ksap6,e_Ksap6,Ksapc3,Ksapc4,Ksapc6,
+ Ksamaj,Ksrms,Ksell,KsPA,Kserrbits,Ksavgap3,Ksclass,KsclasSt,
+ Ksperrbits,Ksseq,Ksxi,Kseta,prim
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end