Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/695/A228 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-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/695/A228 Catalog of dust-echo-like infrared flares (Necker+, 2025)
*================================================================================
*Flaires. A comprehensive catalog of dust-echo-like infrared flares.
* Necker J., Graikou E., Kowalski M., Franckowiak A., Nordin J., Pernice T.,
* van Velzen S., Veres P.M.
* <Astron. Astrophys. 695, A228 (2025)>
* =2025A&A...695A.228N (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Catalog of dust-echo-like infrared flares
integer*4 nr__
parameter (nr__=823) ! Number of records
character*674 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*25 Name (nr__) ! Host Name
real*8 RAdeg (nr__) ! (deg) Right Ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
real*4 HW1mag (nr__) ! ? Host W1 magnitude (Vega)
real*4 HW2mag (nr__) ! ? Host W2 magnitude (Vega)
integer*4 NEWS (nr__) ! ? Index in the NEWS sample
character*20 AllWISE (nr__) ! AllWISE Designation of the Host,
* JHHMMSS.ss+DDMMSS.s
integer*8 PS1 (nr__) ! ? Object identifier of the host
* in Pan-STARRS
integer*8 AllWISEcntr(nr__) ! ? Counter of the host in WISE AllSky
character*24 AllWISEDes (nr__) ! ? Designation of the host in WISE AllSky,
* WISE JHHMMSS.ss+DDMMSS.s
integer*4 NEDLVSIndex(nr__) ! ? Index of the host in NED-LVS
character*25 NEDLVSName (nr__) ! ? Name of the host in NED-LVS
integer*4 ParentSamp (nr__) ! Index in the Flaires parent galaxy sample
real*8 z (nr__) ! ? Redshift
real*8 e_z (nr__) ! ? Uncertainty of the redshift
character*59 r_z (nr__) ! ? The source where the redshift value
* was taken from
real*4 SDSSdist (nr__) ! (arcsec) ? Distance to the associated SDSS object
character*20 SDSSclass (nr__) ! ? BPT class of the associated SDSS object
real*4 TNSdist (nr__) ! (arcsec) ? Distance to the associated TNS object
character*7 TNSobjtype (nr__) ! ? Type of the associated TNS object
character*11 TNSname (nr__) ! ? Name of the associated TNS object
real*8 TNSdate (nr__) ! ? Discovery date [MJD] of the associated
* TNS object
real*8 milliqdist (nr__) ! (arcsec) ? Distance to associated Milliquas object
character*3 milliqtype (nr__) ! ? Broad Type of the associated Milliquas
* object
character*23 MirongName (nr__) ! Name of the host in the MIRONG sample,
* SDSSJHHMMSS.ss+DDMMSS.s
character*11 WTPName (nr__) ! Name of the object in the WISE Transient
* Pipeline TDE sample, WTPNNaaaaaa
real*8 RefTime (nr__) ! (d) Observer frame mean peak time in MJD
real*8 x2W1 (nr__) ! Chi2 value w.r.t. the median spectral flux
* density in W1
integer*4 npointsW1 (nr__) ! [9/21] Number of detections in W1
real*4 FmedW1 (nr__) ! (mJy) Median spectral flux density in W1
real*8 x2W2 (nr__) ! Chi2 value w.r.t. the median spectral flux
* density in W2
integer*4 npointsW2 (nr__) ! [9/21] Number of Detections in W2
real*4 FmedW2 (nr__) ! (mJy) Median spectral flux density in W2
real*4 FbslW1 (nr__) ! (mJy) Measured baseline spectral flux density
* in W1
real*4 e_FbslW1 (nr__) ! (mJy) 1 sigma uncertainty on FbslW1
real*4 FbslW2 (nr__) ! (mJy) Measured baseline spectral flux density
* in W2
real*4 e_FbslW2 (nr__) ! (mJy) 1 sigma uncertainty on FbslW2
real*8 startW1 (nr__) ! (d) MJD of the first epoch of the determined
* excess in W1
real*8 endW1 (nr__) ! (d) MJD of the last epoch of the determined
* excess in W1
integer*4 endedW1 (nr__) ! [0/1] True (1) if the last epoch of the
* excess in W1 is the last available datapoint
real*8 startW2 (nr__) ! (d) MJD of the first epoch of the determined
* excess in W2
real*8 endW2 (nr__) ! (d) MJD of the last epoch of the determined
* excess in W2
integer*4 endedW2 (nr__) ! [0/1] True (1) if the last epoch of the
* excess in W2 is the last available datapoint
real*4 strengthW1 (nr__) ! Dust Echo strength in W1
real*4 strengthW2 (nr__) ! Dust Echo strength in W2
real*4 varW1 (nr__) ! Significance of the extraneous variability
* in W1
real*4 varW2 (nr__) ! Significance of the extraneous variability
* in W2
real*4 MaxFluxW1 (nr__) ! (mW/m2) Maxium flux density in W1
real*4 e_MaxFluxW1(nr__) ! (mW/m2) 1 sigma uncertainty on MaxFluxW1
real*4 MaxFluxW2 (nr__) ! (mW/m2) Maxium flux density in W2
real*4 e_MaxFluxW2(nr__) ! (mW/m2) 1 sigma uncertainty on MaxFluxW2
real*4 FluenceW1 (nr__) ! (mJ/m2) Integrated flux density in W1
real*4 e_FluenceW1(nr__) ! (mJ/m2) Uncertainty on FluenceW1 estimated by
* performing the integration with
* +/- e_MaxFluxW1
real*4 FluenceW2 (nr__) ! (mJ/m2) Integrated flux density in W1
real*4 e_FluenceW2(nr__) ! (mJ/m2) Uncertainty on FluenceW2 estimated by
* performing the integration with
* +/- e_MaxFluxW2
real*4 Sep (nr__) ! (arcsec) Separation of flare data from baseline data
real*4 PeakLbol (nr__) ! (10-7W) ? Peak of the bolometric luminosity
real*8 PeakTime (nr__) ! (d) ? Rest frame time of the peak of the
* bolometric luminosity relative to RefTime
real*4 Ebol (nr__) ! (10-7J) ? Total emitted bolometric energy
real*4 Fluencebol (nr__) ! (mJ/m2) ? Total bolometric fluence as
* Ebol/4{pi}lumdist^2^
C=============================================================================
C Loading file 'tablea1.dat' ! Catalog of dust-echo-like infrared flares
C Format for file interpretation
1 format(
+ A25,1X,F9.5,1X,F9.5,1X,F6.3,1X,F6.3,1X,I8,1X,A20,1X,I18,3X,
+ I19,2X,A24,1X,I7,1X,A25,1X,I8,1X,F9.4,1X,F9.4,1X,A59,1X,F6.3,
+ 1X,A20,1X,F6.3,1X,A7,1X,A11,1X,F9.2,1X,F9.2,1X,A3,1X,A23,1X,
+ A11,1X,F8.2,1X,F9.2,1X,I2,1X,F6.3,1X,F9.2,1X,I2,1X,F6.3,1X,
+ F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F8.2,1X,F8.2,1X,I1,1X,F8.2,1X,
+ F8.2,1X,I1,1X,F6.2,1X,F6.2,1X,F5.2,1X,F5.2,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,F6.2,
+ 1X,E10.4,1X,F8.2,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,823
read(1,'(A674)')ar__
read(ar__,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),HW1mag(i__),HW2mag(i__),
+ NEWS(i__),AllWISE(i__),PS1(i__),AllWISEcntr(i__),
+ AllWISEDes(i__),NEDLVSIndex(i__),NEDLVSName(i__),
+ ParentSamp(i__),z(i__),e_z(i__),r_z(i__),SDSSdist(i__),
+ SDSSclass(i__),TNSdist(i__),TNSobjtype(i__),TNSname(i__),
+ TNSdate(i__),milliqdist(i__),milliqtype(i__),MirongName(i__),
+ WTPName(i__),RefTime(i__),x2W1(i__),npointsW1(i__),
+ FmedW1(i__),x2W2(i__),npointsW2(i__),FmedW2(i__),FbslW1(i__),
+ e_FbslW1(i__),FbslW2(i__),e_FbslW2(i__),startW1(i__),
+ endW1(i__),endedW1(i__),startW2(i__),endW2(i__),endedW2(i__),
+ strengthW1(i__),strengthW2(i__),varW1(i__),varW2(i__),
+ MaxFluxW1(i__),e_MaxFluxW1(i__),MaxFluxW2(i__),
+ e_MaxFluxW2(i__),FluenceW1(i__),e_FluenceW1(i__),
+ FluenceW2(i__),e_FluenceW2(i__),Sep(i__),PeakLbol(i__),
+ PeakTime(i__),Ebol(i__),Fluencebol(i__)
if(ar__(47:52) .EQ. '') HW1mag(i__) = rNULL__
if(ar__(54:59) .EQ. '') HW2mag(i__) = rNULL__
if(ar__(61:68) .EQ. '') NEWS(i__) = iNULL__
if(ar__(91:108) .EQ. '') PS1(i__) = iNULL__
if(ar__(112:130) .EQ. '') AllWISEcntr(i__) = iNULL__
if(ar__(158:164) .EQ. '') NEDLVSIndex(i__) = iNULL__
if(ar__(201:209) .EQ. '') z(i__) = rNULL__
if(ar__(211:219) .EQ. '') e_z(i__) = rNULL__
if(ar__(281:286) .EQ. '') SDSSdist(i__) = rNULL__
if(ar__(309:314) .EQ. '') TNSdist(i__) = rNULL__
if(ar__(336:344) .EQ. '') TNSdate(i__) = rNULL__
if(ar__(346:354) .EQ. '') milliqdist(i__) = rNULL__
if(ar__(634:643) .EQ. '') PeakLbol(i__) = rNULL__
if(ar__(645:652) .EQ. '') PeakTime(i__) = rNULL__
if(ar__(654:663) .EQ. '') Ebol(i__) = rNULL__
if(ar__(665:674) .EQ. '') Fluencebol(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),HW1mag(i__),HW2mag(i__),
+ NEWS(i__),AllWISE(i__),PS1(i__),AllWISEcntr(i__),
+ AllWISEDes(i__),NEDLVSIndex(i__),NEDLVSName(i__),
+ ParentSamp(i__),z(i__),e_z(i__),r_z(i__),SDSSdist(i__),
+ SDSSclass(i__),TNSdist(i__),TNSobjtype(i__),TNSname(i__),
+ TNSdate(i__),milliqdist(i__),milliqtype(i__),MirongName(i__),
+ WTPName(i__),RefTime(i__),x2W1(i__),npointsW1(i__),
+ FmedW1(i__),x2W2(i__),npointsW2(i__),FmedW2(i__),FbslW1(i__),
+ e_FbslW1(i__),FbslW2(i__),e_FbslW2(i__),startW1(i__),
+ endW1(i__),endedW1(i__),startW2(i__),endW2(i__),endedW2(i__),
+ strengthW1(i__),strengthW2(i__),varW1(i__),varW2(i__),
+ MaxFluxW1(i__),e_MaxFluxW1(i__),MaxFluxW2(i__),
+ e_MaxFluxW2(i__),FluenceW1(i__),e_FluenceW1(i__),
+ FluenceW2(i__),e_FluenceW2(i__),Sep(i__),PeakLbol(i__),
+ PeakTime(i__),Ebol(i__),Fluencebol(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end