Conversion of standardized ReadMe file for
catalog 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/ApJ/709/1306 Spectropolarimetry of 60 Cyg and {pi} Aqr (Wisniewski+, 2010)
*================================================================================
*Disk-loss and disk-renewal phases in classical Be stars.
*I. Analysis of long-term spectropolarimetric data.
* Wisniewski J.P., Draper Z.H., Bjorkman K.S., Meade M.R., Bjorkman J.E.,
* Kowalski A.F.
* <Astrophys. J., 709, 1306-1320 (2010)>
* =2010ApJ...709.1306W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Summary of {pi} Aqr data
integer*4 nr__
parameter (nr__=141) ! Number of records
character*57 ar__ ! Full-size record
character*10 Date ! ("YYYY/MM/DD") UT date of the observation
integer*4 JD ! (d) Julian Date of the observation
character*3 Obs ! Facility code (1)
integer*4 f_Obs ! [1,4]? Flag on Obs (2)
real*4 EW_Ha ! (0.1nm) ? H{alpha} equivalent width; in Angstroms
real*4 e_EW_Ha ! (0.1nm) ? Uncertainty in EW.Ha
real*4 VPol ! (%) ? V band polarization
real*4 e_VPol ! (%) ? Uncertainty in VPol
real*4 PA ! (deg) ? Polarization position angle
character*1 State ! [H/L]? Polarization state (3)
*Note (1): Observatory code as follows:
* PBO = Pine Bluff Observatory;
* RIT = Ritter Observatory.
* McD = HPOL polarimetric observations of {pi} Aqr from literature V-band
* polarization measurements (McDavid, Cat. II/220, or J/PASP/106/949).
*Note (2): Flag on observation (only for table 2) as follows:
* 1 = data are from McDavid (1986PASP...98..572M, Cat. II/220);
* 2 = data are from McDavid (1990PASP..102..773M, Cat. II/220);
* 3 = data are from McDavid (1994, Cat. J/PASP/106/949);
* 4 = data are from McDavid (1999PASP..111..494M, Cat. II/220).
*Note (3): A qualitative assessment of whether 60 Cyg (table 1) or {pi} Aqr
* (table 2) is in a polarimetrically high state (H), indicative of the
* presence of strong disk, or in a polarimetrically low state (L), in
* which evidence of the presence of a disk is absent (see Section 3.3.3
* for further discussion of these designations). Note that we have
* labeled the two observations in 1992 August (in table 1) as being
* "high state", primarily because they are isolated from subsequent
* observations by ~5 years.
C=============================================================================
C Loading file 'table2.dat' ! Summary of {pi} Aqr data
C Format for file interpretation
1 format(
+ A10,1X,I7,1X,A3,1X,I1,1X,F6.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.1,
+ 1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'J/ApJ/709/1306/table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,141
read(1,'(A57)')ar__
read(ar__,1)
+ Date,JD,Obs,f_Obs,EW_Ha,e_EW_Ha,VPol,e_VPol,PA,State
if(ar__(24:24) .EQ. '') f_Obs = iNULL__
if(ar__(26:31) .EQ. '') EW_Ha = rNULL__
if(ar__(33:37) .EQ. '') e_EW_Ha = rNULL__
if(ar__(39:43) .EQ. '') VPol = rNULL__
if(ar__(45:49) .EQ. '') e_VPol = rNULL__
if(ar__(51:55) .EQ. '') PA = rNULL__
c ..............Just test output...........
write(6,1)
+ Date,JD,Obs,f_Obs,EW_Ha,e_EW_Ha,VPol,e_VPol,PA,State
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end