Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/598/A108 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-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__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. J/A+A/598/A108 Variability survey of brightest stars (Laur+, 2017) *================================================================================ *Variability survey of brightest stars in selected OB associations. * Laur J., Kolka I., Eenmae T., Tuvikene T., Leedjarv L. * <Astron. Astrophys. 598, A108 (2017)> * =2017A&A...598A.108L (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Observed eclipsing binaries integer*4 nr__ parameter (nr__=80) ! Number of records character*132 ar__ ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg ! (deg) Right Ascension J2000.0 real*8 DEdeg ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field ! Field name integer*4 ID ! Internal ID character*24 Name ! Star name integer*4 RAh ! (h) Right ascension (J2000.0) integer*4 RAm ! (min) Right ascension (J2000.0) real*4 RAs ! (s) Right ascension (J2000.0) character*1 DE_ ! Declination sign (J2000.0) integer*4 DEd ! (deg) Declination (J2000.0) integer*4 DEm ! (arcmin) Declination (J2000.0) integer*4 DEs ! (arcsec) Declination (J2000.0) character*12 Member ! Cluster membership character*13 SPType ! MK Spectral Type (G1) real*4 Vmag ! (mag) Mean magnitude (apparent) in V filter real*4 Amp ! (mag) Amplitude of the light curve real*8 Per ! (d) ?=- Period from this work (1) character*10 VSXinfo ! Variability type in VSX (2) real*8 VSXPer ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Observed {beta} Cephei variables integer*4 nr__1 parameter (nr__1=13) ! Number of records character*132 ar__1 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 ! (deg) Right Ascension J2000.0 real*8 DEdeg_1 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_1 ! Field name integer*4 ID_1 ! Internal ID character*24 Name_1 ! Star name integer*4 RAh_1 ! (h) Right ascension (J2000.0) integer*4 RAm_1 ! (min) Right ascension (J2000.0) real*4 RAs_1 ! (s) Right ascension (J2000.0) character*1 DE__1 ! Declination sign (J2000.0) integer*4 DEd_1 ! (deg) Declination (J2000.0) integer*4 DEm_1 ! (arcmin) Declination (J2000.0) integer*4 DEs_1 ! (arcsec) Declination (J2000.0) character*12 Member_1 ! Cluster membership character*13 SPType_1 ! MK Spectral Type (G1) real*4 Vmag_1 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_1 ! (mag) Amplitude of the light curve real*8 Per_1 ! (d) ?=- Period from this work (1) character*10 VSXinfo_1 ! Variability type in VSX (2) real*8 VSXPer_1 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Observed SPB variables integer*4 nr__2 parameter (nr__2=16) ! Number of records character*132 ar__2 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 ! (deg) Right Ascension J2000.0 real*8 DEdeg_2 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_2 ! Field name integer*4 ID_2 ! Internal ID character*24 Name_2 ! Star name integer*4 RAh_2 ! (h) Right ascension (J2000.0) integer*4 RAm_2 ! (min) Right ascension (J2000.0) real*4 RAs_2 ! (s) Right ascension (J2000.0) character*1 DE__2 ! Declination sign (J2000.0) integer*4 DEd_2 ! (deg) Declination (J2000.0) integer*4 DEm_2 ! (arcmin) Declination (J2000.0) integer*4 DEs_2 ! (arcsec) Declination (J2000.0) character*12 Member_2 ! Cluster membership character*13 SPType_2 ! MK Spectral Type (G1) real*4 Vmag_2 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_2 ! (mag) Amplitude of the light curve real*8 Per_2 ! (d) ?=- Period from this work (1) character*10 VSXinfo_2 ! Variability type in VSX (2) real*8 VSXPer_2 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Observed Cepheid variables integer*4 nr__3 parameter (nr__3=7) ! Number of records character*132 ar__3 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_3 ! (deg) Right Ascension J2000.0 real*8 DEdeg_3 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_3 ! Field name integer*4 ID_3 ! Internal ID character*24 Name_3 ! Star name integer*4 RAh_3 ! (h) Right ascension (J2000.0) integer*4 RAm_3 ! (min) Right ascension (J2000.0) real*4 RAs_3 ! (s) Right ascension (J2000.0) character*1 DE__3 ! Declination sign (J2000.0) integer*4 DEd_3 ! (deg) Declination (J2000.0) integer*4 DEm_3 ! (arcmin) Declination (J2000.0) integer*4 DEs_3 ! (arcsec) Declination (J2000.0) character*12 Member_3 ! Cluster membership character*13 SPType_3 ! MK Spectral Type (G1) real*4 Vmag_3 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_3 ! (mag) Amplitude of the light curve real*8 Per_3 ! (d) ?=- Period from this work (1) character*10 VSXinfo_3 ! Variability type in VSX (2) real*8 VSXPer_3 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! Observed {gamma} Doradus variables integer*4 nr__4 parameter (nr__4=1) ! Number of records character*132 ar__4 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_4 ! (deg) Right Ascension J2000.0 real*8 DEdeg_4 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_4 ! Field name integer*4 ID_4 ! Internal ID character*24 Name_4 ! Star name integer*4 RAh_4 ! (h) Right ascension (J2000.0) integer*4 RAm_4 ! (min) Right ascension (J2000.0) real*4 RAs_4 ! (s) Right ascension (J2000.0) character*1 DE__4 ! Declination sign (J2000.0) integer*4 DEd_4 ! (deg) Declination (J2000.0) integer*4 DEm_4 ! (arcmin) Declination (J2000.0) integer*4 DEs_4 ! (arcsec) Declination (J2000.0) character*12 Member_4 ! Cluster membership character*13 SPType_4 ! MK Spectral Type (G1) real*4 Vmag_4 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_4 ! (mag) Amplitude of the light curve real*8 Per_4 ! (d) ?=- Period from this work (1) character*10 VSXinfo_4 ! Variability type in VSX (2) real*8 VSXPer_4 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! Observed WR variables integer*4 nr__5 parameter (nr__5=3) ! Number of records character*132 ar__5 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_5 ! (deg) Right Ascension J2000.0 real*8 DEdeg_5 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_5 ! Field name integer*4 ID_5 ! Internal ID character*24 Name_5 ! Star name integer*4 RAh_5 ! (h) Right ascension (J2000.0) integer*4 RAm_5 ! (min) Right ascension (J2000.0) real*4 RAs_5 ! (s) Right ascension (J2000.0) character*1 DE__5 ! Declination sign (J2000.0) integer*4 DEd_5 ! (deg) Declination (J2000.0) integer*4 DEm_5 ! (arcmin) Declination (J2000.0) integer*4 DEs_5 ! (arcsec) Declination (J2000.0) character*12 Member_5 ! Cluster membership character*13 SPType_5 ! MK Spectral Type (G1) real*4 Vmag_5 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_5 ! (mag) Amplitude of the light curve real*8 Per_5 ! (d) ?=- Period from this work (1) character*10 VSXinfo_5 ! Variability type in VSX (2) real*8 VSXPer_5 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! Observed {alpha} Cyg variables integer*4 nr__6 parameter (nr__6=30) ! Number of records character*132 ar__6 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_6 ! (deg) Right Ascension J2000.0 real*8 DEdeg_6 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_6 ! Field name integer*4 ID_6 ! Internal ID character*24 Name_6 ! Star name integer*4 RAh_6 ! (h) Right ascension (J2000.0) integer*4 RAm_6 ! (min) Right ascension (J2000.0) real*4 RAs_6 ! (s) Right ascension (J2000.0) character*1 DE__6 ! Declination sign (J2000.0) integer*4 DEd_6 ! (deg) Declination (J2000.0) integer*4 DEm_6 ! (arcmin) Declination (J2000.0) integer*4 DEs_6 ! (arcsec) Declination (J2000.0) character*12 Member_6 ! Cluster membership character*13 SPType_6 ! MK Spectral Type (G1) real*4 Vmag_6 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_6 ! (mag) Amplitude of the light curve real*8 Per_6 ! (d) ?=- Period from this work (1) character*10 VSXinfo_6 ! Variability type in VSX (2) real*8 VSXPer_6 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table8.dat' ! Observed variable Be stars integer*4 nr__7 parameter (nr__7=62) ! Number of records character*132 ar__7 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_7 ! (deg) Right Ascension J2000.0 real*8 DEdeg_7 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_7 ! Field name integer*4 ID_7 ! Internal ID character*24 Name_7 ! Star name integer*4 RAh_7 ! (h) Right ascension (J2000.0) integer*4 RAm_7 ! (min) Right ascension (J2000.0) real*4 RAs_7 ! (s) Right ascension (J2000.0) character*1 DE__7 ! Declination sign (J2000.0) integer*4 DEd_7 ! (deg) Declination (J2000.0) integer*4 DEm_7 ! (arcmin) Declination (J2000.0) integer*4 DEs_7 ! (arcsec) Declination (J2000.0) character*12 Member_7 ! Cluster membership character*13 SPType_7 ! MK Spectral Type (G1) real*4 Vmag_7 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_7 ! (mag) Amplitude of the light curve real*8 Per_7 ! (d) ?=- Period from this work (1) character*10 VSXinfo_7 ! Variability type in VSX (2) real*8 VSXPer_7 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! Observed variable late-type stars integer*4 nr__8 parameter (nr__8=63) ! Number of records character*132 ar__8 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_8 ! (deg) Right Ascension J2000.0 real*8 DEdeg_8 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_8 ! Field name integer*4 ID_8 ! Internal ID character*24 Name_8 ! Star name integer*4 RAh_8 ! (h) Right ascension (J2000.0) integer*4 RAm_8 ! (min) Right ascension (J2000.0) real*4 RAs_8 ! (s) Right ascension (J2000.0) character*1 DE__8 ! Declination sign (J2000.0) integer*4 DEd_8 ! (deg) Declination (J2000.0) integer*4 DEm_8 ! (arcmin) Declination (J2000.0) integer*4 DEs_8 ! (arcsec) Declination (J2000.0) character*12 Member_8 ! Cluster membership character*13 SPType_8 ! MK Spectral Type (G1) real*4 Vmag_8 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_8 ! (mag) Amplitude of the light curve real*8 Per_8 ! (d) ?=- Period from this work (1) character*10 VSXinfo_8 ! Variability type in VSX (2) real*8 VSXPer_8 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table10.dat' ! Observed unclassified variable stars integer*4 nr__9 parameter (nr__9=80) ! Number of records character*132 ar__9 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_9 ! (deg) Right Ascension J2000.0 real*8 DEdeg_9 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_9 ! Field name integer*4 ID_9 ! Internal ID character*24 Name_9 ! Star name integer*4 RAh_9 ! (h) Right ascension (J2000.0) integer*4 RAm_9 ! (min) Right ascension (J2000.0) real*4 RAs_9 ! (s) Right ascension (J2000.0) character*1 DE__9 ! Declination sign (J2000.0) integer*4 DEd_9 ! (deg) Declination (J2000.0) integer*4 DEm_9 ! (arcmin) Declination (J2000.0) integer*4 DEs_9 ! (arcsec) Declination (J2000.0) character*12 Member_9 ! Cluster membership character*13 SPType_9 ! MK Spectral Type (G1) real*4 Vmag_9 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_9 ! (mag) Amplitude of the light curve real*8 Per_9 ! (d) ?=- Period from this work (1) character*10 VSXinfo_9 ! Variability type in VSX (2) real*8 VSXPer_9 ! (d) ?=- Period in VSX (2) *Note (1): "-" means no period found *Note (2): "-" means information missing in the VSX catalogue c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea11.dat' ! Variable stars not in VSX integer*4 nr__10 parameter (nr__10=221) ! Number of records character*116 ar__10 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_10 ! (deg) Right Ascension J2000.0 real*8 DEdeg_10 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*11 Field_10 ! Field name integer*4 ID_10 ! Internal ID character*24 Name_10 ! Star name integer*4 RAh_10 ! (h) Right ascension (J2000.0) integer*4 RAm_10 ! (min) Right ascension (J2000.0) real*4 RAs_10 ! (s) Right ascension (J2000.0) character*1 DE__10 ! Declination sign (J2000.0) integer*4 DEd_10 ! (deg) Declination (J2000.0) integer*4 DEm_10 ! (arcmin) Declination (J2000.0) integer*4 DEs_10 ! (arcsec) Declination (J2000.0) character*12 Member_10 ! Cluster membership character*13 SPType_10 ! MK Spectral Type (G1) real*4 Vmag_10 ! (mag) Mean magnitude (apparent) in V filter real*4 Amp_10 ! (mag) Amplitude of the light curve real*8 Per_10 ! (d) ?=- Period from this work (1) character*4 Type ! Variability type (2) *Note (1): "-" means no period found *Note (2): Variability type as follows: * E = Eclipsing binary * BE = Be star * BCEP = betaCephei type star * SPB = Slowly pulsating B star * ACYG = alpha Cygni type star * CEP = Cepheid * WR = Wolf-Rayet star * SR = late-type star * "-" = unclassified c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table11.dat' ! Light curves of variable stars integer*4 nr__11 parameter (nr__11=204141) ! Number of records character*43 ar__11 ! Full-size record character*11 Field_11 ! Field name integer*4 ID_11 ! Internal ID real*8 Obs_date ! (d) Heliocentric Julian date (HJD-2400000) real*8 mag ! (mag) Differential magnitude in V passband real*4 e_mag ! (mag) Uncertainty in magnitude c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'lc/*' ! Individual light curves integer*4 nr__12 parameter (nr__12=354) ! Number of records character*26 ar__12 ! Full-size record real*8 Obs_date_1 ! (d) Heliocentric Julian date (HJD-2400000) real*8 mag_1 ! (mag) Differential magnitude in V passband real*4 e_mag_1 ! (mag) Uncertainty in magnitude C============================================================================= C Loading file 'table1.dat' ! Observed eclipsing binaries C Format for file interpretation 1 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,80 read(1,'(A132)')ar__ read(ar__,1) + Field,ID,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Member,SPType,Vmag, + Amp,Per,VSXinfo,VSXPer if (idig(ar__(104:111)).EQ.0) Per = rNULL__ if (idig(ar__(123:132)).EQ.0) VSXPer = 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) + Field,ID,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Member,SPType,Vmag, + Amp,Per,VSXinfo,VSXPer write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Observed {beta} Cephei variables C Format for file interpretation 2 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,13 read(1,'(A132)')ar__1 read(ar__1,2) + Field_1,ID_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1, + Member_1,SPType_1,Vmag_1,Amp_1,Per_1,VSXinfo_1,VSXPer_1 if (idig(ar__1(104:111)).EQ.0) Per_1 = rNULL__ if (idig(ar__1(123:132)).EQ.0) VSXPer_1 = rNULL__ RAdeg_1 = rNULL__ DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh_1*15. if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4. if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240. if(DEd_1 .GE. 0) DEdeg_1=DEd_1 if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60. if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600. if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1 c ..............Just test output........... write(6,2) + Field_1,ID_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1, + Member_1,SPType_1,Vmag_1,Amp_1,Per_1,VSXinfo_1,VSXPer_1 write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Observed SPB variables C Format for file interpretation 3 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,16 read(1,'(A132)')ar__2 read(ar__2,3) + Field_2,ID_2,Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2, + Member_2,SPType_2,Vmag_2,Amp_2,Per_2,VSXinfo_2,VSXPer_2 if (idig(ar__2(104:111)).EQ.0) Per_2 = rNULL__ if (idig(ar__2(123:132)).EQ.0) VSXPer_2 = rNULL__ RAdeg_2 = rNULL__ DEdeg_2 = 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 .GT. -180) RAdeg_2=RAh_2*15. if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4. if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240. if(DEd_2 .GE. 0) DEdeg_2=DEd_2 if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60. if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600. if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2 c ..............Just test output........... write(6,3) + Field_2,ID_2,Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2, + Member_2,SPType_2,Vmag_2,Amp_2,Per_2,VSXinfo_2,VSXPer_2 write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Observed Cepheid variables C Format for file interpretation 4 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,7 read(1,'(A132)')ar__3 read(ar__3,4) + Field_3,ID_3,Name_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3, + Member_3,SPType_3,Vmag_3,Amp_3,Per_3,VSXinfo_3,VSXPer_3 if (idig(ar__3(104:111)).EQ.0) Per_3 = rNULL__ if (idig(ar__3(123:132)).EQ.0) VSXPer_3 = rNULL__ RAdeg_3 = rNULL__ DEdeg_3 = 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 .GT. -180) RAdeg_3=RAh_3*15. if(RAm_3 .GT. -180) RAdeg_3=RAdeg_3+RAm_3/4. if(RAs_3 .GT. -180) RAdeg_3=RAdeg_3+RAs_3/240. if(DEd_3 .GE. 0) DEdeg_3=DEd_3 if(DEm_3 .GE. 0) DEdeg_3=DEdeg_3+DEm_3/60. if(DEs_3 .GE. 0) DEdeg_3=DEdeg_3+DEs_3/3600. if(DE__3.EQ.'-'.AND.DEdeg_3.GE.0) DEdeg_3=-DEdeg_3 c ..............Just test output........... write(6,4) + Field_3,ID_3,Name_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3, + Member_3,SPType_3,Vmag_3,Amp_3,Per_3,VSXinfo_3,VSXPer_3 write(6,'(6H Pos: 2F8.4)') RAdeg_3,DEdeg_3 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table5.dat' ! Observed {gamma} Doradus variables C Format for file interpretation 5 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,1 read(1,'(A132)')ar__4 read(ar__4,5) + Field_4,ID_4,Name_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4, + Member_4,SPType_4,Vmag_4,Amp_4,Per_4,VSXinfo_4,VSXPer_4 if (idig(ar__4(104:111)).EQ.0) Per_4 = rNULL__ if (idig(ar__4(123:132)).EQ.0) VSXPer_4 = rNULL__ RAdeg_4 = rNULL__ DEdeg_4 = rNULL__ c Derive coordinates RAdeg_4 and DEdeg_4 from input data c (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown) if(RAh_4 .GT. -180) RAdeg_4=RAh_4*15. if(RAm_4 .GT. -180) RAdeg_4=RAdeg_4+RAm_4/4. if(RAs_4 .GT. -180) RAdeg_4=RAdeg_4+RAs_4/240. if(DEd_4 .GE. 0) DEdeg_4=DEd_4 if(DEm_4 .GE. 0) DEdeg_4=DEdeg_4+DEm_4/60. if(DEs_4 .GE. 0) DEdeg_4=DEdeg_4+DEs_4/3600. if(DE__4.EQ.'-'.AND.DEdeg_4.GE.0) DEdeg_4=-DEdeg_4 c ..............Just test output........... write(6,5) + Field_4,ID_4,Name_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4, + Member_4,SPType_4,Vmag_4,Amp_4,Per_4,VSXinfo_4,VSXPer_4 write(6,'(6H Pos: 2F8.4)') RAdeg_4,DEdeg_4 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! Observed WR variables C Format for file interpretation 6 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,3 read(1,'(A132)')ar__5 read(ar__5,6) + Field_5,ID_5,Name_5,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,DEm_5,DEs_5, + Member_5,SPType_5,Vmag_5,Amp_5,Per_5,VSXinfo_5,VSXPer_5 if (idig(ar__5(104:111)).EQ.0) Per_5 = rNULL__ if (idig(ar__5(123:132)).EQ.0) VSXPer_5 = rNULL__ RAdeg_5 = rNULL__ DEdeg_5 = rNULL__ c Derive coordinates RAdeg_5 and DEdeg_5 from input data c (RAdeg_5 and DEdeg_5 are set to rNULL__ when unknown) if(RAh_5 .GT. -180) RAdeg_5=RAh_5*15. if(RAm_5 .GT. -180) RAdeg_5=RAdeg_5+RAm_5/4. if(RAs_5 .GT. -180) RAdeg_5=RAdeg_5+RAs_5/240. if(DEd_5 .GE. 0) DEdeg_5=DEd_5 if(DEm_5 .GE. 0) DEdeg_5=DEdeg_5+DEm_5/60. if(DEs_5 .GE. 0) DEdeg_5=DEdeg_5+DEs_5/3600. if(DE__5.EQ.'-'.AND.DEdeg_5.GE.0) DEdeg_5=-DEdeg_5 c ..............Just test output........... write(6,6) + Field_5,ID_5,Name_5,RAh_5,RAm_5,RAs_5,DE__5,DEd_5,DEm_5,DEs_5, + Member_5,SPType_5,Vmag_5,Amp_5,Per_5,VSXinfo_5,VSXPer_5 write(6,'(6H Pos: 2F8.4)') RAdeg_5,DEdeg_5 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! Observed {alpha} Cyg variables C Format for file interpretation 7 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,30 read(1,'(A132)')ar__6 read(ar__6,7) + Field_6,ID_6,Name_6,RAh_6,RAm_6,RAs_6,DE__6,DEd_6,DEm_6,DEs_6, + Member_6,SPType_6,Vmag_6,Amp_6,Per_6,VSXinfo_6,VSXPer_6 if (idig(ar__6(104:111)).EQ.0) Per_6 = rNULL__ if (idig(ar__6(123:132)).EQ.0) VSXPer_6 = rNULL__ RAdeg_6 = rNULL__ DEdeg_6 = rNULL__ c Derive coordinates RAdeg_6 and DEdeg_6 from input data c (RAdeg_6 and DEdeg_6 are set to rNULL__ when unknown) if(RAh_6 .GT. -180) RAdeg_6=RAh_6*15. if(RAm_6 .GT. -180) RAdeg_6=RAdeg_6+RAm_6/4. if(RAs_6 .GT. -180) RAdeg_6=RAdeg_6+RAs_6/240. if(DEd_6 .GE. 0) DEdeg_6=DEd_6 if(DEm_6 .GE. 0) DEdeg_6=DEdeg_6+DEm_6/60. if(DEs_6 .GE. 0) DEdeg_6=DEdeg_6+DEs_6/3600. if(DE__6.EQ.'-'.AND.DEdeg_6.GE.0) DEdeg_6=-DEdeg_6 c ..............Just test output........... write(6,7) + Field_6,ID_6,Name_6,RAh_6,RAm_6,RAs_6,DE__6,DEd_6,DEm_6,DEs_6, + Member_6,SPType_6,Vmag_6,Amp_6,Per_6,VSXinfo_6,VSXPer_6 write(6,'(6H Pos: 2F8.4)') RAdeg_6,DEdeg_6 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table8.dat' ! Observed variable Be stars C Format for file interpretation 8 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table8.dat') write(6,*) '....Loading file: table8.dat' do i__=1,62 read(1,'(A132)')ar__7 read(ar__7,8) + Field_7,ID_7,Name_7,RAh_7,RAm_7,RAs_7,DE__7,DEd_7,DEm_7,DEs_7, + Member_7,SPType_7,Vmag_7,Amp_7,Per_7,VSXinfo_7,VSXPer_7 if (idig(ar__7(104:111)).EQ.0) Per_7 = rNULL__ if (idig(ar__7(123:132)).EQ.0) VSXPer_7 = rNULL__ RAdeg_7 = rNULL__ DEdeg_7 = rNULL__ c Derive coordinates RAdeg_7 and DEdeg_7 from input data c (RAdeg_7 and DEdeg_7 are set to rNULL__ when unknown) if(RAh_7 .GT. -180) RAdeg_7=RAh_7*15. if(RAm_7 .GT. -180) RAdeg_7=RAdeg_7+RAm_7/4. if(RAs_7 .GT. -180) RAdeg_7=RAdeg_7+RAs_7/240. if(DEd_7 .GE. 0) DEdeg_7=DEd_7 if(DEm_7 .GE. 0) DEdeg_7=DEdeg_7+DEm_7/60. if(DEs_7 .GE. 0) DEdeg_7=DEdeg_7+DEs_7/3600. if(DE__7.EQ.'-'.AND.DEdeg_7.GE.0) DEdeg_7=-DEdeg_7 c ..............Just test output........... write(6,8) + Field_7,ID_7,Name_7,RAh_7,RAm_7,RAs_7,DE__7,DEd_7,DEm_7,DEs_7, + Member_7,SPType_7,Vmag_7,Amp_7,Per_7,VSXinfo_7,VSXPer_7 write(6,'(6H Pos: 2F8.4)') RAdeg_7,DEdeg_7 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! Observed variable late-type stars C Format for file interpretation 9 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,63 read(1,'(A132)')ar__8 read(ar__8,9) + Field_8,ID_8,Name_8,RAh_8,RAm_8,RAs_8,DE__8,DEd_8,DEm_8,DEs_8, + Member_8,SPType_8,Vmag_8,Amp_8,Per_8,VSXinfo_8,VSXPer_8 if (idig(ar__8(104:111)).EQ.0) Per_8 = rNULL__ if (idig(ar__8(123:132)).EQ.0) VSXPer_8 = rNULL__ RAdeg_8 = rNULL__ DEdeg_8 = rNULL__ c Derive coordinates RAdeg_8 and DEdeg_8 from input data c (RAdeg_8 and DEdeg_8 are set to rNULL__ when unknown) if(RAh_8 .GT. -180) RAdeg_8=RAh_8*15. if(RAm_8 .GT. -180) RAdeg_8=RAdeg_8+RAm_8/4. if(RAs_8 .GT. -180) RAdeg_8=RAdeg_8+RAs_8/240. if(DEd_8 .GE. 0) DEdeg_8=DEd_8 if(DEm_8 .GE. 0) DEdeg_8=DEdeg_8+DEm_8/60. if(DEs_8 .GE. 0) DEdeg_8=DEdeg_8+DEs_8/3600. if(DE__8.EQ.'-'.AND.DEdeg_8.GE.0) DEdeg_8=-DEdeg_8 c ..............Just test output........... write(6,9) + Field_8,ID_8,Name_8,RAh_8,RAm_8,RAs_8,DE__8,DEd_8,DEm_8,DEs_8, + Member_8,SPType_8,Vmag_8,Amp_8,Per_8,VSXinfo_8,VSXPer_8 write(6,'(6H Pos: 2F8.4)') RAdeg_8,DEdeg_8 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table10.dat' ! Observed unclassified variable stars C Format for file interpretation 10 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A10,F10.4) C Effective file loading open(unit=1,status='old',file= +'table10.dat') write(6,*) '....Loading file: table10.dat' do i__=1,80 read(1,'(A132)')ar__9 read(ar__9,10) + Field_9,ID_9,Name_9,RAh_9,RAm_9,RAs_9,DE__9,DEd_9,DEm_9,DEs_9, + Member_9,SPType_9,Vmag_9,Amp_9,Per_9,VSXinfo_9,VSXPer_9 if (idig(ar__9(104:111)).EQ.0) Per_9 = rNULL__ if (idig(ar__9(123:132)).EQ.0) VSXPer_9 = rNULL__ RAdeg_9 = rNULL__ DEdeg_9 = rNULL__ c Derive coordinates RAdeg_9 and DEdeg_9 from input data c (RAdeg_9 and DEdeg_9 are set to rNULL__ when unknown) if(RAh_9 .GT. -180) RAdeg_9=RAh_9*15. if(RAm_9 .GT. -180) RAdeg_9=RAdeg_9+RAm_9/4. if(RAs_9 .GT. -180) RAdeg_9=RAdeg_9+RAs_9/240. if(DEd_9 .GE. 0) DEdeg_9=DEd_9 if(DEm_9 .GE. 0) DEdeg_9=DEdeg_9+DEm_9/60. if(DEs_9 .GE. 0) DEdeg_9=DEdeg_9+DEs_9/3600. if(DE__9.EQ.'-'.AND.DEdeg_9.GE.0) DEdeg_9=-DEdeg_9 c ..............Just test output........... write(6,10) + Field_9,ID_9,Name_9,RAh_9,RAm_9,RAs_9,DE__9,DEd_9,DEm_9,DEs_9, + Member_9,SPType_9,Vmag_9,Amp_9,Per_9,VSXinfo_9,VSXPer_9 write(6,'(6H Pos: 2F8.4)') RAdeg_9,DEdeg_9 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablea11.dat' ! Variable stars not in VSX C Format for file interpretation 11 format( + A11,1X,I4,1X,A24,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X, + A12,1X,A13,1X,F6.3,1X,F5.3,1X,F8.4,1X,A4) C Effective file loading open(unit=1,status='old',file= +'tablea11.dat') write(6,*) '....Loading file: tablea11.dat' do i__=1,221 read(1,'(A116)')ar__10 read(ar__10,11) + Field_10,ID_10,Name_10,RAh_10,RAm_10,RAs_10,DE__10,DEd_10, + DEm_10,DEs_10,Member_10,SPType_10,Vmag_10,Amp_10,Per_10,Type if (idig(ar__10(104:111)).EQ.0) Per_10 = rNULL__ RAdeg_10 = rNULL__ DEdeg_10 = rNULL__ c Derive coordinates RAdeg_10 and DEdeg_10 from input data c (RAdeg_10 and DEdeg_10 are set to rNULL__ when unknown) if(RAh_10 .GT. -180) RAdeg_10=RAh_10*15. if(RAm_10 .GT. -180) RAdeg_10=RAdeg_10+RAm_10/4. if(RAs_10 .GT. -180) RAdeg_10=RAdeg_10+RAs_10/240. if(DEd_10 .GE. 0) DEdeg_10=DEd_10 if(DEm_10 .GE. 0) DEdeg_10=DEdeg_10+DEm_10/60. if(DEs_10 .GE. 0) DEdeg_10=DEdeg_10+DEs_10/3600. if(DE__10.EQ.'-'.AND.DEdeg_10.GE.0) DEdeg_10=-DEdeg_10 c ..............Just test output........... write(6,11) + Field_10,ID_10,Name_10,RAh_10,RAm_10,RAs_10,DE__10,DEd_10, + DEm_10,DEs_10,Member_10,SPType_10,Vmag_10,Amp_10,Per_10,Type write(6,'(6H Pos: 2F8.4)') RAdeg_10,DEdeg_10 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table11.dat' ! Light curves of variable stars C Format for file interpretation 12 format(A11,1X,I4,1X,F11.5,1X,F7.4,1X,F6.4) C Effective file loading open(unit=1,status='old',file= +'table11.dat') write(6,*) '....Loading file: table11.dat' do i__=1,204141 read(1,'(A43)')ar__11 read(ar__11,12)Field_11,ID_11,Obs_date,mag,e_mag c ..............Just test output........... write(6,12)Field_11,ID_11,Obs_date,mag,e_mag c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'lc/*' ! Individual light curves C Format for file interpretation 13 format(F11.5,1X,F7.4,1X,F6.4) C Effective file loading open(unit=1,status='old',file= +'lc/*') write(6,*) '....Loading file: lc/*' do i__=1,354 read(1,'(A26)')ar__12 read(ar__12,13)Obs_date_1,mag_1,e_mag_1 c ..............Just test output........... write(6,13)Obs_date_1,mag_1,e_mag_1 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