Conversion of standardized ReadMe file for
file /./ftp/cats/II/367 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-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 (nr__) ! Unique ID of this merged detection as
* assigned by merge algorithm
* unless PrioOrSec!=0 (SOURCEID)
real*8 RAdeg (nr__) ! (deg) Celestial Right Ascension (J2000)
* (RA2000)
real*8 DEdeg (nr__) ! (deg) [-90/1.5] Celestial Declination (J2000)
* (DEC2000)
character*23 Name (nr__) ! Source name in IAU convention
* (VHS JHHMMSS.ss+DDMMSS.s) (SOURCENAME)
integer*4 cuEventID (nr__) ! [587931] UID of curation event giving
* rise to this record (always "587931")
integer*8 FrameID (nr__) ! UID 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 (5)
real*8 Y_Jpnt (nr__) ! (mag) [-12.1/10.3]? Point source color Y-J
* (using aperMag3) (YMJPNT)
real*8 e_Y_Jpnt (nr__) ! (mag) [4E-4/5999]? Error on point source
* color Y-J (YMJPNTERR)
real*8 J_Hpnt (nr__) ! (mag) [-11.3/9.7]? Point source color J-H
* (using aperMag3) (JMHPNT)
real*8 e_J_Hpnt (nr__) ! (mag) [5E-4/744]? Error on point source
* color J-H (JMHPNTERR)
real*8 H_Kspnt (nr__) ! (mag) [-10.7/11.3]? Point source color H-Ks
* (using aperMag3) (HMKSPNT)
real*8 e_H_Kspnt (nr__) ! (mag) [5E-4/627]? Error on point source
* color H-Ks (HMKSPNTERR)
real*8 J_Kspnt (nr__) ! (mag) [-12/21.4]? Point source color J-Ks
* (using aperMag3) (JMKSPNT)
real*8 e_J_Kspnt (nr__) ! (mag) [5E-4/13621]? Error on point source
* color J-Ks (JMKSPNTERR)
real*8 Y_Jext (nr__) ! (mag) [-12/10.3]? Extended source color Y-J
* (using aperMagNoAperCorr3) (YMJEXT)
real*8 e_Y_Jext (nr__) ! (mag) [4E-4/5999]? Error on extended source
* color Y-J (YMJEXTERR)
real*8 J_Hext (nr__) ! (mag) [-11.2/9.8]? Extended source color J-H
* (using aperMagNoAperCorr3) (JMHEXT)
real*8 e_J_Hext (nr__) ! (mag) [5E-4/744]? Error on extended source
* color J-H (JMHEXTERR)
real*8 H_Ksext (nr__) ! (mag) [-10.6/11.2]? Extended source color
* H-Ks (using aperMagNoAperCorr3)
* (HMKSEXT)
real*8 e_H_Ksext (nr__) ! (mag) [5E-4/627]? Error on extended source
* color H-Ks (HMKSEXTERR)
real*8 J_Ksext (nr__) ! (mag) [-11.8/22.1]? Extended source color
* J-Ks (using aperMagNoAperCorr3)
* (JMKSEXT)
real*8 e_J_Ksext (nr__) ! (mag) [5E-4/13621]? Error on extended source
* color J-Ks (JMKSEXTERR)
real*8 Mstat (nr__) ! [-373/538] Merged N(0,1)
* stellarness-of-profile statistic
* (MERGEDCLASSSTAT) (1)
integer*4 Mclass (nr__) ! [-3/1] Class flag
* (MERGEDCLASS) (2)
real*8 pStar (nr__) ! [0/1] Probability that the source is a
* star
real*8 pGal (nr__) ! [0/1] Probability that the source is a
* galaxy (PGALAXY)
real*8 pNoise (nr__) ! [0/1] Probability that the source is
* noise
real*8 E_B_V (nr__) ! (mag) [4E-4/41] The galactic dust extinction
* value measured from the Schlegel
* (1998ApJ...500..525S) maps (EBV)
real*8 AY (nr__) ! (mag) [5E-4/53] The galactic extinction
* correction in the Y band
real*8 AJ (nr__) ! (mag) [4E-4/38] The galactic extinction
* correction in the J band
real*8 AH (nr__) ! (mag) [2E-4/24] The galactic extinction
* correction in the H band
real*8 AKs (nr__) ! (mag) [1E-4/16] The galactic extinction
* correction in the Ks band
real*8 YMJD (nr__) ! (d) [55160.2/57811.2]? Modified Julian Day
* in Y band
real*8 Ypmag (nr__) ! (mag) [3.6/37]? Extended source Y (Vega) mag
* (Petrosian) (YPETROMAG) (3)
real*8 e_Ypmag (nr__) ! (mag) [0/100]? Error in extended source Y mag
* (Petrosian) (YPETROMAGERR) (3)
real*8 Yap3 (nr__) ! (mag) [4.2/34]? Default point source Y
* aperture corrected Vega mag (2.0" diam)
* (YAPERMAG3) (3)
real*8 e_Yap3 (nr__) ! (mag) [1E-4/100]? Error in default
* point/extended source Y mag
* (YAPERMAG3ERR) (3)
real*8 Yap4 (nr__) ! (mag) [4.1/34]? Point source Y aperture
* corrected Vega mag (2.8" aperture diam)
* (YAPERMAG4) (3)
real*8 e_Yap4 (nr__) ! (mag) [1E-4/100]? Error in point/extended
* source Y mag (2.8 arcsec aperture diam)
* (YAPERMAG4ERR) (3)
real*8 Yap6 (nr__) ! (mag) [3.1/38]? Point source Y aperture
* corrected Vega mag (5.7" aperture diam)
* (YAPERMAG6) (3)
real*8 e_Yap6 (nr__) ! (mag) [1E-4/100]? Error in point/extended
* source Y mag (5.7 arcsec aperture diam)
* (YAPERMAG6ERR) (3)
real*8 Yapc3 (nr__) ! [4.6/35]? Default extended source Y
* aperture Vega mag (2.0 arcsec diam)
* (YAPERMAGNOAPERCORR3)
real*8 Yapc4 (nr__) ! (mag) [4.4/35]? Extended source Y aperture
* Vega mag (2.8 arcsec aperture diameter)
* (YAPERMAGNOAPERCORR4)
real*8 Yapc6 (nr__) ! (mag) [3.2/39]? Extended source Y aperture
* Vega mag (5.7 arcsec aperture diameter)
* (YAPERMAGNOAPERCORR6)
real*8 Yamaj (nr__) ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in Y band
* (YHLCORSMJRADAS)
real*8 Yrms (nr__) ! (pix) [0.7/530]? RMS of axes of ellipse fit
* in Y (YGAUSIG)
real*4 Yell (nr__) ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in Y
real*4 YPA (nr__) ! (deg) [0/180]? Ellipse fit celestial
* orientation in Y
integer*4 Yerrbits (nr__) ! [-1/28] Processing warning/error
* bitwise flags in Y
real*8 Yavgap3 (nr__) ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 Y (YAVERAGECONF)
integer*4 Yclass (nr__) ! [-2/1] Discrete image classification
* flag in Y
real*8 YclasSt (nr__) ! [-156/444]? N(0,1)
* stellarness-of-profile statistic in Y
* (YCLASSSTAT)
integer*4 Yperrbits (nr__) ! [0/12587216] Additional WFAU
* post-processing error bits in Y (4)
integer*4 Yseq (nr__) ! [0/1413053]? The running number of the
* Y detection (YSEQNUM)
real*8 JMJD (nr__) ! (d) [55139/57844.2]? Modified Julian Day in
* J band
real*8 Jpmag (nr__) ! (mag) [-24/38]? Extended source J (Vega) mag
* (Petrosian) (JPETROMAG) (3)
real*8 e_Jpmag (nr__) ! (mag) [0/100]? Error in extended source J mag
* (Petrosian) (JPETROMAGERR) (3)
real*8 Jap3 (nr__) ! (mag) [-21.1/34.1]? Default point source J
* aperture corrected Vega mag (2.0" diam)
* (JAPERMAG3) (3)
real*8 e_Jap3 (nr__) ! (mag) [0/100]? Error in default
* point/extended source J mag
* (JAPERMAG3ERR) (3)
real*8 Jap4 (nr__) ! (mag) [-21.2/39]? Point source J aperture
* corrected Vega mag (2.8" aperture diam)
* (JAPERMAG4) (3)
real*8 e_Jap4 (nr__) ! (mag) [0/100]? Error in point/extended source
* J mag (2.8 arcsec aperture diam)
* (JAPERMAG4ERR) (3)
real*8 Jap6 (nr__) ! (mag) [-21.7/39]? Point source J aperture
* corrected Vega mag (5.7" aperture diam)
* (JAPERMAG6) (3)
real*8 e_Jap6 (nr__) ! (mag) [0/100]? Error in point/extended source
* J mag (5.7 arcsec aperture diam)
* (JAPERMAG6ERR) (3)
real*8 Japc3 (nr__) ! (mag) [-20.5/35]? Default extended source J
* aperture Vega mag (2.0 arcsec diam)
* (JAPERMAGNOAPERCORR3)
real*8 Japc4 (nr__) ! (mag) [-21/39]? Extended source J aperture
* Vega mag (2.8 arcsec aperture diameter)
* (JAPERMAGNOAPERCORR4)
real*8 Japc6 (nr__) ! (mag) [-22/39]? Extended source J aperture
* Vega mag (5.7 arcsec aperture diameter)
* (JAPERMAGNOAPERCORR6)
real*8 Jamaj (nr__) ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in J band
* (JHLCORSMJRADAS)
real*8 Jrms (nr__) ! (pix) [0.7/584]? RMS of axes of ellipse fit
* in J (JGAUSIG)
real*4 Jell (nr__) ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in J
real*4 JPA (nr__) ! (deg) [0/180]? Ellipse fit celestial
* orientation in J
integer*4 Jerrbits (nr__) ! [-1/28]? Processing warning/error
* bitwise flags in J
real*8 Javgap3 (nr__) ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 J (JAVERAGECONF)
integer*4 Jclass (nr__) ! [-2/1]? Discrete image classification
* flag in J
real*8 JclasStat (nr__) ! [-373/458]? N(0,1)
* stellarness-of-profile statistic in J
* (JCLASSSTAT)
integer*4 Jperrbits (nr__) ! [0/12587216]? Additional WFAU
* post-processing error bits in J (4)
integer*4 Jseq (nr__) ! [1/1661208]? The running number of the
* J detection (JSEQNUM)
real*8 Jxi (nr__) ! (arcsec) [-2/2]? Offset of J detection from
* master position (+east/-west)
real*8 Jeta (nr__) ! (arcsec) [-2/2]? Offset of J detection from
* master position (+north/-south)
real*8 HMJD (nr__) ! (d) [55139/56074.2]? Modified Julian Day in
* H band
real*8 Hpmag (nr__) ! (mag) [-24/36]? Extended source H (Vega) mag
* (Petrosian) (HPETROMAG) (3)
real*8 e_Hpmag (nr__) ! (mag) [0/100]? Error in extended source H mag
* (Petrosian) (HPETROMAGERR) (3)
real*8 Hap3 (nr__) ! (mag) [-21.5/29]? Default point source H
* aperture corrected Vega mag (2.0" diam)
* (HAPERMAG3) (3)
real*8 e_Hap3 (nr__) ! (mag) [0/100]? Error in default
* point/extended source H mag
* (HAPERMAG3ERR) (3)
real*8 Hap4 (nr__) ! (mag) [-22/30.1]? Point source H aperture
* corrected Vega mag (2.8" aperture diam)
* (HAPERMAG4) (3)
real*8 e_Hap4 (nr__) ! (mag) [0/100]? Error in point/extended source
* H mag (2.8 arcsec aperture diam)
* (HAPERMAG4ERR) (3)
real*8 Hap6 (nr__) ! (mag) [-22.2/35]? Point source H aperture
* corrected Vega mag (5.7" aperture diam)
* (HAPERMAG6) (3)
real*8 e_Hap6 (nr__) ! (mag) [0/100]? Error in point/extended source
* H mag (5.7 arcsec aperture diam)
* (HAPERMAG6ERR) (3)
real*8 Hapc3 (nr__) ! (mag) [-21/29.1]? Default extended source H
* aperture Vega mag (2.0 arcsec diam)
* (HAPERMAGNOAPERCORR3)
real*8 Hapc4 (nr__) ! (mag) [-21.4/30.3]? Extended source H
* aperture Vega mag
* (2.8 arcsec aperture diameter)
* (HAPERMAGNOAPERCORR4)
real*8 Hapc6 (nr__) ! (mag) [-22.2/35]? Extended source H aperture
* Vega mag (5.7 arcsec aperture diameter)
* (HAPERMAGNOAPERCORR6)
real*8 Hamaj (nr__) ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in H band
* (HHLCORSMJRADAS)
real*8 Hrms (nr__) ! (pix) [0.7/584]? RMS of axes of ellipse fit
* in H (HGAUSIG)
real*4 Hell (nr__) ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in H
real*4 HPA (nr__) ! (deg) [0/180]? Ellipse fit celestial
* orientation in H
integer*4 Herrbits (nr__) ! [-1/28]? Processing warning/error
* bitwise flags in H
real*8 Havgap3 (nr__) ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 H (HAVERAGECONF)
integer*4 Hclass (nr__) ! [-2/1]? Discrete image classification
* flag in H
real*8 HclasSt (nr__) ! [-96/383]? N(0,1)
* stellarness-of-profile statistic in H
* (HCLASSSTAT)
integer*4 Hperrbits (nr__) ! [0/12587216]? Additional WFAU
* post-processing error bits in H (4)
integer*4 Hseq (nr__) ! [1/125514]? The running number of the H
* detection (HSEQNUM)
real*8 Hxi (nr__) ! (arcsec) [-2/2]? Offset of H detection from
* master position (+east/-west)
real*8 Heta (nr__) ! (arcsec) [-2/2]? Offset of H detection from
* master position (+north/-south)
real*8 KsMJD (nr__) ! (d) [55139/57844.2]? Modified Julian Day in
* Ks band
real*8 Kspmag (nr__) ! (mag) [-20.3/34]? Extended source Ks (Vega) mag
* (Petrosian) (KSPETROMAG) (3)
real*8 e_Kspmag (nr__) ! (mag) [0/100]? Error in extended source Ks
* mag (Petrosian) (KSPETROMAGERR) (3)
real*8 Ksap3 (nr__) ! (mag) [-33.5/31]? Default point source Ks
* aperture corrected Vega mag (2.0" diam)
* (KSAPERMAG3) (3)
real*8 e_Ksap3 (nr__) ! (mag) [0/96]? Error in default point/extended
* source Ks mag (KSAPERMAG3ERR) (3)
real*8 Ksap4 (nr__) ! (mag) [-33.4/36]? Point source Ks aperture
* corrected Vega mag (2.8" aperture diam)
* (KSAPERMAG4) (3)
real*8 e_Ksap4 (nr__) ! (mag) [0/100]? Error in point/extended source
* Ks mag (2.8 arcsec aperture diam)
* (KSAPERMAG4ERR) (3)
real*8 Ksap6 (nr__) ! (mag) [-33.3/37]? Point source Ks aperture
* corrected Vega mag (5.7" aperture diam)
* (KSAPERMAG6) (3)
real*8 e_Ksap6 (nr__) ! (mag) [0/100]? Error in point/extended source
* Ks mag (5.7 arcsec aperture diam)
* (KSAPERMAG6ERR) (3)
real*8 Ksapc3 (nr__) ! (mag) [-33.4/32]? Default extended source Ks
* aperture Vega mag (2.0 arcsec diam)
* (KSAPERMAGNOAPERCORR3)
real*8 Ksapc4 (nr__) ! (mag) [-33.3/36]? Extended source Ks aperture
* Vega mag (2.8 arcsec aperture diameter)
* (KSAPERMAGNOAPERCORR4)
real*8 Ksapc6 (nr__) ! (mag) [-33.3/37.1]? Extended source Ks
* aperture Vega mag
* (5.7 arcsec aperture diameter)
* (KSAPERMAGNOAPERCORR6)
real*8 Ksamaj (nr__) ! (arcsec) [0/19]? Seeing corrected half-light,
* semi-major axis in Ks band
* (KSHLCORSMJRADAS)
real*8 Ksrms (nr__) ! (pix) [0.7/549]? RMS of axes of ellipse fit
* in Ks (KSGAUSIG)
real*4 Ksell (nr__) ! [0/1]? 1-b/a, where
* a/b=semi-major/minor axes in Ks
real*4 KsPA (nr__) ! (deg) [0/180]? Ellipse fit celestial
* orientation in Ks
integer*4 Kserrbits (nr__) ! [-1/28]? Processing warning/error
* bitwise flags in Ks
real*8 Ksavgap3 (nr__) ! [0/112]? Average confidence in 2 arcsec
* diam default aperture 3 Ks
* (KSAVERAGECONF)
integer*4 Ksclass (nr__) ! [-2/1]? Discrete image classification
* flag in Ks
real*8 KsclasSt (nr__) ! [-108/385]? N(0,1)
* stellarness-of-profile statistic in Ks
* (KSCLASSSTAT)
integer*4 Ksperrbits (nr__) ! [0/12587216]? Additional WFAU
* post-processing error bits in Ks (4)
integer*4 Ksseq (nr__) ! [1/1455550]? The running number of the
* Ks detection (KSSEQNUM)
real*8 Ksxi (nr__) ! (arcsec) [-2/2]? Offset of Ks detection from
* master position (+east/-west)
real*8 Kseta (nr__) ! (arcsec) [-2/2]? Offset of Ks detection from
* master position (+north/-south)
real*4 prim (nr__) ! [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(i__),RAdeg(i__),DEdeg(i__),Name(i__),cuEventID(i__),
+ FrameID(i__),GLON(i__),GLAT(i__),PriOrSec(i__),Y_Jpnt(i__),
+ e_Y_Jpnt(i__),J_Hpnt(i__),e_J_Hpnt(i__),H_Kspnt(i__),
+ e_H_Kspnt(i__),J_Kspnt(i__),e_J_Kspnt(i__),Y_Jext(i__),
+ e_Y_Jext(i__),J_Hext(i__),e_J_Hext(i__),H_Ksext(i__),
+ e_H_Ksext(i__),J_Ksext(i__),e_J_Ksext(i__),Mstat(i__),
+ Mclass(i__),pStar(i__),pGal(i__),pNoise(i__),E_B_V(i__),
+ AY(i__),AJ(i__),AH(i__),AKs(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__),
+ Yamaj(i__),Yrms(i__),Yell(i__),YPA(i__),Yerrbits(i__),
+ Yavgap3(i__),Yclass(i__),YclasSt(i__),Yperrbits(i__),
+ Yseq(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__),Jamaj(i__),Jrms(i__),
+ Jell(i__),JPA(i__),Jerrbits(i__),Javgap3(i__),Jclass(i__),
+ JclasStat(i__),Jperrbits(i__),Jseq(i__),Jxi(i__),Jeta(i__),
+ HMJD(i__),Hpmag(i__),e_Hpmag(i__),Hap3(i__),e_Hap3(i__),
+ Hap4(i__),e_Hap4(i__),Hap6(i__),e_Hap6(i__),Hapc3(i__),
+ Hapc4(i__),Hapc6(i__),Hamaj(i__),Hrms(i__),Hell(i__),HPA(i__),
+ Herrbits(i__),Havgap3(i__),Hclass(i__),HclasSt(i__),
+ Hperrbits(i__),Hseq(i__),Hxi(i__),Heta(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__),Ksamaj(i__),Ksrms(i__),Ksell(i__),KsPA(i__),
+ Kserrbits(i__),Ksavgap3(i__),Ksclass(i__),KsclasSt(i__),
+ Ksperrbits(i__),Ksseq(i__),Ksxi(i__),Kseta(i__),prim(i__)
if(ar__(115:122) .EQ. '') Y_Jpnt(i__) = rNULL__
if(ar__(124:132) .EQ. '') e_Y_Jpnt(i__) = rNULL__
if(ar__(134:141) .EQ. '') J_Hpnt(i__) = rNULL__
if(ar__(143:150) .EQ. '') e_J_Hpnt(i__) = rNULL__
if(ar__(152:159) .EQ. '') H_Kspnt(i__) = rNULL__
if(ar__(161:168) .EQ. '') e_H_Kspnt(i__) = rNULL__
if(ar__(170:177) .EQ. '') J_Kspnt(i__) = rNULL__
if(ar__(179:188) .EQ. '') e_J_Kspnt(i__) = rNULL__
if(ar__(190:197) .EQ. '') Y_Jext(i__) = rNULL__
if(ar__(199:207) .EQ. '') e_Y_Jext(i__) = rNULL__
if(ar__(209:216) .EQ. '') J_Hext(i__) = rNULL__
if(ar__(218:225) .EQ. '') e_J_Hext(i__) = rNULL__
if(ar__(227:234) .EQ. '') H_Ksext(i__) = rNULL__
if(ar__(236:243) .EQ. '') e_H_Ksext(i__) = rNULL__
if(ar__(245:252) .EQ. '') J_Ksext(i__) = rNULL__
if(ar__(254:263) .EQ. '') e_J_Ksext(i__) = rNULL__
if(ar__(345:356) .EQ. '') YMJD(i__) = rNULL__
if(ar__(358:364) .EQ. '') Ypmag(i__) = rNULL__
if(ar__(366:372) .EQ. '') e_Ypmag(i__) = rNULL__
if(ar__(374:380) .EQ. '') Yap3(i__) = rNULL__
if(ar__(382:388) .EQ. '') e_Yap3(i__) = rNULL__
if(ar__(390:396) .EQ. '') Yap4(i__) = rNULL__
if(ar__(398:404) .EQ. '') e_Yap4(i__) = rNULL__
if(ar__(406:412) .EQ. '') Yap6(i__) = rNULL__
if(ar__(414:420) .EQ. '') e_Yap6(i__) = rNULL__
if(ar__(422:428) .EQ. '') Yapc3(i__) = rNULL__
if(ar__(430:436) .EQ. '') Yapc4(i__) = rNULL__
if(ar__(438:444) .EQ. '') Yapc6(i__) = rNULL__
if(ar__(446:452) .EQ. '') Yamaj(i__) = rNULL__
if(ar__(454:461) .EQ. '') Yrms(i__) = rNULL__
if(ar__(463:468) .EQ. '') Yell(i__) = rNULL__
if(ar__(470:474) .EQ. '') YPA(i__) = rNULL__
if(ar__(479:485) .EQ. '') Yavgap3(i__) = rNULL__
if(ar__(490:497) .EQ. '') YclasSt(i__) = rNULL__
if(ar__(508:514) .EQ. '') Yseq(i__) = iNULL__
if(ar__(516:527) .EQ. '') JMJD(i__) = rNULL__
if(ar__(529:536) .EQ. '') Jpmag(i__) = rNULL__
if(ar__(538:544) .EQ. '') e_Jpmag(i__) = rNULL__
if(ar__(546:553) .EQ. '') Jap3(i__) = rNULL__
if(ar__(555:561) .EQ. '') e_Jap3(i__) = rNULL__
if(ar__(563:570) .EQ. '') Jap4(i__) = rNULL__
if(ar__(572:578) .EQ. '') e_Jap4(i__) = rNULL__
if(ar__(580:587) .EQ. '') Jap6(i__) = rNULL__
if(ar__(589:595) .EQ. '') e_Jap6(i__) = rNULL__
if(ar__(597:604) .EQ. '') Japc3(i__) = rNULL__
if(ar__(606:613) .EQ. '') Japc4(i__) = rNULL__
if(ar__(615:622) .EQ. '') Japc6(i__) = rNULL__
if(ar__(624:630) .EQ. '') Jamaj(i__) = rNULL__
if(ar__(632:639) .EQ. '') Jrms(i__) = rNULL__
if(ar__(641:646) .EQ. '') Jell(i__) = rNULL__
if(ar__(648:652) .EQ. '') JPA(i__) = rNULL__
if(ar__(654:655) .EQ. '') Jerrbits(i__) = iNULL__
if(ar__(657:663) .EQ. '') Javgap3(i__) = rNULL__
if(ar__(665:666) .EQ. '') Jclass(i__) = iNULL__
if(ar__(668:675) .EQ. '') JclasStat(i__) = rNULL__
if(ar__(677:684) .EQ. '') Jperrbits(i__) = iNULL__
if(ar__(686:692) .EQ. '') Jseq(i__) = iNULL__
if(ar__(694:702) .EQ. '') Jxi(i__) = rNULL__
if(ar__(704:712) .EQ. '') Jeta(i__) = rNULL__
if(ar__(714:725) .EQ. '') HMJD(i__) = rNULL__
if(ar__(727:734) .EQ. '') Hpmag(i__) = rNULL__
if(ar__(736:742) .EQ. '') e_Hpmag(i__) = rNULL__
if(ar__(744:751) .EQ. '') Hap3(i__) = rNULL__
if(ar__(753:759) .EQ. '') e_Hap3(i__) = rNULL__
if(ar__(761:768) .EQ. '') Hap4(i__) = rNULL__
if(ar__(770:776) .EQ. '') e_Hap4(i__) = rNULL__
if(ar__(778:785) .EQ. '') Hap6(i__) = rNULL__
if(ar__(787:793) .EQ. '') e_Hap6(i__) = rNULL__
if(ar__(795:802) .EQ. '') Hapc3(i__) = rNULL__
if(ar__(804:811) .EQ. '') Hapc4(i__) = rNULL__
if(ar__(813:820) .EQ. '') Hapc6(i__) = rNULL__
if(ar__(822:828) .EQ. '') Hamaj(i__) = rNULL__
if(ar__(830:837) .EQ. '') Hrms(i__) = rNULL__
if(ar__(839:844) .EQ. '') Hell(i__) = rNULL__
if(ar__(846:850) .EQ. '') HPA(i__) = rNULL__
if(ar__(852:853) .EQ. '') Herrbits(i__) = iNULL__
if(ar__(855:861) .EQ. '') Havgap3(i__) = rNULL__
if(ar__(863:864) .EQ. '') Hclass(i__) = iNULL__
if(ar__(866:872) .EQ. '') HclasSt(i__) = rNULL__
if(ar__(874:881) .EQ. '') Hperrbits(i__) = iNULL__
if(ar__(883:888) .EQ. '') Hseq(i__) = iNULL__
if(ar__(890:898) .EQ. '') Hxi(i__) = rNULL__
if(ar__(900:908) .EQ. '') Heta(i__) = rNULL__
if(ar__(910:921) .EQ. '') KsMJD(i__) = rNULL__
if(ar__(923:930) .EQ. '') Kspmag(i__) = rNULL__
if(ar__(932:938) .EQ. '') e_Kspmag(i__) = rNULL__
if(ar__(940:947) .EQ. '') Ksap3(i__) = rNULL__
if(ar__(949:955) .EQ. '') e_Ksap3(i__) = rNULL__
if(ar__(957:964) .EQ. '') Ksap4(i__) = rNULL__
if(ar__(966:972) .EQ. '') e_Ksap4(i__) = rNULL__
if(ar__(974:981) .EQ. '') Ksap6(i__) = rNULL__
if(ar__(983:989) .EQ. '') e_Ksap6(i__) = rNULL__
if(ar__(991:998) .EQ. '') Ksapc3(i__) = rNULL__
if(ar__(1000:1007) .EQ. '') Ksapc4(i__) = rNULL__
if(ar__(1009:1016) .EQ. '') Ksapc6(i__) = rNULL__
if(ar__(1018:1024) .EQ. '') Ksamaj(i__) = rNULL__
if(ar__(1026:1033) .EQ. '') Ksrms(i__) = rNULL__
if(ar__(1035:1040) .EQ. '') Ksell(i__) = rNULL__
if(ar__(1042:1046) .EQ. '') KsPA(i__) = rNULL__
if(ar__(1048:1049) .EQ. '') Kserrbits(i__) = iNULL__
if(ar__(1051:1057) .EQ. '') Ksavgap3(i__) = rNULL__
if(ar__(1059:1060) .EQ. '') Ksclass(i__) = iNULL__
if(ar__(1062:1069) .EQ. '') KsclasSt(i__) = rNULL__
if(ar__(1071:1078) .EQ. '') Ksperrbits(i__) = iNULL__
if(ar__(1080:1086) .EQ. '') Ksseq(i__) = iNULL__
if(ar__(1088:1096) .EQ. '') Ksxi(i__) = rNULL__
if(ar__(1098:1106) .EQ. '') Kseta(i__) = rNULL__
if(ar__(1108:1110) .EQ. '') prim(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ SrcID(i__),RAdeg(i__),DEdeg(i__),Name(i__),cuEventID(i__),
+ FrameID(i__),GLON(i__),GLAT(i__),PriOrSec(i__),Y_Jpnt(i__),
+ e_Y_Jpnt(i__),J_Hpnt(i__),e_J_Hpnt(i__),H_Kspnt(i__),
+ e_H_Kspnt(i__),J_Kspnt(i__),e_J_Kspnt(i__),Y_Jext(i__),
+ e_Y_Jext(i__),J_Hext(i__),e_J_Hext(i__),H_Ksext(i__),
+ e_H_Ksext(i__),J_Ksext(i__),e_J_Ksext(i__),Mstat(i__),
+ Mclass(i__),pStar(i__),pGal(i__),pNoise(i__),E_B_V(i__),
+ AY(i__),AJ(i__),AH(i__),AKs(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__),
+ Yamaj(i__),Yrms(i__),Yell(i__),YPA(i__),Yerrbits(i__),
+ Yavgap3(i__),Yclass(i__),YclasSt(i__),Yperrbits(i__),
+ Yseq(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__),Jamaj(i__),Jrms(i__),
+ Jell(i__),JPA(i__),Jerrbits(i__),Javgap3(i__),Jclass(i__),
+ JclasStat(i__),Jperrbits(i__),Jseq(i__),Jxi(i__),Jeta(i__),
+ HMJD(i__),Hpmag(i__),e_Hpmag(i__),Hap3(i__),e_Hap3(i__),
+ Hap4(i__),e_Hap4(i__),Hap6(i__),e_Hap6(i__),Hapc3(i__),
+ Hapc4(i__),Hapc6(i__),Hamaj(i__),Hrms(i__),Hell(i__),HPA(i__),
+ Herrbits(i__),Havgap3(i__),Hclass(i__),HclasSt(i__),
+ Hperrbits(i__),Hseq(i__),Hxi(i__),Heta(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__),Ksamaj(i__),Ksrms(i__),Ksell(i__),KsPA(i__),
+ Kserrbits(i__),Ksavgap3(i__),Ksclass(i__),KsclasSt(i__),
+ Ksperrbits(i__),Ksseq(i__),Ksxi(i__),Kseta(i__),prim(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end