Conversion of standardized ReadMe file for
file /./ftp/cats/B/gcvs 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.5, on 2013-Jun-19
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. B/gcvs General Catalogue of Variable Stars (Samus+ 2007-2013)
*================================================================================
*General Catalog of Variable Stars (GCVS database, Version 2012Apr)
* Samus N.N., Durlevich O.V., et al.
* <Institute of Astronomy of Russian Academy of Sciences and Sternberg
* State Astronomical Institute of the Moscow State University >
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcvs_cat.dat' ! The GCVS Catalog (Vol. I-III, version 2013-04-30)
integer*4 nr__
parameter (nr__=47969) ! Number of records
character*174 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension J2000.0
real*8 Dec (nr__) ! (deg) Declination J2000.0
integer*4 VarNum (nr__) ! Numeric designation, made from constellation
* and star names (codes in file "constel.txt")
character*1 m_VarNum (nr__) ! Component identification (1)
character*10 GCVS (nr__) ! Variable star designation (G1)
character*1 n_GCVS (nr__) ! [*] indicates a remark in "gcvs_rem.dat" file
integer*4 RAh (nr__) ! (h) ? Hours RA, equinox J2000.0 (2)
integer*4 RAm (nr__) ! (min) ? Minutes RA, equinox J2000.0 (2)
real*4 RAs (nr__) ! (s) ? Seconds RA, equinox J2000.0 (2)
character*1 DE_ (nr__) ! ? Sign Dec, equinox J2000.0 (2)
integer*4 DEd (nr__) ! (deg) ? Degrees Dec, equinox J2000.0 (2)
integer*4 DEm (nr__) ! (arcmin) ? Minutes Dec, equinox J2000.0 (2)
integer*4 DEs (nr__) ! (arcsec) ? Seconds Dec, equinox J2000.0 (2)
character*1 u_DEs (nr__) ! [:*] position accuracy flags (3)
character*10 VarType (nr__) ! Type of variability (see file "vartype.txt")
character*1 l_magMax (nr__) ! [<>(] Limit or amplitude symbol on magMax (G3)
real*4 magMax (nr__) ! (mag) ? Magnitude at maximum brightness
character*1 u_magMax (nr__) ! Uncertainty flag (:) on magMax
character*1 l_Min1 (nr__) ! [<(] Limit or amplitude symbol on Min1 (G3)
real*4 Min1 (nr__) ! (mag) ? Minimum magnitude or amplitude
character*1 u_Min1 (nr__) ! Uncertainty flag (:) on Min1
character*2 n_Min1 (nr__) ! Alternative photometric system for Min1 (G4)
character*1 l_Min2 (nr__) ! [<(] Limit or amplitude symbol on Min2 (G3)
real*4 Min2 (nr__) ! (mag) ? Secondary minimum magnitude or amplitude
character*1 u_Min2 (nr__) ! Uncertainty flag (:) on Min2
character*2 n_Min2 (nr__) ! Alternative photometric system for Min2 (G4)
character*2 flt (nr__) ! The photometric system for magnitudes (G4)
real*8 Epoch (nr__) ! (d) ? Epoch for maximum light, Julian days (G5)
character*1 u_Epoch (nr__) ! [:+-] Uncertainty flag on Epoch (4)
character*4 Year (nr__) ! Year of outburst for nova or supernova
character*1 u_Year (nr__) ! [:] Uncertainty flag on Year of outburst
character*1 l_Period (nr__) ! [<>(] Code for upper or lower limits (5)
real*8 Period (nr__) ! (d) ? Period of the variable star
character*3 u_Period (nr__) ! [*/N)2: ] Uncertainties on Period (6)
integer*4 M_m_D (nr__) ! (%) ? Rising time or duration of eclipse (G6)
character*1 u_M_m_D (nr__) ! Uncertainty flag (:) on M-m/D
character*1 n_M_m_D (nr__) ! [*] Note for eclipsing variable (G6)
character*17 SpType (nr__) ! MK Spectral type
character*5 Ref1 (nr__) ! Reference to a study of the star (G7)
character*5 Ref2 (nr__) ! Reference to a chart or photograph (G7)
character*1 f_GCVS (nr__) ! [=N+] "N" if the star does not exist (7)
character*10 VarName (nr__) ! Alternative name of the variable (G8)
*Note (1): designations of components of double/multiple stars :
* A,B; 1,2,3 ... (tau 1). Usually only one component of the binary
* is variable, but there exist two cases (CE Cas A,B; EQ Peg A,B) when
* both components are variable and enter the catalogue as separate lines.
*Note (2):
* the position is not listed when the variable is equivalent to another
* catalogued variable (see column VarName)
*Note (3):
* the asterisk (*) means that the coordinates that could not be improved
* were recalculated from old rough coordinates to the equinox J2000.0
* with allowance for precession.
*Note (4): the following codes are used:
* ":" are used if the epoch is uncertain;
* "+" epoch may be later than that reported;
* "-" epoch may be earlier than that reported.
*Note (5):
* ">" or "<" if period is a lower or upper limit, respectively.
* "(" indicates that the period is the mean cycle time of a U Gem
* or recurrent nova (a closing bracket exists in u_Period)
*Note (6): the codes are:
* "*N" if the period may be a multiple of the quoted value;
* "/N" if the period may be a fraction of the quoted value;
* ")" closing bracket of l_Period for U Gem-type or recurrent novae.
*Note (7): Besides "N", the codes are:
* "=" if the star is also known under another variable name
* "+" if the star is in other catalogues (e.g. Hipparcos or Name Lists)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gcvs_rem.dat' ! Remarks on GCVS stars
integer*4 nr__1
parameter (nr__1=14165) ! Number of records
character*94 ar__1 ! Full-size record
character*10 GCVS_1 (nr__1) ! Variable star designation
character*81 Text (nr__1) ! Text of note.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'crossid.dat' ! Cross-identifications (Vol. IV)
integer*4 nr__2
parameter (nr__2=192612) ! Number of records
character*48 ar__2 ! Full-size record
character*30 Name (nr__2) ! Alternative name (catalog or list acronym is
* detailed in the file "abbrev.txt")
character*13 VarName_1 (nr__2) ! Designation in GCVS or NSV (G8)
character*1 u_VarName (nr__2) ! [:?] Flag for doubtful identifications
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'nsv_cat.dat' ! The Suspected Variable stars and Supplement (2011)
integer*4 nr__3
parameter (nr__3=26017) ! Number of records
character*132 ar__3 ! Full-size record
real*8 RA__1 (nr__3) ! (deg) Right Ascension J2000.0
real*8 Dec_1 (nr__3) ! (deg) Declination J2000.0
integer*4 NSV (nr__3) ! NSV number (above 15000 for the Supplement)
character*1 m_NSV (nr__3) ! NSV letter suffix (1)
character*1 u_NSV (nr__3) ! [-] Dubitancy flag if a star's variability
* seems doubtful or erroneous to the compilers
character*1 n_NSV (nr__3) ! [*] indicate a remark in file "nsv.rem.dat"
integer*4 RA1950h (nr__3) ! (h) ? Hours RA, equinox B1950.0 (2)
integer*4 RA1950m (nr__3) ! (min) ? Minutes RA, equinox B1950.0 (2)
real*4 RA1950s (nr__3) ! (s) ? Seconds RA, equinox B1950.0 (2)
character*1 DE1950_ (nr__3) ! ? Sign Dec, equinox B1950.0 (2)
integer*4 DE1950d (nr__3) ! (deg) ? Degrees Dec, equinox B1950.0 (2)
integer*4 DE1950m (nr__3) ! (arcmin) ? Minutes Dec, equinox B1950.0 (2)
integer*4 DE1950s (nr__3) ! (arcsec) ? Seconds Dec, equinox B1950.0 (2)
character*1 u_DE1950s (nr__3) ! [:*] accuracy flags (G2)
integer*4 RAh_1 (nr__3) ! (h) ? Hours RA, equinox J2000.0 (2)
integer*4 RAm_1 (nr__3) ! (min) ? Minutes RA, equinox J2000.0 (2)
real*4 RAs_1 (nr__3) ! (s) ? Seconds RA, equinox J2000.0 (2)
character*1 DE__1 (nr__3) ! ? Sign Dec, equinox J2000.0 (2)
integer*4 DEd_1 (nr__3) ! (deg) ? Degrees Dec, equinox J2000.0 (2)
integer*4 DEm_1 (nr__3) ! (arcmin) ? Minutes Dec, equinox J2000.0 (2)
integer*4 DEs_1 (nr__3) ! (arcsec) ? Seconds Dec, equinox J2000.0 (2)
character*1 u_DEs_1 (nr__3) ! [:*] accuracy flags (G2)
character*6 VarType_1 (nr__3) ! Type of variability (see file "vartype.txt")
character*1 l_magMax_1 (nr__3) ! [><(] Limit or amplitude symbol on magMax (G3)
real*4 magMax_1 (nr__3) ! (mag) ? Magnitude at maximum brightness (3)
character*1 u_magMax_1 (nr__3) ! [:] Uncertainty flag on magMax
character*1 f_magMin (nr__3) ! [(] indicates that magMin is an amplitude
character*1 l_magMin (nr__3) ! [><] Limit symbol on magMin (G3)
real*4 magMin (nr__3) ! (mag) ? Minimum magnitude or amplitude (3)
character*1 u_magMin (nr__3) ! [st:*] Uncertainty or steps on magMin (4)
character*1 n_magMin (nr__3) ! Photometric system for magMin (G4)
character*2 flt_1 (nr__3) ! Photometric system for magnitudes (G4)
character*5 Ref1_1 (nr__3) ! Reference to a study of the star (G7)
character*5 Ref2_1 (nr__3) ! Chart reference (G7)
character*15 Desig (nr__3) ! Designation in paper referenced by Ref1 (5)
character*16 SpType_1 (nr__3) ! Spectral type (6)
character*1 f_NSV (nr__3) ! [=|] '=' for named variables.
character*14 VarName_2 (nr__3) ! Designation in GCVS (G8)
*Note (1):
* The letter "A" is present for one star, NSV 10360A, which is a
* completely different object from NSV 10360.
*Note (2):
* The position is not listed when the variable is equivalent to a GCVS
* star (see column VarName).
*Note (3):
* Magnitudes are reported to hundredths if the observations are
* photoelectric, to tenths or whole magnitudes if they are not.
* If only an amplitude has been measured photoelectrically, then
* the maximum magnitude is generally given to tenths only and the
* minimum is reported to hundredths.
*Note (4):
* "st" indicates that a light amplitude is reported for the minimum
* and is given in steps (instrumental measure).
*Note (5): (empty in the NSV Supplement part):
* identification of the suspected variable in the paper referenced in
* Ref1. DM numbers are given without prefix, the standard naming
* convention of "The Henry Draper Catalogue" being used (BD for zones
* +89 to -22; CD for -23 to -51; CP for -52 to -89). Identifications
* are sometimes made of coordinates, a 6-digit number consisting of
* hours, minutes, and seconds (or tenths of a minute) of time and
* degrees of declination with sign included. An asterisk (*) signifies
* that a designation is given in the remarks to the published catalog.
*Note (6):
* HD spectral types are generally in brackets, e.g. (A0).
* Some special symbols are used:
* AF for A-F stars
* FG for F-G stars
* KM for K-M stars
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'nsv_rem.dat' ! Remarks on NSV (and Supplement) stars
integer*4 nr__4
parameter (nr__4=10768) ! Number of records
character*80 ar__4 ! Full-size record
integer*4 NSV_1 (nr__4) ! NSV number (above 15000 for the Supplement)
character*75 Text_1 (nr__4) ! Text of remark
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'evs_cat.dat' ! Extragalactic Variable Stars. Catalogue (Vol. V)
integer*4 nr__5
parameter (nr__5=10979) ! Number of records
character*160 ar__5 ! Full-size record
real*8 RA__2 (nr__5) ! (deg) Right Ascension B1950.0
real*8 Dec_2 (nr__5) ! (deg) Declination B1950.0
integer*4 VarNum_1 (nr__5) ! [890001/1230042] Variable number (1)
character*12 Name_1 (nr__5) ! Variable star designation (3)
character*1 n_Name (nr__5) ! [*] indicates a remark in file evs_rem.dat
integer*4 RAh_2 (nr__5) ! (h) Right ascension, equinox B1950.0
integer*4 RAm_2 (nr__5) ! (min) (minutes)
real*4 RAs_2 (nr__5) ! (s) (seconds)
character*1 DE__2 (nr__5) ! sign of declination, B1950.0
integer*4 DEd_2 (nr__5) ! (deg) Declination, equinox B1950.0
integer*4 DEm_2 (nr__5) ! (arcmin) (arc minutes)
real*4 DEs_2 (nr__5) ! (arcsec) (arc seconds)
character*1 u_DEs_2 (nr__5) ! [":*!] accuracy flags on position (G2)
character*8 VarType_2 (nr__5) ! Type of variability (see file "vartype.txt")
real*4 magMax_2 (nr__5) ! (mag) ? Magnitude at maximum brightness
character*1 u_magMax_2 (nr__5) ! Uncertainty flag (:) on magMax
character*1 l_magMin_1 (nr__5) ! [<(] Limit or amplitude symbol on magMax (G3)
real*4 magMin_1 (nr__5) ! (mag) ? Magnitude at minimum brightness or amplitude
character*1 u_magMin_1 (nr__5) ! Uncertainty flag (:) on magMin
character*1 flt_2 (nr__5) ! The photometric system for magnitudes (G4)
real*8 Epoch_1 (nr__5) ! (d) ? Epoch for maximum light, Julian days (G5)
character*1 u_Epoch_1 (nr__5) ! Uncertainty flag on Epoch
character*1 l_P (nr__5) ! [(] if mean cycle time of U Gem-type variable
real*8 P (nr__5) ! (d) ? Period (mean for U Gem-type)
character*1 u_P (nr__5) ! Uncertainty flag (:) on P
integer*4 M_m_D_1 (nr__5) ! (%) ? Rising time or duration of eclipse (G6)
character*1 u_M_m_D_1 (nr__5) ! Uncertainty flag (:) on M-m/D
character*1 n_M_m_D_1 (nr__5) ! [*] Note for eclipsing variable (G6)
character*6 SpType_2 (nr__5) ! Spectral type (4)
character*5 Ref1_2 (nr__5) ! Reference for study of the star (G7)
character*5 Ref2_2 (nr__5) ! Reference for chart or photograph (G7)
character*14 Desig1 (nr__5) ! designation in paper Ref1
character*12 Desig2 (nr__5) ! designation in paper Ref2
character*9 VarName_3 (nr__5) ! Designation in GCVS (G8)
character*2 Mm (nr__5) ! [n: ] Non-member note (5)
integer*4 Year_1 (nr__5) ! (yr) ? Year of outburst for novae or supernovae
character*1 u_Year_1 (nr__5) ! [:] Uncertainty flag on YearNova
*Note (1): Number of the variable star.
* This numbering system is made of a 3-digit number identifying the
* galaxy, followed by a 4-digit number identifying the variable in the
* galaxy. The numbering system for galaxies begins with 089, the first
* 88 codes being reserved for constellations (cf. column "VarNum" in
* file "gcvs_cat.dat")
*
* If in the literature for the considered galaxy there exist several
* different variable star numbering systems, we have introduced new
* numbers in the order of increasing right ascension for the equinox
* 1950.0. But if the variable star numbering system for the galaxy is
* practically of general use, we have retained it - if necessary, with
* additional numbers ordered by increasing right ascension.
*Note (3): Variable star Designation.
* The name of the variable in the adopted system of variable star
* designations: name of the galaxy + "V" + variable star No. within the
* galaxy (see column VarNo)
*Note (4):
* If there exist several spectral type determinations for a star, we
* preferred more recent ones, taking into account reliability of each
* determination. We converted symbols "d" (dwarfs) and "g" (giants) into
* luminosity classes V and III. If a nova showed spectral features
* typical of novae during or after the outburst, a symbol (NOVA) is
* given in this column. If the spectrum shows features characteristic of
* U Gem variables, this column contains the symbol (UG). Continuous
* spectrum is designated "cont", the symbol "e" means emission spectrum,
* lower case letters "ea", "eb" in this column stand for "e alpha", "e
* beta". For RR Lyrae stars, spectral types from hydrogen lines are
* given. A plus sign (+) between two spectral type values means that
* spectra of two components of a spectroscopic binary are observed. Two
* spectral type values separated with a minus sign (-) mean the range of
* spectral type variations in the process of brightness changes.
*Note (5):
* Symbols "n" or "n:" for "nonmembers" or "possible nonmembers" of the
* galaxy in question. For the stars that enter the GCVS vols. I-III
* (file gcvs_cat.dat) or the NSV catalogue (file "nsv_cat.dat") and are
* not members of the corresponding galaxy, the present volume gives only
* the running number, coordinates, the symbol "n", and the GCVS or NSV
* name.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'evs_rem.dat' ! Extragalactic Variable Stars. Remarks
integer*4 nr__6
parameter (nr__6=1936) ! Number of records
character*80 ar__6 ! Full-size record
character*12 Name_2 (nr__6) ! Variable star designation
character*67 Text_2 (nr__6) ! Text of remark
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References for all variables (GCVS, NSV, EVS)
integer*4 nr__7
parameter (nr__7=22805) ! Number of records
character*80 ar__7 ! Full-size record
character*5 Ref (nr__7) ! Reference number (G7)
character*74 Text_3 (nr__7) ! Text of reference
C=============================================================================
C Loading file 'gcvs_cat.dat' ! The GCVS Catalog (Vol. I-III, version 2013-04-30)
C Format for file interpretation
1 format(
+ I6,A1,1X,A10,A1,1X,I2,I2,F4.1,A1,I2,I2,I2,A1,1X,A10,1X,A1,
+ F6.3,A1,2X,A1,F6.3,A1,A2,2X,A1,F6.3,A1,A2,2X,A2,1X,F13.5,A1,
+ 1X,A4,A1,1X,A1,F16.10,A3,1X,I2,A1,A1,2X,A17,1X,A5,1X,A5,1X,A1,
+ A10)
C Effective file loading
open(unit=1,file='gcvs_cat.dat', status='old')
write(6,*) '....Loading file: gcvs_cat.dat'
do i__=1,47969
read(1,'(A174)')ar__
read(ar__,1)
+ VarNum(i__),m_VarNum(i__),GCVS(i__),n_GCVS(i__),RAh(i__),
+ RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
+ u_DEs(i__),VarType(i__),l_magMax(i__),magMax(i__),
+ u_magMax(i__),l_Min1(i__),Min1(i__),u_Min1(i__),n_Min1(i__),
+ l_Min2(i__),Min2(i__),u_Min2(i__),n_Min2(i__),flt(i__),
+ Epoch(i__),u_Epoch(i__),Year(i__),u_Year(i__),l_Period(i__),
+ Period(i__),u_Period(i__),M_m_D(i__),u_M_m_D(i__),
+ n_M_m_D(i__),SpType(i__),Ref1(i__),Ref2(i__),f_GCVS(i__),
+ VarName(i__)
if(ar__(21:22) .EQ. '') RAh(i__) = iNULL__
if(ar__(23:24) .EQ. '') RAm(i__) = iNULL__
if(ar__(25:28) .EQ. '') RAs(i__) = rNULL__
if(ar__(30:31) .EQ. '') DEd(i__) = iNULL__
if(ar__(32:33) .EQ. '') DEm(i__) = iNULL__
if(ar__(34:35) .EQ. '') DEs(i__) = iNULL__
if(ar__(50:55) .EQ. '') magMax(i__) = rNULL__
if(ar__(60:65) .EQ. '') Min1(i__) = rNULL__
if(ar__(72:77) .EQ. '') Min2(i__) = rNULL__
if(ar__(86:98) .EQ. '') Epoch(i__) = rNULL__
if(ar__(108:123) .EQ. '') Period(i__) = rNULL__
if(ar__(128:129) .EQ. '') M_m_D(i__) = iNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ VarNum(i__),m_VarNum(i__),GCVS(i__),n_GCVS(i__),RAh(i__),
+ RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),
+ u_DEs(i__),VarType(i__),l_magMax(i__),magMax(i__),
+ u_magMax(i__),l_Min1(i__),Min1(i__),u_Min1(i__),n_Min1(i__),
+ l_Min2(i__),Min2(i__),u_Min2(i__),n_Min2(i__),flt(i__),
+ Epoch(i__),u_Epoch(i__),Year(i__),u_Year(i__),l_Period(i__),
+ Period(i__),u_Period(i__),M_m_D(i__),u_M_m_D(i__),
+ n_M_m_D(i__),SpType(i__),Ref1(i__),Ref2(i__),f_GCVS(i__),
+ VarName(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gcvs_rem.dat' ! Remarks on GCVS stars
C Format for file interpretation
2 format(A10,3X,A81)
C Effective file loading
open(unit=1,file='gcvs_rem.dat', status='old')
write(6,*) '....Loading file: gcvs_rem.dat'
do i__=1,14165
read(1,'(A94)')ar__1
read(ar__1,2)GCVS_1(i__),Text(i__)
c ..............Just test output...........
write(6,2)GCVS_1(i__),Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'crossid.dat' ! Cross-identifications (Vol. IV)
C Format for file interpretation
3 format(A30,3X,A13,1X,A1)
C Effective file loading
open(unit=1,file='crossid.dat', status='old')
write(6,*) '....Loading file: crossid.dat'
do i__=1,192612
read(1,'(A48)')ar__2
read(ar__2,3)Name(i__),VarName_1(i__),u_VarName(i__)
c ..............Just test output...........
write(6,3)Name(i__),VarName_1(i__),u_VarName(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'nsv_cat.dat' ! The Suspected Variable stars and Supplement (2011)
C Format for file interpretation
4 format(
+ I5,A1,A1,A1,1X,I2,I2,F4.1,A1,I2,I2,I2,A1,1X,I2,I2,F4.1,A1,I2,
+ I2,I2,A1,1X,A6,1X,A1,F5.2,A1,1X,A1,A1,F6.3,A1,A1,2X,A2,1X,A5,
+ 1X,A5,1X,A15,1X,A16,A1,A14)
C Effective file loading
open(unit=1,file='nsv_cat.dat', status='old')
write(6,*) '....Loading file: nsv_cat.dat'
do i__=1,26017
read(1,'(A132)')ar__3
read(ar__3,4)
+ NSV(i__),m_NSV(i__),u_NSV(i__),n_NSV(i__),RA1950h(i__),
+ RA1950m(i__),RA1950s(i__),DE1950_(i__),DE1950d(i__),
+ DE1950m(i__),DE1950s(i__),u_DE1950s(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),u_DEs_1(i__),VarType_1(i__),l_magMax_1(i__),
+ magMax_1(i__),u_magMax_1(i__),f_magMin(i__),l_magMin(i__),
+ magMin(i__),u_magMin(i__),n_magMin(i__),flt_1(i__),
+ Ref1_1(i__),Ref2_1(i__),Desig(i__),SpType_1(i__),f_NSV(i__),
+ VarName_2(i__)
if(ar__3(10:11) .EQ. '') RA1950h(i__) = iNULL__
if(ar__3(12:13) .EQ. '') RA1950m(i__) = iNULL__
if(ar__3(14:17) .EQ. '') RA1950s(i__) = rNULL__
if(ar__3(19:20) .EQ. '') DE1950d(i__) = iNULL__
if(ar__3(21:22) .EQ. '') DE1950m(i__) = iNULL__
if(ar__3(23:24) .EQ. '') DE1950s(i__) = iNULL__
if(ar__3(27:28) .EQ. '') RAh_1(i__) = iNULL__
if(ar__3(29:30) .EQ. '') RAm_1(i__) = iNULL__
if(ar__3(31:34) .EQ. '') RAs_1(i__) = rNULL__
if(ar__3(36:37) .EQ. '') DEd_1(i__) = iNULL__
if(ar__3(38:39) .EQ. '') DEm_1(i__) = iNULL__
if(ar__3(40:41) .EQ. '') DEs_1(i__) = iNULL__
if(ar__3(52:56) .EQ. '') magMax_1(i__) = rNULL__
if(ar__3(61:66) .EQ. '') magMin(i__) = rNULL__
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1(i__) = RAh_1(i__)
if(RAh_1(i__) .GE. 0) RA__1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAs_1(i__)/240.
Dec_1(i__) = DEd_1(i__)
if(DEm_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c ..............Just test output...........
write(6,4)
+ NSV(i__),m_NSV(i__),u_NSV(i__),n_NSV(i__),RA1950h(i__),
+ RA1950m(i__),RA1950s(i__),DE1950_(i__),DE1950d(i__),
+ DE1950m(i__),DE1950s(i__),u_DE1950s(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),u_DEs_1(i__),VarType_1(i__),l_magMax_1(i__),
+ magMax_1(i__),u_magMax_1(i__),f_magMin(i__),l_magMin(i__),
+ magMin(i__),u_magMin(i__),n_magMin(i__),flt_1(i__),
+ Ref1_1(i__),Ref2_1(i__),Desig(i__),SpType_1(i__),f_NSV(i__),
+ VarName_2(i__)
write(6,'(6H Pos: 2F8.4)') RA__1(i__),Dec_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'nsv_rem.dat' ! Remarks on NSV (and Supplement) stars
C Format for file interpretation
5 format(I5,A75)
C Effective file loading
open(unit=1,file='nsv_rem.dat', status='old')
write(6,*) '....Loading file: nsv_rem.dat'
do i__=1,10768
read(1,'(A80)')ar__4
read(ar__4,5)NSV_1(i__),Text_1(i__)
c ..............Just test output...........
write(6,5)NSV_1(i__),Text_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'evs_cat.dat' ! Extragalactic Variable Stars. Catalogue (Vol. V)
C Format for file interpretation
6 format(
+ I7,1X,A12,A1,I2,I2,F5.2,A1,I2,I2,F4.1,A1,A8,F5.2,A1,A1,F6.3,
+ A1,1X,A1,F14.6,A1,A1,F13.8,A1,1X,I2,A1,A1,A6,1X,A5,A5,1X,A14,
+ A12,A9,1X,A2,I4,A1)
C Effective file loading
open(unit=1,file='evs_cat.dat', status='old')
write(6,*) '....Loading file: evs_cat.dat'
do i__=1,10979
read(1,'(A160)')ar__5
read(ar__5,6)
+ VarNum_1(i__),Name_1(i__),n_Name(i__),RAh_2(i__),RAm_2(i__),
+ RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
+ u_DEs_2(i__),VarType_2(i__),magMax_2(i__),u_magMax_2(i__),
+ l_magMin_1(i__),magMin_1(i__),u_magMin_1(i__),flt_2(i__),
+ Epoch_1(i__),u_Epoch_1(i__),l_P(i__),P(i__),u_P(i__),
+ M_m_D_1(i__),u_M_m_D_1(i__),n_M_m_D_1(i__),SpType_2(i__),
+ Ref1_2(i__),Ref2_2(i__),Desig1(i__),Desig2(i__),
+ VarName_3(i__),Mm(i__),Year_1(i__),u_Year_1(i__)
if(ar__5(49:53) .EQ. '') magMax_2(i__) = rNULL__
if(ar__5(56:61) .EQ. '') magMin_1(i__) = rNULL__
if(ar__5(65:78) .EQ. '') Epoch_1(i__) = rNULL__
if(ar__5(81:93) .EQ. '') P(i__) = rNULL__
if(ar__5(96:97) .EQ. '') M_m_D_1(i__) = iNULL__
if(ar__5(156:159) .EQ. '') Year_1(i__) = iNULL__
c Derive coordinates RA__2 and Dec_2 from input data
c (RA__2 and Dec_2 are set to rNULL_ when unknown)
RA__2(i__) = RAh_2(i__)
if(RAh_2(i__) .GE. 0) RA__2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GE. 0) RA__2(i__)=RA__2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GE. 0) RA__2(i__)=RA__2(i__)+RAs_2(i__)/240.
Dec_2(i__) = DEd_2(i__)
if(DEm_2(i__) .GE. 0) Dec_2(i__)=Dec_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) Dec_2(i__)=Dec_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.Dec_2(i__).GE.0) Dec_2(i__)=-Dec_2(i__)
c ..............Just test output...........
write(6,6)
+ VarNum_1(i__),Name_1(i__),n_Name(i__),RAh_2(i__),RAm_2(i__),
+ RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),
+ u_DEs_2(i__),VarType_2(i__),magMax_2(i__),u_magMax_2(i__),
+ l_magMin_1(i__),magMin_1(i__),u_magMin_1(i__),flt_2(i__),
+ Epoch_1(i__),u_Epoch_1(i__),l_P(i__),P(i__),u_P(i__),
+ M_m_D_1(i__),u_M_m_D_1(i__),n_M_m_D_1(i__),SpType_2(i__),
+ Ref1_2(i__),Ref2_2(i__),Desig1(i__),Desig2(i__),
+ VarName_3(i__),Mm(i__),Year_1(i__),u_Year_1(i__)
write(6,'(6H Pos: 2F8.4)') RA__2(i__),Dec_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'evs_rem.dat' ! Extragalactic Variable Stars. Remarks
C Format for file interpretation
7 format(A12,1X,A67)
C Effective file loading
open(unit=1,file='evs_rem.dat', status='old')
write(6,*) '....Loading file: evs_rem.dat'
do i__=1,1936
read(1,'(A80)')ar__6
read(ar__6,7)Name_2(i__),Text_2(i__)
c ..............Just test output...........
write(6,7)Name_2(i__),Text_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References for all variables (GCVS, NSV, EVS)
C Format for file interpretation
8 format(A5,1X,A74)
C Effective file loading
open(unit=1,file='refs.dat', status='old')
write(6,*) '....Loading file: refs.dat'
do i__=1,22805
read(1,'(A80)')ar__7
read(ar__7,8)Ref(i__),Text_3(i__)
c ..............Just test output...........
write(6,8)Ref(i__),Text_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end