Conversion of standardized ReadMe file for
file /./ftp/cats/V/86 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-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/86 Interstellar matter in the Galactic Disk (Guarinos J., 1992)
*================================================================================
*Distribution of interstellar matter in the galactic disk from
*visual extinction data
* GUARINOS J.
* <Ph.D. Thesis, Strasbourg Observatory (1992)>
* <in "Astronomy from Large Databases II", Haguenau 14-16 September 1992,
* Ed. A. Heck and F. Murtagh, ESO Conference and Workshop Proceedings No 43,
* ISBN 3-923524-47-1, p. 301>
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ceddu1' ! Catalogue of Extinctions and Distances
Derived from UBV data (parameters)
integer*4 nr__
parameter (nr__=15449) ! Number of records
character*132 ar__ ! Full-size record
C Position composed of: GLON GLAT
character*12 Star (nr__) ! Identifier, either HD or BD/CD/CPD
character*20 MK (nr__) ! The MK spectral type written in
* a standard format
character*1 r_MK (nr__) ! [HJBSK] Source of MK (1)
real*4 Vmag (nr__) ! (mag) Apparent magnitude in UBV system (2)
real*4 B_V (nr__) ! (mag) Colour index in UBV system (2)
real*4 U_B (nr__) ! (mag) ? Colour index in UBV system (3)
real*4 absVmag (nr__) ! (mag) Absolute magnitude (section 3.3.2)
character*4 r_absVmag (nr__) ! Source of absVmag (3)
real*4 v_B_V_0 (nr__) ! (mag) Intrincic colour
character*4 r__B_V_0 (nr__) ! Source of (B-V)0 (4)
real*4 Av (nr__) ! (mag) Interstellar extinction (section 3.3.1)
integer*4 HDist (nr__) ! (pc) Heliocentric distance
integer*4 xDist (nr__) ! (pc) X-component of distance vector
* (Hdist*cos(b)*cos(l))
integer*4 yDist (nr__) ! (pc) Y-component of distance vector
* (Hdist*cos(b)*sin(l))
integer*4 zDist (nr__) ! (pc) Z-component of distance vector
* (Hdist*sin(b))
real*4 GLON (nr__) ! (deg) [0/360[+= Galactic longitude
real*4 GLAT (nr__) ! (deg) Galactic latitude
*Note (1): MK Type source takes the following values:
* B: Third General Catalogue of MK Spectral Classification (Buscombe 1977).
* H: Michigan Spectral Catalogue
* (Houk & Cowley 1975; Houk 1978; Houk 1982; Houk & Smith-Moore 1988).
* J: Catalogue of Selected Spectral Types in the MK System (M. Jaschek 1978).
* K: Catalogue of Jensen (1983).
* S: from SIMBAD for a very small number of stars.
*Note (2): The value is the average of all the non-deviating measurements
* provided by the Photoelectric Photometric Catalogue (Mermilliod 1987)
*Note (3): Source of the table which provided the absolute magnitudes:
* CGMK: Corbally and Garrison (1984)
* GRMK: Grenier et al. (1985), table for magnitude-limited samples
* SKMK: Schmidt-Kaler (1982) MK/Mv table
* SKMQ: Schmidt-Kaler (1982) (B-V)o+Lum.Class/Mv table using the
* photometric method exposed in section 3.3.1 (see
* Garcia, Claria and Levato 1988)
*Note (4): Method of determination of the intrinsic colour:
* MQ: photometric method exposed in section 3.3.1 (see
* Garcia, Claria and Levato 1988)
* SK: from the MK/(B-V)o table from Schmidt-Kaler (1982)
* SK-n: (n either 1 or 2) indicates that a correction of n sub-types
* towards the blue has been performed on the spectral type,
* before using the MK/(B-V)o table from Schmidt-Kaler (1982).
* MQ-1: a correction of 1 sub-type towards the blue has been performed
* on the spectral type, before using the MQ method
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'ceddu2' ! Catalogue of Extinctions and Distances
Derived from UBV data (1)
integer*4 nr__1
parameter (nr__1=15449) ! Number of records
character*132 ar__1 ! Full-size record
C B2000 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg (nr__1) ! (deg) Right Ascension B2000
real*8 DEdeg (nr__1) ! (deg) Declination B2000
C ---------------------------------- ! (position vector(s) in degrees)
character*12 Star_1 (nr__1) ! Identifier, either HD or BD/CD/CPD
character*32 MK_1 (nr__1) ! The full MK spectral type including all
* spectrum peculiarities written in
* a standard format
character*1 r_MK_1 (nr__1) ! [HJBSK] Source of MK (see table ceddu1)
integer*4 RAh (nr__1) ! (h) Right Ascension B2000 (hours)
integer*4 RAm (nr__1) ! (min) Right Ascension B2000 (minutes)
integer*4 RAs (nr__1) ! (s) Right Ascension B2000 (seconds)
character*1 n_RAs (nr__1) ! [*] An '*' exists if RA was originally
* expressed in tenths of minutes
character*1 DE_ (nr__1) ! Declination B2000 (sign)
integer*4 DEd (nr__1) ! (deg) Declination B2000 (degrees)
real*4 DEm (nr__1) ! (arcmin) Declination B2000 (minutes)
real*4 GLON_1 (nr__1) ! (deg) [0/360[+= Galactic longitude
real*4 GLAT_1 (nr__1) ! (deg) Galactic latitude
real*4 Vmag_1 (nr__1) ! (mag) Apparent V magnitude, as in table ceddu1
integer*4 o_Vmag (nr__1) ! Number of measurements used to derive Vmag
real*4 e_Vmag (nr__1) ! (mag) Mean error (rms) of Vmag
real*4 v_B_V (nr__1) ! (mag) Colour index, as in table ceddu1
integer*4 o__B_V (nr__1) ! Number of measurements used to derive (B-V)
real*4 e__B_V (nr__1) ! (mag) Mean error (rms) of (B-V)
real*4 v_U_B (nr__1) ! (mag) ? Colour index (3)
integer*4 o__U_B (nr__1) ! ? Number of observations used to derive (U-B)
real*4 e__U_B (nr__1) ! (mag) ? Mean error (rms) of (U-B)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'asmal1' ! Stars with low reddening Small extinction (1)
integer*4 nr__2
parameter (nr__2=98) ! Number of records
character*115 ar__2 ! Full-size record
character*12 Star_2 (nr__2) ! Identifier, either HD or BD/CD/CPD
character*20 MK_2 (nr__2) ! The MK spectral type written in
* a standard format
character*1 r_MK_2 (nr__2) ! [HJBSK] Source of MK, as in ceddu1 table
real*4 Vmag_2 (nr__2) ! (mag) Apparent magnitude in UBV system
real*4 B_V_1 (nr__2) ! (mag) Colour index in UBV system
real*4 U_B_1 (nr__2) ! (mag) ? Colour index in UBV system
real*4 absVmag_1 (nr__2) ! (mag) ? Absolute magnitude
character*4 r_absVmag_1(nr__2) ! Source of absVmag
real*4 v_B_V_0_1 (nr__2) ! (mag) Intrincic colour
character*4 r__B_V_0_1 (nr__2) ! Source of (B-V)0
real*4 E_B_V (nr__2) ! (mag) Colour excess
integer*4 HDist_1 (nr__2) ! (pc) ? Heliocentric distance
integer*4 xDist_1 (nr__2) ! (pc) ? X-component of distance vector
* (Hdist*cos(b)*cos(l))
integer*4 yDist_1 (nr__2) ! (pc) ? Y-component of distance vector
* (Hdist*cos(b)*sin(l))
integer*4 zDist_1 (nr__2) ! (pc) ? Z-component of distance vector
* (Hdist*sin(b))
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'asmal2' ! Stars with low reddening Small extinction (2)
integer*4 nr__3
parameter (nr__3=98) ! Number of records
character*132 ar__3 ! Full-size record
C B2000 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg_1 (nr__3) ! (deg) Right Ascension B2000
real*8 DEdeg_1 (nr__3) ! (deg) Declination B2000
C ---------------------------------- ! (position vector(s) in degrees)
character*12 Star_3 (nr__3) ! Identifier, either HD or BD/CD/CPD
character*32 MK_3 (nr__3) ! The full MK spectral type including all
* spectrum peculiarities written in
* a standard format
character*1 r_MK_3 (nr__3) ! [HJBSK] Source of MK (see table ceddu1)
integer*4 RAh_1 (nr__3) ! (h) Right Ascension B2000 (hours)
integer*4 RAm_1 (nr__3) ! (min) Right Ascension B2000 (minutes)
integer*4 RAs_1 (nr__3) ! (s) Right Ascension B2000 (seconds)
character*1 n_RAs_1 (nr__3) ! [*] An '*' exists if RA was originally
* expressed in tenths of minutes
character*1 DE__1 (nr__3) ! Declination B2000 (sign)
integer*4 DEd_1 (nr__3) ! (deg) Declination B2000 (degrees)
real*4 DEm_1 (nr__3) ! (arcmin) Declination B2000 (minutes)
real*4 GLON_2 (nr__3) ! (deg) [0/360[+= Galactic longitude
real*4 GLAT_2 (nr__3) ! (deg) Galactic latitude
real*4 Vmag_3 (nr__3) ! (mag) Apparent V magnitude, as in table ceddu1
integer*4 o_Vmag_1 (nr__3) ! Number of measurements used to derive Vmag
real*4 e_Vmag_1 (nr__3) ! (mag) Mean error (rms) of Vmag
real*4 v_B_V_1 (nr__3) ! (mag) Colour index, as in table ceddu1
integer*4 o__B_V_1 (nr__3) ! Number of measurements used to derive (B-V)
real*4 e__B_V_1 (nr__3) ! (mag) Mean error (rms) of (B-V)
real*4 v_U_B_1 (nr__3) ! (mag) ? Colour index (3)
integer*4 o__U_B_1 (nr__3) ! ? Number of observations used to derive (U-B)
real*4 e__U_B_1 (nr__3) ! (mag) ? Mean error (rms) of (U-B)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'avlar1' ! Stars with very Large Extinction (1)
integer*4 nr__4
parameter (nr__4=62) ! Number of records
character*115 ar__4 ! Full-size record
character*12 Star_4 (nr__4) ! Identifier, either HD or BD/CD/CPD
character*20 MK_4 (nr__4) ! The MK spectral type written in
* a standard format
character*1 r_MK_4 (nr__4) ! [HJBSK] Source of MK, as in ceddu1 table
real*4 Vmag_4 (nr__4) ! (mag) Apparent magnitude in UBV system
real*4 B_V_2 (nr__4) ! (mag) Colour index in UBV system
real*4 U_B_2 (nr__4) ! (mag) ? Colour index in UBV system
real*4 absVmag_2 (nr__4) ! (mag) ? Absolute magnitude
character*4 r_absVmag_2(nr__4) ! Source of absVmag
real*4 v_B_V_0_2 (nr__4) ! (mag) Intrincic colour
character*4 r__B_V_0_2 (nr__4) ! Source of (B-V)0
real*4 E_B_V_1 (nr__4) ! (mag) Colour excess
integer*4 HDist_2 (nr__4) ! (pc) ? Heliocentric distance
integer*4 xDist_2 (nr__4) ! (pc) ? X-component of distance vector
* (Hdist*cos(b)*cos(l))
integer*4 yDist_2 (nr__4) ! (pc) ? Y-component of distance vector
* (Hdist*cos(b)*sin(l))
integer*4 zDist_2 (nr__4) ! (pc) ? Z-component of distance vector
* (Hdist*sin(b))
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'avlar2' ! Stars with very Large Extinction (2)
integer*4 nr__5
parameter (nr__5=62) ! Number of records
character*132 ar__5 ! Full-size record
C B2000 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg_2 (nr__5) ! (deg) Right Ascension B2000
real*8 DEdeg_2 (nr__5) ! (deg) Declination B2000
C ---------------------------------- ! (position vector(s) in degrees)
character*12 Star_5 (nr__5) ! Identifier, either HD or BD/CD/CPD
character*32 MK_5 (nr__5) ! The full MK spectral type including all
* spectrum peculiarities written in
* a standard format
character*1 r_MK_5 (nr__5) ! [HJBSK] Source of MK (see table ceddu1)
integer*4 RAh_2 (nr__5) ! (h) Right Ascension B2000 (hours)
integer*4 RAm_2 (nr__5) ! (min) Right Ascension B2000 (minutes)
integer*4 RAs_2 (nr__5) ! (s) Right Ascension B2000 (seconds)
character*1 n_RAs_2 (nr__5) ! [*] An '*' exists if RA was originally
* expressed in tenths of minutes
character*1 DE__2 (nr__5) ! Declination B2000 (sign)
integer*4 DEd_2 (nr__5) ! (deg) Declination B2000 (degrees)
real*4 DEm_2 (nr__5) ! (arcmin) Declination B2000 (minutes)
real*4 GLON_3 (nr__5) ! (deg) [0/360[+= Galactic longitude
real*4 GLAT_3 (nr__5) ! (deg) Galactic latitude
real*4 Vmag_5 (nr__5) ! (mag) Apparent V magnitude, as in table ceddu1
integer*4 o_Vmag_2 (nr__5) ! Number of measurements used to derive Vmag
real*4 e_Vmag_2 (nr__5) ! (mag) Mean error (rms) of Vmag
real*4 v_B_V_1 (nr__5) ! (mag) Colour index, as in table ceddu1
integer*4 o__B_V_1 (nr__5) ! Number of measurements used to derive (B-V)
real*4 e__B_V_1 (nr__5) ! (mag) Mean error (rms) of (B-V)
real*4 v_U_B_1 (nr__5) ! (mag) ? Colour index (3)
integer*4 o__U_B_1 (nr__5) ! ? Number of observations used to derive (U-B)
real*4 e__U_B_1 (nr__5) ! (mag) ? Mean error (rms) of (U-B)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'exneg1' ! Stars with negative extinctions (1)
integer*4 nr__6
parameter (nr__6=390) ! Number of records
character*131 ar__6 ! Full-size record
character*12 Star_6 (nr__6) ! Identifier, either HD or BD/CD/CPD
character*20 MK_6 (nr__6) ! The MK spectral type written in
* a standard format
character*1 r_MK_6 (nr__6) ! [HJBSK] Source of MK, as in ceddu1 table
real*4 Vmag_6 (nr__6) ! (mag) Apparent magnitude in UBV system
real*4 B_V_3 (nr__6) ! (mag) Colour index in UBV system
real*4 U_B_3 (nr__6) ! (mag) ? Colour index in UBV system
real*4 absVmag_3 (nr__6) ! (mag) ? Absolute magnitude
character*4 r_absVmag_3(nr__6) ! Source of absVmag
real*4 v_B_V_0_3 (nr__6) ! (mag) Intrincic colour
character*4 r__B_V_0_3 (nr__6) ! Source of (B-V)0
real*4 E_B_V_1 (nr__6) ! (mag) Colour excess
integer*4 HDist_3 (nr__6) ! (pc) ? Heliocentric distance
integer*4 xDist_3 (nr__6) ! (pc) ? X-component of distance vector
* (Hdist*cos(b)*cos(l))
integer*4 yDist_3 (nr__6) ! (pc) ? Y-component of distance vector
* (Hdist*cos(b)*sin(l))
integer*4 zDist_3 (nr__6) ! (pc) ? Z-component of distance vector
* (Hdist*sin(b))
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'exneg2' ! Stars with negative extinctions (2)
integer*4 nr__7
parameter (nr__7=390) ! Number of records
character*132 ar__7 ! Full-size record
C B2000 position composed of: RAh RAm RAs DE- DEd DEm
real*8 RAdeg_3 (nr__7) ! (deg) Right Ascension B2000
real*8 DEdeg_3 (nr__7) ! (deg) Declination B2000
C ---------------------------------- ! (position vector(s) in degrees)
character*12 Star_7 (nr__7) ! Identifier, either HD or BD/CD/CPD
character*32 MK_7 (nr__7) ! The full MK spectral type including all
* spectrum peculiarities written in
* a standard format
character*1 r_MK_7 (nr__7) ! [HJBSK] Source of MK (see table ceddu1)
integer*4 RAh_3 (nr__7) ! (h) Right Ascension B2000 (hours)
integer*4 RAm_3 (nr__7) ! (min) Right Ascension B2000 (minutes)
integer*4 RAs_3 (nr__7) ! (s) Right Ascension B2000 (seconds)
character*1 n_RAs_3 (nr__7) ! [*] An '*' exists if RA was originally
* expressed in tenths of minutes
character*1 DE__3 (nr__7) ! Declination B2000 (sign)
integer*4 DEd_3 (nr__7) ! (deg) Declination B2000 (degrees)
real*4 DEm_3 (nr__7) ! (arcmin) Declination B2000 (minutes)
real*4 GLON_4 (nr__7) ! (deg) [0/360[+= Galactic longitude
real*4 GLAT_4 (nr__7) ! (deg) Galactic latitude
real*4 Vmag_7 (nr__7) ! (mag) Apparent V magnitude, as in table ceddu1
integer*4 o_Vmag_3 (nr__7) ! Number of measurements used to derive Vmag
real*4 e_Vmag_3 (nr__7) ! (mag) Mean error (rms) of Vmag
real*4 v_B_V_1 (nr__7) ! (mag) Colour index, as in table ceddu1
integer*4 o__B_V_1 (nr__7) ! Number of measurements used to derive (B-V)
real*4 e__B_V_1 (nr__7) ! (mag) Mean error (rms) of (B-V)
real*4 v_U_B_1 (nr__7) ! (mag) ? Colour index (3)
integer*4 o__U_B_1 (nr__7) ! ? Number of observations used to derive (U-B)
real*4 e__U_B_1 (nr__7) ! (mag) ? Mean error (rms) of (U-B)
C=============================================================================
C Loading file 'ceddu1' ! Catalogue of Extinctions and Distances
* Derived from UBV data (parameters)
C Format for file interpretation
1 format(
+ A12,1X,A20,2X,A1,2X,F6.3,1X,F6.3,1X,F6.3,3X,F5.2,1X,A4,3X,
+ F5.2,1X,A4,2X,F5.2,2X,I4,1X,I5,1X,I5,1X,I5,3X,F6.2,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'ceddu1')
write(6,*) '....Loading file: ceddu1'
do i__=1,15449
read(1,'(A132)')ar__
read(ar__,1)
+ Star(i__),MK(i__),r_MK(i__),Vmag(i__),B_V(i__),U_B(i__),
+ absVmag(i__),r_absVmag(i__),v_B_V_0(i__),r__B_V_0(i__),
+ Av(i__),HDist(i__),xDist(i__),yDist(i__),zDist(i__),GLON(i__),
+ GLAT(i__)
if(ar__(53:58) .EQ. '') U_B(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Star(i__),MK(i__),r_MK(i__),Vmag(i__),B_V(i__),U_B(i__),
+ absVmag(i__),r_absVmag(i__),v_B_V_0(i__),r__B_V_0(i__),
+ Av(i__),HDist(i__),xDist(i__),yDist(i__),zDist(i__),GLON(i__),
+ GLAT(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'ceddu2' ! Catalogue of Extinctions and Distances
* Derived from UBV data (1)
C Format for file interpretation
2 format(
+ A12,1X,A32,A1,2X,I2,1X,I2,1X,I2,A1,1X,A1,I2,1X,F4.1,2X,F6.2,
+ 1X,F6.2,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,
+ F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'ceddu2')
write(6,*) '....Loading file: ceddu2'
do i__=1,15449
read(1,'(A132)')ar__1
read(ar__1,2)
+ Star_1(i__),MK_1(i__),r_MK_1(i__),RAh(i__),RAm(i__),RAs(i__),
+ n_RAs(i__),DE_(i__),DEd(i__),DEm(i__),GLON_1(i__),GLAT_1(i__),
+ Vmag_1(i__),o_Vmag(i__),e_Vmag(i__),v_B_V(i__),o__B_V(i__),
+ e__B_V(i__),v_U_B(i__),o__U_B(i__),e__U_B(i__)
if(ar__1(117:122) .EQ. '') v_U_B(i__) = rNULL__
if(ar__1(123:125) .EQ. '') o__U_B(i__) = iNULL__
if(ar__1(127:132) .EQ. '') e__U_B(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(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ Star_1(i__),MK_1(i__),r_MK_1(i__),RAh(i__),RAm(i__),RAs(i__),
+ n_RAs(i__),DE_(i__),DEd(i__),DEm(i__),GLON_1(i__),GLAT_1(i__),
+ Vmag_1(i__),o_Vmag(i__),e_Vmag(i__),v_B_V(i__),o__B_V(i__),
+ e__B_V(i__),v_U_B(i__),o__U_B(i__),e__U_B(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 'asmal1' ! Stars with low reddening Small extinction (1)
C Format for file interpretation
3 format(
+ A12,1X,A20,2X,A1,2X,F6.3,1X,F6.3,1X,F6.3,3X,F5.2,1X,A4,3X,
+ F5.2,1X,A4,2X,F5.2,1X,I5,I6,I6,I6)
C Effective file loading
open(unit=1,status='old',file=
+'asmal1')
write(6,*) '....Loading file: asmal1'
do i__=1,98
read(1,'(A115)')ar__2
read(ar__2,3)
+ Star_2(i__),MK_2(i__),r_MK_2(i__),Vmag_2(i__),B_V_1(i__),
+ U_B_1(i__),absVmag_1(i__),r_absVmag_1(i__),v_B_V_0_1(i__),
+ r__B_V_0_1(i__),E_B_V(i__),HDist_1(i__),xDist_1(i__),
+ yDist_1(i__),zDist_1(i__)
if(ar__2(53:58) .EQ. '') U_B_1(i__) = rNULL__
if(ar__2(62:66) .EQ. '') absVmag_1(i__) = rNULL__
if(ar__2(93:97) .EQ. '') HDist_1(i__) = iNULL__
if(ar__2(98:103) .EQ. '') xDist_1(i__) = iNULL__
if(ar__2(104:109) .EQ. '') yDist_1(i__) = iNULL__
if(ar__2(110:115) .EQ. '') zDist_1(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ Star_2(i__),MK_2(i__),r_MK_2(i__),Vmag_2(i__),B_V_1(i__),
+ U_B_1(i__),absVmag_1(i__),r_absVmag_1(i__),v_B_V_0_1(i__),
+ r__B_V_0_1(i__),E_B_V(i__),HDist_1(i__),xDist_1(i__),
+ yDist_1(i__),zDist_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'asmal2' ! Stars with low reddening Small extinction (2)
C Format for file interpretation
4 format(
+ A12,1X,A32,A1,2X,I2,1X,I2,1X,I2,A1,1X,A1,I2,1X,F4.1,2X,F6.2,
+ 1X,F6.2,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,
+ F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'asmal2')
write(6,*) '....Loading file: asmal2'
do i__=1,98
read(1,'(A132)')ar__3
read(ar__3,4)
+ Star_3(i__),MK_3(i__),r_MK_3(i__),RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),n_RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ GLON_2(i__),GLAT_2(i__),Vmag_3(i__),o_Vmag_1(i__),
+ e_Vmag_1(i__),v_B_V_1(i__),o__B_V_1(i__),e__B_V_1(i__),
+ v_U_B_1(i__),o__U_B_1(i__),e__U_B_1(i__)
if(ar__3(117:122) .EQ. '') v_U_B_1(i__) = rNULL__
if(ar__3(123:125) .EQ. '') o__U_B_1(i__) = iNULL__
if(ar__3(127:132) .EQ. '') e__U_B_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(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,4)
+ Star_3(i__),MK_3(i__),r_MK_3(i__),RAh_1(i__),RAm_1(i__),
+ RAs_1(i__),n_RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ GLON_2(i__),GLAT_2(i__),Vmag_3(i__),o_Vmag_1(i__),
+ e_Vmag_1(i__),v_B_V_1(i__),o__B_V_1(i__),e__B_V_1(i__),
+ v_U_B_1(i__),o__U_B_1(i__),e__U_B_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 'avlar1' ! Stars with very Large Extinction (1)
C Format for file interpretation
5 format(
+ A12,1X,A20,2X,A1,2X,F6.3,1X,F6.3,1X,F6.3,3X,F5.2,1X,A4,3X,
+ F5.2,1X,A4,2X,F5.2,1X,I5,I6,I6,I6)
C Effective file loading
open(unit=1,status='old',file=
+'avlar1')
write(6,*) '....Loading file: avlar1'
do i__=1,62
read(1,'(A115)')ar__4
read(ar__4,5)
+ Star_4(i__),MK_4(i__),r_MK_4(i__),Vmag_4(i__),B_V_2(i__),
+ U_B_2(i__),absVmag_2(i__),r_absVmag_2(i__),v_B_V_0_2(i__),
+ r__B_V_0_2(i__),E_B_V_1(i__),HDist_2(i__),xDist_2(i__),
+ yDist_2(i__),zDist_2(i__)
if(ar__4(53:58) .EQ. '') U_B_2(i__) = rNULL__
if(ar__4(62:66) .EQ. '') absVmag_2(i__) = rNULL__
if(ar__4(93:97) .EQ. '') HDist_2(i__) = iNULL__
if(ar__4(98:103) .EQ. '') xDist_2(i__) = iNULL__
if(ar__4(104:109) .EQ. '') yDist_2(i__) = iNULL__
if(ar__4(110:115) .EQ. '') zDist_2(i__) = iNULL__
c ..............Just test output...........
write(6,5)
+ Star_4(i__),MK_4(i__),r_MK_4(i__),Vmag_4(i__),B_V_2(i__),
+ U_B_2(i__),absVmag_2(i__),r_absVmag_2(i__),v_B_V_0_2(i__),
+ r__B_V_0_2(i__),E_B_V_1(i__),HDist_2(i__),xDist_2(i__),
+ yDist_2(i__),zDist_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'avlar2' ! Stars with very Large Extinction (2)
C Format for file interpretation
6 format(
+ A12,1X,A32,A1,2X,I2,1X,I2,1X,I2,A1,1X,A1,I2,1X,F4.1,2X,F6.2,
+ 1X,F6.2,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,
+ F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'avlar2')
write(6,*) '....Loading file: avlar2'
do i__=1,62
read(1,'(A132)')ar__5
read(ar__5,6)
+ Star_5(i__),MK_5(i__),r_MK_5(i__),RAh_2(i__),RAm_2(i__),
+ RAs_2(i__),n_RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),
+ GLON_3(i__),GLAT_3(i__),Vmag_5(i__),o_Vmag_2(i__),
+ e_Vmag_2(i__),v_B_V_1(i__),o__B_V_1(i__),e__B_V_1(i__),
+ v_U_B_1(i__),o__U_B_1(i__),e__U_B_1(i__)
if(ar__5(117:122) .EQ. '') v_U_B_1(i__) = rNULL__
if(ar__5(123:125) .EQ. '') o__U_B_1(i__) = iNULL__
if(ar__5(127:132) .EQ. '') e__U_B_1(i__) = rNULL__
RAdeg_2(i__) = rNULL__
DEdeg_2(i__) = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c ..............Just test output...........
write(6,6)
+ Star_5(i__),MK_5(i__),r_MK_5(i__),RAh_2(i__),RAm_2(i__),
+ RAs_2(i__),n_RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),
+ GLON_3(i__),GLAT_3(i__),Vmag_5(i__),o_Vmag_2(i__),
+ e_Vmag_2(i__),v_B_V_1(i__),o__B_V_1(i__),e__B_V_1(i__),
+ v_U_B_1(i__),o__U_B_1(i__),e__U_B_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'exneg1' ! Stars with negative extinctions (1)
C Format for file interpretation
7 format(
+ A12,1X,A20,2X,A1,2X,F6.3,1X,F6.3,1X,F6.3,3X,F5.2,1X,A4,3X,
+ F5.2,1X,A4,2X,F5.2,1X,I5,I6,I6,I6)
C Effective file loading
open(unit=1,status='old',file=
+'exneg1')
write(6,*) '....Loading file: exneg1'
do i__=1,390
read(1,'(A131)')ar__6
read(ar__6,7)
+ Star_6(i__),MK_6(i__),r_MK_6(i__),Vmag_6(i__),B_V_3(i__),
+ U_B_3(i__),absVmag_3(i__),r_absVmag_3(i__),v_B_V_0_3(i__),
+ r__B_V_0_3(i__),E_B_V_1(i__),HDist_3(i__),xDist_3(i__),
+ yDist_3(i__),zDist_3(i__)
if(ar__6(53:58) .EQ. '') U_B_3(i__) = rNULL__
if(ar__6(62:66) .EQ. '') absVmag_3(i__) = rNULL__
if(ar__6(93:97) .EQ. '') HDist_3(i__) = iNULL__
if(ar__6(98:103) .EQ. '') xDist_3(i__) = iNULL__
if(ar__6(104:109) .EQ. '') yDist_3(i__) = iNULL__
if(ar__6(110:115) .EQ. '') zDist_3(i__) = iNULL__
c ..............Just test output...........
write(6,7)
+ Star_6(i__),MK_6(i__),r_MK_6(i__),Vmag_6(i__),B_V_3(i__),
+ U_B_3(i__),absVmag_3(i__),r_absVmag_3(i__),v_B_V_0_3(i__),
+ r__B_V_0_3(i__),E_B_V_1(i__),HDist_3(i__),xDist_3(i__),
+ yDist_3(i__),zDist_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'exneg2' ! Stars with negative extinctions (2)
C Format for file interpretation
8 format(
+ A12,1X,A32,A1,2X,I2,1X,I2,1X,I2,A1,1X,A1,I2,1X,F4.1,2X,F6.2,
+ 1X,F6.2,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,F6.3,1X,F6.3,I3,1X,
+ F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'exneg2')
write(6,*) '....Loading file: exneg2'
do i__=1,390
read(1,'(A132)')ar__7
read(ar__7,8)
+ Star_7(i__),MK_7(i__),r_MK_7(i__),RAh_3(i__),RAm_3(i__),
+ RAs_3(i__),n_RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),
+ GLON_4(i__),GLAT_4(i__),Vmag_7(i__),o_Vmag_3(i__),
+ e_Vmag_3(i__),v_B_V_1(i__),o__B_V_1(i__),e__B_V_1(i__),
+ v_U_B_1(i__),o__U_B_1(i__),e__U_B_1(i__)
if(ar__7(117:122) .EQ. '') v_U_B_1(i__) = rNULL__
if(ar__7(123:125) .EQ. '') o__U_B_1(i__) = iNULL__
if(ar__7(127:132) .EQ. '') e__U_B_1(i__) = rNULL__
RAdeg_3(i__) = rNULL__
DEdeg_3(i__) = rNULL__
c Derive coordinates RAdeg_3 and DEdeg_3 from input data
c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
if(RAh_3(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c ..............Just test output...........
write(6,8)
+ Star_7(i__),MK_7(i__),r_MK_7(i__),RAh_3(i__),RAm_3(i__),
+ RAs_3(i__),n_RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),
+ GLON_4(i__),GLAT_4(i__),Vmag_7(i__),o_Vmag_3(i__),
+ e_Vmag_3(i__),v_B_V_1(i__),o__B_V_1(i__),e__B_V_1(i__),
+ v_U_B_1(i__),o__U_B_1(i__),e__U_B_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end