Conversion of standardized ReadMe file for
file /./ftp/cats/I/353 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-14
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. I/353 The Guide Star Catalog, Version 2.4.2 (GSC2.4.2) (STScI, 2020)
*================================================================================
*The Guide Star Catalogue, Version 2.4.2
* Lasker B., Lattanzi M.G., McLean B.J., Bucciarelli B., Drimmel R.,
* Garcia J., Greene G., Guglielmetti F., Hanley C., Hawkins G., Laidler V.G.,
* Loomis C., Meakes M., Mignani R., Morbidelli R., Morrison J., Pannunzio R.,
* Rosenberg A., Sarasso M., Smart R.L., Spagna A., Sturch C.R.,
* Volpicelli A., White R.L., Wolfe D., Zacchei A.
* <Astron. J. 136, 735 (2008) >
* =2008AJ....136..735L
* =2021yCat.1353....0L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gsc242.sam' ! Sample records of GSC 2.4.2
integer*4 nr__
parameter (nr__=1000) ! Number of records
character*945 ar__ ! Full-size record
C Position composed of: RAdeg DEdeg
character*10 GSC2 (nr__) ! GSC2 Public ID (hstID)
real*8 RAdeg (nr__) ! (deg) Right Ascension (ICRS) at Ep=Epoch (ra)
real*8 DEdeg (nr__) ! (deg) Declination (ICRS) at Ep=Epoch (dec)
integer*8 objID (nr__) ! Database ID
integer*4 GSC1 (nr__) ! ? GSC1 ID (gsc1ID)
character*15 GSC2intID (nr__) ! GSC2 Internal ID (gsc2ID)
integer*8 Compass (nr__) ! COMPASS ID (compassgsc2id)
integer*8 GaiaDR1 (nr__) ! ? Gaia DR1 ID (GAIAsourceID)
integer*4 usageFlag (nr__) ! [-1/1] Flag for previously used
* guide star
integer*8 Status (nr__) ! ? Object status codes (sourceStatus) (4)
real*8 Epoch (nr__) ! (yr) Epoch of coordinate
real*8 e_RAdeg (nr__) ! (arcsec) ? Right Ascension error (raErr)
real*8 e_DEdeg (nr__) ! (arcsec) ? Declination error (decErr)
integer*4 r_pos (nr__) ! [2/32] Object coordinates source
* (posSource) (2)
real*8 pmRA (nr__) ! (mas/yr) ? Proper motion along RA (rapm)
real*8 pmDE (nr__) ! (mas/yr) ? Proper motion along DE (decpm)
real*8 e_pmRA (nr__) ! (mas/yr) ? pmRA error (rapmErr)
real*8 e_pmDE (nr__) ! (mas/yr) ? pmDE error (decpmErr)
real*8 dEpoch (nr__) ! (yr) ? Time interval used for PM computation
* (deltaEpoch)
integer*4 r_pm (nr__) ! [0/15]? Proper motion source (pmSource) (2)
real*8 plx (nr__) ! (mas) ? Parallax (parallax)
real*8 e_plx (nr__) ! (mas) ? Parallax error (plxErr)
integer*4 r_plx (nr__) ! [0/15]? Parallax source (plxSource) (2)
real*8 RV (nr__) ! (km/s) ? Radial velocity (radialVelocity)
real*4 e_RV (nr__) ! (km/s) ? Radial velocity error (rvErr)
integer*4 r_RV (nr__) ! [0/15]? Radial Velocity source
* (rvSource) (2)
integer*4 Class (nr__) ! [0/4]? Object classification code
* (classification) (1)
integer*4 r_Class (nr__) ! [0/32]? Classification source
* (classSource) (2)
integer*4 varFlg (nr__) ! [1]? Flag for known variable (variableFlag)
integer*4 r_varFlg (nr__) ! [3/15]? Variability source (varSource) (2)
integer*4 multFlg (nr__) ! [0/1] Flag for known multiple multipleFlag
* (multipleFlag)
integer*4 r_multFlg (nr__) ! [0/3] Multiplicity source (multSource) (2)
real*8 Gmag (nr__) ! (mag) ? Gaia G band magnitude (gaiaGmag)
real*8 e_Gmag (nr__) ! (mag) ? Gaia G band magnitude error (gaiaGmagErr)
integer*4 f_Gmag (nr__) ! [65]?=0 Gaia G band code (gaiaGmagCode) (3)
real*8 OMag (nr__) ! (mag) ? Photographic O band magnitude (OpgMag)
real*8 e_OMag (nr__) ! (mag) ? Photographic O band magnitude error
* (OpgMagErr)
integer*4 f_OMag (nr__) ! [7]?=0 Photographic O band code
* (OpgMagCode) (3)
real*8 Bjmag (nr__) ! (mag) ? Photographic Blue magnitude (JpgMag)
real*8 e_Bjmag (nr__) ! (mag) ? Photographic Blue magnitude error
* (JpgMagErr)
integer*4 f_Bjmag (nr__) ! [0/18]?=0 Photographic Blue bandpass code
* (JpgMagCode) (3)
real*8 VpgMag (nr__) ! (mag) ? Photographic V band magnitude
real*8 e_VpgMag (nr__) ! (mag) ? Photographic V band magnitude error
* (VpgMagErr)
integer*4 f_VpgMag (nr__) ! [1/6]?=0 Photographic V band code
* (VpgMagCode) (3)
real*8 Fpgmag (nr__) ! (mag) ? Photographic Red magnitude
* (FpgMag)
real*8 e_Fpgmag (nr__) ! (mag) ? Photographic Red magnitude error
* (FpgMagErr)
integer*4 f_Fpgmag (nr__) ! [35/36]?=0 Photographic Red bandpass code
* (FpgMagCode) (3)
real*8 Epgmag (nr__) ! (mag) ? Photographic E band magnitude (EpgMag)
real*8 e_Epgmag (nr__) ! (mag) ? Photographic E band magnitude error
* (EpgMagErr)
integer*4 f_Epgmag (nr__) ! [8]?=0 Photographic E band code
* (EpgMagCode) (3)
real*8 Npgmag (nr__) ! (mag) ? Photographic Near-IR magnitude (NpgMag)
real*8 e_Npgmag (nr__) ! (mag) ? Photographic Near-IR magnitude error
* (NpgMagErr)
integer*4 f_Npgmag (nr__) ! [37/38]?=0 Photographic Near-IR bandpass code
* (NpgMagCode) (3)
real*8 Umag (nr__) ! (mag) ? U band magnitude
real*8 e_Umag (nr__) ! (mag) ? U band magnitude error (UmagErr)
integer*4 f_Umag (nr__) ! [19]?=0 U band bandpass code (UmagCode) (3)
integer*4 r_Umag (nr__) ! [3]?=0 U band magnitude source
* (UmagSource) (2)
real*8 Bmag (nr__) ! (mag) ? B band magnitude
real*8 e_Bmag (nr__) ! (mag) ? B band magnitude error (BmagErr)
integer*4 f_Bmag (nr__) ! [3/41]?=0 B band bandpass code (BmagCode) (3)
integer*4 r_Bmag (nr__) ! [3/19]?=0 B band magnitude source
* (BmagSource) (2)
real*8 Vmag (nr__) ! (mag) ? V band magnitude
real*8 e_Vmag (nr__) ! (mag) ? V band magnitude error (VmagErr)
integer*4 f_Vmag (nr__) ! [3/75]?=0 V band bandpass code (VMagCode) (3)
integer*4 r_Vmag (nr__) ! [3/19]?=0 V band magnitude source
* (VmagSource) (2)
real*4 umag_1 (nr__) ! (mag) ? SDSS u band magnitude (SDSSuMag)
real*4 e_umag_1 (nr__) ! (mag) ? SDSS u band magnitude error (SDSSuMagErr)
integer*4 f_umag_1 (nr__) ! [50/74]?=0 SDSS u band code
* (SDSSuMagCode) (3)
integer*4 r_umag_1 (nr__) ! [17/32]?=0 SDSS u band magnitude source
* (SDSSuMagSource) (2)
real*4 gmag_1 (nr__) ! (mag) ? SDSS g band magnitude (SDSSgMag)
real*4 e_gmag_1 (nr__) ! (mag) ? SDSS g band magnitude error (SDSSgMagErr)
integer*4 f_gmag_1 (nr__) ! [51/76]?=0 SDSS g band code
* (SDSSgMagCode) (3)
integer*4 r_gmag (nr__) ! [16/32]?=0 SDSS g band magnitude source
* (SDSSgMagSource) (2)
real*4 rmag (nr__) ! (mag) ? SDSS r band magnitude (SDSSrMag)
real*4 e_rmag (nr__) ! (mag) ? SDSS r band magnitude error (SDSSrMagErr)
integer*4 f_rmag (nr__) ! [52/77]?=0 SDSS r band code
* (SDSSrMagCode) (3)
integer*4 r_rmag (nr__) ! [16/32]?=0 SDSS r band magnitude source
* (SDSSrMagSource) (2)
real*4 imag (nr__) ! (mag) ? SDSS i band magnitude (SDSSiMag)
real*8 e_imag (nr__) ! (mag) ? SDSS i band magnitude error (SDSSiMagErr)
integer*4 f_imag (nr__) ! [53/78]?=0 SDSS i band code
* (SDSSiMagCode) (3)
integer*4 r_imag (nr__) ! [16/32]?=0 SDSS i band magnitude source
* (SDSSiMagSource) (2)
real*4 zmag (nr__) ! (mag) ? SDSS z band magnitude (SDSSzMag)
real*4 e_zmag (nr__) ! (mag) ? SDSS z band magnitude error (SDSSzMagErr)
integer*4 f_zmag (nr__) ! [54/79]?=0 SDSS z band code
* (SDSSzMagCode) (3)
integer*4 r_zmag (nr__) ! [16/31]?=0 SDSS z band magnitude source
* (SDSSzMagSource) (2)
real*8 ymag (nr__) ! (mag) ? PanSTARRS y band magnitude (PS1yMag)
real*8 e_ymag (nr__) ! (mag) ? PanSTARRS y band magnitude error
* (PS1yMagErr)
integer*4 f_ymag (nr__) ! [63]?=0 PanSTARRS y band code
* (PS1yMagCode) (3)
integer*4 r_ymag (nr__) ! [16]?=0 PanSTARRS y band magnitude source
* (PS1ymagSource) (2)
real*8 Jmag (nr__) ! (mag) ? 2MASS J band magnitude (tmassJmag)
real*4 e_Jmag (nr__) ! (mag) ? 2MASS J band magnitude error (tmassJmagErr)
integer*4 f_Jmag (nr__) ! [47]?=0 2MASS J band code (tmassJmagCode) (3)
integer*4 r_Jmag (nr__) ! [9/23]?=0 2MASS J band magnitude source
* (tmassJmagSource) (2)
real*8 Hmag (nr__) ! (mag) ? 2MASS H band magnitude (tmassHmag)
real*4 e_Hmag (nr__) ! (mag) ? 2MASS H band magnitude error (tmassHmagErr)
integer*4 f_Hmag (nr__) ! [48]?=0 2MASS H band code (tmassHmagCode) (3)
integer*4 r_Hmag (nr__) ! [9/23]?=0 2MASS H band magnitude source
* (tmassHmagSource) (2)
real*8 Ksmag (nr__) ! (mag) ? 2MASS Ks band magnitude (tmassKsmag)
real*4 e_Ksmag (nr__) ! (mag) ? 2MASS Ks band magnitude error
* (tmassKsmagErr)
integer*4 f_Ksmag (nr__) ! [49]?=0 2MASS Ks band code
* (tmassKsmagCode) (3)
integer*4 r_Ksmag (nr__) ! [9/23]?=0 2MASS Ks band magnitude source
* (tmassKsmagSource) (2)
real*4 Zmag_1 (nr__) ! (mag) ? VISTA Z band magnitude (VistaZMag)
real*4 e_Zmag_1 (nr__) ! (mag) ? VISTA Z band magnitude error (VistaZMagErr)
integer*4 f_Zmag_1 (nr__) ! [69]?=0 VISTA Z band code (VistaZMagCode) (3)
integer*4 r_Zmag_1 (nr__) ! [23]?=0 VISTA Z band magnitude source
* (VistaZmagSource) (2)
real*8 Ymag_1 (nr__) ! (mag) ? VISTA Y band magnitude (VistaYmag)
real*4 e_Ymag_1 (nr__) ! (mag) ? VISTA Y band magnitude error (VistaYMagErr)
integer*4 f_Ymag_1 (nr__) ! [70]?=0 VISTA Y band code (VistaYMagCode) (3)
integer*4 r_Ymag_1 (nr__) ! [21/23]?=0 VISTA Y band magnitude source
* (VistaYmagSource) (2)
real*4 W1mag (nr__) ! (mag) ? WISE W1 (3.4um) band magnitude (WISEw1Mag)
real*4 e_W1mag (nr__) ! (mag) ? WISE W1 (3.4um) band magnitude error
* (WISEw1MagErr)
integer*4 f_W1mag (nr__) ! [55]?=0 WISE W1 (3.4um) band code
* (WISEw1MagCode) (3)
integer*4 r_W1mag (nr__) ! [10]?=0 WISE W1 (3.4um) band magnitude source
* (WISEw1MagSOurce) (2)
real*4 W2mag (nr__) ! (mag) ? WISE W2 (4.6um) band magnitude (WISEw2Mag)
real*4 e_W2mag (nr__) ! (mag) ? WISE W2 (4.6um) band magnitude error
* (WISEw2MagErr)
integer*4 f_W2mag (nr__) ! [56]?=0 WISE W2 (4.6um) band code
* (WISEw2MagCode) (3)
integer*4 r_W2mag (nr__) ! [10]?=0 WISE W2 (4.6um) band magnitude source
* (WISEw2MagSOurce) (2)
real*4 W3mag (nr__) ! (mag) ? WISE W3 (12um) band magnitude (WISEw3Mag)
real*4 e_W3mag (nr__) ! (mag) ? WISE W3 (12um) band magnitude error
* (WISEw3MagErr)
integer*4 f_W3mag (nr__) ! [57]?=0 WISE W3 (12um) band code
* (WISEw3MagCode) (3)
integer*4 r_W3mag (nr__) ! [10]?=0 WISE W3 (12um) band magnitude source
* (WISEw3MagSOurce) (2)
real*4 W4mag (nr__) ! (mag) ? WISE W4 (22um) band magnitude (WISEw4Mag)
real*4 e_W4mag (nr__) ! (mag) ? WISE W4 (22um) band magnitude error
* (WISEw4MagErr)
integer*4 f_W4mag (nr__) ! [58]?=0 WISE W4 (22um) band code
* (WISEw4MagCode) (3)
integer*4 r_W4mag (nr__) ! [10]?=0 WISE W4 (22um) band magnitude source
* (WISEw4MagSOurce) (2)
real*4 FUVmag (nr__) ! (mag) ? GALEX FUV band magnitude
real*4 e_FUVmag (nr__) ! (mag) ? GALEX FUV band magnitude error (FUVmagErr)
integer*4 f_FUVmag (nr__) ! [80]?=0 GALEX FUV band code (FUVmagCode) (3)
integer*4 r_FUVmag (nr__) ! [11]?=0 GALEX FUV band magnitude source
* (FUVmagSource) (2)
real*4 NUVmag (nr__) ! (mag) ? GALEX NUV band magnitude
real*4 e_NUVmag (nr__) ! (mag) ? GALEX NUV band magnitude error (NUVmagErr)
integer*4 f_NUVmag (nr__) ! [81]?=0 GALEX NUV band code (NUVmagCode) (3)
integer*4 r_NUVmag (nr__) ! [11]?=0 GALEX NUV band magnitude source
* (NUVmagSource) (2)
real*8 a (nr__) ! (pix) ? SemiMajor diameter of object
* (gsc2semiMajorAxis)
real*4 e (nr__) ! (pix) ? Eccentricity of object (gsc2eccentricity)
real*4 PA (nr__) ! (pix) ? Position Angle (North thru East)
* (gsc2positionAngle)
integer*8 GaiaDR2 (nr__) ! ? Gaia DR2 ID (gaiaDr2SourceID)
real*8 RPmag (nr__) ! (mag) ? Gaia Rbp band magnitude (gaiaRbpMag)
real*8 e_RPmag (nr__) ! (mag) ? Gaia Rbp band magnitude (gaiaRbpMagErr)
integer*4 f_RPmag (nr__) ! [66]?=0 Gaia Rbp band code
* (gaiaRbpMagCode) (3)
real*8 BPmag (nr__) ! (mag) ? Gaia Bbp band magnitude (gaiaBbpMag)
real*8 e_BPmag (nr__) ! (mag) ? Gaia Bbp band magnitude (gaiaBbpMagErr)
integer*4 f_BPmag (nr__) ! [67]?=0 Gaia Bbp band code
* (gaiaBbpMagCode) (3)
*Note (1): Classification codes as follows:
* 0 = Star
* 1 = Galaxy
* 2 = Blend
* 3 = Non-Star
* 4 = Unclassified
* 5 = Defect
*Note (2): References as follows:
* 0 = Human, Manual
* 1 = GSC1, DSS1 plates
* 2 = GSC2, DSS2 plates
* 3 = SkyMap, NASA Sky2000v4 catalog
* 4 = USNO-B, USNO DSS2 plates
* 5 = UCAC, USNO CCD Astrograph catalog
* 6 = URAT, USNO Robotic Astrometric Telescope catalog
* 7 = PPMXL, Heidelberg re-reduction of USNOB/2MASS
* 8 = IRAS, NASA Infrared Astronomical Satellite
* 9 = 2MASS, Two Micron All Sky Survey
* 10 = WISE, Wide-Field Infrared Survey Explorer
* 11 = GALEX, Galaxy Evolution Explorer
* 13 = HIPPARCOS, ESA Mission
* 14 = TYCHO, Astrometric Catalog using Star-mappers on Hipparcos mission
* 15 = GAIA, ESA Mission
* 16 = PS1, Panoramic Survey Telescope and Rapid Response System (PanSTARRS)
* 17 = SDSS, Sloan Digital Sky Survey
* 18 = SMSS, SkyMapper Southern Survey
* 19 = APASS, AAVSO Photometric All-Sky Survey
* 20 = VISTA, Visible and Infrared Survey Telescope for Astronomy (ESO)
* 21 = VISTA-VHS, Vista Hemisphere Survey
* 22 = VISTA-UltraVISTA, VISTA deep imaging survey
* 23 = VISTA-VIKING, VISTA Kilo-Degree Infrared Galaxy Survey
* 24 = VISTA-VMC, VISTA Magellanic Survey
* 25 = VISTA-VVV, VISTA Variables in the Via Lactea
* 26 = VISTA-VIDEO, VISTA Deep Extragalactic Observations Survey
* 30 = VST, VLT Survey Telescope (ESO)
* 31 = VST-ATLAS, VST deep survey
* 32 = VST-KIDS, The Kilo-Degree Survey
* 33 = VST-VPAS, The VST Photometric Halpha Survey of the Southern Galactic Plane
*Note (3): coded emulsion / bandpass / filter as follows:
* ----------------------------------------------------------------------
* bcode Name Emulsion/Filter
* ----------------------------------------------------------------------
* 0 = Jpg IIIaJ+GG395 (SERC-J/EJ)
* 1 = V IIaD+W12 (Pal Quick-V)
* 2 = - HST FGS
* 3 = B Johnson
* 4 = V Johnson
* 5 = Fpg IIIaF+RG630 (ESO-R)
* 6 = V495 IIaD+GG495 (Pal QV/AAO XV)
* 7 = O 103aO+no filter (POSS-I Blue)
* 8 = E 103aE+red plexi (POSS-I Red)
* 9 = Fpg IIIaF+RG630
* 10 = - IIaD+GG495+yellow objective (GPO Astrograph)
* 11 = - 103aO+blue objective (Black Birch Astrograph)
* 12 = - 103aO+blue objective Black Birch Astrograph, GSC cal.)
* 13 = - 103aG+GG495+yellow objective (Black Birch Astrograph)
* 14 = - 103aG+GG495+yellow objective (Black Birch Astrograph, GSC)
* 16 = Vpg IIIaJ+GG49
* 18 = Jpg IIIaJ+GG385 (POSS-II Blue)
* 19 = U Johnson
* 20 = R Johnson
* 21 = I Johnson
* 22 = U Cape
* 23 = R Kron
* 24 = I Kron
* 35 = Fpg IIIaF+RG610 (POSS-II Red)
* 36 = Fpg IIIaF+OG590 (SERC-ER/SR, AAO-R/GR)
* 37 = Npg IVN+RG9 (POSS-II IR)
* 38 = Npg IVN+RG715 (SERC-IR)
* 39 = Jpg IIaO+GG385 (ESO-B)
* 41 = BT TYCHO-B
* 42 = VT TYCHO-V
* 43 = - IRAS 12micron
* 44 = - IRAS 25 micron
* 45 = - IRAS 60 micron
* 46 = - IRAS 100 micron
* 47 = J 2MASS
* 48 = H 2MASS
* 49 = K 2MASS
* 50 = u' SDSS
* 51 = g' SDSS
* 52 = r' SDSS
* 53 = i' SDSS
* 54 = z' SDSS
* 55 = W1 WISE
* 56 = W2 WISE
* 57 = W3 WISE
* 58 = W4 WISE
* 59 = g PanSTARRS
* 60 = r PanSTARRS
* 61 = i PanSTARRS
* 62 = z PanSTARRS
* 63 = y PanSTARRS
* 64 = w PanSTARRS
* 65 = G Gaia G
* 66 = RP Gaia Rbp (red)
* 67 = BP Gaia Bbp (blue)
* 68 = Grvs GAIA (rv system)
* 69 = Z VISTA
* 70 = Y VISTA
* 71 = J VISTA
* 72 = H VISTA
* 73 = Ks VISTA
* 74 = u SSMS
* 75 = V SSMS
* 76 = g SSMS
* 77 r SSMS
* 78 = i SSMS
* 79 = z SSMS
* 80 = FUV GALEX
* 81 = NUV GALEX
* ----------------------------------------------------------------------
*Note (4): The sourceStatus field was primarily for indicating plate processing
* information, but augmented to show original catalog if it did not come from
* a plate observation.
* -----------------------------------------------------------------------------
* 1 and 2 = Flag: Number of plates on which this object appeared
* 0 : This may include plate observations that were excluded
* from the GSC2.3.
* 1 : These two digits are read together as an ordinary base
* ten number, eg, 04 or 12.
* -----------------------------------------------------------------------
* 3 = Flag : Centroider type associated with the exported position
* 0 : 0 = Barycenter
* 1 = Circular Gaussian
* 2 = Elliptical Gaussian
* 3 = FPA-applied barycenter
* 4 = multicircular
* 5 = multielliptical
* 6 = FPA + circular
* 7 = FPA + elliptical
* -----------------------------------------------------------------------
* 4 = Flag : Quality of exported J magnitude
* 0 : not present
* 1 : fit
* 2 : extrapolated
* -----------------------------------------------------------------------
* 5 = Flag : Quality of exported F magnitude
* 0 : not present
* 1 : fit
* 2 : extrapolated
* -----------------------------------------------------------------------
* 6 = Flag : Quality of exported V magnitude
* 0 : not present
* 1 : fit
* 2 : extrapolated
* -----------------------------------------------------------------------
* 7 = Flag : Classification|mixed vote
* 0 : mixed vote
* 1 : unanimous vote
* 2 : unanimous defect
* -----------------------------------------------------------------------
* 8 = Flag : Classification voters
* 0 : More than one 15um scan
* 1 : One 15um scan
* 2 : More than one 25um scan
* 3 : One 25um scan
* -----------------------------------------------------------------------
* 9 = Flag : Incomplete processing
* 0 : Completed processing on all plates
* 1 : Object was too big to be cut out (>256x256 pixels)
* on at least one plate
* -----------------------------------------------------------------------
* 10 = Flag : Deblended object
* 0 : Single object on all plates
* 1 : Child (deblended) object on at least one plate
* -----------------------------------------------------------------------
*
* Sources from Catalogs have source code x 100 as source status
* (e.g. sdss = 1700)
C=============================================================================
C Loading file 'gsc242.sam' ! Sample records of GSC 2.4.2
C Format for file interpretation
1 format(
+ A10,1X,F15.11,1X,F15.11,1X,I15,1X,I9,1X,A15,1X,I13,1X,I19,1X,
+ I2,1X,I10,1X,F8.3,1X,F15.7,1X,F15.7,1X,I2,1X,F13.3,1X,F11.3,
+ 1X,F11.3,1X,F11.3,1X,F7.3,1X,I2,1X,F10.4,1X,F8.4,1X,I2,1X,
+ F7.2,1X,F5.2,1X,I2,1X,I1,1X,I2,1X,I1,1X,I2,1X,I1,5X,I1,1X,
+ F9.6,1X,F8.6,1X,I2,1X,F9.6,1X,F8.6,1X,I1,1X,F9.6,1X,F8.6,1X,
+ I2,1X,F9.6,1X,F8.6,1X,I1,1X,F9.6,1X,F8.6,1X,I2,1X,F9.6,1X,
+ F8.6,1X,I1,1X,F9.6,1X,F8.6,1X,I2,1X,F9.6,1X,F8.6,1X,I2,1X,I1,
+ 1X,F9.6,1X,F9.6,1X,I2,1X,I2,1X,F9.6,1X,F9.6,1X,I2,1X,I2,1X,
+ F6.3,1X,F6.3,1X,I2,1X,I2,1X,F6.3,1X,F6.3,1X,I2,1X,I2,1X,F6.3,
+ 1X,F6.3,1X,I2,1X,I2,1X,F6.3,1X,F7.3,1X,I2,1X,I2,1X,F6.3,1X,
+ F6.3,1X,I2,1X,I2,1X,F7.4,1X,F7.4,1X,I2,1X,I2,1X,F7.3,1X,F6.3,
+ 1X,I2,1X,I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,1X,F7.3,1X,F6.3,1X,I2,
+ 1X,I2,1X,F6.3,1X,F6.3,1X,I2,1X,I2,1X,F7.3,1X,F6.3,1X,I2,1X,I2,
+ 1X,F6.3,1X,F5.3,1X,I2,1X,I2,1X,F6.3,1X,F5.3,1X,I2,1X,I2,1X,
+ F6.3,1X,F5.3,1X,I2,1X,I2,1X,F6.3,1X,F5.3,1X,I2,1X,I2,1X,F6.3,
+ 1X,F5.3,1X,I2,1X,I2,1X,F6.3,1X,F5.3,1X,I2,1X,I2,1X,F9.3,1X,
+ F5.3,1X,F5.1,1X,I19,1X,F7.4,1X,F7.4,1X,I2,1X,F7.4,1X,F9.4,1X,
+ I2)
C Effective file loading
open(unit=1,status='old',file=
+'gsc242.sam')
write(6,*) '....Loading file: gsc242.sam'
do i__=1,1000
read(1,'(A945)')ar__
read(ar__,1)
+ GSC2(i__),RAdeg(i__),DEdeg(i__),objID(i__),GSC1(i__),
+ GSC2intID(i__),Compass(i__),GaiaDR1(i__),usageFlag(i__),
+ Status(i__),Epoch(i__),e_RAdeg(i__),e_DEdeg(i__),r_pos(i__),
+ pmRA(i__),pmDE(i__),e_pmRA(i__),e_pmDE(i__),dEpoch(i__),
+ r_pm(i__),plx(i__),e_plx(i__),r_plx(i__),RV(i__),e_RV(i__),
+ r_RV(i__),Class(i__),r_Class(i__),varFlg(i__),r_varFlg(i__),
+ multFlg(i__),r_multFlg(i__),Gmag(i__),e_Gmag(i__),f_Gmag(i__),
+ OMag(i__),e_OMag(i__),f_OMag(i__),Bjmag(i__),e_Bjmag(i__),
+ f_Bjmag(i__),VpgMag(i__),e_VpgMag(i__),f_VpgMag(i__),
+ Fpgmag(i__),e_Fpgmag(i__),f_Fpgmag(i__),Epgmag(i__),
+ e_Epgmag(i__),f_Epgmag(i__),Npgmag(i__),e_Npgmag(i__),
+ f_Npgmag(i__),Umag(i__),e_Umag(i__),f_Umag(i__),r_Umag(i__),
+ Bmag(i__),e_Bmag(i__),f_Bmag(i__),r_Bmag(i__),Vmag(i__),
+ e_Vmag(i__),f_Vmag(i__),r_Vmag(i__),umag_1(i__),e_umag_1(i__),
+ f_umag_1(i__),r_umag_1(i__),gmag_1(i__),e_gmag_1(i__),
+ f_gmag_1(i__),r_gmag(i__),rmag(i__),e_rmag(i__),f_rmag(i__),
+ r_rmag(i__),imag(i__),e_imag(i__),f_imag(i__),r_imag(i__),
+ zmag(i__),e_zmag(i__),f_zmag(i__),r_zmag(i__),ymag(i__),
+ e_ymag(i__),f_ymag(i__),r_ymag(i__),Jmag(i__),e_Jmag(i__),
+ f_Jmag(i__),r_Jmag(i__),Hmag(i__),e_Hmag(i__),f_Hmag(i__),
+ r_Hmag(i__),Ksmag(i__),e_Ksmag(i__),f_Ksmag(i__),r_Ksmag(i__),
+ Zmag_1(i__),e_Zmag_1(i__),f_Zmag_1(i__),r_Zmag_1(i__),
+ Ymag_1(i__),e_Ymag_1(i__),f_Ymag_1(i__),r_Ymag_1(i__),
+ W1mag(i__),e_W1mag(i__),f_W1mag(i__),r_W1mag(i__),W2mag(i__),
+ e_W2mag(i__),f_W2mag(i__),r_W2mag(i__),W3mag(i__),
+ e_W3mag(i__),f_W3mag(i__),r_W3mag(i__),W4mag(i__),
+ e_W4mag(i__),f_W4mag(i__),r_W4mag(i__),FUVmag(i__),
+ e_FUVmag(i__),f_FUVmag(i__),r_FUVmag(i__),NUVmag(i__),
+ e_NUVmag(i__),f_NUVmag(i__),r_NUVmag(i__),a(i__),e(i__),
+ PA(i__),GaiaDR2(i__),RPmag(i__),e_RPmag(i__),f_RPmag(i__),
+ BPmag(i__),e_BPmag(i__),f_BPmag(i__)
if(ar__(60:68) .EQ. '') GSC1(i__) = iNULL__
if(ar__(100:118) .EQ. '') GaiaDR1(i__) = iNULL__
if(ar__(123:132) .EQ. '') Status(i__) = iNULL__
if(ar__(143:157) .EQ. '') e_RAdeg(i__) = rNULL__
if(ar__(159:173) .EQ. '') e_DEdeg(i__) = rNULL__
if(ar__(178:190) .EQ. '') pmRA(i__) = rNULL__
if(ar__(192:202) .EQ. '') pmDE(i__) = rNULL__
if(ar__(204:214) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__(216:226) .EQ. '') e_pmDE(i__) = rNULL__
if(ar__(228:234) .EQ. '') dEpoch(i__) = rNULL__
if(ar__(236:237) .EQ. '') r_pm(i__) = iNULL__
if(ar__(239:248) .EQ. '') plx(i__) = rNULL__
if(ar__(250:257) .EQ. '') e_plx(i__) = rNULL__
if(ar__(259:260) .EQ. '') r_plx(i__) = iNULL__
if(ar__(262:268) .EQ. '') RV(i__) = rNULL__
if(ar__(270:274) .EQ. '') e_RV(i__) = rNULL__
if(ar__(276:277) .EQ. '') r_RV(i__) = iNULL__
if(ar__(279:279) .EQ. '') Class(i__) = iNULL__
if(ar__(281:282) .EQ. '') r_Class(i__) = iNULL__
if(ar__(284:284) .EQ. '') varFlg(i__) = iNULL__
if(ar__(286:287) .EQ. '') r_varFlg(i__) = iNULL__
if(ar__(297:305) .EQ. '') Gmag(i__) = rNULL__
if(ar__(307:314) .EQ. '') e_Gmag(i__) = rNULL__
if(ar__(319:327) .EQ. '') OMag(i__) = rNULL__
if(ar__(329:336) .EQ. '') e_OMag(i__) = rNULL__
if(ar__(340:348) .EQ. '') Bjmag(i__) = rNULL__
if(ar__(350:357) .EQ. '') e_Bjmag(i__) = rNULL__
if(ar__(362:370) .EQ. '') VpgMag(i__) = rNULL__
if(ar__(372:379) .EQ. '') e_VpgMag(i__) = rNULL__
if(ar__(383:391) .EQ. '') Fpgmag(i__) = rNULL__
if(ar__(393:400) .EQ. '') e_Fpgmag(i__) = rNULL__
if(ar__(405:413) .EQ. '') Epgmag(i__) = rNULL__
if(ar__(415:422) .EQ. '') e_Epgmag(i__) = rNULL__
if(ar__(426:434) .EQ. '') Npgmag(i__) = rNULL__
if(ar__(436:443) .EQ. '') e_Npgmag(i__) = rNULL__
if(ar__(448:456) .EQ. '') Umag(i__) = rNULL__
if(ar__(458:465) .EQ. '') e_Umag(i__) = rNULL__
if(ar__(472:480) .EQ. '') Bmag(i__) = rNULL__
if(ar__(482:490) .EQ. '') e_Bmag(i__) = rNULL__
if(ar__(498:506) .EQ. '') Vmag(i__) = rNULL__
if(ar__(508:516) .EQ. '') e_Vmag(i__) = rNULL__
if(ar__(524:529) .EQ. '') umag_1(i__) = rNULL__
if(ar__(531:536) .EQ. '') e_umag_1(i__) = rNULL__
if(ar__(544:549) .EQ. '') gmag_1(i__) = rNULL__
if(ar__(551:556) .EQ. '') e_gmag_1(i__) = rNULL__
if(ar__(564:569) .EQ. '') rmag(i__) = rNULL__
if(ar__(571:576) .EQ. '') e_rmag(i__) = rNULL__
if(ar__(584:589) .EQ. '') imag(i__) = rNULL__
if(ar__(591:597) .EQ. '') e_imag(i__) = rNULL__
if(ar__(605:610) .EQ. '') zmag(i__) = rNULL__
if(ar__(612:617) .EQ. '') e_zmag(i__) = rNULL__
if(ar__(625:631) .EQ. '') ymag(i__) = rNULL__
if(ar__(633:639) .EQ. '') e_ymag(i__) = rNULL__
if(ar__(647:653) .EQ. '') Jmag(i__) = rNULL__
if(ar__(655:660) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(668:674) .EQ. '') Hmag(i__) = rNULL__
if(ar__(676:681) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(689:695) .EQ. '') Ksmag(i__) = rNULL__
if(ar__(697:702) .EQ. '') e_Ksmag(i__) = rNULL__
if(ar__(710:715) .EQ. '') Zmag_1(i__) = rNULL__
if(ar__(717:722) .EQ. '') e_Zmag_1(i__) = rNULL__
if(ar__(730:736) .EQ. '') Ymag_1(i__) = rNULL__
if(ar__(738:743) .EQ. '') e_Ymag_1(i__) = rNULL__
if(ar__(751:756) .EQ. '') W1mag(i__) = rNULL__
if(ar__(758:762) .EQ. '') e_W1mag(i__) = rNULL__
if(ar__(770:775) .EQ. '') W2mag(i__) = rNULL__
if(ar__(777:781) .EQ. '') e_W2mag(i__) = rNULL__
if(ar__(789:794) .EQ. '') W3mag(i__) = rNULL__
if(ar__(796:800) .EQ. '') e_W3mag(i__) = rNULL__
if(ar__(808:813) .EQ. '') W4mag(i__) = rNULL__
if(ar__(815:819) .EQ. '') e_W4mag(i__) = rNULL__
if(ar__(827:832) .EQ. '') FUVmag(i__) = rNULL__
if(ar__(834:838) .EQ. '') e_FUVmag(i__) = rNULL__
if(ar__(846:851) .EQ. '') NUVmag(i__) = rNULL__
if(ar__(853:857) .EQ. '') e_NUVmag(i__) = rNULL__
if(ar__(865:873) .EQ. '') a(i__) = rNULL__
if(ar__(875:879) .EQ. '') e(i__) = rNULL__
if(ar__(881:885) .EQ. '') PA(i__) = rNULL__
if(ar__(887:905) .EQ. '') GaiaDR2(i__) = iNULL__
if(ar__(907:913) .EQ. '') RPmag(i__) = rNULL__
if(ar__(915:921) .EQ. '') e_RPmag(i__) = rNULL__
if(ar__(926:932) .EQ. '') BPmag(i__) = rNULL__
if(ar__(934:942) .EQ. '') e_BPmag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ GSC2(i__),RAdeg(i__),DEdeg(i__),objID(i__),GSC1(i__),
+ GSC2intID(i__),Compass(i__),GaiaDR1(i__),usageFlag(i__),
+ Status(i__),Epoch(i__),e_RAdeg(i__),e_DEdeg(i__),r_pos(i__),
+ pmRA(i__),pmDE(i__),e_pmRA(i__),e_pmDE(i__),dEpoch(i__),
+ r_pm(i__),plx(i__),e_plx(i__),r_plx(i__),RV(i__),e_RV(i__),
+ r_RV(i__),Class(i__),r_Class(i__),varFlg(i__),r_varFlg(i__),
+ multFlg(i__),r_multFlg(i__),Gmag(i__),e_Gmag(i__),f_Gmag(i__),
+ OMag(i__),e_OMag(i__),f_OMag(i__),Bjmag(i__),e_Bjmag(i__),
+ f_Bjmag(i__),VpgMag(i__),e_VpgMag(i__),f_VpgMag(i__),
+ Fpgmag(i__),e_Fpgmag(i__),f_Fpgmag(i__),Epgmag(i__),
+ e_Epgmag(i__),f_Epgmag(i__),Npgmag(i__),e_Npgmag(i__),
+ f_Npgmag(i__),Umag(i__),e_Umag(i__),f_Umag(i__),r_Umag(i__),
+ Bmag(i__),e_Bmag(i__),f_Bmag(i__),r_Bmag(i__),Vmag(i__),
+ e_Vmag(i__),f_Vmag(i__),r_Vmag(i__),umag_1(i__),e_umag_1(i__),
+ f_umag_1(i__),r_umag_1(i__),gmag_1(i__),e_gmag_1(i__),
+ f_gmag_1(i__),r_gmag(i__),rmag(i__),e_rmag(i__),f_rmag(i__),
+ r_rmag(i__),imag(i__),e_imag(i__),f_imag(i__),r_imag(i__),
+ zmag(i__),e_zmag(i__),f_zmag(i__),r_zmag(i__),ymag(i__),
+ e_ymag(i__),f_ymag(i__),r_ymag(i__),Jmag(i__),e_Jmag(i__),
+ f_Jmag(i__),r_Jmag(i__),Hmag(i__),e_Hmag(i__),f_Hmag(i__),
+ r_Hmag(i__),Ksmag(i__),e_Ksmag(i__),f_Ksmag(i__),r_Ksmag(i__),
+ Zmag_1(i__),e_Zmag_1(i__),f_Zmag_1(i__),r_Zmag_1(i__),
+ Ymag_1(i__),e_Ymag_1(i__),f_Ymag_1(i__),r_Ymag_1(i__),
+ W1mag(i__),e_W1mag(i__),f_W1mag(i__),r_W1mag(i__),W2mag(i__),
+ e_W2mag(i__),f_W2mag(i__),r_W2mag(i__),W3mag(i__),
+ e_W3mag(i__),f_W3mag(i__),r_W3mag(i__),W4mag(i__),
+ e_W4mag(i__),f_W4mag(i__),r_W4mag(i__),FUVmag(i__),
+ e_FUVmag(i__),f_FUVmag(i__),r_FUVmag(i__),NUVmag(i__),
+ e_NUVmag(i__),f_NUVmag(i__),r_NUVmag(i__),a(i__),e(i__),
+ PA(i__),GaiaDR2(i__),RPmag(i__),e_RPmag(i__),f_RPmag(i__),
+ BPmag(i__),e_BPmag(i__),f_BPmag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end