Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/458/245 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-Aug-11
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/458/245 Obs. of 2 Unidentified Gamma-ray Sources (La Palombara+, 2006)
*================================================================================
*XMM-Newton and ESO observations of the two Unidentified Gamma-ray Sources
*3EG J0616-3310 and 3EG J1249-8330.
* La Palombara N., Mignani R.P., Hatziminaoglou E., Schirmer M., Bignami G.F.,
* Caraveo P.
* <Astron. Astrophys. 458, 245 (2006)>
* =2006A&A...458..245L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Main characteristics of the 8 EPIC observations of
3EG J0616-3310 (Obs=1-4) & 3EG J1249-8330 (Obs=5-8)
integer*4 nr__
parameter (nr__=8) ! Number of records
character*70 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Obs (nr__) ! EPIC field sequential reference number
integer*4 Rev (nr__) ! Revolution number of XMM
character*10 ObsDate (nr__) ! ("YYYY-MM-DD") Observation date
character*8 ObsTime (nr__) ! ("h:m:s") Observation time (hh:mm:ss)
integer*4 RAh (nr__) ! (h) Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Right ascension (J2000)
real*4 RAs (nr__) ! (s) Right ascension (J2000)
character*1 DE_ (nr__) ! Declination sign (J2000)
integer*4 DEd (nr__) ! (deg) Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Declination (J2000)
real*4 PN (nr__) ! (ks) Exposure time for PN instrument
real*4 MOS1 (nr__) ! (ks) Exposure time for MOS1 instrument
real*4 MOS2 (nr__) ! (ks) Exposure time for MOS2 instrument
real*4 NH (nr__) ! (10+20cm-2) Hydrogen column density
integer*4 Det (nr__) ! Number of detected sources
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Main characteristics of the X-ray sources with no
candidate optical counterpart and log(FX/Fopt)>~2
integer*4 nr__1
parameter (nr__1=9) ! Number of records
character*96 ar__1 ! Full-size record
integer*4 Obs_1 (nr__1) ! EPIC field sequential reference number
character*16 XMMU (nr__1) ! XMMU designation (JHHMMSS.s+DDMMSS)
character*23 Bands (nr__1) ! Detection bands
character*2 l_kTbt (nr__1) ! [<>= ] Limit flag on kTbt
real*4 kTbt (nr__1) ! (keV) Range of possible kT temperature for
* thermal bremsstrahlung spectrum
integer*4 kTbt2 (nr__1) ! (keV) ? Range of possible kT temperature for
* thermal bremsstrahlung spectrum
character*2 l_kTbb (nr__1) ! [<>= ] Limit flag on kTbb
real*4 kTbb (nr__1) ! (keV) Range of possible kT temperature for
* thermal blackbody spectrum
real*4 kTbb2 (nr__1) ! (keV) ? Range of possible kT temperature for
* thermal blackbody spectrum
character*2 l_Gamma (nr__1) ! [<>= ] Limit flag on Gamma
real*4 Gamma (nr__1) ! (keV) Range of photon index for a power law
* spectrum
real*4 Gamma2 (nr__1) ! (keV) ? Range of photon index for a power law
* spectrum
real*4 Flux (nr__1) ! (10-17W/m2) Source flux in 0.3-10keV energy range
real*4 e_Flux (nr__1) ! (10-17W/m2) rms uncertainty on Flux
real*4 log_FX_Fopt(nr__1) ! ([-]) Lower limit of X-ray-to-optical flux ratio
real*4 e_log_FX_Fopt(nr__1) ! ([-]) rms uncertainty on log(FX/Fopt)
C=============================================================================
C Loading file 'table1.dat' ! Main characteristics of the 8 EPIC observations of
* 3EG J0616-3310 (Obs=1-4) & 3EG J1249-8330 (Obs=5-8)
C Format for file interpretation
1 format(
+ I1,1X,I3,1X,A10,1X,A8,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,
+ F4.1,1X,F3.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,8
read(1,'(A70)')ar__
read(ar__,1)
+ Obs(i__),Rev(i__),ObsDate(i__),ObsTime(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),PN(i__),
+ MOS1(i__),MOS2(i__),NH(i__),Det(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Obs(i__),Rev(i__),ObsDate(i__),ObsTime(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),PN(i__),
+ MOS1(i__),MOS2(i__),NH(i__),Det(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Main characteristics of the X-ray sources with no
* candidate optical counterpart and log(FX/Fopt)>~2
C Format for file interpretation
2 format(
+ I1,6X,A16,1X,A23,1X,A2,F3.1,1X,I1,1X,A2,F3.1,1X,F3.1,1X,A2,
+ F3.1,1X,F3.1,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,9
read(1,'(A96)')ar__1
read(ar__1,2)
+ Obs_1(i__),XMMU(i__),Bands(i__),l_kTbt(i__),kTbt(i__),
+ kTbt2(i__),l_kTbb(i__),kTbb(i__),kTbb2(i__),l_Gamma(i__),
+ Gamma(i__),Gamma2(i__),Flux(i__),e_Flux(i__),log_FX_Fopt(i__),
+ e_log_FX_Fopt(i__)
if(ar__1(55:55) .EQ. '') kTbt2(i__) = iNULL__
if(ar__1(63:65) .EQ. '') kTbb2(i__) = rNULL__
if(ar__1(73:75) .EQ. '') Gamma2(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Obs_1(i__),XMMU(i__),Bands(i__),l_kTbt(i__),kTbt(i__),
+ kTbt2(i__),l_kTbb(i__),kTbb(i__),kTbb2(i__),l_Gamma(i__),
+ Gamma(i__),Gamma2(i__),Flux(i__),e_Flux(i__),log_FX_Fopt(i__),
+ e_log_FX_Fopt(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end