Conversion of standardized ReadMe file for
file /./ftp/cats/J/AN/325/740 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. J/AN/325/740 Astrophysical supplements to ASCC-2.5 (Kharchenko+, 2004)
*================================================================================
*Astrophysical supplements to the ASCC-2.5.
*II. Membership probabilities in 520 Galactic open cluster sky areas.
* Kharchenko N.V., Piskunov A.E., Roeser S., Schilbach E., Scholz R.-D.
* <Astron. Nachr., 325, 740-748 (2004)>
* =2004AN....325..740K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'csoca.dat' ! The Catalog of Stars in the Open Cluster Areas
integer*4 nr__
parameter (nr__=171319) ! Number of records
character*227 ar__ ! Full-size record
C J2000.0 position composed of: RAhour DEdeg (Epoch=1991.25)
real*8 RAhour (nr__) ! (h) *Right Ascension J2000.0, Epoch=1991.25
real*8 DEdeg (nr__) ! (deg) Declination J2000.0, Epoch=1991.25
real*4 e_RAhour (nr__) ! (mas) Standard error in RA*cos(DE)
real*4 e_DEdeg (nr__) ! (mas) Standard error in DEdeg
real*8 Plx (nr__) ! (mas) ?=9999.99 Trigonometric parallax
real*8 e_Plx (nr__) ! (mas) ?=999.99 Standard error in Plx
real*8 pmRA (nr__) ! (mas/yr) Proper Motion in RA*cos(DEdeg)
real*8 pmDE (nr__) ! (mas/yr) Proper Motion in DE
real*8 e_pmRA (nr__) ! (mas/yr) Standard error in pmRA
real*8 e_pmDE (nr__) ! (mas/yr) Standard error in pmDE
real*8 Bmag (nr__) ! (mag) ?=99.999 B magnitude in Johnson system
real*8 Vmag (nr__) ! (mag) ?=99.999 V magnitude in Johnson system
real*8 e_Bmag (nr__) ! (mag) ?=9.999 Standard error on B magnitude
real*8 e_Vmag (nr__) ! (mag) ?=9.999 Standard error on V magnitude
real*8 Scat (nr__) ! (mag) ?=9.999 Scatter on magnitude
character*1 v1 (nr__) ! *[GN ] Known variability from GCVS/NSV
character*1 v2 (nr__) ! *[UVW ] Variability from Tycho-1
character*1 v3 (nr__) ! *[CDMPRU ] Variability type
character*1 v4 (nr__) ! *[VYIXR ] Variability from CMC11
character*2 d12 (nr__) ! *[A-S ] CCDM component identifier (Cat. <I/274>)
character*1 d3 (nr__) ! *[A-S ] Component identifier
character*1 d4 (nr__) ! *[DRSYZ ] Duplicity from Tycho-1
character*1 d5 (nr__) ! *[CGOVX ] Double/Multiple Systems flag
character*1 d6 (nr__) ! [D ] Duplicity flag from PPM
character*20 SpType (nr__) ! MK Spectral type (from <I/280> or <III/231>)
integer*4 HIP (nr__) ! ?=0 (I/239) Hipparcos number
integer*4 HD (nr__) ! ?=0 HD (III/135) number
integer*4 ASCC (nr__) ! ASCC-2.5 (I/280) number
integer*4 CRVAD (nr__) ! ?=0 CRVAD (III/239) number
real*8 RV (nr__) ! (km/s) ?=999.99 Average Radial Velocity
real*8 e_RV (nr__) ! (km/s) ?=0.00 Mean standard error in RV
integer*4 Seq (nr__) ! [1-520] Cluster sequential number, in
* cluster.dat file
real*8 rad (nr__) ! (deg) Distance from the cluster center
real*8 Pkin (nr__) ! [0,1] Membership probability from proper motion
real*8 Pph (nr__) ! [0,1] Membership probability from Photometry
integer*4 Psp (nr__) ! *[0,1] Membership probability from position
*Note on RAhour:
* RA = RAhour+24, if RAhour < 0 h
* RA = RAhour-24, if RAhour > 24 h
*Note on v1:
* Byte 292 from Tycho-1 file tyc_main.dat
* This flag is set when the variability is known, as follows:
* G = variable known in General Catalog of Variable Stars
* N = variable known in New Suspected Variables catalog
*Note on v2:
* Byte 294 from Tycho-1 file tyc_main.dat
* This flag has the following meaning:
* U = apparent variability in the Tycho data; may be due to duplicity
* V = strong evidence of intrinsic variability
* W = suspected intrinsic variability
*Note on v3:
* Byte 322 from Hipparcos file hip_main.dat
* Hipparcos-defined type of variability (a blank entry signifies that the
* entry could not be classified as variable or constant), as follows:
* C = no variability detected ("constant")
* D = duplicity-induced variability
* M = possibly micro-variable (amplitude < 0.03mag)
* P = periodic variable
* R = V-I colour index was revised due to variability analysis
* U = unsolved variable which does not fall in the other categories
*Note on v4:
* Flag from CMC11 catalogue
* V : variable star
* Y or I : if the CAMC did not measure the mag, an approximate value
* from other sources is given to a precision of 0.1 mag
* X or R : star fainter than 12 with unknown mag is given as 12.
*Note on d12:
* Bytes 149, 150 from Tycho-2 file catalog.dat.
* The CCDM component identifiers for double or multiple Hipparcos stars
* contributing to this Tycho-2 entry.
*Note on d3:
* Byte 41 from Hipparcos file hip_dm_com.dat
*Note on d4:
* Byte 296 from Tycho-1 file tyc_main.dat
* Unresolved duplicity status from Tycho data analysis, as follows:
* D = duplicity clearly indicated (B_Tmag and V_Tmag refer to combined light)
* R = duplicity weakly indicated, combined with indication of variability
* S = duplicity suspected
* Y = investigation for duplicity carried out on Tycho data, no indication of
* duplicity was found
* Z = investigation for duplicity not carried out
* blank = Hipparcos star not observed by Tycho
*Note on d5:
* Byte 347 from Hipparcos file hip_main.dat
* Indicates that further details are given in the Hipparcos Double and
* Multiple Systems Annex, as follows:
* C = solutions for the components
* G = acceleration or higher order terms
* O = orbital solutions
* V = variability-induced movers (apparent motion arises from variability)
* X = stochastic solution (probably astrometric binaries with short period)
*Note on Psp:
* Psp = 1 if rad <= cluster radius
* Psp = 0 if rad > cluster radius
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'cluster.dat' ! Cluster list
integer*4 nr__1
parameter (nr__1=520) ! Number of records
character*53 ar__1 ! Full-size record
C J2000.0 position composed of: RAhour DEdeg
integer*4 Seq_1 (nr__1) ! [1-520] Cluster sequential number
character*33 Name (nr__1) ! NGC, IC or other common designation
real*8 RAhour_1 (nr__1) ! (h) Right Ascension J2000.0
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000.0
C=============================================================================
C Loading file 'csoca.dat' ! The Catalog of Stars in the Open Cluster Areas
C Format for file interpretation
1 format(
+ F13.9,F13.8,F6.1,F6.1,F9.2,F9.2,F9.2,F9.2,F7.2,F7.2,F7.3,F7.3,
+ F7.3,F7.3,F7.3,1X,A1,A1,A1,A1,A2,A1,A1,A1,A1,1X,A20,I7,I7,I8,
+ I6,F9.2,F7.2,1X,I4,F7.4,F7.4,F7.4,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'csoca.dat')
write(6,*) '....Loading file: csoca.dat'
do i__=1,171319
read(1,'(A227)')ar__
read(ar__,1)
+ RAhour(i__),DEdeg(i__),e_RAhour(i__),e_DEdeg(i__),Plx(i__),
+ e_Plx(i__),pmRA(i__),pmDE(i__),e_pmRA(i__),e_pmDE(i__),
+ Bmag(i__),Vmag(i__),e_Bmag(i__),e_Vmag(i__),Scat(i__),v1(i__),
+ v2(i__),v3(i__),v4(i__),d12(i__),d3(i__),d4(i__),d5(i__),
+ d6(i__),SpType(i__),HIP(i__),HD(i__),ASCC(i__),CRVAD(i__),
+ RV(i__),e_RV(i__),Seq(i__),rad(i__),Pkin(i__),Pph(i__),
+ Psp(i__)
c ..............Just test output...........
write(6,1)
+ RAhour(i__),DEdeg(i__),e_RAhour(i__),e_DEdeg(i__),Plx(i__),
+ e_Plx(i__),pmRA(i__),pmDE(i__),e_pmRA(i__),e_pmDE(i__),
+ Bmag(i__),Vmag(i__),e_Bmag(i__),e_Vmag(i__),Scat(i__),v1(i__),
+ v2(i__),v3(i__),v4(i__),d12(i__),d3(i__),d4(i__),d5(i__),
+ d6(i__),SpType(i__),HIP(i__),HD(i__),ASCC(i__),CRVAD(i__),
+ RV(i__),e_RV(i__),Seq(i__),rad(i__),Pkin(i__),Pph(i__),
+ Psp(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'cluster.dat' ! Cluster list
C Format for file interpretation
2 format(I4,1X,A33,1X,F7.3,F7.2)
C Effective file loading
open(unit=1,status='old',file=
+'cluster.dat')
write(6,*) '....Loading file: cluster.dat'
do i__=1,520
read(1,'(A53)')ar__1
read(ar__1,2)Seq_1(i__),Name(i__),RAhour_1(i__),DEdeg_1(i__)
c ..............Just test output...........
write(6,2)Seq_1(i__),Name(i__),RAhour_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end