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