Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/883/94 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.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 (nr__) ! ("Y/M/D") Date of observation (UT) at exposure start
character*1 f_Date (nr__) ! [*c] Flag on Date (1)
real*8 MJD (nr__) ! (d) [58183.5/58622.6] Modified Julian Date
* at exposure start
integer*4 delT (nr__) ! (d) [72/503] Time since ATLAS first detection of
* the optical flare (MJD=58110.2),
* in rounded rest-frame days
character*22 Tel (nr__) ! Telescope
character*6 Inst (nr__) ! Instrument
character*1 u_Exp (nr__) ! [~] Uncertainty flag on Exp
integer*4 Exp (nr__) ! (s) [300/3600] Exposure time
character*1 f_Exp (nr__) ! [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 (nr__1) ! (d) [56734.5/58630.6] Modified Julian Date
* of Observation (JD-2400000.5)
real*8 Phase (nr__1) ! (d) [-1349.42/510.5] Rest-frame days since
* 2017 December 23 (MJD=58110.2)
character*1 l_mag (nr__1) ! [l] 3-sigma non-detection limit flag on mag
real*4 mag (nr__1) ! (mag) [14.8/19.5] Apparent magnitude in Filt
real*4 e_mag (nr__1) ! (mag) [0.014/0.4]?=0 Uncertainty in mag (G1)
character*7 Tel_1 (nr__1) ! Telescope (ASAS-SN or ATLAS)
character*1 Filt (nr__1) ! Filter (g, V or c, o)
character*11 q_mag (nr__1) ! 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 (nr__2) ! (d) [56734.5/58630.6] Modified Julian Date
* of Observation (JD-2400000.5)
real*8 Phase_1 (nr__2) ! (d) [-1349.42/510.5] Phase in days since
* 2017 December 23 (MJD=58110.2)
real*8 EW (nr__2) ! (0.1nm) [0.014/205]? Equivalent width in Angstrom units
real*4 e_EW (nr__2) ! (0.1nm) [0.014/27.2]? Lower uncertainty on EW
real*4 E_EW_1 (nr__2) ! (0.1nm) [0.2/13.4]? Upper uncertainty on EW
character*8 Line (nr__2) ! H{alpha} or H{beta}
real*8 Fx (nr__2) ! [0.3/1298]? 0.5-10keV X-ray flux
real*4 e_Fx (nr__2) ! [0/24]? Lower uncertainty on Fx
real*4 E_Fx_1 (nr__2) ! [0.1/16.5]? Upper uncertainty on Fx
character*10 Tel_2 (nr__2) ! Telescope (ASAS-SN, ATLAS, Swift-UVOT,
* Swift-XRT, NICER or XMM)
character*1 l_mag_1 (nr__2) ! [l] 3-sigma non-detection limit flag on mag
real*4 mag_1 (nr__2) ! (mag) [13.89/19.5]? Apparent magnitude in Filt
real*4 e_mag_1 (nr__2) ! (mag) [0.014/0.4]? Uncertainty on mag (G1)
character*4 Filt_1 (nr__2) ! 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 (nr__3) ! (0.1nm) [3198.2/10706] Observed wavelength
real*4 Flux (nr__3) ! (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(i__),f_Date(i__),MJD(i__),delT(i__),Tel(i__),Inst(i__),
+ u_Exp(i__),Exp(i__),f_Exp(i__)
c ..............Just test output...........
write(6,1)
+ Date(i__),f_Date(i__),MJD(i__),delT(i__),Tel(i__),Inst(i__),
+ u_Exp(i__),Exp(i__),f_Exp(i__)
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(i__),Phase(i__),l_mag(i__),mag(i__),e_mag(i__),
+ Tel_1(i__),Filt(i__),q_mag(i__)
c ..............Just test output...........
write(6,2)
+ MJD_1(i__),Phase(i__),l_mag(i__),mag(i__),e_mag(i__),
+ Tel_1(i__),Filt(i__),q_mag(i__)
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(i__),Phase_1(i__),EW(i__),e_EW(i__),E_EW_1(i__),
+ Line(i__),Fx(i__),e_Fx(i__),E_Fx_1(i__),Tel_2(i__),
+ l_mag_1(i__),mag_1(i__),e_mag_1(i__),Filt_1(i__)
if(ar__2(21:27) .EQ. '') EW(i__) = rNULL__
if(ar__2(29:34) .EQ. '') e_EW(i__) = rNULL__
if(ar__2(36:41) .EQ. '') E_EW_1(i__) = rNULL__
if(ar__2(52:59) .EQ. '') Fx(i__) = rNULL__
if(ar__2(61:66) .EQ. '') e_Fx(i__) = rNULL__
if(ar__2(68:73) .EQ. '') E_Fx_1(i__) = rNULL__
if(ar__2(88:93) .EQ. '') mag_1(i__) = rNULL__
if(ar__2(95:99) .EQ. '') e_mag_1(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ MJD_2(i__),Phase_1(i__),EW(i__),e_EW(i__),E_EW_1(i__),
+ Line(i__),Fx(i__),e_Fx(i__),E_Fx_1(i__),Tel_2(i__),
+ l_mag_1(i__),mag_1(i__),e_mag_1(i__),Filt_1(i__)
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(i__),Flux(i__)
c ..............Just test output...........
write(6,4)lambda(i__),Flux(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end