Conversion of standardized ReadMe file for
file /./ftp/cats/I/353 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-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 ! GSC2 Public ID (hstID)
real*8 RAdeg ! (deg) Right Ascension (ICRS) at Ep=Epoch (ra)
real*8 DEdeg ! (deg) Declination (ICRS) at Ep=Epoch (dec)
integer*8 objID ! Database ID
integer*4 GSC1 ! ? GSC1 ID (gsc1ID)
character*15 GSC2intID ! GSC2 Internal ID (gsc2ID)
integer*8 Compass ! COMPASS ID (compassgsc2id)
integer*8 GaiaDR1 ! ? Gaia DR1 ID (GAIAsourceID)
integer*4 usageFlag ! [-1/1] Flag for previously used
* guide star
integer*8 Status ! ? Object status codes (sourceStatus) (4)
real*8 Epoch ! (yr) Epoch of coordinate
real*8 e_RAdeg ! (arcsec) ? Right Ascension error (raErr)
real*8 e_DEdeg ! (arcsec) ? Declination error (decErr)
integer*4 r_pos ! [2/32] Object coordinates source
* (posSource) (2)
real*8 pmRA ! (mas/yr) ? Proper motion along RA (rapm)
real*8 pmDE ! (mas/yr) ? Proper motion along DE (decpm)
real*8 e_pmRA ! (mas/yr) ? pmRA error (rapmErr)
real*8 e_pmDE ! (mas/yr) ? pmDE error (decpmErr)
real*8 dEpoch ! (yr) ? Time interval used for PM computation
* (deltaEpoch)
integer*4 r_pm ! [0/15]? Proper motion source (pmSource) (2)
real*8 plx ! (mas) ? Parallax (parallax)
real*8 e_plx ! (mas) ? Parallax error (plxErr)
integer*4 r_plx ! [0/15]? Parallax source (plxSource) (2)
real*8 RV ! (km/s) ? Radial velocity (radialVelocity)
real*4 e_RV ! (km/s) ? Radial velocity error (rvErr)
integer*4 r_RV ! [0/15]? Radial Velocity source
* (rvSource) (2)
integer*4 Class ! [0/4]? Object classification code
* (classification) (1)
integer*4 r_Class ! [0/32]? Classification source
* (classSource) (2)
integer*4 varFlg ! [1]? Flag for known variable (variableFlag)
integer*4 r_varFlg ! [3/15]? Variability source (varSource) (2)
integer*4 multFlg ! [0/1] Flag for known multiple multipleFlag
* (multipleFlag)
integer*4 r_multFlg ! [0/3] Multiplicity source (multSource) (2)
real*8 Gmag ! (mag) ? Gaia G band magnitude (gaiaGmag)
real*8 e_Gmag ! (mag) ? Gaia G band magnitude error (gaiaGmagErr)
integer*4 f_Gmag ! [65]?=0 Gaia G band code (gaiaGmagCode) (3)
real*8 OMag ! (mag) ? Photographic O band magnitude (OpgMag)
real*8 e_OMag ! (mag) ? Photographic O band magnitude error
* (OpgMagErr)
integer*4 f_OMag ! [7]?=0 Photographic O band code
* (OpgMagCode) (3)
real*8 Bjmag ! (mag) ? Photographic Blue magnitude (JpgMag)
real*8 e_Bjmag ! (mag) ? Photographic Blue magnitude error
* (JpgMagErr)
integer*4 f_Bjmag ! [0/18]?=0 Photographic Blue bandpass code
* (JpgMagCode) (3)
real*8 VpgMag ! (mag) ? Photographic V band magnitude
real*8 e_VpgMag ! (mag) ? Photographic V band magnitude error
* (VpgMagErr)
integer*4 f_VpgMag ! [1/6]?=0 Photographic V band code
* (VpgMagCode) (3)
real*8 Fpgmag ! (mag) ? Photographic Red magnitude
* (FpgMag)
real*8 e_Fpgmag ! (mag) ? Photographic Red magnitude error
* (FpgMagErr)
integer*4 f_Fpgmag ! [35/36]?=0 Photographic Red bandpass code
* (FpgMagCode) (3)
real*8 Epgmag ! (mag) ? Photographic E band magnitude (EpgMag)
real*8 e_Epgmag ! (mag) ? Photographic E band magnitude error
* (EpgMagErr)
integer*4 f_Epgmag ! [8]?=0 Photographic E band code
* (EpgMagCode) (3)
real*8 Npgmag ! (mag) ? Photographic Near-IR magnitude (NpgMag)
real*8 e_Npgmag ! (mag) ? Photographic Near-IR magnitude error
* (NpgMagErr)
integer*4 f_Npgmag ! [37/38]?=0 Photographic Near-IR bandpass code
* (NpgMagCode) (3)
real*8 Umag ! (mag) ? U band magnitude
real*8 e_Umag ! (mag) ? U band magnitude error (UmagErr)
integer*4 f_Umag ! [19]?=0 U band bandpass code (UmagCode) (3)
integer*4 r_Umag ! [3]?=0 U band magnitude source
* (UmagSource) (2)
real*8 Bmag ! (mag) ? B band magnitude
real*8 e_Bmag ! (mag) ? B band magnitude error (BmagErr)
integer*4 f_Bmag ! [3/41]?=0 B band bandpass code (BmagCode) (3)
integer*4 r_Bmag ! [3/19]?=0 B band magnitude source
* (BmagSource) (2)
real*8 Vmag ! (mag) ? V band magnitude
real*8 e_Vmag ! (mag) ? V band magnitude error (VmagErr)
integer*4 f_Vmag ! [3/75]?=0 V band bandpass code (VMagCode) (3)
integer*4 r_Vmag ! [3/19]?=0 V band magnitude source
* (VmagSource) (2)
real*4 umag_1 ! (mag) ? SDSS u band magnitude (SDSSuMag)
real*4 e_umag_1 ! (mag) ? SDSS u band magnitude error (SDSSuMagErr)
integer*4 f_umag_1 ! [50/74]?=0 SDSS u band code
* (SDSSuMagCode) (3)
integer*4 r_umag_1 ! [17/32]?=0 SDSS u band magnitude source
* (SDSSuMagSource) (2)
real*4 gmag_1 ! (mag) ? SDSS g band magnitude (SDSSgMag)
real*4 e_gmag_1 ! (mag) ? SDSS g band magnitude error (SDSSgMagErr)
integer*4 f_gmag_1 ! [51/76]?=0 SDSS g band code
* (SDSSgMagCode) (3)
integer*4 r_gmag ! [16/32]?=0 SDSS g band magnitude source
* (SDSSgMagSource) (2)
real*4 rmag ! (mag) ? SDSS r band magnitude (SDSSrMag)
real*4 e_rmag ! (mag) ? SDSS r band magnitude error (SDSSrMagErr)
integer*4 f_rmag ! [52/77]?=0 SDSS r band code
* (SDSSrMagCode) (3)
integer*4 r_rmag ! [16/32]?=0 SDSS r band magnitude source
* (SDSSrMagSource) (2)
real*4 imag ! (mag) ? SDSS i band magnitude (SDSSiMag)
real*8 e_imag ! (mag) ? SDSS i band magnitude error (SDSSiMagErr)
integer*4 f_imag ! [53/78]?=0 SDSS i band code
* (SDSSiMagCode) (3)
integer*4 r_imag ! [16/32]?=0 SDSS i band magnitude source
* (SDSSiMagSource) (2)
real*4 zmag ! (mag) ? SDSS z band magnitude (SDSSzMag)
real*4 e_zmag ! (mag) ? SDSS z band magnitude error (SDSSzMagErr)
integer*4 f_zmag ! [54/79]?=0 SDSS z band code
* (SDSSzMagCode) (3)
integer*4 r_zmag ! [16/31]?=0 SDSS z band magnitude source
* (SDSSzMagSource) (2)
real*8 ymag ! (mag) ? PanSTARRS y band magnitude (PS1yMag)
real*8 e_ymag ! (mag) ? PanSTARRS y band magnitude error
* (PS1yMagErr)
integer*4 f_ymag ! [63]?=0 PanSTARRS y band code
* (PS1yMagCode) (3)
integer*4 r_ymag ! [16]?=0 PanSTARRS y band magnitude source
* (PS1ymagSource) (2)
real*8 Jmag ! (mag) ? 2MASS J band magnitude (tmassJmag)
real*4 e_Jmag ! (mag) ? 2MASS J band magnitude error (tmassJmagErr)
integer*4 f_Jmag ! [47]?=0 2MASS J band code (tmassJmagCode) (3)
integer*4 r_Jmag ! [9/23]?=0 2MASS J band magnitude source
* (tmassJmagSource) (2)
real*8 Hmag ! (mag) ? 2MASS H band magnitude (tmassHmag)
real*4 e_Hmag ! (mag) ? 2MASS H band magnitude error (tmassHmagErr)
integer*4 f_Hmag ! [48]?=0 2MASS H band code (tmassHmagCode) (3)
integer*4 r_Hmag ! [9/23]?=0 2MASS H band magnitude source
* (tmassHmagSource) (2)
real*8 Ksmag ! (mag) ? 2MASS Ks band magnitude (tmassKsmag)
real*4 e_Ksmag ! (mag) ? 2MASS Ks band magnitude error
* (tmassKsmagErr)
integer*4 f_Ksmag ! [49]?=0 2MASS Ks band code
* (tmassKsmagCode) (3)
integer*4 r_Ksmag ! [9/23]?=0 2MASS Ks band magnitude source
* (tmassKsmagSource) (2)
real*4 Zmag_1 ! (mag) ? VISTA Z band magnitude (VistaZMag)
real*4 e_Zmag_1 ! (mag) ? VISTA Z band magnitude error (VistaZMagErr)
integer*4 f_Zmag_1 ! [69]?=0 VISTA Z band code (VistaZMagCode) (3)
integer*4 r_Zmag_1 ! [23]?=0 VISTA Z band magnitude source
* (VistaZmagSource) (2)
real*8 Ymag_1 ! (mag) ? VISTA Y band magnitude (VistaYmag)
real*4 e_Ymag_1 ! (mag) ? VISTA Y band magnitude error (VistaYMagErr)
integer*4 f_Ymag_1 ! [70]?=0 VISTA Y band code (VistaYMagCode) (3)
integer*4 r_Ymag_1 ! [21/23]?=0 VISTA Y band magnitude source
* (VistaYmagSource) (2)
real*4 W1mag ! (mag) ? WISE W1 (3.4um) band magnitude (WISEw1Mag)
real*4 e_W1mag ! (mag) ? WISE W1 (3.4um) band magnitude error
* (WISEw1MagErr)
integer*4 f_W1mag ! [55]?=0 WISE W1 (3.4um) band code
* (WISEw1MagCode) (3)
integer*4 r_W1mag ! [10]?=0 WISE W1 (3.4um) band magnitude source
* (WISEw1MagSOurce) (2)
real*4 W2mag ! (mag) ? WISE W2 (4.6um) band magnitude (WISEw2Mag)
real*4 e_W2mag ! (mag) ? WISE W2 (4.6um) band magnitude error
* (WISEw2MagErr)
integer*4 f_W2mag ! [56]?=0 WISE W2 (4.6um) band code
* (WISEw2MagCode) (3)
integer*4 r_W2mag ! [10]?=0 WISE W2 (4.6um) band magnitude source
* (WISEw2MagSOurce) (2)
real*4 W3mag ! (mag) ? WISE W3 (12um) band magnitude (WISEw3Mag)
real*4 e_W3mag ! (mag) ? WISE W3 (12um) band magnitude error
* (WISEw3MagErr)
integer*4 f_W3mag ! [57]?=0 WISE W3 (12um) band code
* (WISEw3MagCode) (3)
integer*4 r_W3mag ! [10]?=0 WISE W3 (12um) band magnitude source
* (WISEw3MagSOurce) (2)
real*4 W4mag ! (mag) ? WISE W4 (22um) band magnitude (WISEw4Mag)
real*4 e_W4mag ! (mag) ? WISE W4 (22um) band magnitude error
* (WISEw4MagErr)
integer*4 f_W4mag ! [58]?=0 WISE W4 (22um) band code
* (WISEw4MagCode) (3)
integer*4 r_W4mag ! [10]?=0 WISE W4 (22um) band magnitude source
* (WISEw4MagSOurce) (2)
real*4 FUVmag ! (mag) ? GALEX FUV band magnitude
real*4 e_FUVmag ! (mag) ? GALEX FUV band magnitude error (FUVmagErr)
integer*4 f_FUVmag ! [80]?=0 GALEX FUV band code (FUVmagCode) (3)
integer*4 r_FUVmag ! [11]?=0 GALEX FUV band magnitude source
* (FUVmagSource) (2)
real*4 NUVmag ! (mag) ? GALEX NUV band magnitude
real*4 e_NUVmag ! (mag) ? GALEX NUV band magnitude error (NUVmagErr)
integer*4 f_NUVmag ! [81]?=0 GALEX NUV band code (NUVmagCode) (3)
integer*4 r_NUVmag ! [11]?=0 GALEX NUV band magnitude source
* (NUVmagSource) (2)
real*8 a ! (pix) ? SemiMajor diameter of object
* (gsc2semiMajorAxis)
real*4 e ! (pix) ? Eccentricity of object (gsc2eccentricity)
real*4 PA ! (pix) ? Position Angle (North thru East)
* (gsc2positionAngle)
integer*8 GaiaDR2 ! ? Gaia DR2 ID (gaiaDr2SourceID)
real*8 RPmag ! (mag) ? Gaia Rbp band magnitude (gaiaRbpMag)
real*8 e_RPmag ! (mag) ? Gaia Rbp band magnitude (gaiaRbpMagErr)
integer*4 f_RPmag ! [66]?=0 Gaia Rbp band code
* (gaiaRbpMagCode) (3)
real*8 BPmag ! (mag) ? Gaia Bbp band magnitude (gaiaBbpMag)
real*8 e_BPmag ! (mag) ? Gaia Bbp band magnitude (gaiaBbpMagErr)
integer*4 f_BPmag ! [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,RAdeg,DEdeg,objID,GSC1,GSC2intID,Compass,GaiaDR1,
+ usageFlag,Status,Epoch,e_RAdeg,e_DEdeg,r_pos,pmRA,pmDE,e_pmRA,
+ e_pmDE,dEpoch,r_pm,plx,e_plx,r_plx,RV,e_RV,r_RV,Class,r_Class,
+ varFlg,r_varFlg,multFlg,r_multFlg,Gmag,e_Gmag,f_Gmag,OMag,
+ e_OMag,f_OMag,Bjmag,e_Bjmag,f_Bjmag,VpgMag,e_VpgMag,f_VpgMag,
+ Fpgmag,e_Fpgmag,f_Fpgmag,Epgmag,e_Epgmag,f_Epgmag,Npgmag,
+ e_Npgmag,f_Npgmag,Umag,e_Umag,f_Umag,r_Umag,Bmag,e_Bmag,
+ f_Bmag,r_Bmag,Vmag,e_Vmag,f_Vmag,r_Vmag,umag_1,e_umag_1,
+ f_umag_1,r_umag_1,gmag_1,e_gmag_1,f_gmag_1,r_gmag,rmag,e_rmag,
+ f_rmag,r_rmag,imag,e_imag,f_imag,r_imag,zmag,e_zmag,f_zmag,
+ r_zmag,ymag,e_ymag,f_ymag,r_ymag,Jmag,e_Jmag,f_Jmag,r_Jmag,
+ Hmag,e_Hmag,f_Hmag,r_Hmag,Ksmag,e_Ksmag,f_Ksmag,r_Ksmag,
+ Zmag_1,e_Zmag_1,f_Zmag_1,r_Zmag_1,Ymag_1,e_Ymag_1,f_Ymag_1,
+ r_Ymag_1,W1mag,e_W1mag,f_W1mag,r_W1mag,W2mag,e_W2mag,f_W2mag,
+ r_W2mag,W3mag,e_W3mag,f_W3mag,r_W3mag,W4mag,e_W4mag,f_W4mag,
+ r_W4mag,FUVmag,e_FUVmag,f_FUVmag,r_FUVmag,NUVmag,e_NUVmag,
+ f_NUVmag,r_NUVmag,a,e,PA,GaiaDR2,RPmag,e_RPmag,f_RPmag,BPmag,
+ e_BPmag,f_BPmag
if(ar__(60:68) .EQ. '') GSC1 = iNULL__
if(ar__(100:118) .EQ. '') GaiaDR1 = iNULL__
if(ar__(123:132) .EQ. '') Status = iNULL__
if(ar__(143:157) .EQ. '') e_RAdeg = rNULL__
if(ar__(159:173) .EQ. '') e_DEdeg = rNULL__
if(ar__(178:190) .EQ. '') pmRA = rNULL__
if(ar__(192:202) .EQ. '') pmDE = rNULL__
if(ar__(204:214) .EQ. '') e_pmRA = rNULL__
if(ar__(216:226) .EQ. '') e_pmDE = rNULL__
if(ar__(228:234) .EQ. '') dEpoch = rNULL__
if(ar__(236:237) .EQ. '') r_pm = iNULL__
if(ar__(239:248) .EQ. '') plx = rNULL__
if(ar__(250:257) .EQ. '') e_plx = rNULL__
if(ar__(259:260) .EQ. '') r_plx = iNULL__
if(ar__(262:268) .EQ. '') RV = rNULL__
if(ar__(270:274) .EQ. '') e_RV = rNULL__
if(ar__(276:277) .EQ. '') r_RV = iNULL__
if(ar__(279:279) .EQ. '') Class = iNULL__
if(ar__(281:282) .EQ. '') r_Class = iNULL__
if(ar__(284:284) .EQ. '') varFlg = iNULL__
if(ar__(286:287) .EQ. '') r_varFlg = iNULL__
if(ar__(297:305) .EQ. '') Gmag = rNULL__
if(ar__(307:314) .EQ. '') e_Gmag = rNULL__
if(ar__(319:327) .EQ. '') OMag = rNULL__
if(ar__(329:336) .EQ. '') e_OMag = rNULL__
if(ar__(340:348) .EQ. '') Bjmag = rNULL__
if(ar__(350:357) .EQ. '') e_Bjmag = rNULL__
if(ar__(362:370) .EQ. '') VpgMag = rNULL__
if(ar__(372:379) .EQ. '') e_VpgMag = rNULL__
if(ar__(383:391) .EQ. '') Fpgmag = rNULL__
if(ar__(393:400) .EQ. '') e_Fpgmag = rNULL__
if(ar__(405:413) .EQ. '') Epgmag = rNULL__
if(ar__(415:422) .EQ. '') e_Epgmag = rNULL__
if(ar__(426:434) .EQ. '') Npgmag = rNULL__
if(ar__(436:443) .EQ. '') e_Npgmag = rNULL__
if(ar__(448:456) .EQ. '') Umag = rNULL__
if(ar__(458:465) .EQ. '') e_Umag = rNULL__
if(ar__(472:480) .EQ. '') Bmag = rNULL__
if(ar__(482:490) .EQ. '') e_Bmag = rNULL__
if(ar__(498:506) .EQ. '') Vmag = rNULL__
if(ar__(508:516) .EQ. '') e_Vmag = rNULL__
if(ar__(524:529) .EQ. '') umag_1 = rNULL__
if(ar__(531:536) .EQ. '') e_umag_1 = rNULL__
if(ar__(544:549) .EQ. '') gmag_1 = rNULL__
if(ar__(551:556) .EQ. '') e_gmag_1 = rNULL__
if(ar__(564:569) .EQ. '') rmag = rNULL__
if(ar__(571:576) .EQ. '') e_rmag = rNULL__
if(ar__(584:589) .EQ. '') imag = rNULL__
if(ar__(591:597) .EQ. '') e_imag = rNULL__
if(ar__(605:610) .EQ. '') zmag = rNULL__
if(ar__(612:617) .EQ. '') e_zmag = rNULL__
if(ar__(625:631) .EQ. '') ymag = rNULL__
if(ar__(633:639) .EQ. '') e_ymag = rNULL__
if(ar__(647:653) .EQ. '') Jmag = rNULL__
if(ar__(655:660) .EQ. '') e_Jmag = rNULL__
if(ar__(668:674) .EQ. '') Hmag = rNULL__
if(ar__(676:681) .EQ. '') e_Hmag = rNULL__
if(ar__(689:695) .EQ. '') Ksmag = rNULL__
if(ar__(697:702) .EQ. '') e_Ksmag = rNULL__
if(ar__(710:715) .EQ. '') Zmag_1 = rNULL__
if(ar__(717:722) .EQ. '') e_Zmag_1 = rNULL__
if(ar__(730:736) .EQ. '') Ymag_1 = rNULL__
if(ar__(738:743) .EQ. '') e_Ymag_1 = rNULL__
if(ar__(751:756) .EQ. '') W1mag = rNULL__
if(ar__(758:762) .EQ. '') e_W1mag = rNULL__
if(ar__(770:775) .EQ. '') W2mag = rNULL__
if(ar__(777:781) .EQ. '') e_W2mag = rNULL__
if(ar__(789:794) .EQ. '') W3mag = rNULL__
if(ar__(796:800) .EQ. '') e_W3mag = rNULL__
if(ar__(808:813) .EQ. '') W4mag = rNULL__
if(ar__(815:819) .EQ. '') e_W4mag = rNULL__
if(ar__(827:832) .EQ. '') FUVmag = rNULL__
if(ar__(834:838) .EQ. '') e_FUVmag = rNULL__
if(ar__(846:851) .EQ. '') NUVmag = rNULL__
if(ar__(853:857) .EQ. '') e_NUVmag = rNULL__
if(ar__(865:873) .EQ. '') a = rNULL__
if(ar__(875:879) .EQ. '') e = rNULL__
if(ar__(881:885) .EQ. '') PA = rNULL__
if(ar__(887:905) .EQ. '') GaiaDR2 = iNULL__
if(ar__(907:913) .EQ. '') RPmag = rNULL__
if(ar__(915:921) .EQ. '') e_RPmag = rNULL__
if(ar__(926:932) .EQ. '') BPmag = rNULL__
if(ar__(934:942) .EQ. '') e_BPmag = rNULL__
c ..............Just test output...........
write(6,1)
+ GSC2,RAdeg,DEdeg,objID,GSC1,GSC2intID,Compass,GaiaDR1,
+ usageFlag,Status,Epoch,e_RAdeg,e_DEdeg,r_pos,pmRA,pmDE,e_pmRA,
+ e_pmDE,dEpoch,r_pm,plx,e_plx,r_plx,RV,e_RV,r_RV,Class,r_Class,
+ varFlg,r_varFlg,multFlg,r_multFlg,Gmag,e_Gmag,f_Gmag,OMag,
+ e_OMag,f_OMag,Bjmag,e_Bjmag,f_Bjmag,VpgMag,e_VpgMag,f_VpgMag,
+ Fpgmag,e_Fpgmag,f_Fpgmag,Epgmag,e_Epgmag,f_Epgmag,Npgmag,
+ e_Npgmag,f_Npgmag,Umag,e_Umag,f_Umag,r_Umag,Bmag,e_Bmag,
+ f_Bmag,r_Bmag,Vmag,e_Vmag,f_Vmag,r_Vmag,umag_1,e_umag_1,
+ f_umag_1,r_umag_1,gmag_1,e_gmag_1,f_gmag_1,r_gmag,rmag,e_rmag,
+ f_rmag,r_rmag,imag,e_imag,f_imag,r_imag,zmag,e_zmag,f_zmag,
+ r_zmag,ymag,e_ymag,f_ymag,r_ymag,Jmag,e_Jmag,f_Jmag,r_Jmag,
+ Hmag,e_Hmag,f_Hmag,r_Hmag,Ksmag,e_Ksmag,f_Ksmag,r_Ksmag,
+ Zmag_1,e_Zmag_1,f_Zmag_1,r_Zmag_1,Ymag_1,e_Ymag_1,f_Ymag_1,
+ r_Ymag_1,W1mag,e_W1mag,f_W1mag,r_W1mag,W2mag,e_W2mag,f_W2mag,
+ r_W2mag,W3mag,e_W3mag,f_W3mag,r_W3mag,W4mag,e_W4mag,f_W4mag,
+ r_W4mag,FUVmag,e_FUVmag,f_FUVmag,r_FUVmag,NUVmag,e_NUVmag,
+ f_NUVmag,r_NUVmag,a,e,PA,GaiaDR2,RPmag,e_RPmag,f_RPmag,BPmag,
+ e_BPmag,f_BPmag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end