Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/371/79 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.5, on 2013-May-26
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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/A+A/371/79 NGC 5548 UBVRI phot. and Halpha, Hbeta fluxes (Dietrich+, 2001)
*================================================================================
*A spectroscopic and photometric study of short-timescale variability
*in NGC 5548
* Dietrich M., Bender C.F., Bergmann D.J., Bills T.E., Bochkarev N.G.,
* Burenkov A.N., Gaskell C.M., Gutzmer D.D., Grove R., Hiller M.E.,
* Huchra J.P., Klimek E.S., Lund C., Merkulova N.I., Pebley S., Poulsen M.A.,
* Pronik V.I., Sergeev S.G., Sergeev E.A., Shapovalova A.I., Vlasyuk V.V.,
* Wilkes B.
* <Astron. Astrophys. 371, 79 (2001)>
* =2001A&A...371...79D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1a.dat' ! log of the photometric observations. The flux of
the broad band fluxes is given in units
of [10^-15erg/s/cm^2^/{AA}]
integer*4 nr__
parameter (nr__=178) ! Number of records
character*81 ar__ ! Full-size record
character*10 Date (nr__) ! ("YYYY/MM/DD") Observation date
real*8 HJD (nr__) ! (d) Heliocentric Julian date of observation
character*1 Code (nr__) ! [DQ] Code (G1)
real*4 Umag (nr__) ! (10-17W/m2/nm) ? U magnitude
real*4 e_Umag (nr__) ! (10-17W/m2/nm) ? rms uncertainty on Umag
real*4 Bmag (nr__) ! (10-17W/m2/nm) ? B magnitude
real*4 e_Bmag (nr__) ! (10-17W/m2/nm) ? rms uncertainty on Bmag
real*4 Vmag (nr__) ! (10-17W/m2/nm) V magnitude
real*4 e_Vmag (nr__) ! (10-17W/m2/nm) rms uncertainty on Vmag
real*4 Rmag (nr__) ! (10-17W/m2/nm) ? R magnitude
real*4 e_Rmag (nr__) ! (10-17W/m2/nm) ? rms uncertainty on Rmag
real*4 Imag (nr__) ! (10-17W/m2/nm) ? I magnitude
real*4 e_Imag (nr__) ! (10-17W/m2/nm) ? rms uncertainty on Imag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2a.dat' ! log of the spectroscopic observations
integer*4 nr__1
parameter (nr__1=84) ! Number of records
character*83 ar__1 ! Full-size record
character*10 Date_1 (nr__1) ! ("YYYY/MM/DD") Observation date
real*8 Middle (nr__1) ! (d) Julian date of middle of observation
character*2 Code_1 (nr__1) ! Code (G1)
integer*4 Tint (nr__1) ! (s) Integration time
integer*4 linf (nr__1) ! (0.1nm) Lower limit of spectroscopic range
integer*4 lsup (nr__1) ! (0.1nm) Upper limit of spectroscopic range
real*4 Aper1 (nr__1) ! (arcsec) Aperture
real*4 Aper2 (nr__1) ! (arcsec) Aperture
integer*4 PA (nr__1) ! (deg) Position angle
real*4 Seeing (nr__1) ! (arcsec) Seeing
real*4 Res (nr__1) ! (0.1nm) ? Resolution
integer*4 Res2 (nr__1) ! (0.1nm) ? Upper limit of resolution when interval
character*10 Filename (nr__1) ! Filename of the observation
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3a.dat' ! Emission line and continuum flux measurements.
integer*4 nr__2
parameter (nr__2=56) ! Number of records
character*58 ar__2 ! Full-size record
character*10 Date_2 (nr__2) ! ("YYYY/MM/DD") Observation date
real*8 HJD_1 (nr__2) ! (d) Heliocentric Julian date of observation
character*2 Code_2 (nr__2) ! Code (G1)
real*4 Fl5100 (nr__2) ! (10-17W/m2/nm) ? Continuum flux (1)
real*4 e_Fl5100 (nr__2) ! (10-17W/m2/nm) ? rms uncertainty on Fl5100 (1)
real*4 FHbeta (nr__2) ! (10-16W/m2) ? H{beta} emission line flux (2)
real*4 e_FHbeta (nr__2) ! (10-16W/m2) ? rms uncertainty on FHbeta (2)
real*4 FHalpha (nr__2) ! (10-16W/m2) ? H{alpha} emission line flux (2)
real*4 e_FHalpha (nr__2) ! (10-16W/m2) ? rms uncertainty on FHalpha (2)
*Note (1): Wavelength range 5185{AA}-5195{AA}. In units of 10^-15erg/s/cm^2^/{AA}
*Note (2): In units of 10^-13erg/s/cm^2^
C=============================================================================
C Loading file 'table1a.dat' ! log of the photometric observations. The flux of
* the broad band fluxes is given in units
* of [10^-15erg/s/cm^2^/{AA}]
C Format for file interpretation
1 format(
+ A10,1X,F12.4,1X,A1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,file='table1a.dat', status='old')
write(6,*) '....Loading file: table1a.dat'
do i__=1,178
read(1,'(A81)')ar__
read(ar__,1)
+ Date(i__),HJD(i__),Code(i__),Umag(i__),e_Umag(i__),Bmag(i__),
+ e_Bmag(i__),Vmag(i__),e_Vmag(i__),Rmag(i__),e_Rmag(i__),
+ Imag(i__),e_Imag(i__)
if(ar__(27:31) .EQ. '') Umag(i__) = rNULL__
if(ar__(33:36) .EQ. '') e_Umag(i__) = rNULL__
if(ar__(38:42) .EQ. '') Bmag(i__) = rNULL__
if(ar__(44:47) .EQ. '') e_Bmag(i__) = rNULL__
if(ar__(60:64) .EQ. '') Rmag(i__) = rNULL__
if(ar__(66:69) .EQ. '') e_Rmag(i__) = rNULL__
if(ar__(71:75) .EQ. '') Imag(i__) = rNULL__
if(ar__(77:80) .EQ. '') e_Imag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Date(i__),HJD(i__),Code(i__),Umag(i__),e_Umag(i__),Bmag(i__),
+ e_Bmag(i__),Vmag(i__),e_Vmag(i__),Rmag(i__),e_Rmag(i__),
+ Imag(i__),e_Imag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2a.dat' ! log of the spectroscopic observations
C Format for file interpretation
2 format(
+ A10,1X,F13.5,1X,A2,1X,I4,1X,I4,1X,I4,1X,F4.2,3X,F5.2,1X,I3,2X,
+ F3.1,1X,F4.2,1X,I1,2X,A10)
C Effective file loading
open(unit=1,file='table2a.dat', status='old')
write(6,*) '....Loading file: table2a.dat'
do i__=1,84
read(1,'(A83)')ar__1
read(ar__1,2)
+ Date_1(i__),Middle(i__),Code_1(i__),Tint(i__),linf(i__),
+ lsup(i__),Aper1(i__),Aper2(i__),PA(i__),Seeing(i__),Res(i__),
+ Res2(i__),Filename(i__)
if(ar__1(66:69) .EQ. '') Res(i__) = rNULL__
if(ar__1(71:71) .EQ. '') Res2(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Date_1(i__),Middle(i__),Code_1(i__),Tint(i__),linf(i__),
+ lsup(i__),Aper1(i__),Aper2(i__),PA(i__),Seeing(i__),Res(i__),
+ Res2(i__),Filename(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3a.dat' ! Emission line and continuum flux measurements.
C Format for file interpretation
3 format(
+ A10,1X,F11.3,1X,A2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,file='table3a.dat', status='old')
write(6,*) '....Loading file: table3a.dat'
do i__=1,56
read(1,'(A58)')ar__2
read(ar__2,3)
+ Date_2(i__),HJD_1(i__),Code_2(i__),Fl5100(i__),e_Fl5100(i__),
+ FHbeta(i__),e_FHbeta(i__),FHalpha(i__),e_FHalpha(i__)
if(ar__2(27:31) .EQ. '') Fl5100(i__) = rNULL__
if(ar__2(33:36) .EQ. '') e_Fl5100(i__) = rNULL__
if(ar__2(38:42) .EQ. '') FHbeta(i__) = rNULL__
if(ar__2(44:47) .EQ. '') e_FHbeta(i__) = rNULL__
if(ar__2(49:53) .EQ. '') FHalpha(i__) = rNULL__
if(ar__2(55:58) .EQ. '') e_FHalpha(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Date_2(i__),HJD_1(i__),Code_2(i__),Fl5100(i__),e_Fl5100(i__),
+ FHbeta(i__),e_FHbeta(i__),FHalpha(i__),e_FHalpha(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end