Conversion of standardized ReadMe file for
file /./ftp/cats/J/PASP/109/849 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-18
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. J/PASP/109/849 Probing the LHS Catalog (Gizis+ 1997)
*================================================================================
*Probing the LHS Catalog I. New Nearby Stars and the Coolest Subdwarf
* Gizis J.E., Reid, I.N.
* <Publ. Astron. Soc. Pac. 109, 849 (1997)>
* =1997PASP..109..849G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Candidate nearby stars: basic data (table 1)
integer*4 nr__
parameter (nr__=71) ! Number of records
character*73 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*13 Name (nr__) ! *Star name
character*1 Bin (nr__) ! *[ps] Primary or secondary binary comp.
character*1 Note1 (nr__) ! *[bcde] Individual Note
integer*4 RAh (nr__) ! (h) Right ascension, J2000
integer*4 RAm (nr__) ! (min) Right ascension, J2000
real*4 RAs (nr__) ! (s) Right ascension, J2000
character*1 DE_ (nr__) ! Declination sign
integer*4 DEd (nr__) ! (deg) Declination, J2000
integer*4 DEm (nr__) ! (arcmin) Declination, J2000
integer*4 DEs (nr__) ! (arcsec) Declination, J2000
integer*4 Src (nr__) ! *[5] Source of positional data
real*4 MV (nr__) ! (mag) Absolute V magnitude
real*4 r (nr__) ! (pc) *Distance
integer*4 er (nr__) ! (%) Distance uncertainty
integer*4 wPi (nr__) ! (%) *? Weight for trigonometric parallax
integer*4 wSp (nr__) ! (%) ? Weight for spectroscopic parallax
character*4 SpType (nr__) ! Spectral type
real*4 Vmag (nr__) ! (mag) Apparent V magnitude
*Note on Name:
* Gl - Gliese (1969, Veroff. Astr. Rechen-Instituts, Heidelberg, Nr. 22)
* GJ - Gliese & Jahreiss =1979A&AS...38..423G
* Wo - Woolley et al (1970, R. Obs. Ann., No. 5) =Cat. <V/32>
* LHS - Luyten Half-Second Catalogue =Cat. <I/87>
* LTT - Luyten Two-Tenths Catalogue
* LP - Luyten Palomar proper-motion catalogue
* G - Giclas et al (1971, Lowell Proper Motion Survey) =Cat. <I/79>
* Steph - Stephenson (1986, AJ, 91, 137) =Cat. <III/123>
* (Nomenclature note: should be named "StKM")
*Note on Bin:
* Components in a binary system are identified as
* p (primary) or s (secondary, tertiary).
*Note on Note1:
* b G 035-027 and LP 469-118: Photometry based upon Fleming
* et al. (1988) observed V_{system} = 14.20 and assumed
* Delta V = 0.42
* c Steph 2.908: This is the 809th star in Stephenson (1986).
* Photometry is from Weis (1991).
* d BD+28 3698: Jahreiss notes that this object is presumably a dK
* binary with Delta m = 0.8, see Heintz (1993). The distance
* was estimated by Jahreiss.
* e LHS 3705: This is also Cou 2234. Estimated distance assumes
* the system is a binary with Delta m = 0.0 as determined by
* Couteau (1985).
*Note on Src:
* Source of positional data:
*
* 5- Digital Sky Survey
* Codes 1-4 were used in Reid et al. (1995),Hawley et al. (1996),
* Hawley+ (1997)
*Note on r:
* Distance and estimated uncertainty (%) from the authors' data,
* using a weighted combination of trigonometric and spectroscopic
* parallaxes.
*Note on wPi:
* An entry of -1 indicates that Jahreiss's distance estimates were adopted.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Candidate nearby stars: bandstrengths (table 2)
integer*4 nr__1
parameter (nr__1=71) ! Number of records
character*77 ar__1 ! Full-size record
character*13 Name_1 (nr__1) ! Star name
real*4 TiO1 (nr__1) ! TiO1 index: F(6718-6723)/F(6703-6708)
real*4 TiO2 (nr__1) ! TiO2 index: F(7058-7061)/F(7043-7046)
real*4 TiO3 (nr__1) ! TiO3 index: F(7092-7097)/F(7079-7084)
real*4 TiO4 (nr__1) ! TiO4 index: F(7130-7135)/F(7115-7120)
real*4 TiO5 (nr__1) ! TiO5 index: F(7126-7135)/F(7042-7046)
real*4 CaH1 (nr__1) ! CaH1 index: F(6380-6390)/F
* (6345-6355,6410-6420))
real*4 CaH2 (nr__1) ! CaH2 index: F(6814-6846)/F(7042-7046)
real*4 CaH3 (nr__1) ! CaH3 index: F(6960-6990)/F(7042-7046)
real*4 CaOH (nr__1) ! ? CaOH index: F(6230-6240)/F(6345-6354)
real*4 Halpha (nr__1) ! (0.1nm) ? H-alpha Eq. Width
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Candidate nearby stars: Space Velocities
integer*4 nr__2
parameter (nr__2=71) ! Number of records
character*56 ar__2 ! Full-size record
character*13 Name_2 (nr__2) ! Star name
real*8 pmRA (nr__2) ! (arcsec/yr) Proper motion in right ascension,
* from CNS3
real*8 pmDE (nr__2) ! (arcsec/yr) Proper motion in declination, from CNS3
real*4 RV (nr__2) ! (km/s) Radial velocity
integer*4 U (nr__2) ! (km/s) U space velocity
integer*4 V (nr__2) ! (km/s) V space velocity
integer*4 W (nr__2) ! (km/s) W space velocity
real*4 MV_1 (nr__2) ! (mag) Absolute V magnitude (from table1.dat)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Faint LHS stars
integer*4 nr__3
parameter (nr__3=40) ! Number of records
character*58 ar__3 ! Full-size record
character*8 Name_3 (nr__3) ! *Star name
character*1 Note1_1 (nr__3) ! *[a] Individual Note
real*4 mr (nr__3) ! r magnitude from Luyten (1979)
real*4 TiO5_1 (nr__3) ! TiO5 index: F(7126-7135)/F(7042-7046)
real*4 CaH1_1 (nr__3) ! CaH1 index: F(6380-6390)/F
* (6345-6355,6410-6420))
real*4 CaH2_1 (nr__3) ! CaH2 index: F(6814-6846)/F(7042-7046)
real*4 CaH3_1 (nr__3) ! CaH3 index: F(6960-6990)/F(7042-7046)
real*4 Halpha_1 (nr__3) ! (0.1nm) ? H alpha Equivalent Width
character*7 SpType_1 (nr__3) ! Spectral Type
integer*4 RV_1 (nr__3) ! (km/s) ? Radial Velocity
*Note on Name:
* LHS - Luyten Half-Second Catalogue =Cat. <I/87>
*Note on Note1:
* a Probably misidentified.
C=============================================================================
C Loading file 'table1.dat' ! Candidate nearby stars: basic data (table 1)
C Format for file interpretation
1 format(
+ 1X,A13,1X,A1,1X,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,
+ 1X,I1,1X,F5.2,1X,F4.1,1X,I2,1X,I2,1X,I3,1X,A4,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,71
read(1,'(A73)')ar__
read(ar__,1)
+ Name(i__),Bin(i__),Note1(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Src(i__),MV(i__),r(i__),
+ er(i__),wPi(i__),wSp(i__),SpType(i__),Vmag(i__)
if(ar__(57:58) .EQ. '') wPi(i__) = iNULL__
if(ar__(60:62) .EQ. '') wSp(i__) = iNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),Bin(i__),Note1(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),Src(i__),MV(i__),r(i__),
+ er(i__),wPi(i__),wSp(i__),SpType(i__),Vmag(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Candidate nearby stars: bandstrengths (table 2)
C Format for file interpretation
2 format(
+ 2X,A13,3X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,71
read(1,'(A77)')ar__1
read(ar__1,2)
+ Name_1(i__),TiO1(i__),TiO2(i__),TiO3(i__),TiO4(i__),TiO5(i__),
+ CaH1(i__),CaH2(i__),CaH3(i__),CaOH(i__),Halpha(i__)
if(ar__1(67:71) .EQ. '') CaOH(i__) = rNULL__
if(ar__1(73:77) .EQ. '') Halpha(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),TiO1(i__),TiO2(i__),TiO3(i__),TiO4(i__),TiO5(i__),
+ CaH1(i__),CaH2(i__),CaH3(i__),CaOH(i__),Halpha(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Candidate nearby stars: Space Velocities
C Format for file interpretation
3 format(1X,A13,1X,F7.4,1X,F7.4,1X,F6.1,1X,I3,1X,I4,1X,I3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,71
read(1,'(A56)')ar__2
read(ar__2,3)
+ Name_2(i__),pmRA(i__),pmDE(i__),RV(i__),U(i__),V(i__),W(i__),
+ MV_1(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),pmRA(i__),pmDE(i__),RV(i__),U(i__),V(i__),W(i__),
+ MV_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Faint LHS stars
C Format for file interpretation
4 format(
+ 1X,A8,2X,A1,1X,F4.1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F3.1,
+ 1X,A7,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,40
read(1,'(A58)')ar__3
read(ar__3,4)
+ Name_3(i__),Note1_1(i__),mr(i__),TiO5_1(i__),CaH1_1(i__),
+ CaH2_1(i__),CaH3_1(i__),Halpha_1(i__),SpType_1(i__),RV_1(i__)
if(ar__3(43:45) .EQ. '') Halpha_1(i__) = rNULL__
if(ar__3(55:58) .EQ. '') RV_1(i__) = iNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),Note1_1(i__),mr(i__),TiO5_1(i__),CaH1_1(i__),
+ CaH2_1(i__),CaH3_1(i__),Halpha_1(i__),SpType_1(i__),RV_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end