Conversion of standardized ReadMe file for
file /./ftp/cats/V/106 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-16
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/106 Catalogue of X-ray binaries (Liu+ 2001)
*================================================================================
*A catalogue of high- and low- mass X-ray binaries
* Liu Q.Z., van Paradijs J., van den Heuvel E.P.J.
* <Astron. Astrophys. Suppl. Ser. 147, 25 (2000)>
* =2000A&AS..147...25L
* <Astron. Astrophys. Suppl. Ser. 368, 1021 (2001)>
* =2001A&A...368.1021L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hmxbcat.dat' ! High-mass X-ray binaries catalog
integer*4 nr__
parameter (nr__=130) ! Number of records
character*359 ar__ ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000.0
real*8 DEdeg (nr__) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*16 Name (nr__) ! Source Name (1)
character*6 Type (nr__) ! [ABTUDGPVZ (?)] X-ray type (2)
integer*4 RAh (nr__) ! (h) Right ascension (J2000.0) (9)
integer*4 RAm (nr__) ! (min) Right ascension (J2000.0) (9)
real*4 RAs (nr__) ! (s) ? Right ascension (J2000.0) (9)
character*1 u_RAs (nr__) ! Uncertainty flag on RAs
character*1 DE_ (nr__) ! Declination sign
integer*4 DEd (nr__) ! (deg) Declination (J2000.0) (9)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000.0) (9)
real*4 DEs (nr__) ! (arcsec) ? Declination (J2000.0) (9)
character*1 u_DEs (nr__) ! Uncertainty flag on DEs
real*4 GLON (nr__) ! (deg) Galactic longitude
real*8 GLAT (nr__) ! (deg) Galactic latitude
character*4 Pos (nr__) ! Type of observation from which the source
* position has been derived (3)
real*8 e_Pos (nr__) ! (arcsec) ? Accuracy of the position (4)
character*17 Opt (nr__) ! Optical counterpart (5)
character*15 r_Opt (nr__) ! References for Optical counterpart (finding
* charts) (in [lh]mxbrefs.dat files)
real*4 Vmag (nr__) ! (mag) ? V magnitude (6)
real*4 Vmagl (nr__) ! (mag) ? Lower magnitude when interval
character*1 u_Vmag (nr__) ! Uncertainty flag on Vmag
real*4 B_V (nr__) ! (mag) ? B-V colour index
character*1 u_B_V (nr__) ! Uncertainty flag on B-V
real*4 B_Vl (nr__) ! (mag) ? Lower B-V colour index when interval
real*4 U_B (nr__) ! (mag) ? U-B colour index
character*1 l_E_B_V (nr__) ! [~<>] Limit flag on E(B-V)
real*4 E_B_V (nr__) ! (mag) ? Interstellar reddening
character*2 l_E_B_V_2 (nr__) ! [>= ] Limit flag on E(B-V)2
real*4 E_B_V_2 (nr__) ! (mag) ? Second value of E(B-V) when interval
character*1 u_E_B_V (nr__) ! [:P] Uncertainty flag on E(B-V)
character*16 r_Vmag (nr__) ! References for UBV photometry,
* in [lh]mxbrefs.dat files
character*1 l_Fx (nr__) ! Limit flag on Fx
real*8 Fx (nr__) ! (uJy) ? Average X-ray flux (2-10keV)
real*8 Fxu (nr__) ! (uJy) ? Upper limit of X-ray flux when interval
character*14 Range (nr__) ! Energy range for flux if not 2-20keV
character*18 r_Fx (nr__) ! References for X-ray flux,
* in [lh]mxbrefs.dat files
real*8 Porb (nr__) ! (d) ? Orbital period
real*4 Porb2 (nr__) ! (d) ? Second value of Porb when interval
character*1 u_Porb (nr__) ! [?:d] Uncertainty flag on Porb
real*8 Ppulse (nr__) ! (s) ? Pulse period for X-ray pulsars
character*1 u_Ppulse (nr__) ! [?] Uncertainty flag on Ppulse
character*15 r_Ppulse (nr__) ! References for Pulse period,
* in [lh]mxbrefs.dat files
character*25 Cat (nr__) ! Survey catalogues and experiments in which the
* source was listed and detected, respectively (7)
character*12 SpType (nr__) ! Spectral type
character*13 Name2 (nr__) ! Other name
character*1 u_Name2 (nr__) ! [?] Uncertainty flag on Name2
character*15 Name3 (nr__) ! Other name (8)
*Note (1): Source name, with rough information on its sky location according
* to the convention hhmm+/-ddd. Here hh and mm indicate the hours and
* minutes of right ascension, ddd the declination in units of 0.1 degree
* (in a small number of cases, the coordinates shown in the name are
* given with more, or fewer, digits). However, for a ROSAT source the
* name is always given in the form of hhmm.m+/-ddmm. The prefix J
* indicates a name based on J2000 coordinates. Otherwise, 1950
* coordinates were used in the name.
*Note (2): Type:
* P: X-ray pulsars
* T: transient X-ray source
* U: Ultra-soft X-ray spectrum These sources include black-hole
* candidates; some 'extreme ultra-soft' (EUS) source may be white
* dwarf (WD) on whose surface steady nuclear burning takes place.
* V: variable
* A: atoll source
* B: X-ray burst source
* D: "dipping" low-mass X-ray binary
* G: globular-cluster X-ray source
* Z: Z-type
*Note (3): o: optical; x: X-ray; r: radio; IR: infrared
*Note (4): We give an indication of the accuracy of this position, in the
* form of equivalent (90 percent confidence level) error radii, but in
* several cases this can only be considered an approximation (e.g. when
* the error box is not circular). When no accuracy is quoted, it is
* about one arcsecond or better.
*Note (5): Name of the optical counterpart.
* An asterisk followed by a number or letter refers to star number used
* in the finding chart; "star" refers a star in the finding chart that
* has not been assigned a number or letter.
*Note (6): Individual notes for magnitudes:
* In hmxb.dat:
* J0529.8-6556: B=14.5,R=13.3
* J0532.4-6535: J=11.9,H=11.1, K=10.7
* 2030+407 : I=20.0
* In lmxb.dat:
* 1550-631 Bmag >=20.
* 1552-645 Bmag = 17.5
* 1636-550 Imag ~ 18.2
* 1756-365 Kmag = 15.1
* J1810.8-2609: B=21.5, R=19.5
* 1811-171 Kmag = 12
* 1820-303 Bmag = 18.7
* 1837+049 Bmag = 19.2
* 1915+105 Imag = 23.4
*Note (7): Survey catalogues and experiments abbreviations:
* A: Ariel V sky survey
* AS: ASCA
* B: BeppoSAX
* C: Compton {gamma}-ray Observatory
* E: Einstein Observatory
* Exo: Exosat
* G: Ginga
* Gr: Granat
* H: HEAO A-1 sky survey
* Ha: Hakucho
* I: Indian X-ray Astronomy Experiment (IXAE)
* K: Kvant
* M: Mit OSO-7 sky survey
* OAO: Orbiting Astronomical Observatory
* R: ROSAT
* S: SAS 3
* SL: Space Lab
* T: Tenma
* U: Uhuru sky survey
* V: Vela-5 and -6 satellites
* X: Rossi XTE
*Note (8): In hmxb.dat:
* For 1258-613, No. 2 in Mason et al., 1978MNRAS.184P..45M
* For 1538-522, No. 12 in Schwartz et al., 1978Natur.275..517S
*Note (9): The printed table gives positions in J2000 or B1950,
* depending whether the source name stars by a 'J' or not.
* (All positions were converted to J2000 in the electronic version at CDS)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lmxbcat.dat' ! Low-mass X-ray binaries catalog
integer*4 nr__1
parameter (nr__1=150) ! Number of records
character*359 ar__1 ! Full-size record
C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*16 Name_1 (nr__1) ! Source Name (1)
character*6 Type_1 (nr__1) ! [ABTUDGPVZ (?)] X-ray type (2)
integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000.0) (9)
integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000.0) (9)
real*4 RAs_1 (nr__1) ! (s) ? Right ascension (J2000.0) (9)
character*1 u_RAs_1 (nr__1) ! Uncertainty flag on RAs
character*1 DE__1 (nr__1) ! Declination sign
integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000.0) (9)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000.0) (9)
real*4 DEs_1 (nr__1) ! (arcsec) ? Declination (J2000.0) (9)
character*1 u_DEs_1 (nr__1) ! Uncertainty flag on DEs
real*4 GLON_1 (nr__1) ! (deg) Galactic longitude
real*8 GLAT_1 (nr__1) ! (deg) Galactic latitude
character*4 Pos_1 (nr__1) ! Type of observation from which the source
* position has been derived (3)
real*8 e_Pos_1 (nr__1) ! (arcsec) ? Accuracy of the position (4)
character*17 Opt_1 (nr__1) ! Optical counterpart (5)
character*15 r_Opt_1 (nr__1) ! References for Optical counterpart (finding
* charts) (in [lh]mxbrefs.dat files)
real*4 Vmag_1 (nr__1) ! (mag) ? V magnitude (6)
real*4 Vmagl_1 (nr__1) ! (mag) ? Lower magnitude when interval
character*1 u_Vmag_1 (nr__1) ! Uncertainty flag on Vmag
real*4 B_V_1 (nr__1) ! (mag) ? B-V colour index
character*1 u_B_V_1 (nr__1) ! Uncertainty flag on B-V
real*4 B_Vl_1 (nr__1) ! (mag) ? Lower B-V colour index when interval
real*4 U_B_1 (nr__1) ! (mag) ? U-B colour index
character*1 l_E_B_V_1 (nr__1) ! [~<>] Limit flag on E(B-V)
real*4 E_B_V_1 (nr__1) ! (mag) ? Interstellar reddening
character*2 l_E_B_V_2_1(nr__1) ! [>= ] Limit flag on E(B-V)2
real*4 E_B_V_2_1 (nr__1) ! (mag) ? Second value of E(B-V) when interval
character*1 u_E_B_V_1 (nr__1) ! [:P] Uncertainty flag on E(B-V)
character*16 r_Vmag_1 (nr__1) ! References for UBV photometry,
* in [lh]mxbrefs.dat files
character*1 l_Fx_1 (nr__1) ! Limit flag on Fx
real*8 Fx_1 (nr__1) ! (uJy) ? Average X-ray flux (2-10keV)
real*8 Fxu_1 (nr__1) ! (uJy) ? Upper limit of X-ray flux when interval
character*14 Range_1 (nr__1) ! Energy range for flux if not 2-20keV
character*18 r_Fx_1 (nr__1) ! References for X-ray flux,
* in [lh]mxbrefs.dat files
real*8 Porb_1 (nr__1) ! (d) ? Orbital period
real*4 Porb2_1 (nr__1) ! (d) ? Second value of Porb when interval
character*1 u_Porb_1 (nr__1) ! [?:d] Uncertainty flag on Porb
real*8 Ppulse_1 (nr__1) ! (s) ? Pulse period for X-ray pulsars
character*1 u_Ppulse_1 (nr__1) ! [?] Uncertainty flag on Ppulse
character*15 r_Ppulse_1 (nr__1) ! References for Pulse period,
* in [lh]mxbrefs.dat files
character*25 Cat_1 (nr__1) ! Survey catalogues and experiments in which the
* source was listed and detected, respectively (7)
character*12 SpType_1 (nr__1) ! Spectral type
character*13 Name2_1 (nr__1) ! Other name
character*1 u_Name2_1 (nr__1) ! [?] Uncertainty flag on Name2
character*15 Name3_1 (nr__1) ! Other name (8)
*Note (1): Source name, with rough information on its sky location according
* to the convention hhmm+/-ddd. Here hh and mm indicate the hours and
* minutes of right ascension, ddd the declination in units of 0.1 degree
* (in a small number of cases, the coordinates shown in the name are
* given with more, or fewer, digits). However, for a ROSAT source the
* name is always given in the form of hhmm.m+/-ddmm. The prefix J
* indicates a name based on J2000 coordinates. Otherwise, 1950
* coordinates were used in the name.
*Note (2): Type:
* P: X-ray pulsars
* T: transient X-ray source
* U: Ultra-soft X-ray spectrum These sources include black-hole
* candidates; some 'extreme ultra-soft' (EUS) source may be white
* dwarf (WD) on whose surface steady nuclear burning takes place.
* V: variable
* A: atoll source
* B: X-ray burst source
* D: "dipping" low-mass X-ray binary
* G: globular-cluster X-ray source
* Z: Z-type
*Note (3): o: optical; x: X-ray; r: radio; IR: infrared
*Note (4): We give an indication of the accuracy of this position, in the
* form of equivalent (90 percent confidence level) error radii, but in
* several cases this can only be considered an approximation (e.g. when
* the error box is not circular). When no accuracy is quoted, it is
* about one arcsecond or better.
*Note (5): Name of the optical counterpart.
* An asterisk followed by a number or letter refers to star number used
* in the finding chart; "star" refers a star in the finding chart that
* has not been assigned a number or letter.
*Note (6): Individual notes for magnitudes:
* In hmxb.dat:
* J0529.8-6556: B=14.5,R=13.3
* J0532.4-6535: J=11.9,H=11.1, K=10.7
* 2030+407 : I=20.0
* In lmxb.dat:
* 1550-631 Bmag >=20.
* 1552-645 Bmag = 17.5
* 1636-550 Imag ~ 18.2
* 1756-365 Kmag = 15.1
* J1810.8-2609: B=21.5, R=19.5
* 1811-171 Kmag = 12
* 1820-303 Bmag = 18.7
* 1837+049 Bmag = 19.2
* 1915+105 Imag = 23.4
*Note (7): Survey catalogues and experiments abbreviations:
* A: Ariel V sky survey
* AS: ASCA
* B: BeppoSAX
* C: Compton {gamma}-ray Observatory
* E: Einstein Observatory
* Exo: Exosat
* G: Ginga
* Gr: Granat
* H: HEAO A-1 sky survey
* Ha: Hakucho
* I: Indian X-ray Astronomy Experiment (IXAE)
* K: Kvant
* M: Mit OSO-7 sky survey
* OAO: Orbiting Astronomical Observatory
* R: ROSAT
* S: SAS 3
* SL: Space Lab
* T: Tenma
* U: Uhuru sky survey
* V: Vela-5 and -6 satellites
* X: Rossi XTE
*Note (8): In hmxb.dat:
* For 1258-613, No. 2 in Mason et al., 1978MNRAS.184P..45M
* For 1538-522, No. 12 in Schwartz et al., 1978Natur.275..517S
*Note (9): The printed table gives positions in J2000 or B1950,
* depending whether the source name stars by a 'J' or not.
* (All positions were converted to J2000 in the electronic version at CDS)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hmxbnote.dat' ! Individual notes (HMXB)
integer*4 nr__2
parameter (nr__2=425) ! Number of records
character*80 ar__2 ! Full-size record
character*16 Name_2 (nr__2) ! Source name
character*63 Note (nr__2) ! Text of the note (1)
*Note (1): The following abbreviations have been used:
* o, optical; x, X-ray; r, radio; IR, infrared.
* A reference on the source position is given under 'pos.'.
* The references code are from [lh]mxbrefs.dat files
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lmxbnote.dat' ! Individual notes (LMXB)
integer*4 nr__3
parameter (nr__3=780) ! Number of records
character*80 ar__3 ! Full-size record
character*16 Name_3 (nr__3) ! Source name
character*63 Note_1 (nr__3) ! Text of the note (1)
*Note (1): The following abbreviations have been used:
* o, optical; x, X-ray; r, radio; IR, infrared.
* A reference on the source position is given under 'pos.'.
* The references code are from [lh]mxbrefs.dat files
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'hmxbrefs.dat' ! References (HMXB)
integer*4 nr__4
parameter (nr__4=814) ! Number of records
character*176 ar__4 ! Full-size record
integer*4 Ref (nr__4) ! Reference number
character*1 m_Ref (nr__4) ! [abc] Multiplicity index on Ref for
* associated reference
character*19 BibCode (nr__4) ! BibCode
character*37 Aut (nr__4) ! Author's name
character*113 Com (nr__4) ! Comments
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lmxbrefs.dat' ! References (LMXB)
integer*4 nr__5
parameter (nr__5=1472) ! Number of records
character*176 ar__5 ! Full-size record
integer*4 Ref_1 (nr__5) ! Reference number
character*1 m_Ref_1 (nr__5) ! [abc] Multiplicity index on Ref for
* associated reference
character*19 BibCode_1 (nr__5) ! BibCode
character*37 Aut_1 (nr__5) ! Author's name
character*113 Com_1 (nr__5) ! Comments
C=============================================================================
C Loading file 'hmxbcat.dat' ! High-mass X-ray binaries catalog
C Format for file interpretation
1 format(
+ A16,1X,A6,1X,I2,1X,I2,1X,F6.3,A1,A1,I2,1X,I2,1X,F5.2,A1,1X,
+ F6.2,1X,F8.4,1X,A4,1X,F7.2,1X,A17,2X,A15,2X,F6.3,1X,F5.2,A1,
+ 1X,F6.3,A1,1X,F4.2,1X,F5.2,1X,A1,F6.3,1X,A2,F4.2,A1,2X,A16,2X,
+ A1,F11.3,1X,F8.2,1X,A14,2X,A18,2X,F10.6,1X,F4.1,A1,1X,F11.5,
+ A1,2X,A15,2X,A25,1X,A12,1X,A13,A1,1X,A15)
C Effective file loading
open(unit=1,status='old',file=
+'hmxbcat.dat')
write(6,*) '....Loading file: hmxbcat.dat'
do i__=1,130
read(1,'(A359)')ar__
read(ar__,1)
+ Name(i__),Type(i__),RAh(i__),RAm(i__),RAs(i__),u_RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),u_DEs(i__),GLON(i__),
+ GLAT(i__),Pos(i__),e_Pos(i__),Opt(i__),r_Opt(i__),Vmag(i__),
+ Vmagl(i__),u_Vmag(i__),B_V(i__),u_B_V(i__),B_Vl(i__),U_B(i__),
+ l_E_B_V(i__),E_B_V(i__),l_E_B_V_2(i__),E_B_V_2(i__),
+ u_E_B_V(i__),r_Vmag(i__),l_Fx(i__),Fx(i__),Fxu(i__),
+ Range(i__),r_Fx(i__),Porb(i__),Porb2(i__),u_Porb(i__),
+ Ppulse(i__),u_Ppulse(i__),r_Ppulse(i__),Cat(i__),SpType(i__),
+ Name2(i__),u_Name2(i__),Name3(i__)
if(ar__(31:36) .EQ. '') RAs(i__) = rNULL__
if(ar__(45:49) .EQ. '') DEs(i__) = rNULL__
if(ar__(73:79) .EQ. '') e_Pos(i__) = rNULL__
if(ar__(117:122) .EQ. '') Vmag(i__) = rNULL__
if(ar__(124:128) .EQ. '') Vmagl(i__) = rNULL__
if(ar__(131:136) .EQ. '') B_V(i__) = rNULL__
if(ar__(139:142) .EQ. '') B_Vl(i__) = rNULL__
if(ar__(144:148) .EQ. '') U_B(i__) = rNULL__
if(ar__(151:156) .EQ. '') E_B_V(i__) = rNULL__
if(ar__(160:163) .EQ. '') E_B_V_2(i__) = rNULL__
if(ar__(186:196) .EQ. '') Fx(i__) = rNULL__
if(ar__(198:205) .EQ. '') Fxu(i__) = rNULL__
if(ar__(243:252) .EQ. '') Porb(i__) = rNULL__
if(ar__(254:257) .EQ. '') Porb2(i__) = rNULL__
if(ar__(260:270) .EQ. '') Ppulse(i__) = rNULL__
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__),Type(i__),RAh(i__),RAm(i__),RAs(i__),u_RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),u_DEs(i__),GLON(i__),
+ GLAT(i__),Pos(i__),e_Pos(i__),Opt(i__),r_Opt(i__),Vmag(i__),
+ Vmagl(i__),u_Vmag(i__),B_V(i__),u_B_V(i__),B_Vl(i__),U_B(i__),
+ l_E_B_V(i__),E_B_V(i__),l_E_B_V_2(i__),E_B_V_2(i__),
+ u_E_B_V(i__),r_Vmag(i__),l_Fx(i__),Fx(i__),Fxu(i__),
+ Range(i__),r_Fx(i__),Porb(i__),Porb2(i__),u_Porb(i__),
+ Ppulse(i__),u_Ppulse(i__),r_Ppulse(i__),Cat(i__),SpType(i__),
+ Name2(i__),u_Name2(i__),Name3(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 'lmxbcat.dat' ! Low-mass X-ray binaries catalog
C Format for file interpretation
2 format(
+ A16,1X,A6,1X,I2,1X,I2,1X,F6.3,A1,A1,I2,1X,I2,1X,F5.2,A1,1X,
+ F6.2,1X,F8.4,1X,A4,1X,F7.2,1X,A17,2X,A15,2X,F6.3,1X,F5.2,A1,
+ 1X,F6.3,A1,1X,F4.2,1X,F5.2,1X,A1,F6.3,1X,A2,F4.2,A1,2X,A16,2X,
+ A1,F11.3,1X,F8.2,1X,A14,2X,A18,2X,F10.6,1X,F4.1,A1,1X,F11.5,
+ A1,2X,A15,2X,A25,1X,A12,1X,A13,A1,1X,A15)
C Effective file loading
open(unit=1,status='old',file=
+'lmxbcat.dat')
write(6,*) '....Loading file: lmxbcat.dat'
do i__=1,150
read(1,'(A359)')ar__1
read(ar__1,2)
+ Name_1(i__),Type_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ u_RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
+ u_DEs_1(i__),GLON_1(i__),GLAT_1(i__),Pos_1(i__),e_Pos_1(i__),
+ Opt_1(i__),r_Opt_1(i__),Vmag_1(i__),Vmagl_1(i__),
+ u_Vmag_1(i__),B_V_1(i__),u_B_V_1(i__),B_Vl_1(i__),U_B_1(i__),
+ l_E_B_V_1(i__),E_B_V_1(i__),l_E_B_V_2_1(i__),E_B_V_2_1(i__),
+ u_E_B_V_1(i__),r_Vmag_1(i__),l_Fx_1(i__),Fx_1(i__),Fxu_1(i__),
+ Range_1(i__),r_Fx_1(i__),Porb_1(i__),Porb2_1(i__),
+ u_Porb_1(i__),Ppulse_1(i__),u_Ppulse_1(i__),r_Ppulse_1(i__),
+ Cat_1(i__),SpType_1(i__),Name2_1(i__),u_Name2_1(i__),
+ Name3_1(i__)
if(ar__1(31:36) .EQ. '') RAs_1(i__) = rNULL__
if(ar__1(45:49) .EQ. '') DEs_1(i__) = rNULL__
if(ar__1(73:79) .EQ. '') e_Pos_1(i__) = rNULL__
if(ar__1(117:122) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__1(124:128) .EQ. '') Vmagl_1(i__) = rNULL__
if(ar__1(131:136) .EQ. '') B_V_1(i__) = rNULL__
if(ar__1(139:142) .EQ. '') B_Vl_1(i__) = rNULL__
if(ar__1(144:148) .EQ. '') U_B_1(i__) = rNULL__
if(ar__1(151:156) .EQ. '') E_B_V_1(i__) = rNULL__
if(ar__1(160:163) .EQ. '') E_B_V_2_1(i__) = rNULL__
if(ar__1(186:196) .EQ. '') Fx_1(i__) = rNULL__
if(ar__1(198:205) .EQ. '') Fxu_1(i__) = rNULL__
if(ar__1(243:252) .EQ. '') Porb_1(i__) = rNULL__
if(ar__1(254:257) .EQ. '') Porb2_1(i__) = rNULL__
if(ar__1(260:270) .EQ. '') Ppulse_1(i__) = rNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Type_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ u_RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),
+ u_DEs_1(i__),GLON_1(i__),GLAT_1(i__),Pos_1(i__),e_Pos_1(i__),
+ Opt_1(i__),r_Opt_1(i__),Vmag_1(i__),Vmagl_1(i__),
+ u_Vmag_1(i__),B_V_1(i__),u_B_V_1(i__),B_Vl_1(i__),U_B_1(i__),
+ l_E_B_V_1(i__),E_B_V_1(i__),l_E_B_V_2_1(i__),E_B_V_2_1(i__),
+ u_E_B_V_1(i__),r_Vmag_1(i__),l_Fx_1(i__),Fx_1(i__),Fxu_1(i__),
+ Range_1(i__),r_Fx_1(i__),Porb_1(i__),Porb2_1(i__),
+ u_Porb_1(i__),Ppulse_1(i__),u_Ppulse_1(i__),r_Ppulse_1(i__),
+ Cat_1(i__),SpType_1(i__),Name2_1(i__),u_Name2_1(i__),
+ Name3_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hmxbnote.dat' ! Individual notes (HMXB)
C Format for file interpretation
3 format(A16,1X,A63)
C Effective file loading
open(unit=1,status='old',file=
+'hmxbnote.dat')
write(6,*) '....Loading file: hmxbnote.dat'
do i__=1,425
read(1,'(A80)')ar__2
read(ar__2,3)Name_2(i__),Note(i__)
c ..............Just test output...........
write(6,3)Name_2(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lmxbnote.dat' ! Individual notes (LMXB)
C Format for file interpretation
4 format(A16,1X,A63)
C Effective file loading
open(unit=1,status='old',file=
+'lmxbnote.dat')
write(6,*) '....Loading file: lmxbnote.dat'
do i__=1,780
read(1,'(A80)')ar__3
read(ar__3,4)Name_3(i__),Note_1(i__)
c ..............Just test output...........
write(6,4)Name_3(i__),Note_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'hmxbrefs.dat' ! References (HMXB)
C Format for file interpretation
5 format(I4,A1,1X,A19,1X,A37,A113)
C Effective file loading
open(unit=1,status='old',file=
+'hmxbrefs.dat')
write(6,*) '....Loading file: hmxbrefs.dat'
do i__=1,814
read(1,'(A176)')ar__4
read(ar__4,5)
+ Ref(i__),m_Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,5)
+ Ref(i__),m_Ref(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lmxbrefs.dat' ! References (LMXB)
C Format for file interpretation
6 format(I4,A1,1X,A19,1X,A37,A113)
C Effective file loading
open(unit=1,status='old',file=
+'lmxbrefs.dat')
write(6,*) '....Loading file: lmxbrefs.dat'
do i__=1,1472
read(1,'(A176)')ar__5
read(ar__5,6)
+ Ref_1(i__),m_Ref_1(i__),BibCode_1(i__),Aut_1(i__),Com_1(i__)
c ..............Just test output...........
write(6,6)
+ Ref_1(i__),m_Ref_1(i__),BibCode_1(i__),Aut_1(i__),Com_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end