Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/591/A21 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 2026-Jun-08
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/591/A21 Optical flux behaviour of Fermi blazars (Marchesini+, 2016)
*================================================================================
*Optical flux behaviour of a sample of Fermi blazars.
* Marchesini E.J., Andruchow I., Cellone S.A., Combi J.A., Zibecchi L.,
* Marti J., Romero G.E., Munoz-Arjonilla A.J., Luque-Escamilla P.,
* Sanchez-Sutil J.R.
* <Astron. Astrophys. 591, A21 (2016)>
* =2016A&A...591A..21M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sample.dat' ! Sample
integer*4 nr__
parameter (nr__=9) ! Number of records
character*61 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*15 Name ! Blazar name
character*12 Type ! Type of source
integer*4 RAh ! (h) Right ascension (J2000)
integer*4 RAm ! (min) Right ascension (J2000)
real*4 RAs ! (s) Right ascension (J2000)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) Declination (J2000)
integer*4 DEm ! (arcmin) Declination (J2000)
integer*4 DEs ! (arcsec) Declination (J2000)
real*4 Vmag ! (mag) Visual magnitude (as found in NED)
real*4 z ! Redshift (as found in NED)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Statistical results for the differential light
curves
integer*4 nr__1
parameter (nr__1=30) ! Number of records
character*91 ar__1 ! Full-size record
character*15 Name_1 ! Blazar name
character*10 Date ! Observation UT date or "All nights"
real*4 sigmaV ! (mag) V-band observational error (1)
real*4 sigmaR ! (mag) R-band observational error (1)
character*2 VarV ! [V NV] V-band variability flag
character*2 VarR ! [V NV] R-band variability flag
real*8 FV ! V-band confidence parameter from the F-test (2)
real*8 FR ! R-band confidence parameter from the F-test (2)
real*4 GammaV ! V-band gamma corrective factor (3)
real*4 GammaR ! R-band gamma corrective factor (3)
integer*4 NV ! Number of points in the V light curves
integer*4 NR ! Number of points in the R light curves
real*4 FtV ! V-band F^t^ critical value (2)
real*4 FtR ! R-band F^t^ critical value (2)
*Note (1): Observational error, Sigma, obtained from the standard deviation of
* the control-comparison differential light-curve.
*Note (2): For each light curve, we estimated the F parameter, defined as the
* ratio of the variance of the object-comparison light curve ({sigma}^2^_o-c_)
* to the variance of the control-comparison light curve ({sigma}^2^_k-c_).
* If F=({sigma}_o-c_)^2^/({sigma}_k-c_)^2^>=F_n_^{alpha}^, the object is said
* to be variable, being F_n_^{alpha}^ a critical value. This critical value is
* calculated for a set of n=N-1 degrees of freedom, where N is the number of
* points in the light-curve, while {alpha} is chosen to determine the desired
* level of confidence. If {alpha}=0.01, then the F test has a 99%
* confidence level.
*Note (3): gamma corrective factor as described in "Howell S.B., Warnock
* III A. & Mitchell K.J., 1988AJ.....95..247H".
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'phot.dat' ! BV differential photometry for the 9 BL Lacs
integer*4 nr__2
parameter (nr__2=481) ! Number of records
character*71 ar__2 ! Full-size record
character*15 Name_2 ! Blazar name
character*1 Band ! [RV] Observed band
integer*4 Seq ! Sequential number of observations within
* the source and the band
real*8 HJD ! (d) Heliocentric Julian Date
real*8 dHJD ! (d) Differential Heliocentric Julian Date
* (counting from the last integer date)
real*4 dmag ! (mag) Differential magnitude (1)
real*4 e_dmag ! (mag) rms uncertainty on dmag
real*4 dmagC ! (mag) Differential control magnitude (2)
real*4 e_dmagC ! (mag) ? rms uncertainty on dmagC
*Note (1): i.e. the instrumental magnitude of the science object minus the
* instrumental magnitude of the comparison star, both on the same frame.
*Note (2): i.e. the instrumental magnitude of the control object minus the
* instrumental magnitude of the comparison star, both on the same frame.
C=============================================================================
C Loading file 'sample.dat' ! Sample
C Format for file interpretation
1 format(
+ A15,1X,A12,2X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F4.1,
+ 1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'sample.dat')
write(6,*) '....Loading file: sample.dat'
do i__=1,9
read(1,'(A61)')ar__
read(ar__,1)Name,Type,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,z
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)Name,Type,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,z
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Statistical results for the differential light
* curves
C Format for file interpretation
2 format(
+ A15,1X,A10,2X,F5.3,1X,F5.3,1X,A2,1X,A2,2X,F7.2,2X,F7.2,1X,
+ F4.2,1X,F4.2,1X,I2,1X,I2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,30
read(1,'(A91)')ar__1
read(ar__1,2)
+ Name_1,Date,sigmaV,sigmaR,VarV,VarR,FV,FR,GammaV,GammaR,NV,NR,
+ FtV,FtR
c ..............Just test output...........
write(6,2)
+ Name_1,Date,sigmaV,sigmaR,VarV,VarR,FV,FR,GammaV,GammaR,NV,NR,
+ FtV,FtR
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'phot.dat' ! BV differential photometry for the 9 BL Lacs
C Format for file interpretation
3 format(
+ A15,1X,A1,1X,I2,1X,F15.7,1X,F8.6,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'phot.dat')
write(6,*) '....Loading file: phot.dat'
do i__=1,481
read(1,'(A71)')ar__2
read(ar__2,3)
+ Name_2,Band,Seq,HJD,dHJD,dmag,e_dmag,dmagC,e_dmagC
if(ar__2(67:71) .EQ. '') e_dmagC = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,Band,Seq,HJD,dHJD,dmag,e_dmag,dmagC,e_dmagC
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end