Conversion of standardized ReadMe file for
file /./ftp/cats/V/137D 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 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. V/137D Extended Hipparcos Compilation (XHIP) (Anderson+, 2012)
*================================================================================
*XHIP: An Extended Hipparcos Compilation
* Anderson E., Francis C.
* <Astron. Letters 38 (2012)>
* =2012AstL...38..331A (http://arxiv.org/abs/1108.4971)
* =2012PAZh...38..374A
* =2012yCat.5137....0A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'main.dat' ! Astrometry, spectrography, space motions, exoplanets
integer*4 nr__
parameter (nr__=117955) ! Number of records
character*417 ar__ ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=J1991.25)
integer*4 HIP ! Hipparcos identifier
character*6 Comp ! Component(s) (1)
character*41 Classes ! SIMBAD classifications (comma separated) (2)
integer*4 Gr ! ? Cluster or Association membership(s) (3)
real*8 RAdeg ! (deg) Right ascension (ICRS, Epoch=J1991.25)
real*8 DEdeg ! (deg) Declination (ICRS, Epoch=J1991.25)
real*4 Plx ! (mas) Trigonometric parallax
real*8 pmRA ! (mas/yr) Proper motion in RA*cos(DEdeg)
real*8 pmDE ! (mas/yr) Proper motion in Declination
real*4 e_RAdeg ! (mas) Standard error on RA*cos(DEdeg)
real*4 e_DEdeg ! (mas) Standard error on DE
real*4 e_Plx ! (mas) Standard error on Plx
real*4 e_pmRA ! (mas/yr) Standard error on pmRA
real*4 e_pmDE ! (mas/yr) Standard error on pmDE
integer*4 r_HIP ! [1/2] Reference for Hipparcos astrometry (4)
integer*4 r_pm ! [1/3] Reference for proper motion (5)
real*8 GLon ! (deg) Galactic longitude (6)
real*8 GLat ! (deg) Galactic latitude (6)
real*8 Dist ! (pc) ? Heliocentric distance (7)
real*4 e_Dist ! (%) ? Distance error expressed as percentage (8)
real*8 pmGLon ! (mas/yr) Proper motion in GLon*cos(GLat)
real*8 pmGLat ! (mas/yr) Proper motion in GLat
real*4 X ! (pc) ? Heliocentric distance towards Gal. center
real*8 Y ! (pc) ? Heliocentric distance towards Gal. rotation
real*4 Z ! (pc) ? Heliocentric distance towards N. Gal. Pole
integer*4 RGal ! (pc) ? Galactocentric distance (9)
real*4 vT ! (km/s) ? Transverse velocity
character*26 SpType ! Spectral type (MK, HD, or other)
integer*4 Tc ! ]0/140[? Temperature class codified (10)
integer*4 Lc ! [1/6]? Luminosity class codified (11)
real*8 RV ! (km/s) ? Radial velocity
real*4 e_RV ! (km/s) ? Standard error on RV (12)
character*1 q_RV ! [ABCD] Quality flag on RV (13)
real*4 v_Fe_H_ ! ([Sun]) ? Iron abundance
real*4 e__Fe_H_ ! ([Sun]) ? Standard error on [Fe/H] (14)
character*1 q__Fe_H_ ! [ABC] Quality flag on [Fe/H] (15)
real*4 age ! (Gyr) ? Age, in billions of years
real*4 b_age ! (Gyr) ? Lower confidence limit on age
real*4 B_age_1 ! (Gyr) ? Upper confidence limit on age
real*4 U ! (km/s) ? Heliocentric velocity towards Gal. center
real*4 V ! (km/s) ? Heliocentric velocity towards Gal. rotation
real*4 W ! (km/s) ? Heliocentric velocity towards N. Gal. Pole
real*4 UVW ! (km/s) ? Total heliocentric velocity
real*8 Dmin ! (pc) ? Minimum distance at solar encounter (16)
integer*4 Tmin ! (kyr) ? Timing of Dmin, in thousands of years (17)
real*4 e ! ? Total orbital eccentricity (18)
real*4 phi ! (deg) [-180/180]? Pericenter position angle (18)
integer*4 a ! (pc) ? Semi-major axis of orbit (18)
integer*4 b ! (pc) ? Semi-minor axis of orbit (18)
integer*4 c ! (pc) ? Focus-to-center distance of orbit (18)
integer*4 L ! (deg) ? Semilatus rectum of orbit (18)
integer*4 Rmin ! (pc) ? Orbital radius at pericenter (18)
integer*4 Rmax ! (pc) ? Orbital radius at apocenter (18)
integer*4 Npl ! ? Number of exoplanets (known in April 2012)
character*5 Mpl ! Exoplanet discovery method(s) (19)
*Note (1): from CCDM (Cat. I/274) or WDS (Cat. B/wds)
*Note (2): object types, see
* http://simbad.u-strasbg.fr/simbad/sim-display?data=otypes
*Note (3): Group IDs are in biblio.dat (GrpName)
* see also: groups.dat and Cat. B/Ocl
*Note (4): flag indicates astrometry source:
* 1 = HIP (Cat. I/239) for 3,504 problematic cases in HIP2
* 2 = HIP2 (Cat. I/311) for all remaining cases
*Note (5): flag indicates the following:
* 1 = From HIP2 (Cat. I/311)
* 2 = From Tycho-2 (Cat. I/259)
* 3 = HIP2 & Tycho-2 combined for improved error bounds
*Note (6): transformation from ICRS defined by:
* NGP RAdeg = 192.85948,
* NGP DEdeg = +27.12825,
* GLon of ascending node @ celestial equator = 32.93192 deg
*Note (7): null where parallax errors <20%
* Dist = 1000/Plx * (1+1.2*(e_Plx/Plx)^2^)
*Note (8): null where Dist is fitted to cluster
*Note (9): adopts Sgr A* (X,Y,Z) as (7400,-7.2,-6.0) pc
*Note (10): Temperature classes are:
* O=10 B=20 A=30 F=40 G=50 K=60 M=70
* L=80 T=90 S=100 C=110 R=120 N=130
* Subclasses 0-9 are additionally summed.
*Note (11): I=1 II=2 III=3 IV=4 V=5 VI=6
*Note (12): 999 = error not available though RV is available
*Note (13): quality flag means:
* A = Most likely to be within stated error bounds
* B = May have small systematic errors
* C = May have larger systematic errors
* D = No error available or other serious problems
*Note (14): 9.99 = error not available though [Fe/H] is available
*Note (15): quality flag means:
* A = Calibrated
* B = Calibrated, but no error bounds assigned
* C = Uncalibrated and no error bounds assigned
*Note (16): based on linear approximation to space motion
*Note (17): Past times are negative; future times positive
*Note (18): LSR (U,V,W) = (14.1,14.6,6.9) km/s are
* adopted. Note (9) above with regard to Rgal also applies.
*Note (19): method is codified:
* RA = "radial velocity or astrometric methods"
* I = "imaging"
* T = "timing"
* X = "transit"
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'photo.dat' ! Photometry
integer*4 nr__1
parameter (nr__1=117955) ! Number of records
character*235 ar__1 ! Full-size record
integer*4 HIP_1 ! (pc) Hipparcos identifier
real*8 Hpmag ! (mag) Median magnitude in Hipparcos system
real*4 e_Hpmag ! (mag) Standard error on Hpmag
character*1 m_Hpmag ! [A-Z*-] Reference flag for Hpmag (20)
real*4 Hpmax ! (mag) ? Hpmag at maximum (5th percentile)
real*4 Hpmin ! (mag) ? Hpmag at minimum (95th percentile)
real*4 Per ! (d) ? Variability period
character*1 Hvar ! (mag) [CDMPRU] Variability type (21)
real*4 Umag ! (mag) ? Apparent magnitude in Johnson U
real*4 Bmag ! (mag) ? Apparent magnitude in Johnson B
real*4 Vmag ! (mag) Apparent magnitude in Johnson V
real*4 Rmag ! (mag) ? Apparent magnitude in R
real*4 Imag ! (mag) ? Apparent magnitude in I
real*4 Jmag ! (mag) ? J selected default magnitude (22)
real*4 Hmag ! (mag) ? H selected default magnitude (22)
real*4 Kmag ! (mag) ? K selected default magnitude (22)
real*4 e_Jmag ! (mag) ? J total magnitude uncertainty (22)
real*4 e_Hmag ! (mag) ? H total magnitude uncertainty (22)
real*4 e_Kmag ! (mag) ? K total magnitude uncertainty (22)
character*17 v2MASS ! 2MASS source designation (22)
character*3 q2M ! JHK photometric quality flag (22)
real*4 B_V ! (mag) ? Johnson B-V color index
real*4 V_I ! (mag) ? Color index in Cousins' system
real*4 e_B_V ! (mag) ? Standard error on B-V
real*4 e_V_I ! (mag) ? Standard error on V-I
real*4 HpMag_1 ! (mag) ? Absolute Magnitude of Hpmag (23)
real*4 UMag_1 ! (mag) ? Absolute Magnitude in Johnson U (23)
real*4 BMag_1 ! (mag) ? Absolute Magnitude in Johnson B (23)
real*4 VMag_1 ! (mag) ? Absolute Magnitude in Johnson V (23)
real*4 RMag_1 ! (mag) ? Absolute Magnitude in R (23)
real*4 IMag_1 ! (mag) ? Absolute Magnitude in I (23)
real*4 JMag_1 ! (mag) ? Absolute Magnitude J (23)
real*4 HMag_1 ! (mag) ? Absolute Magnitude H (23)
real*4 KMag_1 ! (mag) ? Absolute Magnitude K (23)
real*8 Lum ! (Lsun) ? Stellar luminosity (23)
real*4 magmin ! (mag) ? Apparent magnitude V at Tmin (main.dat)
*Note (20): 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 (21): Hipparcos-defined type of variability
* blank = 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 (22): See Cat. II/246
*Note (23): does not factor absorption
* by the interstellar medium; may also contain integrated components.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'biblio.dat' ! References
integer*4 nr__2
parameter (nr__2=117955) ! Number of records
character*376 ar__2 ! Full-size record
integer*4 HIP_2 ! Hipparcos identifier
integer*4 HD ! ? Henry Draper catalog identifier
character*3 Cst ! Constellation membership
integer*4 Atlas ! Millennium Star Atlas page number
character*10 Coords ! RA, DE in compact format; n = + ; s = -
character*48 Name ! Star name(s)
character*106 GrpName ! Cluster or Association name(s)
character*16 CompID ! Reference for Comp (main.dat)
character*25 rSpType ! Reference for SpType (main.dat) (24)
character*15 r_RV ! Reference for RV (main.dat)
character*127 r__Fe_H_ ! Reference(s) for [Fe/H] (main.dat)
*Note (24): BibCodes "1993BICDS..43....5T" & "2003AJ....125..359W" are
* augmented (after a colon :) with the content of the r_Sp or r_SpType
* columns from those sources respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! Reference key for r_RV & r_[Fe/H]
integer*4 nr__3
parameter (nr__3=674) ! Number of records
character*23 ar__3 ! Full-size record
integer*4 KeyCode ! [1,703] 3-digit shorthand code (25)
character*19 BibCode ! Corresponding 19-digit bibliographic code (26)
*Note (25): Range of key code:
* * [001-176]: Calibrated [Fe/H] sources
* * [201-495]: Uncalibrated [Fe/H] sources
* * [501-703]: Radial velocity sources
*Note (26): Radial velocities obtained through SIMBAD having no source
* reference are indicated with the pseudo-code "2010.........SIMBAD"
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'groups.dat' ! Statistics on open clusters and stellar associations
integer*4 nr__4
parameter (nr__4=87) ! Number of records
character*318 ar__4 ! Full-size record
C Position composed of: RAdeg DEdeg
character*21 GrpName_1 ! Group name
character*21 AltName ! Alternative group name
integer*4 N ! Number of Hiipparcos candidate stars
real*4 n_1 ! [3.5/16] Separation parameter (27)
character*13 Stream ! Kinematic group (28)
real*4 RAdeg_1 ! (deg) Group mean right ascension (ICRS)
real*4 DEdeg_1 ! (deg) Group mean declination (ICRS)
real*4 Dist_1 ! (pc) Group mean parallax distance
real*4 pmRA_1 ! (mas/yr) Group mean proper motion in RA*cos(DEdeg)
real*4 pmDE_1 ! (mas/yr) Mean proper motion in declination
real*4 RV_1 ! (km/s) Group mean radial velocity
real*4 s_RAdeg ! (deg) Group width on RA*cos(DEdeg) (1-{sigma})
real*4 s_DEdeg ! (deg) Group width on DE (1-{sigma})
real*4 s_Dist ! (pc) Group width on Dist (1-{sigma})
real*4 s_pmRA ! (mas/yr) Group width on pmRA (1-{sigma})
real*4 s_pmDE ! (mas/yr) Group width on pmDE (1-{sigma})
real*4 s_RV ! (km/s) Group width on radial velocity
real*4 e_RAdeg_1 ! (deg) Standard error on RA*cos(DEdeg)
real*4 e_DEdeg_1 ! (deg) Standard error on DE
real*4 e_Dist_1 ! (pc) Standard error on Dist
real*4 e_pmRA_1 ! (mas/yr) Standard error on pmRA
real*4 e_pmDE_1 ! (mas/yr) Standard error on pmDE
real*4 e_RV_1 ! (km/s) Standard error on radial velocity
real*4 rad ! (pc) Physical radius (1-{sigma})
real*4 arad ! (deg) Angular radius (1-{sigma})
real*4 X_1 ! (pc) Heliocentric distance towards Gal. center
real*4 Y_1 ! (pc) Heliocentric distance towards Gal. rotation
real*4 Z_1 ! (pc) Heliocentric distance towards N. Gal. Pole
real*4 U_1 ! (km/s) Heliocentric velocity towards Gal. center
real*4 V_1 ! (km/s) Heliocentric velocity towards Gal. rotation
real*4 W_1 ! (km/s) Heliocentric velocity towards N. Gal. Pole
real*4 s_X ! (pc) Group width on X (1-{sigma})
real*4 s_Y ! (pc) Group width on Y (1-{sigma})
real*4 s_Z ! (pc) Group width on Z (1-{sigma})
real*4 s_U ! (km/s) Group width on U (1-{sigma})
real*4 s_V ! (km/s) Group width on V (1-{sigma})
real*4 s_W ! (km/s) Group width on W (1-{sigma})
real*4 e_X ! (pc) Standard error on X
real*4 e_Y ! (pc) Standard error on Y
real*4 e_Z ! (pc) Standard error on Z
real*4 e_U ! (km/s) Standard error on U
real*4 e_V ! (km/s) Standard error on V
real*4 e_W ! (km/s) Standard error on W
real*4 v_Fe_H__1 ! ([Sun]) ? Iron abundance
real*4 e__Fe_H__1 ! ([Sun]) ? Standard error on [Fe/H]
integer*4 n__Fe_H_ ! ? Number of Hipparcos member stars with [Fe/H]
real*4 E_B_V_1 ! (mag) ? Color excess in B-V, CDS Cat. B/ocl
real*4 E_B_V_2 ! (mag) ? Color excess in B-V, CDS Cat. J/A+A/477/165
*Note (27): the separation parameter n is a measure of the concentration of
* a cluster, and of its separation from the surrounding star field (low
* values show poor dynamical separation from surrounding stars). A value
* of 5 is usual for a cluster, while 3.5 is typical for an
* association.
*Note (28): see 2012MNRAS.422.1283F (arXiv:1202.1375)
C=============================================================================
C Loading file 'main.dat' ! Astrometry, spectrography, space motions, exoplanets
C Format for file interpretation
1 format(
+ I6,1X,A6,1X,A41,1X,I1,1X,F12.8,1X,F12.8,1X,F6.2,1X,F8.2,1X,
+ F8.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F6.2,1X,F5.2,1X,I1,1X,I1,1X,
+ F12.8,1X,F12.8,1X,F7.2,1X,F4.1,1X,F8.2,1X,F8.2,1X,F6.1,1X,
+ F7.1,1X,F6.1,1X,I4,1X,F5.1,1X,A26,1X,I3,1X,I1,1X,F7.2,1X,F6.2,
+ 1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F6.1,
+ 1X,F6.1,1X,F6.1,1X,F5.1,1X,F7.2,1X,I5,1X,F6.4,1X,F6.1,1X,I7,
+ 1X,I5,1X,I7,1X,I5,1X,I4,1X,I7,1X,I1,1X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'main.dat')
write(6,*) '....Loading file: main.dat'
do i__=1,117955
read(1,'(A417)')ar__
read(ar__,1)
+ HIP,Comp,Classes,Gr,RAdeg,DEdeg,Plx,pmRA,pmDE,e_RAdeg,e_DEdeg,
+ e_Plx,e_pmRA,e_pmDE,r_HIP,r_pm,GLon,GLat,Dist,e_Dist,pmGLon,
+ pmGLat,X,Y,Z,RGal,vT,SpType,Tc,Lc,RV,e_RV,q_RV,v_Fe_H_,
+ e__Fe_H_,q__Fe_H_,age,b_age,B_age_1,U,V,W,UVW,Dmin,Tmin,e,phi,
+ a,b,c,L,Rmin,Rmax,Npl,Mpl
if(ar__(57:57) .EQ. '') Gr = iNULL__
if(ar__(173:179) .EQ. '') Dist = rNULL__
if(ar__(181:184) .EQ. '') e_Dist = rNULL__
if(ar__(204:209) .EQ. '') X = rNULL__
if(ar__(211:217) .EQ. '') Y = rNULL__
if(ar__(219:224) .EQ. '') Z = rNULL__
if(ar__(226:229) .EQ. '') RGal = iNULL__
if(ar__(231:235) .EQ. '') vT = rNULL__
if(ar__(264:266) .EQ. '') Tc = iNULL__
if(ar__(268:268) .EQ. '') Lc = iNULL__
if(ar__(270:276) .EQ. '') RV = rNULL__
if(ar__(278:283) .EQ. '') e_RV = rNULL__
if(ar__(287:291) .EQ. '') v_Fe_H_ = rNULL__
if(ar__(293:296) .EQ. '') e__Fe_H_ = rNULL__
if(ar__(300:303) .EQ. '') age = rNULL__
if(ar__(305:308) .EQ. '') b_age = rNULL__
if(ar__(310:313) .EQ. '') B_age_1 = rNULL__
if(ar__(315:320) .EQ. '') U = rNULL__
if(ar__(322:327) .EQ. '') V = rNULL__
if(ar__(329:334) .EQ. '') W = rNULL__
if(ar__(336:340) .EQ. '') UVW = rNULL__
if(ar__(342:348) .EQ. '') Dmin = rNULL__
if(ar__(350:354) .EQ. '') Tmin = iNULL__
if(ar__(356:361) .EQ. '') e = rNULL__
if(ar__(363:368) .EQ. '') phi = rNULL__
if(ar__(370:376) .EQ. '') a = iNULL__
if(ar__(378:382) .EQ. '') b = iNULL__
if(ar__(384:390) .EQ. '') c = iNULL__
if(ar__(392:396) .EQ. '') L = iNULL__
if(ar__(398:401) .EQ. '') Rmin = iNULL__
if(ar__(403:409) .EQ. '') Rmax = iNULL__
if(ar__(411:411) .EQ. '') Npl = iNULL__
c ..............Just test output...........
write(6,1)
+ HIP,Comp,Classes,Gr,RAdeg,DEdeg,Plx,pmRA,pmDE,e_RAdeg,e_DEdeg,
+ e_Plx,e_pmRA,e_pmDE,r_HIP,r_pm,GLon,GLat,Dist,e_Dist,pmGLon,
+ pmGLat,X,Y,Z,RGal,vT,SpType,Tc,Lc,RV,e_RV,q_RV,v_Fe_H_,
+ e__Fe_H_,q__Fe_H_,age,b_age,B_age_1,U,V,W,UVW,Dmin,Tmin,e,phi,
+ a,b,c,L,Rmin,Rmax,Npl,Mpl
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'photo.dat' ! Photometry
C Format for file interpretation
2 format(
+ I6,1X,F7.4,1X,F6.4,1X,A1,1X,F5.2,1X,F5.2,1X,F6.2,1X,A1,1X,
+ F6.3,1X,F6.3,1X,F5.2,1X,F5.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,A17,1X,A3,1X,F6.3,1X,F5.2,1X,F5.3,
+ 1X,F4.2,1X,F6.3,1X,F6.3,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F6.2,1X,F9.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'photo.dat')
write(6,*) '....Loading file: photo.dat'
do i__=1,117955
read(1,'(A235)')ar__1
read(ar__1,2)
+ HIP_1,Hpmag,e_Hpmag,m_Hpmag,Hpmax,Hpmin,Per,Hvar,Umag,Bmag,
+ Vmag,Rmag,Imag,Jmag,Hmag,Kmag,e_Jmag,e_Hmag,e_Kmag,v2MASS,q2M,
+ B_V,V_I,e_B_V,e_V_I,HpMag_1,UMag_1,BMag_1,VMag_1,RMag_1,
+ IMag_1,JMag_1,HMag_1,KMag_1,Lum,magmin
if(ar__1(25:29) .EQ. '') Hpmax = rNULL__
if(ar__1(31:35) .EQ. '') Hpmin = rNULL__
if(ar__1(37:42) .EQ. '') Per = rNULL__
if(ar__1(46:51) .EQ. '') Umag = rNULL__
if(ar__1(53:58) .EQ. '') Bmag = rNULL__
if(ar__1(66:70) .EQ. '') Rmag = rNULL__
if(ar__1(72:77) .EQ. '') Imag = rNULL__
if(ar__1(79:84) .EQ. '') Jmag = rNULL__
if(ar__1(86:91) .EQ. '') Hmag = rNULL__
if(ar__1(93:98) .EQ. '') Kmag = rNULL__
if(ar__1(100:104) .EQ. '') e_Jmag = rNULL__
if(ar__1(106:110) .EQ. '') e_Hmag = rNULL__
if(ar__1(112:116) .EQ. '') e_Kmag = rNULL__
if(ar__1(140:145) .EQ. '') B_V = rNULL__
if(ar__1(147:151) .EQ. '') V_I = rNULL__
if(ar__1(153:157) .EQ. '') e_B_V = rNULL__
if(ar__1(159:162) .EQ. '') e_V_I = rNULL__
if(ar__1(164:169) .EQ. '') HpMag_1 = rNULL__
if(ar__1(171:176) .EQ. '') UMag_1 = rNULL__
if(ar__1(178:182) .EQ. '') BMag_1 = rNULL__
if(ar__1(184:188) .EQ. '') VMag_1 = rNULL__
if(ar__1(190:194) .EQ. '') RMag_1 = rNULL__
if(ar__1(196:200) .EQ. '') IMag_1 = rNULL__
if(ar__1(202:206) .EQ. '') JMag_1 = rNULL__
if(ar__1(208:212) .EQ. '') HMag_1 = rNULL__
if(ar__1(214:219) .EQ. '') KMag_1 = rNULL__
if(ar__1(221:229) .EQ. '') Lum = rNULL__
if(ar__1(231:235) .EQ. '') magmin = rNULL__
c ..............Just test output...........
write(6,2)
+ HIP_1,Hpmag,e_Hpmag,m_Hpmag,Hpmax,Hpmin,Per,Hvar,Umag,Bmag,
+ Vmag,Rmag,Imag,Jmag,Hmag,Kmag,e_Jmag,e_Hmag,e_Kmag,v2MASS,q2M,
+ B_V,V_I,e_B_V,e_V_I,HpMag_1,UMag_1,BMag_1,VMag_1,RMag_1,
+ IMag_1,JMag_1,HMag_1,KMag_1,Lum,magmin
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'biblio.dat' ! References
C Format for file interpretation
3 format(
+ I6,1X,I6,1X,A3,1X,I4,1X,A10,1X,A48,1X,A106,1X,A16,1X,A25,1X,
+ A15,1X,A127)
C Effective file loading
open(unit=1,status='old',file=
+'biblio.dat')
write(6,*) '....Loading file: biblio.dat'
do i__=1,117955
read(1,'(A376)')ar__2
read(ar__2,3)
+ HIP_2,HD,Cst,Atlas,Coords,Name,GrpName,CompID,rSpType,r_RV,
+ r__Fe_H_
if(ar__2(8:13) .EQ. '') HD = iNULL__
c ..............Just test output...........
write(6,3)
+ HIP_2,HD,Cst,Atlas,Coords,Name,GrpName,CompID,rSpType,r_RV,
+ r__Fe_H_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! Reference key for r_RV & r_[Fe/H]
C Format for file interpretation
4 format(I3,1X,A19)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,674
read(1,'(A23)')ar__3
read(ar__3,4)KeyCode,BibCode
c ..............Just test output...........
write(6,4)KeyCode,BibCode
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'groups.dat' ! Statistics on open clusters and stellar associations
C Format for file interpretation
5 format(
+ A21,1X,A21,1X,I3,1X,F4.1,1X,A13,1X,F6.2,1X,F6.2,1X,F6.1,1X,
+ F6.1,1X,F6.1,1X,F5.1,1X,F6.2,1X,F6.2,1X,F6.1,1X,F5.1,1X,F5.1,
+ 1X,F4.1,1X,F5.2,1X,F5.2,1X,F6.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,
+ F5.1,1X,F5.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,F5.1,1X,F5.1,1X,F5.1,
+ 1X,F5.1,1X,F5.1,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F5.1,1X,
+ F5.1,1X,F5.1,1X,F3.1,1X,F3.1,1X,F3.1,1X,F5.2,1X,F4.2,1X,I3,1X,
+ F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'groups.dat')
write(6,*) '....Loading file: groups.dat'
do i__=1,87
read(1,'(A318)')ar__4
read(ar__4,5)
+ GrpName_1,AltName,N,n_1,Stream,RAdeg_1,DEdeg_1,Dist_1,pmRA_1,
+ pmDE_1,RV_1,s_RAdeg,s_DEdeg,s_Dist,s_pmRA,s_pmDE,s_RV,
+ e_RAdeg_1,e_DEdeg_1,e_Dist_1,e_pmRA_1,e_pmDE_1,e_RV_1,rad,
+ arad,X_1,Y_1,Z_1,U_1,V_1,W_1,s_X,s_Y,s_Z,s_U,s_V,s_W,e_X,e_Y,
+ e_Z,e_U,e_V,e_W,v_Fe_H__1,e__Fe_H__1,n__Fe_H_,E_B_V_1,E_B_V_2
if(ar__4(295:299) .EQ. '') v_Fe_H__1 = rNULL__
if(ar__4(301:304) .EQ. '') e__Fe_H__1 = rNULL__
if(ar__4(306:308) .EQ. '') n__Fe_H_ = iNULL__
if(ar__4(310:313) .EQ. '') E_B_V_1 = rNULL__
if(ar__4(315:318) .EQ. '') E_B_V_2 = rNULL__
c ..............Just test output...........
write(6,5)
+ GrpName_1,AltName,N,n_1,Stream,RAdeg_1,DEdeg_1,Dist_1,pmRA_1,
+ pmDE_1,RV_1,s_RAdeg,s_DEdeg,s_Dist,s_pmRA,s_pmDE,s_RV,
+ e_RAdeg_1,e_DEdeg_1,e_Dist_1,e_pmRA_1,e_pmDE_1,e_RV_1,rad,
+ arad,X_1,Y_1,Z_1,U_1,V_1,W_1,s_X,s_Y,s_Z,s_U,s_V,s_W,e_X,e_Y,
+ e_Z,e_U,e_V,e_W,v_Fe_H__1,e__Fe_H__1,n__Fe_H_,E_B_V_1,E_B_V_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end