Conversion of standardized ReadMe file for
catalog into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2024-Apr-20
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/239 The Hipparcos and Tycho Catalogues (ESA 1997)
*================================================================================
*The Hipparcos and Tycho Catalogues
* ESA 1997
* <ESA, 1997, The Hipparcos Catalogue, ESA SP-1200>
* <ESA, 1997, The Tycho Catalogue, ESA SP-1200>
* =1997HIP...C......0E
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_main.dat' ! The Hipparcos Main Catalogue
integer*4 nr__
parameter (nr__=118218) ! Number of records
character*450 ar__ ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=J1991.25)
character*1 Catalog (nr__) ! [H] Catalogue (H=Hipparcos) (H0)
integer*4 HIP (nr__) ! Identifier (HIP number) (H1)
character*1 Proxy (nr__) ! *[HT] Proximity flag (H2)
character*11 RAhms (nr__) ! Right ascension in h m s, ICRS (J1991.25) (H3)
character*11 DEdms (nr__) ! Declination in deg ' ", ICRS (J1991.25) (H4)
real*4 Vmag (nr__) ! (mag) ? Magnitude in Johnson V (H5)
integer*4 VarFlag (nr__) ! *[1,3]? Coarse variability flag (H6)
character*1 r_Vmag (nr__) ! *[GHT] Source of magnitude (H7)
real*8 RAdeg (nr__) ! (deg) *? alpha, degrees (ICRS, Epoch=J1991.25) (H8)
real*8 DEdeg (nr__) ! (deg) *? delta, degrees (ICRS, Epoch=J1991.25) (H9)
character*1 AstroRef (nr__) ! *[*+A-Z] Reference flag for astrometry (H10)
real*8 Plx (nr__) ! (mas) ? Trigonometric parallax (H11)
real*8 pmRA (nr__) ! (mas/yr) *? Proper motion mu_alpha.cos(delta), ICRS(H12)
real*8 pmDE (nr__) ! (mas/yr) *? Proper motion mu_delta, ICRS (H13)
real*4 e_RAdeg (nr__) ! (mas) *? Standard error in RA*cos(DEdeg) (H14)
real*4 e_DEdeg (nr__) ! (mas) *? Standard error in DE (H15)
real*4 e_Plx (nr__) ! (mas) ? Standard error in Plx (H16)
real*4 e_pmRA (nr__) ! (mas/yr) ? Standard error in pmRA (H17)
real*4 e_pmDE (nr__) ! (mas/yr) ? Standard error in pmDE (H18)
real*4 DE_RA (nr__) ! [-1/1]? Correlation, DE/RA*cos(delta) (H19)
real*4 Plx_RA (nr__) ! [-1/1]? Correlation, Plx/RA*cos(delta) (H20)
real*4 Plx_DE (nr__) ! [-1/1]? Correlation, Plx/DE (H21)
real*4 pmRA_RA (nr__) ! [-1/1]? Correlation, pmRA/RA*cos(delta) (H22)
real*4 pmRA_DE (nr__) ! [-1/1]? Correlation, pmRA/DE (H23)
real*4 pmRA_Plx (nr__) ! [-1/1]? Correlation, pmRA/Plx (H24)
real*4 pmDE_RA (nr__) ! [-1/1]? Correlation, pmDE/RA*cos(delta) (H25)
real*4 pmDE_DE (nr__) ! [-1/1]? Correlation, pmDE/DE (H26)
real*4 pmDE_Plx (nr__) ! [-1/1]? Correlation, pmDE/Plx (H27)
real*4 pmDE_pmRA (nr__) ! [-1/1]? Correlation, pmDE/pmRA (H28)
integer*4 F1 (nr__) ! (%) ? Percentage of rejected data (H29)
real*4 F2 (nr__) ! *? Goodness-of-fit parameter (H30)
real*4 BTmag (nr__) ! (mag) ? Mean BT magnitude (H32)
real*4 e_BTmag (nr__) ! (mag) ? Standard error on BTmag (H33)
real*4 VTmag (nr__) ! (mag) ? Mean VT magnitude (H34)
real*4 e_VTmag (nr__) ! (mag) ? Standard error on VTmag (H35)
character*1 m_BTmag (nr__) ! *[A-Z*-] Reference flag for BT and VTmag (H36)
real*4 B_V (nr__) ! (mag) ? Johnson B-V colour (H37)
real*4 e_B_V (nr__) ! (mag) ? Standard error on B-V (H38)
character*1 r_B_V (nr__) ! [GT] Source of B-V from Ground or Tycho (H39)
real*4 V_I (nr__) ! (mag) ? Colour index in Cousins' system (H40)
real*4 e_V_I (nr__) ! (mag) ? Standard error on V-I (H41)
character*1 r_V_I (nr__) ! *[A-T] Source of V-I (H42)
character*1 CombMag (nr__) ! [*] Flag for combined Vmag, B-V, V-I (H43)
real*8 Hpmag (nr__) ! (mag) *? Median magnitude in Hipparcos system (H44)
real*4 e_Hpmag (nr__) ! (mag) *? Standard error on Hpmag (H45)
real*4 Hpscat (nr__) ! (mag) ? Scatter on Hpmag (H46)
integer*4 o_Hpmag (nr__) ! ? Number of observations for Hpmag (H47)
character*1 m_Hpmag (nr__) ! *[A-Z*-] Reference flag for Hpmag (H48)
real*4 Hpmax (nr__) ! (mag) ? Hpmag at maximum (5th percentile) (H49)
real*4 HPmin (nr__) ! (mag) ? Hpmag at minimum (95th percentile) (H50)
real*8 Period (nr__) ! (d) ? Variability period (days) (H51)
character*1 HvarType (nr__) ! *[CDMPRU]? variability type (H52)
character*1 moreVar (nr__) ! *[12] Additional data about variability (H53)
character*1 morePhoto (nr__) ! [ABC] Light curve Annex (H54)
character*10 CCDM (nr__) ! CCDM identifier (H55)
character*1 n_CCDM (nr__) ! *[HIM] Historical status flag (H56)
integer*4 Nsys (nr__) ! ? Number of entries with same CCDM (H57)
integer*4 Ncomp (nr__) ! ? Number of components in this entry (H58)
character*1 MultFlag (nr__) ! *[CGOVX] Double/Multiple Systems flag (H59)
character*1 Source (nr__) ! *[PFILS] Astrometric source flag (H60)
character*1 Qual (nr__) ! *[ABCDS] Solution quality (H61)
character*2 m_HIP (nr__) ! Component identifiers (H62)
integer*4 theta (nr__) ! (deg) ? Position angle between components (H63)
real*8 rho (nr__) ! (arcsec) ? Angular separation between components (H64)
real*4 e_rho (nr__) ! (arcsec) ? Standard error on rho (H65)
real*4 dHp (nr__) ! (mag) ? Magnitude difference of components (H66)
real*4 e_dHp (nr__) ! (mag) ? Standard error on dHp (H67)
character*1 Survey (nr__) ! [S] Flag indicating a Survey Star (H68)
character*1 Chart (nr__) ! *[DG] Identification Chart (H69)
character*1 Notes (nr__) ! *[DGPWXYZ] Existence of notes (H70)
integer*4 HD (nr__) ! [1/359083]? HD number <III/135> (H71)
character*10 BD (nr__) ! Bonner DM <I/119>, <I/122> (H72)
character*10 CoD (nr__) ! Cordoba Durchmusterung (DM) <I/114> (H73)
character*10 CPD (nr__) ! Cape Photographic DM <I/108> (H74)
real*4 v_V_I_red (nr__) ! (mag) V-I used for reductions (H75)
character*12 SpType (nr__) ! Spectral type (H76)
character*1 r_SpType (nr__) ! *[1234GKSX]? Source of spectral type (H77)
*Note on Proxy: this flag provides a coarse indication of the presence
* of nearby objects within 10arcsec of the given entry.
* If non-blank, it indicates that
* 'H' there is one or more distinct Hipparcos Catalogue entries,
* or distinct components of system from h_dm_com.dat
* 'T' there is one or more distinct Tycho entries
* If 'H' and 'T' apply, 'H' is adopted.
* The 'T' flag implies either an inconsistency between the Hipparcos
* and Tycho catalogues, or a deficiency in one or both of the
* catalogues.
*Note on RAdeg, DEdeg: right ascension and declination are
* expressed in degrees for epoch J1991.25 (JD2448349.0625 (TT)) in the
* ICRS (International Celestial Reference System, consistent with
* J2000) reference system.
* There are 263 cases where these fields are missing (no astrometric
* solution could be found)
*Note on pmRA, pmDE:
* The proper motions refer to the ICRS and to the epoch J1991.25.
*Note on e_RAdeg, e_DEdeg:
* The standard errors refer to the epoch J1991.25, and represent a
* minimum of the error on the position. The actual standard error
* on the positions is increasing for epochs increasingly differing
* from the nominal J1991.25 epoch.
*Note on VarFlag: the values are
* 1: < 0.06mag ; 2: 0.06-0.6mag ; 3: >0.6mag
*Note on r_Vmag: the source is
* G = ground-based, H=HIP, T=Tycho
*Note on AstroRef: this flag indicates that the astrometric parameters in H3-4
* and H8-30 refer to:
* A to Z: the letter indicates the component of a double or multiple system
* *: the photocentre of a double or multiple system
* +: the centre of mass
*Note on F2: values exceeding +3 indicate a bad fit to the data.
*Note on m_BTmag: this flag indicates the component or combined photometry:
* A to Z : the letter indicates the component measured in Tycho
* (non-single star)
* * : the photometry refers to all components of the Hipparcos entry
* - : single-pointing triple or quadruple system
*Note on r_V-I: the origin of the V-I colour, in summary:
* 'A' for an observation of V-I in Cousins' system;
* 'B' to 'K' when V-I derived from measurements in other
* bands/photoelectric systems
* 'L' to 'P' when V-I derived from Hipparcos and Star Mapper photometry
* 'Q' for long-period variables
* 'R' to 'T' when colours are unknown
*Note on Hpmag, e_Hpmag:
* the Hipparcos magnitude could not be determined for 14 stars.
*Note on m_Hpmag: this flag indicates for double or multiple entries:
* A to Z : the letter indicates the specified component measured
* * : combined Hpmag of a double system, corrected for attenuation
* - : combined Hpmag of a multiple system, not corrected for attenuation
*Note on HvarType: Hipparcos-defined type of variability (a blank entry
* signifies that the entry could not be classified as variable or constant):
* C : no variability detected ("constant")
* D : duplicity-induced variability
* M : possibly micro-variable (amplitude < 0.03mag)
* P : periodic variable
* R : V-I colour index was revised due to variability analysis
* U : unsolved variable which does not fall in the other categories
*Note on moreVar: more data about periodic variability are provided
*Note on n_CCDM: the flag takes the following values:
* H : determined multiple by Hipparcos, previously unknown
* I : system previously identified as multiple in HIC <I/196> (annex1)
* M : miscellaneous (system identified after publication of HIC)
*Note on MultFlag: indicates that further details are given in the Double
* and Multiple Systems Annex:
* C : solutions for the components
* G : acceleration or higher order terms
* O : orbital solutions
* V : variability-induced movers (apparent motion arises from variability)
* X : stochastic solution (probably astrometric binaries with short period)
*Note on Source: qualifies the source of the astrometric parameters H8-30
* with a 'C' in MultFlag:
* P : primary target of a 2- or 3-pointing system
* F : secondary or tertiary of a 2- or 3-pointing 'fixed' system
* (common parallax and proper motions)
* I : secondary or tertiary of a 2- or 3-pointing 'independent' system
* (no constraints on parallax or proper motions)
* L : secondary or tertiary of a 2- or 3-pointing 'linear' system
* (common parallax)
* S : astrometric parameters from 'single-star merging' process.
*Note on Qual: Reliability of the double or multiple star solution:
* A=good, B=fair, C=poor, D=uncertain, S=suspected non-single
*Note on Chart: the chart was produced:
* D : from the STScI Digitized Sky Survey
* G : from the Guide Star Catalog
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on r_SpType: the flag indicates the source, as:
* 1 : Michigan catalogue for the HD stars, vol. 1 (Houk+, 1975) <III/31>
* 2 : Michigan catalogue for the HD stars, vol. 2 (Houk, 1978) <III/51>
* 3 : Michigan Catalogue for the HD stars, vol. 3 (Houk, 1982) <III/80>
* 4 : Michigan Catalogue for the HD stars, vol. 4 (Houk+, 1988) <III/133>
* G : updated after publication of the HIC <I/196>
* K : General Catalog of Variable Stars, 4th Ed. (Kholopov+ 1988) <II/214>
* S : SIMBAD data-base http://cdsweb.u-strasbg.fr/Simbad.html
* X : Miscellaneous
* A blank entry has no corresponding information.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h_dm_com.dat' ! Double and Multiples: Component solutions -COMP
integer*4 nr__1
parameter (nr__1=24588) ! Number of records
character*238 ar__1 ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=J1991.25)
character*10 CCDM_1 (nr__1) ! CCDM number (DC1)
integer*4 S (nr__1) ! *Solution identifier (DC2)
character*1 Type (nr__1) ! *[FIL]Type of solution (DC3)
character*1 Source_1 (nr__1) ! *[CFN] Solution source (DC4)
character*1 Qual_1 (nr__1) ! *[ABCD] Solution quality (DC5)
character*1 Notes_1 (nr__1) ! *[DGPWXYZ] Existence of notes (DC6)
integer*4 Nsys_1 (nr__1) ! Number of solutions for the system (DCM1)
integer*4 Ncomp_1 (nr__1) ! Number of components in this solution (DCM2)
integer*4 Nparm (nr__1) ! Number of free parameters in solution (DCM3)
integer*4 Ncorr (nr__1) ! Number of correlation records (DCM4)
integer*4 seq (nr__1) ! Sequential component number (DCM6)
character*1 comp_id (nr__1) ! Component identifier (DC7)
integer*4 HIP_1 (nr__1) ! HIP number (DC8)
real*4 Hp (nr__1) ! (mag) Magnitude of component (DC9)
real*4 e_Hp (nr__1) ! (mag) Standard error of Hp magnitude (DC10)
real*4 BT (nr__1) ! (mag) ? Magnitude of component, BT (DC11)
real*4 e_BT (nr__1) ! (mag) ? Standard error of BT (DC12)
real*4 VT (nr__1) ! (mag) ? Magnitude of component, VT (DC13)
real*4 e_VT (nr__1) ! (mag) ? Standard error of VT (DC14)
real*8 RAdeg_1 (nr__1) ! (deg) alpha, degrees (ICRS, Epoch=J1991.25) (DC15)
real*8 DEdeg_1 (nr__1) ! (deg) delta, degrees (ICRS, Epoch=J1991.25) (DC16)
real*8 Plx_1 (nr__1) ! (mas) Trigonometric parallax (DC17)
real*8 pmRA_1 (nr__1) ! (mas/yr) Proper motion in mu_alpha.cos(delta)ICRS(DC18)
real*8 pmDE_1 (nr__1) ! (mas/yr) Proper motion in mu_delta in ICRS (DC19)
real*4 e_RAdeg_1 (nr__1) ! (mas) Standard error in RA*cos(DEdeg) (DC20)
real*4 e_DEdeg_1 (nr__1) ! (mas) Standard error in DE (DC21)
real*4 e_Plx_1 (nr__1) ! (mas) Standard error in Plx (DC22)
real*4 e_pmRA_1 (nr__1) ! (mas/yr) Standard error in pmRA (DC23)
real*4 e_pmDE_1 (nr__1) ! (mas/yr) Standard error in pmDE (DC24)
character*1 ref (nr__1) ! Reference component for following data (DC25)
real*8 theta_1 (nr__1) ! (deg) *? Position angle (DC26)
real*8 rho_1 (nr__1) ! (arcsec) ? Separation from reference component (DC27)
real*8 d_theta (nr__1) ! (deg/yr) ? Rate of change of theta (DC28)
real*4 d_rho (nr__1) ! (arcsec/yr) ? Rate of change of separation (DC29)
integer*4 seq_ref (nr__1) ! *Sequential record number (DCM7)
character*25 flag (nr__1) ! *Status flags for parameters (DCM8)
*Note on S: a digit identifies different solutions pertaining to the
* same CCDM number.
*Note on Type: Summary of double or multiple star solution:
* F: fixed double or multiple system
* (identical proper motions and parallaxes)
* I: individual parallaxes and linear (relative) motion
* (possible optical double star)
* L: linear double or multiple system
* (may have different proper motions but assumed to have same parallax)
*Note on Source: the source of the solution is given by this flag:
* C: combined FAST and NDAC solution
* F: solution taken from the FAST Consortium only
* N: solution taken from the NDAC Consortium only
*Note on Qual: Reliability of the double or multiple star solution:
* A=good B=fair C=poor D=uncertain
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on theta: position angle relative to reference component
*Note on seq_ref: Sequential record number for the reference component, this
* field is set to zero if DC25 is blank.
*Note on flag: Status flags for Hp, RAdeg, DEdeg, Plx, pmRA, pmDE.
* 1 = estimated,
* 0 = constrained to the value of the first component.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'h_dm_cor.dat' ! Double and Multiples: Component solutions -CORR
integer*4 nr__2
parameter (nr__2=12591) ! Number of records
character*238 ar__2 ! Full-size record
character*10 CCDM_2 (nr__2) ! CCDM number (DC1)
integer*4 S_1 (nr__2) ! *Solution identifier (DC2)
character*1 Type_1 (nr__2) ! *[FIL]Type of solution (DC3)
character*1 Source_2 (nr__2) ! *[CFN] Solution source (DC4)
character*1 Qual_2 (nr__2) ! *[ABCD] Solution quality (DC5)
character*1 Notes_2 (nr__2) ! *[DGPWXYZ] Existence of notes (DC6)
integer*4 Nsys_2 (nr__2) ! Number of solutions for the system (DCM1)
integer*4 Ncomp_2 (nr__2) ! Number of components in this solution (DCM2)
integer*4 Nparm_1 (nr__2) ! Number of free parameters in solution (DCM3)
integer*4 Ncorr_1 (nr__2) ! Number of correlation records (DCM4)
integer*4 seq_1 (nr__2) ! Sequential component number (DCM6)
integer*4 corr (66,nr__2) ! *[-99/999]? Correlation coefficients (DCM7)
*Note on S: a digit identifies different solutions pertaining to the
* same CCDM number.
*Note on Type: Summary of double or multiple star solution:
* F: fixed double or multiple system
* (identical proper motions and parallaxes)
* I: individual parallaxes and linear (relative) motion
* (possible optical double star)
* L: linear double or multiple system
* (may have different proper motions but assumed to have same parallax)
*Note on Source: the source of the solution is given by this flag:
* C: combined FAST and NDAC solution
* F: solution taken from the FAST Consortium only
* N: solution taken from the NDAC Consortium only
*Note on Qual: Reliability of the double or multiple star solution:
* A=good B=fair C=poor D=uncertain
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on corr: the correlation records contain the correlation coefficients among
* the maximum set of 6N_C parameters. The number of correlation coefficients
* is 6N_C(6N_C-1)/2=66, 153 and 276 for N_C=2,3,4. A correlation record
* contains up to 66 coefficients; thus respectively 1,3 and 5 records are
* needed for a double, triple or quadruple star. The correlation coefficients
* are coded as integers between -99 and 999 using the arcsine
* transformation. The order of the correlation coefficients is indicated by
* the following table:
* Hp1 RA1 DE1 Plx1 pmRA1 pmDE1 Hp2 RA2 DE2 Plx2 pmRA2 pmDE2 Hp3 ...
* Hp1 1 r1 r2 r4 r7 r11 r16 r22 r29 r37 r46 r56 r67 ...
* RA1 r1 1 r3 r5 r8 r12 r17 r23 r30 r38 r47 r57 r68 ...
* DE1 r2 r3 1 r6 r9 r13 r18 r24 r31 r39 r48 r58 r69 ...
* Plx1 r4 r5 r6 1 r10 r14 r19 r25 r32 r40 r49 r59 r70 ...
* pmRA1 r7 r8 r9 r10 1 r15 r20 r26 r33 r41 r50 r60 r71 ...
* pmDE1 r11 r12 r13 r14 r15 1 r21 r27 r34 r42 r51 r61 r72 ...
* Hp2 r16 r17 r18 r19 r20 r21 1 r28 r35 r43 r52 r62 r73 ...
* RA2 r22 r23 r24 r25 r26 r27 r28 1 r36 r44 r53 r63 r74 ...
* DE2 r29 r30 r31 r32 r33 r34 r35 r36 1 r45 r54 r64 r75 ...
* Plx2 r37 r38 r39 r40 r41 r42 r43 r44 r45 1 r55 r65 r76 ...
* pmRA2 r46 r47 r48 r49 r50 r51 r52 r53 r54 r55 1 r66 r77 ...
* pmDE2 r56 r57 r58 r59 r60 r61 r62 r63 r64 r65 r66 1 r78 ...
* Hp3 r67 r68 r69 r70 r71 r72 r73 r74 r75 r76 r77 r78 1 ...
* ... ... ... ... ... ... ... ... ... ... ... ... ... ... ...
* These correlation coefficients are written as 66I3
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_dm_g.dat' ! Double and Multiples: Acceleration solutions
integer*4 nr__3
parameter (nr__3=2622) ! Number of records
character*195 ar__3 ! Full-size record
integer*4 HIP_2 (nr__3) ! Identifier (HIP number) (DG1)
real*8 gRA (nr__3) ! (mas/yr2) *Acceleration gRA = d(pmRA)/dt (DG2)
real*8 gDE (nr__3) ! (mas/yr2) *Acceleration gDE = d(pmDE)/dt (DG3)
real*8 e_gRA (nr__3) ! (mas/yr2) Standard error of gRA (DG4)
real*8 e_gDE (nr__3) ! (mas/yr2) Standard error of gDE (DG5)
real*4 Fg (nr__3) ! *Significance of the g terms (DG6)
real*8 dgRA (nr__3) ! (mas/yr3) *? dgRA = d2(pmRA)/dt2 (DG7)
real*8 dgDE (nr__3) ! (mas/yr3) *? dgDE = d2(pmDE)/dt2 (DG8)
real*8 e_dgRA (nr__3) ! (mas/yr3) ? Standard error of dgRA (DG9)
real*8 e_dgDE (nr__3) ! (mas/yr3) ? Standard error of dgDE (DG10)
real*4 Fdg (nr__3) ! *? Signifance of the dg terms (DG11)
character*1 Notes_3 (nr__3) ! *[DGPWXYZ] Existence of notes (DG12)
integer*4 num (nr__3) ! *[7,9] Number of astrometric parameters (DGM1)
integer*4 corr_1 (36,nr__3) ! *[-99/999]? Correlation coefficients (DGM2)
*Note on gRA: component in right ascension of the apparent acceleration
* of the photocentre at epoch J1991.25.
*Note on gDE: component in declination of the apparent acceleration
* of the photocentre at epoch J1991.25.
*Note on Fg: the quadratic model is only adopted if the g terms
* are significant (Fg>3.44)
*Note on dgRA: component in right ascension of the rate of change of the
* apparent acceleration of the photocentre at epoch J1991.25.
*Note on dgDE: component in declination of the rate of change of the
* apparent acceleration of the photocentre at epoch J1991.25.
*Note on Fdg: The cubic model is only adopted if the dg terms are
* significant (Fdg>3.44)
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on num: for a quadratic model of the photocentric motion this
* parameter is 7, for a cubic model it is 9.
*Note on corr: the complete set of n(n-1)*0.5 correlation coefficients
* (where n=7 for a quadratic and n=9 for a cubic model of the
* photocentric motion) is given in the order indicated by the
* following table:
* RA Dec Plx pmRA pmDE gRA gDE dgRA dgDE
* RA 1 r1 r2 r4 r7 r11 r16 r22 r29
* Dec r1 1 r3 r5 r8 r12 r17 r23 r30
* Plx r2 r3 1 r6 r9 r13 r18 r24 r31
* pmRA r4 r5 r6 1 r10 r14 r19 r25 r32
* pmDE r7 r8 r9 r10 1 r15 r20 r26 r33
* gRA r11 r12 r13 r14 r15 1 r21 r27 r34
* gDE r16 r17 r18 r19 r20 r21 1 r28 r35
* dgRA r22 r23 r24 r25 r26 r27 r28 1 r36
* dgDE r29 r30 r31 r32 r33 r34 r35 r36 1
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_dm_o.dat' ! Double and Multiples: Orbital solutions
integer*4 nr__4
parameter (nr__4=235) ! Number of records
character*337 ar__4 ! Full-size record
integer*4 HIP_3 (nr__4) ! Identifier (HIP) (D01)
real*8 P (nr__4) ! (d) Orbital period (DO2)
real*8 T (nr__4) ! (d) *Time of periastron passage (DO3)
real*8 a0 (nr__4) ! (mas) Semi-major axis of photocentric orbit (DO4)
real*4 ecc (nr__4) ! [0,1] Eccentricity (DO5)
real*4 w (nr__4) ! (deg) *[0,360] Argument of periastron (DO6)
real*4 i (nr__4) ! (deg) *[0,180] Inclination (DO7)
real*4 Omega (nr__4) ! (deg) *[0,360] Position angle of the node (DO8)
real*8 e_P (nr__4) ! (d) ? Standard error of P (DO9)
real*8 e_T (nr__4) ! (d) ? Standard error of T (DO10)
real*4 e_a0 (nr__4) ! (mas) ? Standard error of a0 (DO11)
real*4 e_ecc (nr__4) ! ? Standard error of ecc (DO12)
real*4 e_w (nr__4) ! (deg) ? Standard error of w (DO13)
real*4 e_i (nr__4) ! (deg) ? Standard error of i (DO14)
real*4 e_Omega (nr__4) ! (deg) ? Standard error of Omega (DO15)
integer*4 dmRef (nr__4) ! *? reference to the literature (DO16)
character*1 Notes_4 (nr__4) ! *[DGPWXYZ] Existence of notes (DO17)
character*12 flag_1 (nr__4) ! *Status flags for the parameters (DOM1)
integer*4 corr_2 (66,nr__4) ! *[-99/999]?=450 Correlation coefficients (DOM2)
*Note on T: this is the date when the photocentre is closest to the
* centre of mass in the orbital plane. It is equivalent to the
* closest approach of the stellar components.
*Note on w: this is the angle in the orbital plane from the line of
* nodes to the major axis, measured from the nodal point (DO8) to
* the periastron in the direction of motion.
*Note on i: the inclination of the orbital plane to the tangent plane
* of the sky. Taken to be in the first quadrant if the apparent
* motion is direct (counter-clockwise) and in the second quadrant
* for retrograde (clockwise) apparent motion.
*Note on Omega: This is the position angle (measured counter-clockwise
* as seen on the sky from the +Dec direction) of the line of nodes,
* or the intersection of the orbital and tangent planes. If the
* radial motion of the components is known from spectroscopic studies,
* then Omega should give the position angle of the ascending node,
* at which the primary star crosses the tangent plane while receding
* from the observer. In the absence of spectroscopic information
* Omega refers to the mode with the smallest positive position angle.
*Note on dmRef: A reference number in this field points to references in the
* printed catalogue, explained in dmsa_o.doc file.
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on flag: status flags for the 12 astrometric and orbital
* parameters taken in the order indicated below.
* 1= estimated; 0= not estimated.
*Note on corr: the correlation coefficients in bytes 140-337 are given
* in the table below. Correlation coefficients which are undefined
* (corresponding to blanks in DO9-15, and a status flag=0 in
* DOB127) are set to zero. All correlation coefficients are coded
* as integers between -99 and 999 using the arcsine transformation.
* Undefined values are coded as 450.
* RA Dec Plx pmRA pmDE P T a0 e w i Omega
* RA 1 r1 r2 r4 r7 r11 r16 r22 r29 r37 r46 r56
* Dec r1 1 r3 r5 r8 r12 r17 r23 r30 r36 r47 r57
* Plx r2 r3 1 r6 r9 r13 r18 r24 r31 r37 r46 r58
* pmRA r4 r5 r6 1 r10 r14 r19 r25 r32 r38 r47 r59
* pmDE r7 r8 r9 r10 1 r15 r20 r26 r33 r39 r48 r60
* P r11 r12 r13 r14 r15 1 r21 r27 r34 r40 r49 r61
* T r16 r17 r18 r19 r20 r21 1 r28 r35 r41 r50 r62
* a0 r22 r23 r24 r25 r26 r27 r28 1 r36 r42 r51 r63
* e r29 r30 r31 r32 r33 r34 r35 r36 1 r43 r52 r64
* w r37 r38 r39 r40 r41 r42 r43 r44 r45 1 r53 r65
* i r46 r47 r48 r49 r50 r51 r52 r53 r54 r55 1 r66
* Omega r56 r57 r59 r59 r60 r61 r62 r63 r64 r65 r66 1
* These correlation coefficients are written as 66I3
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_dm_v.dat' ! Double and Multiples: VIM solutions
integer*4 nr__5
parameter (nr__5=288) ! Number of records
character*144 ar__5 ! Full-size record
integer*4 HIP_4 (nr__5) ! Identifier (HIP) (DV1)
real*4 Hp_ref (nr__5) ! (mag) *Reference magnitude (DV2)
real*8 DRA (nr__5) ! (mas) VIM element in right ascension (DV3)
real*8 DDE (nr__5) ! (mas) VIM element in declination (DV4)
real*8 e_DRA (nr__5) ! (mas) Standard error of DRA (DV5)
real*8 e_DDE (nr__5) ! (mas) Standard error of DDE (DV6)
real*4 FD (nr__5) ! *Significance of the VIM elements (DV7)
real*4 theta_C (nr__5) ! (deg) *Position angle of the constant component (DV8)
real*4 e_theta_C (nr__5) ! (deg) Standard error of theta_C (DV9)
real*4 minSep (nr__5) ! (mas) Lower limit for separation of binary (DV10)
real*4 dvar (nr__5) ! (mas) *Displacement of photocentre (DV11)
character*1 Notes_5 (nr__5) ! *[DGPWXYZ] Existence of notes (DV12)
integer*4 corr_3 (21,nr__5) ! *[-99/999]Correlation coefficients (DVM1)
*Note on Hp_ref: the reference magnitude is freely chosen and defines the
* reference point for the object. The positional parameters RA and Dec (H8,
* H9 in hip_main.dat) and the VIM elements DV3,DV4 depend on the chosen
* Hp_ref. DV7-11 do not depend on Hp_ref.
*Note on FD: The VIM solution is only accepted if FD>2.15
*Note on theta_C: position angle of the constant component of the binary with
* respect to the variable component, measured counterclockwise as seen on
* the sky from the +Dec direction.
*Note on dvar: displacement of photocentre between minimum and maximum
* luminosity of the system. Indicates the size of the VIM effect.
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on corr: the correlation coefficients in bytes 82-144 are given in the
* following sequence:
* RA Dec Plx pmRA pmDE DRA DDE
* RA 1 r1 r2 r4 r7 r11 r16
* Dec r1 1 r3 r5 r8 r12 r17
* Plx r2 r3 1 r6 r9 r13 r18
* pmRA r4 r5 r6 1 r10 r14 r19
* pmDE r7 r8 r9 r10 1 r15 r20
* DRA r11 r12 r13 r14 r15 1 r21
* DDE r16 r17 r18 r19 r20 r21 1
* These correlation coefficients are written as 21I3
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_dm_x.dat' ! Double and Multiples: Stochastic solutions
integer*4 nr__6
parameter (nr__6=1561) ! Number of records
character*22 ar__6 ! Full-size record
integer*4 HIP_5 (nr__6) ! Identifier (HIP) (DX1)
real*4 epsilon (nr__6) ! (mas) Cosmic error, epsilon (DX2)
real*4 e_epsilon (nr__6) ! (mas) Standard error of epsilon (DX3)
character*1 Notes_6 (nr__6) ! *[DGPWXYZ] Existence of notes (DX4)
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_va_1.dat' ! Variability Annex: Periodic variables
integer*4 nr__7
parameter (nr__7=2712) ! Number of records
character*142 ar__7 ! Full-size record
integer*4 HIP_6 (nr__7) ! Identifier (HIP) (P1)
character*1 NewVar (nr__7) ! [*] Flag if new variable (P2)
character*12 SpType_1 (nr__7) ! Spectral type (P3)
character*1 HvarType_1 (nr__7) ! *[CDMPRU] Variability type (1-letter) (P4)
character*6 VarType (nr__7) ! *Variability type as in GCVS/NSV (P5)
character*1 n_VarType (nr__7) ! [*] Flag if newly classified by Hipparcos (P6)
real*4 maxMag (nr__7) ! (mag) Magnitude at max from curve fitting (P7)
character*1 l_minMag (nr__7) ! *[>] Limit flag (>) (P8)
real*4 minMag (nr__7) ! (mag) Magnitude at min from curve fitting (P9)
real*4 log_sA_A (nr__7) ! ? log_10(sigma_A / A) (P10)
real*8 Period_1 (nr__7) ! (d) ? Mean period in days (P11)
real*4 log_sP (nr__7) ! ([d]) ? log_10(sigma_P) (P12)
real*8 Ep_2440000 (nr__7) ! (d) ? Epoch (JD-2440000) of zero phase (P13)
integer*4 q_Ep_2440000(nr__7) ! *[0/5]? Precision flag (P14)
character*1 morePhoto_1(nr__7) ! [ABC] Light curve Annex (P15)
character*1 Notes_7 (nr__7) ! *[DGPWXYZ] Existence of notes (P16)
character*12 VarName (nr__7) ! Variable star name (P17)
real*8 period_2 (nr__7) ! (d) ? Period from literature (P18)
real*8 epoch (nr__7) ! (d) ? Epoch from literature (P19)
real*4 max (nr__7) ! (mag) ? Magnitude at max from literature (P20)
real*4 min (nr__7) ! (mag) ? Magnitude at min from literature (P21)
character*1 Band (nr__7) ! *[UBVKIRPYb] Photometric band (P22)
character*1 refFlag (nr__7) ! [R] Reference in printed catalogue (P23)
*Note on HvarType: Hipparcos-defined type of variability:
* C : no variability detected ("constant")
* D : duplicity-induced variability
* M : possibly micro-variable (amplitude < 0.03mag)
* P : periodic variable
* R : V-I colour index was revised due to variability analysis
* U : unsolved variable which does not fall in the other categories
*Note on VarType: this is the 6-letter type defined
* in the General Catalog of Variable Stars <II/214>
* (see also the "Types of Variability" section below)
*Note on l_minMag: the flag (>) indicates that the true magnitude at minimum
* luminosity is likely to be larger than the value of minMag given in P9.
*Note on q_Ep-2440000: the value is 1-log10(sigma_epoch), i.e.
* 0 for an accuracy of about 10 days,
* 1 for an accuracy of about 1 day,
* 2 for an accuracy of about 0.1day,
* 3 for an accuracy of about 0.01day,
* 4 for an accuracy of about 0.001day,
* 5 for an accuracy of about 0.0001day.
*Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on Band: U,B,V,K,I,R refer to Johnson system, or closely related
* systems; P refers to photographic magnitudes, Y and b for the
* Stroemgren y and b bands.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hip_va_2.dat' ! Variability Annex: Unsolved variables
integer*4 nr__8
parameter (nr__8=5542) ! Number of records
character*142 ar__8 ! Full-size record
integer*4 HIP_7 (nr__8) ! Identifier (HIP number) (U1)
character*1 NewVar_1 (nr__8) ! [*] Flag if new variable (U2)
character*12 SpType_2 (nr__8) ! Spectral type (U3)
character*1 HvarType_2 (nr__8) ! *[CDMPRU] Variability type (1-letter) (U4)
character*6 VarType_1 (nr__8) ! *Variability type (6-letter) as in GCVS/NSV(U5)
character*1 n_VarType_1(nr__8) ! [*] Flag if newly classified by Hipparcos (U6)
real*4 maxMag_1 (nr__8) ! (mag) *Magnitude at max from curve fitting (U7)
character*1 l_minMag_1 (nr__8) ! [>] Limit flag (>) (U8)
real*4 minMag_1 (nr__8) ! (mag) *Magnitude at min from curve fitting (U9)
real*4 med_Hp (nr__8) ! (mag) Median Hp (U10)
real*8 A (nr__8) ! (mag) Intrinsic variability amplitude (U11)
real*4 e_A (nr__8) ! (mag) Standard error of A (U12)
character*1 morePhoto_2(nr__8) ! [ABC] Light curve Annex (U15)
character*1 Notes_8 (nr__8) ! *[DGPWXYZ] Notes (U16)
character*12 VarName_1 (nr__8) ! Variable star name (U17)
real*8 period_3 (nr__8) ! (d) ? Period from literature (U18)
real*8 Ep_2440000_1(nr__8) ! (d) ? Epoch (JD-2440000) from literature (U19)
real*4 max_1 (nr__8) ! (mag) ? Magnitude at max from literature (U20)
real*4 min_1 (nr__8) ! (mag) ? Magnitude at min from literature (U21)
character*1 Band_1 (nr__8) ! *[UBVKIRPYb] Photometric band (U22)
character*1 refFlag_1 (nr__8) ! [R] Reference in printed catalogue (U23)
*Note on HvarType: Hipparcos-defined type of variability:
* C : no variability detected ("constant")
* D : duplicity-induced variability
* M : possibly micro-variable (amplitude < 0.03mag)
* P : periodic variable
* R : V-I colour index was revised due to variability analysis
* U : unsolved variable which does not fall in the other categories
*Note on VarType: this is the 6-letter type defined
* in the General Catalog of Variable Stars <II/214>
* (see also the "Types of Variability" section below)
*Note on maxMag: this is given if significant
*Note on minMag: this is given if significant
*Note on Notes: Note on Notes: the flag has the following meaning:
* D : double and multiple systems note only (note in hd_notes.doc file)
* G : general note only (note in hg_notes.doc file)
* P : photometric notes only (note in hp_notes.doc file)
* W : D + P
* X : D + G
* Y : G + P
* Z : D + G + P
*Note on Band: U,B,V,K,I,R refer to Johnson system, or closely related
* systems; P refers to photographic magnitudes, Y and b for the
* Stroemgren y and b bands.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'solar_ha.dat' ! Solar System Annex: Astrometric catalogue
integer*4 nr__9
parameter (nr__9=5609) ! Number of records
character*64 ar__9 ! Full-size record
C Position composed of: RAdeg DEdeg
integer*4 ID (nr__9) ! Object number (SHA1)
real*8 RAdeg_2 (nr__9) ! (deg) Reference point RA (ICRS system) (SHA2)
real*8 DEdeg_2 (nr__9) ! (deg) Reference point Dec (ICRS system) (SHA3)
real*8 Ep_2440000_2(nr__9) ! (d) *Measurement epoch (SHA4)
real*4 delay (nr__9) ! (s) *Light delay time (SHA5)
real*8 theta_2 (nr__9) ! (deg) *Position angle, theta (SHA6)
real*4 e_lambda (nr__9) ! (mas) Estimated standard error of abscissa (SHA7)
integer*4 flag_2 (nr__9) ! *[1,2]FAST or NDAC flag (SHA8)
*Note on Ep-2440000: the measurement epoch is specified in JD with respect to
* JD(TT)2440000.0 and is corrected to the geocentre.
*Note on delay: this gives the applied light time delay in the geocentric
* direction of the observed object between the satellite and the Earth.
*Note on theta: The position angle of the slit coordinate direction w is
* reckoned positive from North through East.
*Note on flag: If the transit corresponds to an NDAC record the flag is 1.
* For a FAST record the flag is 2.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'solar_hp.dat' ! Solar System Annex: Photometric catalogue
integer*4 nr__10
parameter (nr__10=2639) ! Number of records
character*63 ar__10 ! Full-size record
integer*4 ID_1 (nr__10) ! Object number (SHP1)
real*8 Ep_2440000_3(nr__10) ! (d) *Measurement epoch (SHP2)
real*8 Hp_dc (nr__10) ! (mag) Magnitude from unmodulated signal (SHP3)
real*4 e_Hp_dc (nr__10) ! (mag) Standard error on Hp_dc (SHP4)
real*8 Hp_ac (nr__10) ! (mag) Magnitude from modulated signal (SHP5)
real*4 e_Hp_ac (nr__10) ! (mag) Standard error on Hp_ac (SHP6)
real*4 r (nr__10) ! (AU) Distance: Sun-asteroid (SHP7)
real*4 Delta (nr__10) ! (AU) Distance: satellite-asteroid (SHP8)
real*4 alpha (nr__10) ! (deg) *Solar phase angle (SHP9)
*Note on Ep-2440000: the measurement epoch is specified in JD with respect
* to JD(TT)2440000.0 and is not corrected to the geocentre (i.e. the
* light-time delay between the satellite and the earth is also neglected.
*Note on alpha: the solar phase angle is the angle between the Sun and the
* satellite as viewed from the asteroid.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'solar_t.dat' ! Solar System Annex: Tycho astrometry/photometry
integer*4 nr__11
parameter (nr__11=291) ! Number of records
character*95 ar__11 ! Full-size record
C Position composed of: RAdeg DEdeg
integer*4 ID_2 (nr__11) ! Object number (ST1)
real*8 Ep_2440000_4(nr__11) ! (d) *Measurement epoch (ST2)
real*8 RAdeg_3 (nr__11) ! (deg) Right ascension (ICRS system) (ST3)
real*8 DEdeg_3 (nr__11) ! (deg) Declination (ICRS system) (ST4)
real*4 BTmag_1 (nr__11) ! (mag) ? Mean BT magnitude (ST5)
real*4 VTmag_1 (nr__11) ! (mag) ? Mean VT magnitude (ST6)
integer*4 Flag_3 (nr__11) ! *[1,2] Transit flag (ST7)
real*4 e_RAdeg_2 (nr__11) ! (mas) Standard error on RA (ST8)
real*4 e_DEdeg_2 (nr__11) ! (mas) Standard error on Dec (ST9)
real*4 DE_RA_1 (nr__11) ! Correlation, Dec/RA*cos(delta) (ST10)
real*4 theta_3 (nr__11) ! (deg) Position angle of slit (direction w) (ST11)
integer*4 sign_z (nr__11) ! *[-1/1] Inclined slit flag (ST12)
real*4 e_incl (nr__11) ! (mas) Standard error on slit pos. (inclined) (ST13)
real*4 e_vert (nr__11) ! (mas) Standard error on slit pos. (vertical) (ST14)
*Note on Ep-2440000: the measurement epoch is specified in JD with respect
* to JD(TT)2440000.0 and is not corrected to the geocentre (i.e. the
* light-time delay between the satellite and the earth is also neglected.
*Note on Flag: if flag =1 then only one crossing of the field of view has
* been detected or retained. If Flag n>1, then n candidate observed
* transits have been detected or retained for that predicted
* observation. ST1, ST2, ST7, ST11, ST12 are identical in those cases.
*Note on sign_z: if the transit occurred in the upper branch of the
* inclined slits, the flag is +1, if the transit occurred in the
* lower branch of the inclined slits, the flag is -1.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hd_notes.doc' ! Hipparcos notes: Double and multiple systems
integer*4 nr__12
parameter (nr__12=2622) ! Number of records
character*97 ar__12 ! Full-size record
integer*4 HIP_8 (nr__12) ! HIP number
character*1 Note1 (nr__12) ! [DG] Double and multiple or General note
character*1 Note2 (nr__12) ! [DP] Double and multiple or Photometric note
integer*4 Ntot (nr__12) ! Total number of lines for the HIP object
integer*4 Nline (nr__12) ! Running line number in range [1,Ntot]
character*80 Text (nr__12) ! Text of note.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hg_notes.doc' ! Hipparcos notes: General notes
integer*4 nr__13
parameter (nr__13=3898) ! Number of records
character*97 ar__13 ! Full-size record
integer*4 HIP_9 (nr__13) ! HIP number
character*1 Note1_1 (nr__13) ! [DG] Double and multiple or General note
character*1 Note2_1 (nr__13) ! [DP] Double and multiple or Photometric note
integer*4 Ntot_1 (nr__13) ! Total number of lines for the HIP object
integer*4 Nline_1 (nr__13) ! Running line number in range [1,Ntot]
character*80 Text_1 (nr__13) ! Text of note.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hp_notes.doc' ! Hipparcos notes: Photometric notes
integer*4 nr__14
parameter (nr__14=2444) ! Number of records
character*97 ar__14 ! Full-size record
integer*4 HIP_10 (nr__14) ! HIP number
character*1 Note1_2 (nr__14) ! [DG] Double and multiple or General note
character*1 Note2_2 (nr__14) ! [DP] Double and multiple or Photometric note
integer*4 Ntot_2 (nr__14) ! Total number of lines for the HIP object
integer*4 Nline_2 (nr__14) ! Running line number in range [1,Ntot]
character*80 Text_2 (nr__14) ! Text of note.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hp_refs.doc' ! References Hipparcos stars
integer*4 nr__15
parameter (nr__15=33769) ! Number of records
character*19 ar__15 ! Full-size record
integer*4 HIP_11 (nr__15) ! HIP number
integer*4 Ntot_3 (nr__15) ! Total number of lines for the HIP object
integer*4 Nline_3 (nr__15) ! Running line number in range [1,Ntot]
real*4 nRef (nr__15) ! Reference number (explained in hp_auth.doc)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hp_auth.doc' ! References of hp_notes.doc
integer*4 nr__16
parameter (nr__16=4335) ! Number of records
character*80 ar__16 ! Full-size record
real*4 nRef_1 (nr__16) ! Reference number from hp_refs.doc
character*70 Text_3 (nr__16) ! Text of reference
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'dmsa_o.doc' ! References of hip_dm_o.dat
integer*4 nr__17
parameter (nr__17=118) ! Number of records
character*80 ar__17 ! Full-size record
integer*4 dmRef_1 (nr__17) ! Reference number (hip_dm_o.dat file)
integer*4 Ntot_4 (nr__17) ! Total number of lines for the reference
integer*4 Nline_4 (nr__17) ! Running line number in range [1,Ntot]
character*72 Text_4 (nr__17) ! Text of reference
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tyc_main.dat' ! The main part of Tycho Catalogue
integer*4 nr__18
parameter (nr__18=1058332) ! Number of records
character*350 ar__18 ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=J1991.25)
character*1 Catalog_1 (nr__18) ! [T] Catalogue (T = Tycho) (T0)
character*12 TYC (nr__18) ! *TYC1-3 (TYC number) (T1)
character*1 Proxy_1 (nr__18) ! [HT]? Proximity flag (T2)
character*11 RAhms_1 (nr__18) ! Right ascension in h m s, ICRS (J1991.25) (T3)
character*11 DEdms_1 (nr__18) ! Declination in deg ' ", ICRS (J1991.25) (T4)
real*4 Vmag_1 (nr__18) ! (mag) ? Magnitude in Johnson V (T5)
character*1 r_Vmag_1 (nr__18) ! *[BDTV] Source of magnitude (T7)
real*8 RAdeg_4 (nr__18) ! (deg) *alpha, degrees (ICRS, Epoch=J1991.25) (T8)
real*8 DEdeg_4 (nr__18) ! (deg) *delta, degrees (ICRS, Epoch=J1991.25) (T9)
character*1 AstroRef_1 (nr__18) ! *[X]? Reference flag for astrometry (T10)
real*8 Plx_2 (nr__18) ! (mas) *? Trigonometric parallax (T11)
real*8 pmRA_2 (nr__18) ! (mas/yr) *? Proper motion mu_alpha.cos(delta), ICRS(T12)
real*8 pmDE_2 (nr__18) ! (mas/yr) *? Proper motion mu_delta, ICRS (T13)
real*4 e_RAdeg_3 (nr__18) ! (mas) *? Standard error in RA*cos(delta) (T14)
real*4 e_DEdeg_3 (nr__18) ! (mas) *? Standard error in DE (T15)
real*4 e_Plx_2 (nr__18) ! (mas) *? Standard error in Plx (T16)
real*4 e_pmRA_2 (nr__18) ! (mas/yr) *? Standard error in pmRA (T17)
real*4 e_pmDE_2 (nr__18) ! (mas/yr) *? Standard error in pmDE (T18)
real*4 DE_RA_2 (nr__18) ! [-1/1]? Correlation, DE/RA*cos(delta) (T19)
real*4 Plx_RA_1 (nr__18) ! [-1/1]? Correlation, Plx/RA*cos(delta) (T20)
real*4 Plx_DE_1 (nr__18) ! [-1/1]? Correlation, Plx/DE (T21)
real*4 pmRA_RA_1 (nr__18) ! [-1/1]? Correlation, pmRA/RA*cos(delta) (T22)
real*4 pmRA_DE_1 (nr__18) ! [-1/1]? Correlation, pmRA/DE (T23)
real*4 pmRA_Plx_1 (nr__18) ! [-1/1]? Correlation, pmRA/Plx (T24)
real*4 pmDE_RA_1 (nr__18) ! [-1/1]? Correlation, pmDE/RA*cos(delta) (T25)
real*4 pmDE_DE_1 (nr__18) ! [-1/1]? Correlation, pmDE/DE (T26)
real*4 pmDE_Plx_1 (nr__18) ! [-1/1]? Correlation, pmDE/Plx (T27)
real*4 pmDE_pmRA_1(nr__18) ! [-1/1]? Correlation, pmDE/pmRA (T28)
integer*4 Nastro (nr__18) ! ? Number of transits for astrometry (T29)
real*4 F2_1 (nr__18) ! *? Goodness-of-fit parameter (T30)
integer*4 HIP_12 (nr__18) ! ? Hipparcos HIP number (T31)
real*4 BTmag_2 (nr__18) ! (mag) ? Mean BT magnitude (T32)
real*4 e_BTmag_1 (nr__18) ! (mag) ? Standard error in BTmag (T33)
real*4 VTmag_2 (nr__18) ! (mag) ? Mean VT magnitude (T34)
real*4 e_VTmag_1 (nr__18) ! (mag) ? Standard error in VTmag (T35)
character*1 r_BTmag (nr__18) ! *[DMNT] Source of photometry (T36)
real*4 B_V_1 (nr__18) ! (mag) ? Johnson B-V colour (T37)
real*4 e_B_V_1 (nr__18) ! (mag) ? Standard error on B-V (T38)
integer*4 Q (nr__18) ! *? Astrometric quality flag, Q (T40)
real*4 Fs (nr__18) ! ? Signal-to-noise ratio of the star image(T41)
character*1 Source_3 (nr__18) ! *[HPR] Source of astrometric data (T42)
integer*4 Nphoto (nr__18) ! ? Number of transits for photometry (T43)
real*4 VTscat (nr__18) ! (mag) ? Estimate of VTmag scatter (T44)
real*4 VTmax (nr__18) ! (mag) ? VTmag at maximum (15th percentile) (T45)
real*4 VTmin (nr__18) ! (mag) ? VTmag at minimum (85th percentile) (T46)
character*1 Var (nr__18) ! *[GN]? Known variability from GCVS/NSV (T47)
character*1 VarFlag_1 (nr__18) ! *[UVW]? Variability from Tycho (T48)
character*1 MultFlag_1 (nr__18) ! *[DRSYZ]? Duplicity from Tycho (T49)
character*1 morePhoto_3(nr__18) ! [AB] Epoch photometry in Annex A or B (T50)
character*2 m_HIP_1 (nr__18) ! CCDM component identifier (T51)
integer*4 PPM (nr__18) ! *[1/789676]? PPM and Supplement (T52)
integer*4 HD_1 (nr__18) ! [1/359083]? HD cat. <III/135> (T53)
character*10 BD_1 (nr__18) ! Bonner DM <I/119>, <I/122> (T54)
character*10 CoD_1 (nr__18) ! Cordoba DM <I/114> (T55)
character*10 CPD_1 (nr__18) ! Cape Photographic DM <I/108> (T56)
character*1 Remark (nr__18) ! *[JKLM] Notes (T57)
*Note on TYC: the designation of an object in the Tycho Catalogue uses the
* Guide Star Catalog numbering system (a region number (TYC1) and a number
* within the region (TYC2)) followed by a Tycho specific component number
* (TYC3).
*Note on r_Vmag: if non blank, the field has the following meaning:
* B : no VTmag available, therefore BTmag was adopted
* D : derived from approximate BTmag and VTmag (r_BTmag field T36 is 'D')
* T : derived from approximate VTmag (r_BTmag field T36 is 'T')
* V : no BTmag available, therefore VTmag was adopted in Vmag
*Note on RAdeg, DEdeg: right ascension and declination are
* expressed in degrees for epoch J1991.25 (JD2448349.0625 (TT)) in the
* ICRS (International Celestial Reference System, close to
* J2000) reference system.
*Note on AstroRef: 'X' indicates a dubious astrometric reference star
* in the context of the Tycho catalogue.
*Note on Plx, pmRA, pmDE, e_RAdeg, e_DEdeg, e_Plx, e_pmRA, e_pmDE:
* For these fields the second decimal digit is always blank
*Note on F2: values exceeding +2.5 to +3 indicate a bad fit to the data.
*Note on r_BTmag: the source flag has the following meaning:
* D : approximate BTmag, VTmag obtained for resolved double stars
* M : BTmag and VTmag are median values rather than de-censored means
* N : BTmag and VTmag are de-censored means
* T : BTmag is not given, and VTmag is an estimate; these magnitudes
* are systematically too bright by up to 1mag.
* A blank indicates an Hipparcos star not observed by Tycho (T42 = H).
*Note on Q: the astrometric quality flag is defined by the following table,
* where sigma(max) is the largest of the 5 astrometric standard errors:
* ---------------------------------------------------------
* Q sigma(max) Astrometric quality
* ---------------------------------------------------------
* 1 < 5 very high
* 2 5- 10 very high
* 3 10- 25 high
* 4 25- 50 high
* 5 50-150 medium
* 6 < 150 perhaps non-single
* 7 < 150 low
* 8 < 150 perhaps non-stellar
* 9 ~ 200 low, position derived from TICR
* ('R' in Source field T42)
* blank --- unassigned, 'H' in Source field T42
* (Hipparcos star not observed by Tycho)
* ---------------------------------------------------------
*Note on Source:
* H : Hipparcos star not observed by Tycho
* P : only the position was determined (no proper motion, no parallax)
* R : the position is derived from TICR catalog <I/197>
*Note on Var: this flag is set when the variability is known:
* G : variable known in General Catalog of Variable Stars <II/214>
* N : variable known in New Suspected Variables catalog <II/140>
*Note on VarFlag: this flag has the following meaning:
* U : apparent variability in the Tycho data; may be due to duplicity
* V : strong evidence of intrinsic variability
* W : suspected intrinsic variability
*Note on MultFlag: unresolved duplicity status from Tycho data analysis:
* D : duplicity clearly indicated (BTmag and VTmag refer to combined light)
* R : duplicity weakly indicated, combined with indication of variability
* S : duplicity suspected
* Y : investigation for duplicity carried out on Tycho data,
* no indication of duplicity was found
* Z : investigation for duplicity not carried out
* blank : Hipparcos star not observed by Tycho
*Note on PPM: from the 3 parts of the PPM catalogue
* North <I/146>, South <I/193>, and 90000Supplement <I/208>
*Note on Remark: if non-blank, the remark has the following meaning:
* J : disagreement with position of magnitude in GSC1.1 catalog <I/220>
* K : dubious Tycho parallax (Plx)
* L : dubious Tycho proper motion (disagrees with PPM catalogue)
* M : very uncertain Tycho magnitude (standard error larger than 0.3mag)
C=============================================================================
C Loading file 'hip_main.dat' ! The Hipparcos Main Catalogue
C Format for file interpretation
1 format(
+ A1,7X,I6,1X,A1,1X,A11,1X,A11,1X,F5.2,1X,I1,1X,A1,1X,F12.8,1X,
+ F12.8,1X,A1,1X,F7.2,1X,F8.2,1X,F8.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F5.2,8X,F6.3,1X,
+ F5.3,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1,1X,F4.2,1X,
+ F4.2,1X,A1,1X,A1,1X,F7.4,1X,F6.4,1X,F5.3,1X,I3,1X,A1,1X,F5.2,
+ 1X,F5.2,1X,F7.2,1X,A1,1X,A1,1X,A1,1X,A10,1X,A1,1X,I2,1X,I2,1X,
+ A1,1X,A1,1X,A1,1X,A2,1X,I3,1X,F7.3,1X,F5.3,1X,F5.2,1X,F4.2,1X,
+ A1,1X,A1,1X,A1,1X,I6,1X,A10,1X,A10,1X,A10,1X,F4.2,1X,A12,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'hip_main.dat')
write(6,*) '....Loading file: hip_main.dat'
do i__=1,118218
read(1,'(A450)')ar__
read(ar__,1)
+ Catalog(i__),HIP(i__),Proxy(i__),RAhms(i__),DEdms(i__),
+ Vmag(i__),VarFlag(i__),r_Vmag(i__),RAdeg(i__),DEdeg(i__),
+ AstroRef(i__),Plx(i__),pmRA(i__),pmDE(i__),e_RAdeg(i__),
+ e_DEdeg(i__),e_Plx(i__),e_pmRA(i__),e_pmDE(i__),DE_RA(i__),
+ Plx_RA(i__),Plx_DE(i__),pmRA_RA(i__),pmRA_DE(i__),
+ pmRA_Plx(i__),pmDE_RA(i__),pmDE_DE(i__),pmDE_Plx(i__),
+ pmDE_pmRA(i__),F1(i__),F2(i__),BTmag(i__),e_BTmag(i__),
+ VTmag(i__),e_VTmag(i__),m_BTmag(i__),B_V(i__),e_B_V(i__),
+ r_B_V(i__),V_I(i__),e_V_I(i__),r_V_I(i__),CombMag(i__),
+ Hpmag(i__),e_Hpmag(i__),Hpscat(i__),o_Hpmag(i__),m_Hpmag(i__),
+ Hpmax(i__),HPmin(i__),Period(i__),HvarType(i__),moreVar(i__),
+ morePhoto(i__),CCDM(i__),n_CCDM(i__),Nsys(i__),Ncomp(i__),
+ MultFlag(i__),Source(i__),Qual(i__),m_HIP(i__),theta(i__),
+ rho(i__),e_rho(i__),dHp(i__),e_dHp(i__),Survey(i__),
+ Chart(i__),Notes(i__),HD(i__),BD(i__),CoD(i__),CPD(i__),
+ v_V_I_red(i__),SpType(i__),r_SpType(i__)
if(ar__(42:46) .EQ. '') Vmag(i__) = rNULL__
if(ar__(48:48) .EQ. '') VarFlag(i__) = iNULL__
if(ar__(52:63) .EQ. '') RAdeg(i__) = rNULL__
if(ar__(65:76) .EQ. '') DEdeg(i__) = rNULL__
if(ar__(80:86) .EQ. '') Plx(i__) = rNULL__
if(ar__(88:95) .EQ. '') pmRA(i__) = rNULL__
if(ar__(97:104) .EQ. '') pmDE(i__) = rNULL__
if(ar__(106:111) .EQ. '') e_RAdeg(i__) = rNULL__
if(ar__(113:118) .EQ. '') e_DEdeg(i__) = rNULL__
if(ar__(120:125) .EQ. '') e_Plx(i__) = rNULL__
if(ar__(127:132) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__(134:139) .EQ. '') e_pmDE(i__) = rNULL__
if(ar__(141:145) .EQ. '') DE_RA(i__) = rNULL__
if(ar__(147:151) .EQ. '') Plx_RA(i__) = rNULL__
if(ar__(153:157) .EQ. '') Plx_DE(i__) = rNULL__
if(ar__(159:163) .EQ. '') pmRA_RA(i__) = rNULL__
if(ar__(165:169) .EQ. '') pmRA_DE(i__) = rNULL__
if(ar__(171:175) .EQ. '') pmRA_Plx(i__) = rNULL__
if(ar__(177:181) .EQ. '') pmDE_RA(i__) = rNULL__
if(ar__(183:187) .EQ. '') pmDE_DE(i__) = rNULL__
if(ar__(189:193) .EQ. '') pmDE_Plx(i__) = rNULL__
if(ar__(195:199) .EQ. '') pmDE_pmRA(i__) = rNULL__
if(ar__(201:203) .EQ. '') F1(i__) = iNULL__
if(ar__(205:209) .EQ. '') F2(i__) = rNULL__
if(ar__(218:223) .EQ. '') BTmag(i__) = rNULL__
if(ar__(225:229) .EQ. '') e_BTmag(i__) = rNULL__
if(ar__(231:236) .EQ. '') VTmag(i__) = rNULL__
if(ar__(238:242) .EQ. '') e_VTmag(i__) = rNULL__
if(ar__(246:251) .EQ. '') B_V(i__) = rNULL__
if(ar__(253:257) .EQ. '') e_B_V(i__) = rNULL__
if(ar__(261:264) .EQ. '') V_I(i__) = rNULL__
if(ar__(266:269) .EQ. '') e_V_I(i__) = rNULL__
if(ar__(275:281) .EQ. '') Hpmag(i__) = rNULL__
if(ar__(283:288) .EQ. '') e_Hpmag(i__) = rNULL__
if(ar__(290:294) .EQ. '') Hpscat(i__) = rNULL__
if(ar__(296:298) .EQ. '') o_Hpmag(i__) = iNULL__
if(ar__(302:306) .EQ. '') Hpmax(i__) = rNULL__
if(ar__(308:312) .EQ. '') HPmin(i__) = rNULL__
if(ar__(314:320) .EQ. '') Period(i__) = rNULL__
if(ar__(341:342) .EQ. '') Nsys(i__) = iNULL__
if(ar__(344:345) .EQ. '') Ncomp(i__) = iNULL__
if(ar__(356:358) .EQ. '') theta(i__) = iNULL__
if(ar__(360:366) .EQ. '') rho(i__) = rNULL__
if(ar__(368:372) .EQ. '') e_rho(i__) = rNULL__
if(ar__(374:378) .EQ. '') dHp(i__) = rNULL__
if(ar__(380:383) .EQ. '') e_dHp(i__) = rNULL__
if(ar__(391:396) .EQ. '') HD(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Catalog(i__),HIP(i__),Proxy(i__),RAhms(i__),DEdms(i__),
+ Vmag(i__),VarFlag(i__),r_Vmag(i__),RAdeg(i__),DEdeg(i__),
+ AstroRef(i__),Plx(i__),pmRA(i__),pmDE(i__),e_RAdeg(i__),
+ e_DEdeg(i__),e_Plx(i__),e_pmRA(i__),e_pmDE(i__),DE_RA(i__),
+ Plx_RA(i__),Plx_DE(i__),pmRA_RA(i__),pmRA_DE(i__),
+ pmRA_Plx(i__),pmDE_RA(i__),pmDE_DE(i__),pmDE_Plx(i__),
+ pmDE_pmRA(i__),F1(i__),F2(i__),BTmag(i__),e_BTmag(i__),
+ VTmag(i__),e_VTmag(i__),m_BTmag(i__),B_V(i__),e_B_V(i__),
+ r_B_V(i__),V_I(i__),e_V_I(i__),r_V_I(i__),CombMag(i__),
+ Hpmag(i__),e_Hpmag(i__),Hpscat(i__),o_Hpmag(i__),m_Hpmag(i__),
+ Hpmax(i__),HPmin(i__),Period(i__),HvarType(i__),moreVar(i__),
+ morePhoto(i__),CCDM(i__),n_CCDM(i__),Nsys(i__),Ncomp(i__),
+ MultFlag(i__),Source(i__),Qual(i__),m_HIP(i__),theta(i__),
+ rho(i__),e_rho(i__),dHp(i__),e_dHp(i__),Survey(i__),
+ Chart(i__),Notes(i__),HD(i__),BD(i__),CoD(i__),CPD(i__),
+ v_V_I_red(i__),SpType(i__),r_SpType(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h_dm_com.dat' ! Double and Multiples: Component solutions -COMP
C Format for file interpretation
2 format(
+ A10,1X,I1,1X,A1,1X,A1,1X,A1,1X,A1,1X,I1,1X,I2,1X,I2,1X,I2,6X,
+ I2,1X,A1,1X,I6,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F12.8,1X,F12.8,1X,F7.2,1X,F8.2,1X,F8.2,1X,F6.2,1X,
+ F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,A1,1X,F7.3,1X,F8.3,1X,F8.3,1X,
+ F6.3,1X,I2,1X,A25)
C Effective file loading
open(unit=1,status='old',file=
+'h_dm_com.dat')
write(6,*) '....Loading file: h_dm_com.dat'
do i__=1,24588
read(1,'(A238)')ar__1
read(ar__1,2)
+ CCDM_1(i__),S(i__),Type(i__),Source_1(i__),Qual_1(i__),
+ Notes_1(i__),Nsys_1(i__),Ncomp_1(i__),Nparm(i__),Ncorr(i__),
+ seq(i__),comp_id(i__),HIP_1(i__),Hp(i__),e_Hp(i__),BT(i__),
+ e_BT(i__),VT(i__),e_VT(i__),RAdeg_1(i__),DEdeg_1(i__),
+ Plx_1(i__),pmRA_1(i__),pmDE_1(i__),e_RAdeg_1(i__),
+ e_DEdeg_1(i__),e_Plx_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),
+ ref(i__),theta_1(i__),rho_1(i__),d_theta(i__),d_rho(i__),
+ seq_ref(i__),flag(i__)
if(ar__1(63:68) .EQ. '') BT(i__) = rNULL__
if(ar__1(70:74) .EQ. '') e_BT(i__) = rNULL__
if(ar__1(76:81) .EQ. '') VT(i__) = rNULL__
if(ar__1(83:87) .EQ. '') e_VT(i__) = rNULL__
if(ar__1(178:184) .EQ. '') theta_1(i__) = rNULL__
if(ar__1(186:193) .EQ. '') rho_1(i__) = rNULL__
if(ar__1(195:202) .EQ. '') d_theta(i__) = rNULL__
if(ar__1(204:209) .EQ. '') d_rho(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ CCDM_1(i__),S(i__),Type(i__),Source_1(i__),Qual_1(i__),
+ Notes_1(i__),Nsys_1(i__),Ncomp_1(i__),Nparm(i__),Ncorr(i__),
+ seq(i__),comp_id(i__),HIP_1(i__),Hp(i__),e_Hp(i__),BT(i__),
+ e_BT(i__),VT(i__),e_VT(i__),RAdeg_1(i__),DEdeg_1(i__),
+ Plx_1(i__),pmRA_1(i__),pmDE_1(i__),e_RAdeg_1(i__),
+ e_DEdeg_1(i__),e_Plx_1(i__),e_pmRA_1(i__),e_pmDE_1(i__),
+ ref(i__),theta_1(i__),rho_1(i__),d_theta(i__),d_rho(i__),
+ seq_ref(i__),flag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'h_dm_cor.dat' ! Double and Multiples: Component solutions -CORR
C Format for file interpretation
3 format(
+ A10,1X,I1,1X,A1,1X,A1,1X,A1,1X,A1,1X,I1,1X,I2,1X,I2,1X,I2,6X,
+ I2,1X,66I3)
C Effective file loading
open(unit=1,status='old',file=
+'h_dm_cor.dat')
write(6,*) '....Loading file: h_dm_cor.dat'
do i__=1,12591
read(1,'(A238)')ar__2
read(ar__2,3)
+ CCDM_2(i__),S_1(i__),Type_1(i__),Source_2(i__),Qual_2(i__),
+ Notes_2(i__),Nsys_2(i__),Ncomp_2(i__),Nparm_1(i__),
+ Ncorr_1(i__),seq_1(i__),(corr(j__,i__),j__=1,66)
if(ar__2(41:238) .EQ. '') corr(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ CCDM_2(i__),S_1(i__),Type_1(i__),Source_2(i__),Qual_2(i__),
+ Notes_2(i__),Nsys_2(i__),Ncomp_2(i__),Nparm_1(i__),
+ Ncorr_1(i__),seq_1(i__),(corr(j__,i__),j__=1,66)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hip_dm_g.dat' ! Double and Multiples: Acceleration solutions
C Format for file interpretation
4 format(
+ I6,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,F7.2,1X,
+ F7.2,1X,F7.2,1X,F5.2,1X,A1,1X,I1,1X,36I3)
C Effective file loading
open(unit=1,status='old',file=
+'hip_dm_g.dat')
write(6,*) '....Loading file: hip_dm_g.dat'
do i__=1,2622
read(1,'(A195)')ar__3
read(ar__3,4)
+ HIP_2(i__),gRA(i__),gDE(i__),e_gRA(i__),e_gDE(i__),Fg(i__),
+ dgRA(i__),dgDE(i__),e_dgRA(i__),e_dgDE(i__),Fdg(i__),
+ Notes_3(i__),num(i__),(corr_1(j__,i__),j__=1,36)
if(ar__3(46:52) .EQ. '') dgRA(i__) = rNULL__
if(ar__3(54:60) .EQ. '') dgDE(i__) = rNULL__
if(ar__3(62:68) .EQ. '') e_dgRA(i__) = rNULL__
if(ar__3(70:76) .EQ. '') e_dgDE(i__) = rNULL__
if(ar__3(78:82) .EQ. '') Fdg(i__) = rNULL__
if(ar__3(88:195) .EQ. '') corr_1(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ HIP_2(i__),gRA(i__),gDE(i__),e_gRA(i__),e_gDE(i__),Fg(i__),
+ dgRA(i__),dgDE(i__),e_dgRA(i__),e_dgDE(i__),Fdg(i__),
+ Notes_3(i__),num(i__),(corr_1(j__,i__),j__=1,36)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hip_dm_o.dat' ! Double and Multiples: Orbital solutions
C Format for file interpretation
5 format(
+ I6,1X,F10.4,1X,F11.4,1X,F8.2,1X,F6.4,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F8.4,1X,F9.4,1X,F5.2,1X,F6.4,1X,F6.2,1X,F6.2,1X,F6.2,1X,I3,
+ 1X,A1,1X,A12,1X,66I3)
C Effective file loading
open(unit=1,status='old',file=
+'hip_dm_o.dat')
write(6,*) '....Loading file: hip_dm_o.dat'
do i__=1,235
read(1,'(A337)')ar__4
read(ar__4,5)
+ HIP_3(i__),P(i__),T(i__),a0(i__),ecc(i__),w(i__),i(i__),
+ Omega(i__),e_P(i__),e_T(i__),e_a0(i__),e_ecc(i__),e_w(i__),
+ e_i(i__),e_Omega(i__),dmRef(i__),Notes_4(i__),flag_1(i__),
+ (corr_2(j__,i__),j__=1,66)
if(ar__4(68:75) .EQ. '') e_P(i__) = rNULL__
if(ar__4(77:85) .EQ. '') e_T(i__) = rNULL__
if(ar__4(87:91) .EQ. '') e_a0(i__) = rNULL__
if(ar__4(93:98) .EQ. '') e_ecc(i__) = rNULL__
if(ar__4(100:105) .EQ. '') e_w(i__) = rNULL__
if(ar__4(107:112) .EQ. '') e_i(i__) = rNULL__
if(ar__4(114:119) .EQ. '') e_Omega(i__) = rNULL__
if(ar__4(121:123) .EQ. '') dmRef(i__) = iNULL__
c ..............Just test output...........
write(6,5)
+ HIP_3(i__),P(i__),T(i__),a0(i__),ecc(i__),w(i__),i(i__),
+ Omega(i__),e_P(i__),e_T(i__),e_a0(i__),e_ecc(i__),e_w(i__),
+ e_i(i__),e_Omega(i__),dmRef(i__),Notes_4(i__),flag_1(i__),
+ (corr_2(j__,i__),j__=1,66)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hip_dm_v.dat' ! Double and Multiples: VIM solutions
C Format for file interpretation
6 format(
+ I6,1X,F5.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,F5.2,1X,F6.2,1X,
+ F6.2,1X,F6.1,1X,F6.1,1X,A1,1X,21I3)
C Effective file loading
open(unit=1,status='old',file=
+'hip_dm_v.dat')
write(6,*) '....Loading file: hip_dm_v.dat'
do i__=1,288
read(1,'(A144)')ar__5
read(ar__5,6)
+ HIP_4(i__),Hp_ref(i__),DRA(i__),DDE(i__),e_DRA(i__),
+ e_DDE(i__),FD(i__),theta_C(i__),e_theta_C(i__),minSep(i__),
+ dvar(i__),Notes_5(i__),(corr_3(j__,i__),j__=1,21)
c ..............Just test output...........
write(6,6)
+ HIP_4(i__),Hp_ref(i__),DRA(i__),DDE(i__),e_DRA(i__),
+ e_DDE(i__),FD(i__),theta_C(i__),e_theta_C(i__),minSep(i__),
+ dvar(i__),Notes_5(i__),(corr_3(j__,i__),j__=1,21)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hip_dm_x.dat' ! Double and Multiples: Stochastic solutions
C Format for file interpretation
7 format(I6,1X,F6.2,1X,F6.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'hip_dm_x.dat')
write(6,*) '....Loading file: hip_dm_x.dat'
do i__=1,1561
read(1,'(A22)')ar__6
read(ar__6,7)
+ HIP_5(i__),epsilon(i__),e_epsilon(i__),Notes_6(i__)
c ..............Just test output...........
write(6,7)
+ HIP_5(i__),epsilon(i__),e_epsilon(i__),Notes_6(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hip_va_1.dat' ! Variability Annex: Periodic variables
C Format for file interpretation
8 format(
+ I6,1X,A1,1X,A12,1X,A1,1X,A6,1X,A1,1X,F6.3,1X,A1,1X,F6.3,1X,
+ F6.1,1X,F12.7,1X,F6.1,1X,F9.4,1X,I1,1X,A1,1X,A1,1X,A12,1X,
+ F10.5,1X,F10.2,1X,F5.2,1X,F5.2,1X,A1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'hip_va_1.dat')
write(6,*) '....Loading file: hip_va_1.dat'
do i__=1,2712
read(1,'(A142)')ar__7
read(ar__7,8)
+ HIP_6(i__),NewVar(i__),SpType_1(i__),HvarType_1(i__),
+ VarType(i__),n_VarType(i__),maxMag(i__),l_minMag(i__),
+ minMag(i__),log_sA_A(i__),Period_1(i__),log_sP(i__),
+ Ep_2440000(i__),q_Ep_2440000(i__),morePhoto_1(i__),
+ Notes_7(i__),VarName(i__),period_2(i__),epoch(i__),max(i__),
+ min(i__),Band(i__),refFlag(i__)
if(ar__7(50:55) .EQ. '') log_sA_A(i__) = rNULL__
if(ar__7(57:68) .EQ. '') Period_1(i__) = rNULL__
if(ar__7(70:75) .EQ. '') log_sP(i__) = rNULL__
if(ar__7(77:85) .EQ. '') Ep_2440000(i__) = rNULL__
if(ar__7(87:87) .EQ. '') q_Ep_2440000(i__) = iNULL__
if(ar__7(106:115) .EQ. '') period_2(i__) = rNULL__
if(ar__7(117:126) .EQ. '') epoch(i__) = rNULL__
if(ar__7(128:132) .EQ. '') max(i__) = rNULL__
if(ar__7(134:138) .EQ. '') min(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ HIP_6(i__),NewVar(i__),SpType_1(i__),HvarType_1(i__),
+ VarType(i__),n_VarType(i__),maxMag(i__),l_minMag(i__),
+ minMag(i__),log_sA_A(i__),Period_1(i__),log_sP(i__),
+ Ep_2440000(i__),q_Ep_2440000(i__),morePhoto_1(i__),
+ Notes_7(i__),VarName(i__),period_2(i__),epoch(i__),max(i__),
+ min(i__),Band(i__),refFlag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hip_va_2.dat' ! Variability Annex: Unsolved variables
C Format for file interpretation
9 format(
+ I6,1X,A1,1X,A12,1X,A1,1X,A6,1X,A1,1X,F6.3,1X,A1,1X,F6.3,1X,
+ F6.3,1X,F12.3,1X,F6.3,13X,A1,1X,A1,1X,A12,1X,F10.5,1X,F10.2,
+ 1X,F5.2,1X,F5.2,1X,A1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'hip_va_2.dat')
write(6,*) '....Loading file: hip_va_2.dat'
do i__=1,5542
read(1,'(A142)')ar__8
read(ar__8,9)
+ HIP_7(i__),NewVar_1(i__),SpType_2(i__),HvarType_2(i__),
+ VarType_1(i__),n_VarType_1(i__),maxMag_1(i__),l_minMag_1(i__),
+ minMag_1(i__),med_Hp(i__),A(i__),e_A(i__),morePhoto_2(i__),
+ Notes_8(i__),VarName_1(i__),period_3(i__),Ep_2440000_1(i__),
+ max_1(i__),min_1(i__),Band_1(i__),refFlag_1(i__)
if(ar__8(106:115) .EQ. '') period_3(i__) = rNULL__
if(ar__8(117:126) .EQ. '') Ep_2440000_1(i__) = rNULL__
if(ar__8(128:132) .EQ. '') max_1(i__) = rNULL__
if(ar__8(134:138) .EQ. '') min_1(i__) = rNULL__
c ..............Just test output...........
write(6,9)
+ HIP_7(i__),NewVar_1(i__),SpType_2(i__),HvarType_2(i__),
+ VarType_1(i__),n_VarType_1(i__),maxMag_1(i__),l_minMag_1(i__),
+ minMag_1(i__),med_Hp(i__),A(i__),e_A(i__),morePhoto_2(i__),
+ Notes_8(i__),VarName_1(i__),period_3(i__),Ep_2440000_1(i__),
+ max_1(i__),min_1(i__),Band_1(i__),refFlag_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'solar_ha.dat' ! Solar System Annex: Astrometric catalogue
C Format for file interpretation
10 format(
+ I3,1X,F11.7,1X,F11.7,1X,F13.7,1X,F5.2,1X,F7.3,1X,F6.2,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'solar_ha.dat')
write(6,*) '....Loading file: solar_ha.dat'
do i__=1,5609
read(1,'(A64)')ar__9
read(ar__9,10)
+ ID(i__),RAdeg_2(i__),DEdeg_2(i__),Ep_2440000_2(i__),
+ delay(i__),theta_2(i__),e_lambda(i__),flag_2(i__)
c ..............Just test output...........
write(6,10)
+ ID(i__),RAdeg_2(i__),DEdeg_2(i__),Ep_2440000_2(i__),
+ delay(i__),theta_2(i__),e_lambda(i__),flag_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'solar_hp.dat' ! Solar System Annex: Photometric catalogue
C Format for file interpretation
11 format(
+ I3,1X,F11.5,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F5.3,1X,F5.3,
+ 1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'solar_hp.dat')
write(6,*) '....Loading file: solar_hp.dat'
do i__=1,2639
read(1,'(A63)')ar__10
read(ar__10,11)
+ ID_1(i__),Ep_2440000_3(i__),Hp_dc(i__),e_Hp_dc(i__),
+ Hp_ac(i__),e_Hp_ac(i__),r(i__),Delta(i__),alpha(i__)
c ..............Just test output...........
write(6,11)
+ ID_1(i__),Ep_2440000_3(i__),Hp_dc(i__),e_Hp_dc(i__),
+ Hp_ac(i__),e_Hp_ac(i__),r(i__),Delta(i__),alpha(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'solar_t.dat' ! Solar System Annex: Tycho astrometry/photometry
C Format for file interpretation
12 format(
+ I3,1X,F13.7,1X,F11.7,1X,F11.7,1X,F5.2,1X,F5.2,1X,I1,1X,F5.1,
+ 1X,F5.1,1X,F5.2,1X,F6.2,1X,I2,1X,F5.1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'solar_t.dat')
write(6,*) '....Loading file: solar_t.dat'
do i__=1,291
read(1,'(A95)')ar__11
read(ar__11,12)
+ ID_2(i__),Ep_2440000_4(i__),RAdeg_3(i__),DEdeg_3(i__),
+ BTmag_1(i__),VTmag_1(i__),Flag_3(i__),e_RAdeg_2(i__),
+ e_DEdeg_2(i__),DE_RA_1(i__),theta_3(i__),sign_z(i__),
+ e_incl(i__),e_vert(i__)
if(ar__11(43:47) .EQ. '') BTmag_1(i__) = rNULL__
if(ar__11(49:53) .EQ. '') VTmag_1(i__) = rNULL__
c ..............Just test output...........
write(6,12)
+ ID_2(i__),Ep_2440000_4(i__),RAdeg_3(i__),DEdeg_3(i__),
+ BTmag_1(i__),VTmag_1(i__),Flag_3(i__),e_RAdeg_2(i__),
+ e_DEdeg_2(i__),DE_RA_1(i__),theta_3(i__),sign_z(i__),
+ e_incl(i__),e_vert(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hd_notes.doc' ! Hipparcos notes: Double and multiple systems
C Format for file interpretation
13 format(I6,1X,A1,1X,A1,1X,I2,1X,I2,1X,A80)
C Effective file loading
open(unit=1,status='old',file=
+'hd_notes.doc')
write(6,*) '....Loading file: hd_notes.doc'
do i__=1,2622
read(1,'(A97)')ar__12
read(ar__12,13)
+ HIP_8(i__),Note1(i__),Note2(i__),Ntot(i__),Nline(i__),
+ Text(i__)
c ..............Just test output...........
write(6,13)
+ HIP_8(i__),Note1(i__),Note2(i__),Ntot(i__),Nline(i__),
+ Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hg_notes.doc' ! Hipparcos notes: General notes
C Format for file interpretation
14 format(I6,1X,A1,1X,A1,1X,I2,1X,I2,1X,A80)
C Effective file loading
open(unit=1,status='old',file=
+'hg_notes.doc')
write(6,*) '....Loading file: hg_notes.doc'
do i__=1,3898
read(1,'(A97)')ar__13
read(ar__13,14)
+ HIP_9(i__),Note1_1(i__),Note2_1(i__),Ntot_1(i__),Nline_1(i__),
+ Text_1(i__)
c ..............Just test output...........
write(6,14)
+ HIP_9(i__),Note1_1(i__),Note2_1(i__),Ntot_1(i__),Nline_1(i__),
+ Text_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hp_notes.doc' ! Hipparcos notes: Photometric notes
C Format for file interpretation
15 format(I6,1X,A1,1X,A1,1X,I2,1X,I2,1X,A80)
C Effective file loading
open(unit=1,status='old',file=
+'hp_notes.doc')
write(6,*) '....Loading file: hp_notes.doc'
do i__=1,2444
read(1,'(A97)')ar__14
read(ar__14,15)
+ HIP_10(i__),Note1_2(i__),Note2_2(i__),Ntot_2(i__),
+ Nline_2(i__),Text_2(i__)
c ..............Just test output...........
write(6,15)
+ HIP_10(i__),Note1_2(i__),Note2_2(i__),Ntot_2(i__),
+ Nline_2(i__),Text_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hp_refs.doc' ! References Hipparcos stars
C Format for file interpretation
16 format(I6,1X,I2,1X,I2,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'hp_refs.doc')
write(6,*) '....Loading file: hp_refs.doc'
do i__=1,33769
read(1,'(A19)')ar__15
read(ar__15,16)HIP_11(i__),Ntot_3(i__),Nline_3(i__),nRef(i__)
c ..............Just test output...........
write(6,16)HIP_11(i__),Ntot_3(i__),Nline_3(i__),nRef(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hp_auth.doc' ! References of hp_notes.doc
C Format for file interpretation
17 format(F6.3,1X,A70)
C Effective file loading
open(unit=1,status='old',file=
+'hp_auth.doc')
write(6,*) '....Loading file: hp_auth.doc'
do i__=1,4335
read(1,'(A80)')ar__16
read(ar__16,17)nRef_1(i__),Text_3(i__)
c ..............Just test output...........
write(6,17)nRef_1(i__),Text_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'dmsa_o.doc' ! References of hip_dm_o.dat
C Format for file interpretation
18 format(I3,1X,I1,1X,I1,1X,A72)
C Effective file loading
open(unit=1,status='old',file=
+'dmsa_o.doc')
write(6,*) '....Loading file: dmsa_o.doc'
do i__=1,118
read(1,'(A80)')ar__17
read(ar__17,18)dmRef_1(i__),Ntot_4(i__),Nline_4(i__),Text_4(i__)
c ..............Just test output...........
write(6,18)dmRef_1(i__),Ntot_4(i__),Nline_4(i__),Text_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tyc_main.dat' ! The main part of Tycho Catalogue
C Format for file interpretation
19 format(
+ A1,1X,A12,1X,A1,1X,A11,1X,A11,1X,F5.2,3X,A1,1X,F12.8,1X,F12.8,
+ 1X,A1,1X,F7.2,1X,F8.2,1X,F8.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,I3,1X,F5.2,1X,I6,1X,F6.3,1X,
+ F5.3,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F5.3,3X,I1,1X,F4.1,1X,
+ A1,1X,I3,1X,F5.3,1X,F5.2,1X,F5.2,1X,A1,1X,A1,1X,A1,1X,A1,1X,
+ A2,1X,I6,1X,I6,1X,A10,1X,A10,1X,A10,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tyc_main.dat')
write(6,*) '....Loading file: tyc_main.dat'
do i__=1,1058332
read(1,'(A350)')ar__18
read(ar__18,19)
+ Catalog_1(i__),TYC(i__),Proxy_1(i__),RAhms_1(i__),
+ DEdms_1(i__),Vmag_1(i__),r_Vmag_1(i__),RAdeg_4(i__),
+ DEdeg_4(i__),AstroRef_1(i__),Plx_2(i__),pmRA_2(i__),
+ pmDE_2(i__),e_RAdeg_3(i__),e_DEdeg_3(i__),e_Plx_2(i__),
+ e_pmRA_2(i__),e_pmDE_2(i__),DE_RA_2(i__),Plx_RA_1(i__),
+ Plx_DE_1(i__),pmRA_RA_1(i__),pmRA_DE_1(i__),pmRA_Plx_1(i__),
+ pmDE_RA_1(i__),pmDE_DE_1(i__),pmDE_Plx_1(i__),
+ pmDE_pmRA_1(i__),Nastro(i__),F2_1(i__),HIP_12(i__),
+ BTmag_2(i__),e_BTmag_1(i__),VTmag_2(i__),e_VTmag_1(i__),
+ r_BTmag(i__),B_V_1(i__),e_B_V_1(i__),Q(i__),Fs(i__),
+ Source_3(i__),Nphoto(i__),VTscat(i__),VTmax(i__),VTmin(i__),
+ Var(i__),VarFlag_1(i__),MultFlag_1(i__),morePhoto_3(i__),
+ m_HIP_1(i__),PPM(i__),HD_1(i__),BD_1(i__),CoD_1(i__),
+ CPD_1(i__),Remark(i__)
if(ar__18(42:46) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__18(80:86) .EQ. '') Plx_2(i__) = rNULL__
if(ar__18(88:95) .EQ. '') pmRA_2(i__) = rNULL__
if(ar__18(97:104) .EQ. '') pmDE_2(i__) = rNULL__
if(ar__18(106:111) .EQ. '') e_RAdeg_3(i__) = rNULL__
if(ar__18(113:118) .EQ. '') e_DEdeg_3(i__) = rNULL__
if(ar__18(120:125) .EQ. '') e_Plx_2(i__) = rNULL__
if(ar__18(127:132) .EQ. '') e_pmRA_2(i__) = rNULL__
if(ar__18(134:139) .EQ. '') e_pmDE_2(i__) = rNULL__
if(ar__18(141:145) .EQ. '') DE_RA_2(i__) = rNULL__
if(ar__18(147:151) .EQ. '') Plx_RA_1(i__) = rNULL__
if(ar__18(153:157) .EQ. '') Plx_DE_1(i__) = rNULL__
if(ar__18(159:163) .EQ. '') pmRA_RA_1(i__) = rNULL__
if(ar__18(165:169) .EQ. '') pmRA_DE_1(i__) = rNULL__
if(ar__18(171:175) .EQ. '') pmRA_Plx_1(i__) = rNULL__
if(ar__18(177:181) .EQ. '') pmDE_RA_1(i__) = rNULL__
if(ar__18(183:187) .EQ. '') pmDE_DE_1(i__) = rNULL__
if(ar__18(189:193) .EQ. '') pmDE_Plx_1(i__) = rNULL__
if(ar__18(195:199) .EQ. '') pmDE_pmRA_1(i__) = rNULL__
if(ar__18(201:203) .EQ. '') Nastro(i__) = iNULL__
if(ar__18(205:209) .EQ. '') F2_1(i__) = rNULL__
if(ar__18(211:216) .EQ. '') HIP_12(i__) = iNULL__
if(ar__18(218:223) .EQ. '') BTmag_2(i__) = rNULL__
if(ar__18(225:229) .EQ. '') e_BTmag_1(i__) = rNULL__
if(ar__18(231:236) .EQ. '') VTmag_2(i__) = rNULL__
if(ar__18(238:242) .EQ. '') e_VTmag_1(i__) = rNULL__
if(ar__18(246:251) .EQ. '') B_V_1(i__) = rNULL__
if(ar__18(253:257) .EQ. '') e_B_V_1(i__) = rNULL__
if(ar__18(261:261) .EQ. '') Q(i__) = iNULL__
if(ar__18(263:266) .EQ. '') Fs(i__) = rNULL__
if(ar__18(270:272) .EQ. '') Nphoto(i__) = iNULL__
if(ar__18(274:278) .EQ. '') VTscat(i__) = rNULL__
if(ar__18(280:284) .EQ. '') VTmax(i__) = rNULL__
if(ar__18(286:290) .EQ. '') VTmin(i__) = rNULL__
if(ar__18(303:308) .EQ. '') PPM(i__) = iNULL__
if(ar__18(310:315) .EQ. '') HD_1(i__) = iNULL__
c ..............Just test output...........
write(6,19)
+ Catalog_1(i__),TYC(i__),Proxy_1(i__),RAhms_1(i__),
+ DEdms_1(i__),Vmag_1(i__),r_Vmag_1(i__),RAdeg_4(i__),
+ DEdeg_4(i__),AstroRef_1(i__),Plx_2(i__),pmRA_2(i__),
+ pmDE_2(i__),e_RAdeg_3(i__),e_DEdeg_3(i__),e_Plx_2(i__),
+ e_pmRA_2(i__),e_pmDE_2(i__),DE_RA_2(i__),Plx_RA_1(i__),
+ Plx_DE_1(i__),pmRA_RA_1(i__),pmRA_DE_1(i__),pmRA_Plx_1(i__),
+ pmDE_RA_1(i__),pmDE_DE_1(i__),pmDE_Plx_1(i__),
+ pmDE_pmRA_1(i__),Nastro(i__),F2_1(i__),HIP_12(i__),
+ BTmag_2(i__),e_BTmag_1(i__),VTmag_2(i__),e_VTmag_1(i__),
+ r_BTmag(i__),B_V_1(i__),e_B_V_1(i__),Q(i__),Fs(i__),
+ Source_3(i__),Nphoto(i__),VTscat(i__),VTmax(i__),VTmin(i__),
+ Var(i__),VarFlag_1(i__),MultFlag_1(i__),morePhoto_3(i__),
+ m_HIP_1(i__),PPM(i__),HD_1(i__),BD_1(i__),CoD_1(i__),
+ CPD_1(i__),Remark(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end