Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/883/94 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-Sep-10
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/883/94 Optical and UV follow-up of 1ES 1927+654 (Trakhtenbrot+, 2019)
*================================================================================
*1ES 1927+654: an AGN caught changing look on a timescale of months.
* Trakhtenbrot B., Arcavi I., MacLeod C.L., Ricci C., Kara E., Graham M.L.,
* Stern D., Harrison F.A., Burke J., Hiramatsu D., Hosseinzadeh G.,
* Howell D.A., Smartt S.J., Rest A., Prieto J.L., Shappee B.J.,
* Holoien T.W.-S., Bersier D., Filippenko A.V., Brink T.G., Zheng W., Li R.,
* Remillard R.A., Loewenstein M.
* <Astrophys. J., 883, 94 (2019)>
* =2019ApJ...883...94T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Log of spectroscopic observations
integer*4 nr__
parameter (nr__=28) ! Number of records
character*61 ar__ ! Full-size record
character*10 Date ! ("Y/M/D") Date of observation (UT) at exposure start
character*1 f_Date ! [*c] Flag on Date (1)
real*8 MJD ! (d) [58183.5/58622.6] Modified Julian Date
* at exposure start
integer*4 delT ! (d) [72/503] Time since ATLAS first detection of
* the optical flare (MJD=58110.2),
* in rounded rest-frame days
character*22 Tel ! Telescope
character*6 Inst ! Instrument
character*1 u_Exp ! [~] Uncertainty flag on Exp
integer*4 Exp ! (s) [300/3600] Exposure time
character*1 f_Exp ! [d] Flag on Exp (1)
*Note (1): Flag as follows:
* * = spectra taken under different weather and/or instrument conditions
* compared to a standard star, thus affecting the validity of their
* continuum shape.
* c = A low-resolution spectrum, not shown in Figure 3.
* d = The Shane/Kast exposure times were 3x500s for the red arm, and
* 1x1560s for the blue.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fig1.dat' ! Long-term optical light curve of 1ES1927+654
integer*4 nr__1
parameter (nr__1=714) ! Number of records
character*56 ar__1 ! Full-size record
real*8 MJD_1 ! (d) [56734.5/58630.6] Modified Julian Date
* of Observation (JD-2400000.5)
real*8 Phase ! (d) [-1349.42/510.5] Rest-frame days since
* 2017 December 23 (MJD=58110.2)
character*1 l_mag ! [l] 3-sigma non-detection limit flag on mag
real*4 mag ! (mag) [14.8/19.5] Apparent magnitude in Filt
real*4 e_mag ! (mag) [0.014/0.4]?=0 Uncertainty in mag (G1)
character*7 Tel_1 ! Telescope (ASAS-SN or ATLAS)
character*1 Filt ! Filter (g, V or c, o)
character*11 q_mag ! Quality for mag: "bad weather"
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fig2.dat' ! Light curve and equivalent-width evolution of
1ES 1927+654 during the flare
integer*4 nr__2
parameter (nr__2=2909) ! Number of records
character*104 ar__2 ! Full-size record
real*8 MJD_2 ! (d) [56734.5/58630.6] Modified Julian Date
* of Observation (JD-2400000.5)
real*8 Phase_1 ! (d) [-1349.42/510.5] Phase in days since
* 2017 December 23 (MJD=58110.2)
real*8 EW ! (0.1nm) [0.014/205]? Equivalent width in Angstrom units
real*4 e_EW ! (0.1nm) [0.014/27.2]? Lower uncertainty on EW
real*4 E_EW_1 ! (0.1nm) [0.2/13.4]? Upper uncertainty on EW
character*8 Line ! H{alpha} or H{beta}
real*8 Fx ! [0.3/1298]? 0.5-10keV X-ray flux
real*4 e_Fx ! [0/24]? Lower uncertainty on Fx
real*4 E_Fx_1 ! [0.1/16.5]? Upper uncertainty on Fx
character*10 Tel_2 ! Telescope (ASAS-SN, ATLAS, Swift-UVOT,
* Swift-XRT, NICER or XMM)
character*1 l_mag_1 ! [l] 3-sigma non-detection limit flag on mag
real*4 mag_1 ! (mag) [13.89/19.5]? Apparent magnitude in Filt
real*4 e_mag_1 ! (mag) [0.014/0.4]? Uncertainty on mag (G1)
character*4 Filt_1 ! Filter (g, V; c, o or U, UVM2, UVW1, UVW2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'fig3/*' ! Optical spectra of 1ES 1927+654 in ASCII format
for 26 observation dates
integer*4 nr__3
parameter (nr__3=26) ! Number of records
character*20 ar__3 ! Full-size record
real*8 lambda ! (0.1nm) [3198.2/10706] Observed wavelength
real*4 Flux ! (cW/m2/nm) [1.5e-17/7e-14] Flux density,
* in erg/cm^2^/s/{AA} units
C=============================================================================
C Loading file 'table1.dat' ! Log of spectroscopic observations
C Format for file interpretation
1 format(A10,1X,A1,1X,F7.1,1X,I3,1X,A22,1X,A6,1X,A1,I4,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,28
read(1,'(A61)')ar__
read(ar__,1)Date,f_Date,MJD,delT,Tel,Inst,u_Exp,Exp,f_Exp
c ..............Just test output...........
write(6,1)Date,f_Date,MJD,delT,Tel,Inst,u_Exp,Exp,f_Exp
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fig1.dat' ! Long-term optical light curve of 1ES1927+654
C Format for file interpretation
2 format(F9.3,1X,F9.3,1X,A1,1X,F6.3,1X,F5.3,1X,A7,1X,A1,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'fig1.dat')
write(6,*) '....Loading file: fig1.dat'
do i__=1,714
read(1,'(A56)')ar__1
read(ar__1,2)MJD_1,Phase,l_mag,mag,e_mag,Tel_1,Filt,q_mag
c ..............Just test output...........
write(6,2)MJD_1,Phase,l_mag,mag,e_mag,Tel_1,Filt,q_mag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fig2.dat' ! Light curve and equivalent-width evolution of
* 1ES 1927+654 during the flare
C Format for file interpretation
3 format(
+ F9.3,1X,F9.3,1X,F7.3,1X,F6.3,1X,F6.3,1X,A8,1X,F8.3,1X,F6.3,1X,
+ F6.3,1X,A10,1X,A1,1X,F6.3,1X,F5.3,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'fig2.dat')
write(6,*) '....Loading file: fig2.dat'
do i__=1,2909
read(1,'(A104)')ar__2
read(ar__2,3)
+ MJD_2,Phase_1,EW,e_EW,E_EW_1,Line,Fx,e_Fx,E_Fx_1,Tel_2,
+ l_mag_1,mag_1,e_mag_1,Filt_1
if(ar__2(21:27) .EQ. '') EW = rNULL__
if(ar__2(29:34) .EQ. '') e_EW = rNULL__
if(ar__2(36:41) .EQ. '') E_EW_1 = rNULL__
if(ar__2(52:59) .EQ. '') Fx = rNULL__
if(ar__2(61:66) .EQ. '') e_Fx = rNULL__
if(ar__2(68:73) .EQ. '') E_Fx_1 = rNULL__
if(ar__2(88:93) .EQ. '') mag_1 = rNULL__
if(ar__2(95:99) .EQ. '') e_mag_1 = rNULL__
c ..............Just test output...........
write(6,3)
+ MJD_2,Phase_1,EW,e_EW,E_EW_1,Line,Fx,e_Fx,E_Fx_1,Tel_2,
+ l_mag_1,mag_1,e_mag_1,Filt_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'fig3/*' ! Optical spectra of 1ES 1927+654 in ASCII format
* for 26 observation dates
C Format for file interpretation
4 format(F8.2,1X,E11.4)
C Effective file loading
open(unit=1,status='old',file=
+'fig3/*')
write(6,*) '....Loading file: fig3/*'
do i__=1,26
read(1,'(A20)')ar__3
read(ar__3,4)lambda,Flux
c ..............Just test output...........
write(6,4)lambda,Flux
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end