Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/597/A47 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-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/A+A/597/A47 CARMENES input catalogue of M dwarfs II (Cortes-Contreras+ 2017)
*================================================================================
*CARMENES input catalogue of M dwarfs.
*II. High-resolution imaging with FastCam.
* Cortes-Contreras M., Bejar V.J.S., Caballero J.A., Gauza B., Montes D.,
* Alonso-Floriano F.J., Jeffers S.V., Morales J.C., Reiners A., Ribas I.,
* Schoefer P., Quirrenbach A., Amado P.J., Mundt R., Seifert W.
* <Astron. Astrophys. 597, A47 (2017)>
* =2017A&A...597A..47C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Log of observed stars: list of 490 observed
M-dwarf targets
integer*4 nr__
parameter (nr__=490) ! Number of records
character*267 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq ! Identification number
character*11 Karmn ! Our Carmencita identifier, JHHMMm+DDdA (1)
character*21 Name ! Common name
integer*4 RAh ! (h) Right ascension (J2000) (2)
integer*4 RAm ! (min) Right ascension (J2000) (2)
real*4 RAs ! (s) Right ascension (J2000) (2)
character*1 DE_ ! Declination sign (J2000) (2)
integer*4 DEd ! (deg) Declination (J2000) (2)
integer*4 DEm ! (arcmin) Declination (J2000) (2)
real*4 DEs ! (arcsec) Declination (J2000) (2)
real*4 Jmag ! (mag) 2MASS J-band magnitude (2)
character*6 SpType ! MK spectral type
character*5 r_SpType ! Reference of the spectral type (3)
real*4 Dist ! (pc) ?=- Distance
character*9 r_Dist ! Reference of the distance (4)
character*11 Obs_date1 ! FastCam observation date
character*7 ExpTime1 ! Exposure time for Obs.date1 (in s)
character*11 Obs_date2 ! FastCam observation date
character*7 ExpTime2 ! Exposure time for Obs.date2 (in s)
character*11 Obs_date3 ! FastCam observation date
character*7 ExpTime3 ! Exposure time for Obs.date3 (in s)
character*11 Obs_date4 ! FastCam observation date
character*7 ExpTime4 ! Exposure time for Obs.date4 (in s)
character*11 Obs_date5 ! FastCam observation date
character*7 ExpTime5 ! Exposure time for Obs.date5 (in s)
character*11 Obs_date6 ! FastCam observation date
character*7 ExpTime6 ! Exposure time for Obs.date6 (in s)
character*11 Obs_date7 ! FastCam observation date
character*7 ExpTime7 ! Exposure time for Obs.date7 (in s)
character*11 Obs_date8 ! FastCam observation date
character*7 ExpTime8 ! Exposure time for Obs.date8 (in s)
*Note (1): Carmencita identifier (Quirrenbach et al. 2015, 18th Cambridge
* Workshop on Cool Stars, Stellar Systems, and the Sun, 18, 897;
* Alonso-Floriano et al., Cat. J/A+A/577/A128).
*Note (2): from 2MASX (Skrutskie et al., 2006, Cat. VII/233)
*Note (3): References for spectral type as follows:
* AF15 : Alonso-Floriano et al., 2015, Cat. J/A+A/577/A128
* Cru03: Cruz et al., 2003, Cat. J/AJ/126/2421
* Dea12: Deacon et al., 2012ApJ...757..100D
* Gig10: Gigoyan et al., 2010, Cat. J/other/Ap/53.123
* Gra03: Gray et al., 2003, Cat. J/AJ/126/2048
* Gra06: Gray et al., 2006, Cat. J/AJ/132/161
* Klutsch et al. priv. comm.
* Koe10: Koen et al., 2010, Cat. J/MNRAS/403/1949
* Lep13: Lepine et al., 2013, Cat. J/AJ/145/102
* Mon01: Montes et al., 2001, Cat. J/MNRAS/328/45
* New14: Newton et al., 2014, Cat. J/AJ/147/20
* PMSU : (Palomar/Michigan State University survey catalogue of nearby stars)
* Reid et al., 1995AJ....110.1838R and Hawley et al., 1996AJ....112.2799H,
* Cat. III/198, Gizis et al. 2002, Cat. J/AJ/123/3356
* Ria06: Riaz et al., 2006, Cat. J/AJ/132/866
* Sch05: Scholz et al., 2005, Cat. J/A+A/442/211
* Sim15: Simon-Diaz et al., 2015ApJ...799..169S
* ZS04 : Zuckerman & Song, 2004, Cat, J/other/ARA+A/42.685
*Note (4): References for Distance as follows:
* Ben00: Benedict et al., 2000AJ....120.1106B
* Cru03: Cruz et al., 2003, Cat. J/AJ/126/2421
* Daw05: Dawson et al., 2005, Cat. J/PASP/117/1
* Dit14: Dittmann et al., 2014ApJ...784..156D
* GC09 : Gatewood & Coban, 2009AJ....137..402G
* HD80 : Harrington & Dahn, 1980AJ.....85..454H
* Hen06: Henry et al., 2006AJ....132.2360H
* HT98 : Hershey & Taff, 1998AJ....116.1440H
* HIP2 : van Leeuwen 2007, Cat. I/311
* Ire08: Ireland et al., 2008ApJ...678..463I
* Jen52: Jenkins, 1952, General catalogue of trigonometric stellar parallaxes,
* Yale University Observatory, USA
* Jen63: Jenkins, 1963, General catalogue of trigonometric stellar parallaxes,
* Yale University Observatory, USA
* Jen09: Jenkins et al., 2009, Cat. J/ApJ/704/975
* Lep13: Lepine et al., 2013, Cat. J/AJ/145/102
* New14: Newton et al., 2014, Cat. J/AJ/147/20
* PMSU : (Palomar/Michigan State University survey catalogue of nearby stars)
* Reid et al., 1995AJ....110.1838R and Hawley et al., 1996AJ....112.2799H,
* Cat. III/198, Gizis et al. 2002, Cat. J/AJ/123/3356
* Rei02 : Reid et al., 2002, Cat. J/AJ/123/2806
* Ried10: Reidel et al., 2010AJ....140..897R
* Ria06 : Riaz et al., 2006, Cat. J/AJ/132/866
* Ried14: Riedel et al., 2014, Cat. J/AJ/147/85
* Sub09 : Subasavage et al., 2009AJ....137.4547S
* vAl95 : van Altena et al., 1995, General catalogue of trigonometric
* stellar parallaxes, Yale University Observatory, USA, Cat. I/238
* Wein16: Weinberger et al., 2016AJ....152...24W
*
* An "a" or "b" preceding the reference indicates that no measure for this
* component was found but we used instead the measure of the A or B companion,
* respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! ADS standard stars
integer*4 nr__1
parameter (nr__1=68) ! Number of records
character*94 ar__1 ! Full-size record
integer*4 ADS ! ADS number (Aitken, 1932ADS...C......0A)
character*6 Comp ! Component
character*1 n_ADS ! [b] Note on ADS 4241 (A-B) (1)
real*8 rhoL ! (arcsec) Previously published value of the angular
* separation
real*4 e_rhoL ! (arcsec) ? rms uncertainty on rhoL
real*4 thetaL ! (deg) Previously published value of the
* position angle
real*4 e_thetaL ! (deg) ? rms uncertainty on thetaL
character*9 EpL ! (yr) Epoch of observation of the published angular
* separation and position angle values
* (JYYYY.yyy)
character*6 r_rhoL ! Reference of the published values (2)
real*4 rho ! (arcsec) Measured angular separation with FastCam
real*4 e_rho ! (arcsec) rms uncertainty on rho
real*4 theta ! (deg) Measured position angle with FastCam
real*4 e_theta ! (deg) rms uncertainty on theta
character*9 Ep ! (yr) Epoch of observation with FastCam (JYYYY.yyy)
*Note (1): DeltaI=0.35mag for A-B and DeltaI=5.39mag for AB-C (this work).
*Note (2): References as follows:
* Ben14 = Benavides, 2014, El Observador de Estrellas Dobles, 12, 21
* Des11 = Desidera et al., 2011, Cat. J/A+A/533/A90
* Har11 = Hartkopf et al., 2011, Cat. J/AJ/142/56
* Mas12 = Mason et al., 2012, Cat. J/AJ/143/124
* Mas13 = Mason et al., 2013, Cat. J/AJ/146/56
* Sca11 = Scardia et al., 2011AN....332..508S
* Sca13 = Scardia et al., 2013MNRAS.434.2803S
* Schl13 = Schlimmer, 2013, Journal of Double Star Observations, 9, 230
* Sim15 = Simon-Diaz et al., 2015ApJ...799..169S
* Thor11 = Thorel et al., 2011, Observations & Travaux 78, 20
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! List of observed stars with confirmed visual
(unbound) companions
integer*4 nr__2
parameter (nr__2=69) ! Number of records
character*33 ar__2 ! Full-size record
character*11 Karmn_1 ! Carmencita identifier (JHHMMm+DDdA)
integer*4 CompNo ! Companion identification number
real*4 rho_1 ! (arcsec) measured angular separation with FastCam
integer*4 theta_1 ! (deg) measured position angle with FastCam
character*9 Ep_1 ! (yr) epoch of observation with FastCam (JYYYY.yyy)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea4.dat' ! Astrometric properties of the physically and
likely bound binaries in the sample
integer*4 nr__3
parameter (nr__3=192) ! Number of records
character*99 ar__3 ! Full-size record
character*1 Type ! [PL] Physically or Likely bound systems
character*11 Karmn_2 ! Carmencita identifier with component
character*5 m_Karmn ! Multiplicity index on Karmn
character*1 n_Karmn ! Note (1)
character*11 WDS ! Washington Double Star catalogue name or New
character*7 Disc ! Discoverer code (Mason et al.,
* 2001AJ....122.3466M, Cat. B/wds)
character*4 Ref ! Reference when no WDS name (2)
character*1 l_rho ! [~] Limit flag on rho
real*4 rho_2 ! (arcsec) ?=- Measured angular separation with FastCam
real*4 e_rho_1 ! (arcsec) ? rms uncertainty on rho
character*1 l_theta ! [~] Limit flag on theta
real*8 theta_2 ! (deg) ?=- Measured position angle with FastCam
real*4 e_theta_1 ! (deg) ? rms uncertainty on theta
real*4 DeltaI ! (mag) ?=- Measured magnitude difference in the I-band
* with FastCam
real*4 e_DeltaI ! (mag) ? rms uncertainty on DeltaI
character*9 Ep_2 ! (yr) Epoch of observation with FastCam (JYYYY.yyy)
character*1 GSC1_2 ! [FT-] Multiplicity flag from version 1.2 of
* the Guide Star Catalog (Morrison et al.,
* 2001AJ....121.1752M, Cat. I/254) (3)
*Note (1): Individual notes as follows:
* b = rho and theta measures in the second epoch correspond to observations with
* the CARMENES acquisition and guiding camera, which were carried out only
* for confirming physical association.
* c = Similar periods between Behall & Harrington (1976PASP...88..204B) and this
* work suggest that we are resolving the same pair.
* d = The B component is a double lined spectroscopic binary with a period
* shorter than 215.0d (Shkolnik et al., 2010ApJ...716.1522S) and is not
* the third component resolved here.
* e = We considered that this resolved binary is the spectroscopic binary
* identified by Bonfils et al. (2013, Cat. J/A+A/549/A109), according to
* the radial-velocity amplitude and period estimation of the pair.
* Nevertheless, we can not affirm whether the system is double or triple.
* f = Close visual binary in Bowler et al. (2015, Cat. J/ApJS/216/7) in one
* epoch. Our multi-epoch images confirm physical binding.
* g = It is also a spectroscopic binary.
* h = Faint pair for which we could not measure rho and theta with precision.
* The companion, also detected by Tokovinin et al. (2015,
* Cat. J/AJ/150/50), seems to be physically related and consistent over
* time with the companion identified by Montagnier et al.
* (2006A&A...460L..19M) at 0.068 arcsec.
* i = Also observed and identified as single in Law et al.
* (2008MNRAS.384..150L), probably due to the crossing of the companion
* behind the primary at the observing epoch, in agreement with the
* 13 a period estimated in this work.
*Note (2): References as follows:
* WD15 = Ward-Duong et al. (2015, Cat. J/MNRAS/449/2618)
* A15 = Ansdell et al. (2015, Cat. J/ApJ/798/41)
*Note (3): Guide Star Catalog multiplicity flag as follows:
* F = false
* T = true
* - = no data
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea5.dat' ! Derived parameters for confirmed physical pairs
integer*4 nr__4
parameter (nr__4=97) ! Number of records
character*95 ar__4 ! Full-size record
character*11 Karmn_3 ! Carmencita identifier
character*1 n_Karmn_1 ! [Y] Y for stars associated to young stellar
* populations (see Table 3 in the article)
character*6 Comp_1 ! Stellar component (A for the primary, B for the
* secondary and so on)
real*4 DeltaI_1 ! (mag) ?=- Difference in the I-band between components
real*4 e_DeltaI_1 ! (mag) ? rms uncertainty on DeltaI
real*4 Imag1 ! (mag) ?=- Individual I magnitude of the primary
real*4 Imag2 ! (mag) ?=- Individual I magnitude of the secondary
character*6 SpType_1 ! Spectral type of the system
character*7 SpType1 ! Individual spectral type of the primary
character*4 SpType2 ! Individual spectral type of the secondary
real*4 M1 ! (Msun) ?=- Individual mass of the primary
real*4 e_M1 ! (Msun) ? rms uncertainty on M1
character*1 l_M2 ! Limit flag on M2
real*4 M2 ! (Msun) ?=- Individual mass of the secondary
real*4 e_M2 ! (Msun) ? rms uncertainty on M2
real*8 Porb ! (a) ?=- Estimated orbital period (1)
character*1 n_Porb ! [c] Note on Porb (2)
*Note (1): The periods given are a lower limit, as we equal our maximum
* separation measured to the semimajor axis.
*Note (2): c: The B component of the system is too faint to estimate its mass
* with our MI-mass relation and the period was calculated using the stellar mass
* limit (0.07M_{sun}_). Hence, the period of this system should be considered as
* a lower approximation.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea6.dat' ! Known binaries at angular separations larger
than 5 arcsec
integer*4 nr__5
parameter (nr__5=81) ! Number of records
character*99 ar__5 ! Full-size record
character*10 WDS_1 ! Washington Double Star catalogue discoverer
* code (Mason et al.,
* 2001AJ....122.3466M, Cat. B/wds)
character*19 Name1 ! Carmencita name or common name of the primary
character*16 SpType1_1 ! spectral type of the primary
character*23 Name2 ! Carmencita name or common name of the secondary
character*18 SpType2_1 ! spectral type of the secondary
real*8 Sep ! (arcsec) published angular separation
character*1 Notes ! Notes (1)
*Note (1): Notes as follows:
* a = C is background (GRB34)
* b = AB is separated by 2.5 arcsec (WOR4)
* c = "B-G" are background (LMP3)
* d = Simbad indicates that J03574-011 is a spectroscopic binary but we did
* not find any reference
* e = Triple system
* f = Primary is SB2 (Llamas, 2014, MSc thesis, Universidad Complutense de
* Madrid, Spain)
* g = Primary is an astrometric binary separated by 0.07 arcseconds
* (Strand, 1977AJ.....82....9S)
* h = The primary is a SB2 (Schoefer, 2015, MSc thesis, Institut fur Astrophysik
* Goettingen, Germany)
* i = It was also observed with the CAMELOT low resolution imager at the
* Observatorio del Teide (Tenerife) in September 2015 with the BVIgri
* filters in order to obtain more photometric information but we could
* not avoid the saturation of the primary
* j = BC is separated by 0.8 arcsec
* k = Primary is a SB (Simbad)
* l = Background source at 9.62 arcsec (TNN6)
* m = WDS VBS41 at 4 arcsec and 208 deg was not detected in this work nor in
* Davison et al. (2015AJ....149..106D). It could be an unrelated companion.
* The WDS "AC" designation refers to the pair in the table
* n = Hierarchical quadruple
* o = Primary is a SB1 (Schofer 2015, MSc thesis, Institut fur Astrophysik
* Goettingen, Germany). Two other WDS entries under
* the same discoverer code (STF1321) are not physically bound components
* p = SB2 (Bonfils et al., 2013, Cat. J/A+A/549/A109). Hawley et al.
* (1996AJ....112.2799H, Cat. III/198) listed in Table 1.(b) a companion at
* 5.0-arcsec 1.2mag fainter in V.
* 2MASS resolved a source 6.5 mag fainter in the J-band at 5.9-arcsec and
* 179deg (quality flag: AUU). Neither Law et al. (2008MNRAS.384..150L)
* nor Dieterich et al. (2012, Cat. J/AJ/144/64 - with NICMOS
* onboard Hubble detected it. We believe that Hawley et al.
* (1996AJ....112.2799H, Cat. III/198) made reference to the background
* star 2MASS J12141817+0037297
* q = Other WDS entries under BU 1442 and STG 6 are unrelated sources
* r = Quintuple system
* s = Simbad mixes up the true primary GJ 3999 ("L 845-016"), the true secondary
* 2MASS J17174454-1148261, and the background star GJ 4000 B
* ("L 845-015", see Table A.3).
* We preserve the current (wrong) nomenclature
* t = C is background (HEL3)
* u = Other WDS entries under RAO 23 are unrelated sources
* v = Secondary is SB1 (Schoefer, 2015, MSc thesis, Institut fur Astrophysik
* Goettingen, Germany)
* w = AB is separated by 0.05 arcsec (BWL56)
* x = Bab is separated by 0.26 arcsec (BWL59). AC is background (BWL59)
* y = A and B are SB1 (Delfosse et al., 1999A&A...344..897D).
* Other WDS entries under LMP 24 are unrelated sources
* z = Secondary is SB2 (Schoefer, 2015, MSc thesis, Institut fur Astrophysik
* Goettingen, Germany)
C=============================================================================
C Loading file 'tablea1.dat' ! Log of observed stars: list of 490 observed
* M-dwarf targets
C Format for file interpretation
1 format(
+ I3,1X,A11,2X,A21,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,
+ 1X,F6.3,1X,A6,2X,A5,1X,F6.3,1X,A9,1X,A11,1X,A7,2X,A11,1X,A7,
+ 2X,A11,1X,A7,2X,A11,1X,A7,2X,A11,1X,A7,2X,A11,1X,A7,2X,A11,1X,
+ A7,2X,A11,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,490
read(1,'(A267)')ar__
read(ar__,1)
+ Seq,Karmn,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Jmag,SpType,
+ r_SpType,Dist,r_Dist,Obs_date1,ExpTime1,Obs_date2,ExpTime2,
+ Obs_date3,ExpTime3,Obs_date4,ExpTime4,Obs_date5,ExpTime5,
+ Obs_date6,ExpTime6,Obs_date7,ExpTime7,Obs_date8,ExpTime8
if (idig(ar__(85:90)).EQ.0) Dist = rNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Seq,Karmn,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Jmag,SpType,
+ r_SpType,Dist,r_Dist,Obs_date1,ExpTime1,Obs_date2,ExpTime2,
+ Obs_date3,ExpTime3,Obs_date4,ExpTime4,Obs_date5,ExpTime5,
+ Obs_date6,ExpTime6,Obs_date7,ExpTime7,Obs_date8,ExpTime8
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! ADS standard stars
C Format for file interpretation
2 format(
+ I5,1X,A6,A1,1X,F7.4,1X,F6.4,1X,F6.2,1X,F4.2,3X,A9,1X,A6,1X,
+ F6.3,1X,F5.3,1X,F6.2,1X,F4.2,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,68
read(1,'(A94)')ar__1
read(ar__1,2)
+ ADS,Comp,n_ADS,rhoL,e_rhoL,thetaL,e_thetaL,EpL,r_rhoL,rho,
+ e_rho,theta,e_theta,Ep
if(ar__1(23:28) .EQ. '') e_rhoL = rNULL__
if(ar__1(37:40) .EQ. '') e_thetaL = rNULL__
c ..............Just test output...........
write(6,2)
+ ADS,Comp,n_ADS,rhoL,e_rhoL,thetaL,e_thetaL,EpL,r_rhoL,rho,
+ e_rho,theta,e_theta,Ep
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea3.dat' ! List of observed stars with confirmed visual
* (unbound) companions
C Format for file interpretation
3 format(A11,1X,I1,1X,F5.2,1X,I3,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,69
read(1,'(A33)')ar__2
read(ar__2,3)Karmn_1,CompNo,rho_1,theta_1,Ep_1
c ..............Just test output...........
write(6,3)Karmn_1,CompNo,rho_1,theta_1,Ep_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea4.dat' ! Astrometric properties of the physically and
* likely bound binaries in the sample
C Format for file interpretation
4 format(
+ A1,1X,A11,A5,2X,A1,1X,A11,1X,A7,1X,A4,3X,A1,F6.3,1X,F5.3,1X,
+ A1,F7.3,1X,F5.2,1X,F4.2,1X,F4.2,1X,A9,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea4.dat')
write(6,*) '....Loading file: tablea4.dat'
do i__=1,192
read(1,'(A99)')ar__3
read(ar__3,4)
+ Type,Karmn_2,m_Karmn,n_Karmn,WDS,Disc,Ref,l_rho,rho_2,e_rho_1,
+ l_theta,theta_2,e_theta_1,DeltaI,e_DeltaI,Ep_2,GSC1_2
if (idig(ar__3(51:56)).EQ.0) rho_2 = rNULL__
if(ar__3(58:62) .EQ. '') e_rho_1 = rNULL__
if (idig(ar__3(65:71)).EQ.0) theta_2 = rNULL__
if(ar__3(73:77) .EQ. '') e_theta_1 = rNULL__
if (idig(ar__3(79:82)).EQ.0) DeltaI = rNULL__
if(ar__3(84:87) .EQ. '') e_DeltaI = rNULL__
c ..............Just test output...........
write(6,4)
+ Type,Karmn_2,m_Karmn,n_Karmn,WDS,Disc,Ref,l_rho,rho_2,e_rho_1,
+ l_theta,theta_2,e_theta_1,DeltaI,e_DeltaI,Ep_2,GSC1_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea5.dat' ! Derived parameters for confirmed physical pairs
C Format for file interpretation
5 format(
+ A11,2X,A1,1X,A6,2X,F4.2,1X,F4.2,1X,F5.2,1X,F6.3,1X,A6,1X,A7,
+ 1X,A4,1X,F4.2,1X,F4.2,1X,A1,F4.2,1X,F4.2,1X,F7.1,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea5.dat')
write(6,*) '....Loading file: tablea5.dat'
do i__=1,97
read(1,'(A95)')ar__4
read(ar__4,5)
+ Karmn_3,n_Karmn_1,Comp_1,DeltaI_1,e_DeltaI_1,Imag1,Imag2,
+ SpType_1,SpType1,SpType2,M1,e_M1,l_M2,M2,e_M2,Porb,n_Porb
if (idig(ar__4(24:27)).EQ.0) DeltaI_1 = rNULL__
if(ar__4(29:32) .EQ. '') e_DeltaI_1 = rNULL__
if (idig(ar__4(34:38)).EQ.0) Imag1 = rNULL__
if (idig(ar__4(40:45)).EQ.0) Imag2 = rNULL__
if (idig(ar__4(67:70)).EQ.0) M1 = rNULL__
if(ar__4(72:75) .EQ. '') e_M1 = rNULL__
if (idig(ar__4(78:81)).EQ.0) M2 = rNULL__
if(ar__4(83:86) .EQ. '') e_M2 = rNULL__
if (idig(ar__4(88:94)).EQ.0) Porb = rNULL__
c ..............Just test output...........
write(6,5)
+ Karmn_3,n_Karmn_1,Comp_1,DeltaI_1,e_DeltaI_1,Imag1,Imag2,
+ SpType_1,SpType1,SpType2,M1,e_M1,l_M2,M2,e_M2,Porb,n_Porb
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea6.dat' ! Known binaries at angular separations larger
* than 5 arcsec
C Format for file interpretation
6 format(A10,1X,A19,1X,A16,A23,1X,A18,1X,F7.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea6.dat')
write(6,*) '....Loading file: tablea6.dat'
do i__=1,81
read(1,'(A99)')ar__5
read(ar__5,6)WDS_1,Name1,SpType1_1,Name2,SpType2_1,Sep,Notes
c ..............Just test output...........
write(6,6)WDS_1,Name1,SpType1_1,Name2,SpType2_1,Sep,Notes
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end