Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/569/A95 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-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/569/A95 SDSS optically selected BL Lac candidates (Kuegler+, 2014)
*================================================================================
*Properties of optically selected BL Lacertae candidates from the SDSS.
* Kuegler S.D., Nilsson K., Heidt J., Esser J., Schultz T.
* <Astron. Astrophys., 569, A95-95 (2014)>
* =2014A&A...569A..95K
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Host galaxy results for the bulge model
integer*4 nr__
parameter (nr__=66) ! Number of records
character*93 ar__ ! Full-size record
character*18 SDSS ! Target SDSS designation (HHMMSS.ss+DDMMSS.s)
character*3 Tel ! Telescope used: CA=Calar Alto 2.2m,
* NOT=Nordic Optical Telescope, NTT at ESO
real*4 z ! [0.04/0.71]?=- Redshift, listed only when
* the redshift determination is secure
real*4 rmag_c ! (mag) ?=- SDSS r' magnitude of the core m_c_
real*4 e_rmag_c ! (mag) ? rms uncertainty in rmag.c
real*4 rmag_g ! (mag) r' magnitude of the host galaxy m_g_
real*4 e_rmag_g ! (mag) rms uncertainty in rmag.g
real*4 re ! (arcsec) [0.5/8.4] Effective radius r_e_
real*4 e_re ! (arcsec) rms uncertainty on re
real*4 ell_g ! [0/1] Ellipticity of host galaxy {epsilon}_g_
integer*4 pa_g ! (deg) Position angle of the host galaxy {theta}_g_
character*1 T_g ! [BU] Host galaxy type: B=bulge, U=undefined
real*4 fc ! [0/1] Fraction of the core flux in total flux
real*4 e_fc ! rms uncertainty on fc
character*1 u_fc ! [)] Uncertainty flag on fc
real*4 rMAG ! (mag) ?=- Absolute r' magnitude M_r'_
real*4 Re_1 ! (kpc) ?=- Effective radius in kpc
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Main properties objects discussed in Heidt & Nilsson
(2011, J/A+A/529/A162) and the present work
integer*4 nr__1
parameter (nr__1=195) ! Number of records
character*90 ar__1 ! Full-size record
character*18 SDSS_1 ! Target SDSS designation (HHMMSS.ss+DDMMSS.s)
character*1 l_z ! Limit flag on z
real*4 z_1 ! [0.04/5.1]?=- SDSS spectroscopic redshift
character*1 u_z ! [u] SDSS z flag (u for uncertain)
character*1 l_aox ! Limit flag on aox
real*4 aox ! Spectral index Optical/X-ray {alpha}_ox_ (1)
character*1 l_aro ! Limit flag on aro
real*4 aro ! ?=- Spectral index Radio/Optical {alpha}_ro_ (1)
character*1 l_Pol ! Limit flag on Pol
real*4 Pol ! (%) ?=- Polarisation
real*4 e_Pol ! (%) ?=- rms uncertainty on Pol
real*4 VL ! (mag) [0.03/0.15] Variability limit
character*1 f_Vamp ! [-] Variability amplitude is negative?
real*4 Vamp ! (mag) Variability amplitude
character*1 l_fc ! Limit flag on fc
real*4 fc_1 ! Core (AGN) fraction over the total flux
* fraction as derived from the host galaxy fits
real*4 e_fc_1 ! ?=- rms uncertainty on fc
real*4 lognuP ! ([Hz]) ?=- Peak frequency determined from SED fits
character*6 f_lognuP ! Flag for uncertain lognuP with respect to
* the SED fits (2)
*Note (1): All these properties are from Collinge et al. (2005AJ....129.2542C,
* Cat. J/AJ/129/2542).
*Note (2): The flags indicate whether the peak estimate is uncertain due to the
* shallowness of the radio data:
* radio = if only one or no radio point available
* galaxy = due to the presence of a strong host galaxy if the host galaxy
* fraction of the total flux exceeds 0.5
C=============================================================================
C Loading file 'table3.dat' ! Host galaxy results for the bulge model
C Format for file interpretation
1 format(
+ A18,1X,A3,1X,F6.4,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F3.1,1X,
+ F3.1,1X,F4.2,1X,I3,1X,A1,1X,F4.2,1X,F4.2,A1,1X,F6.2,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,66
read(1,'(A93)')ar__
read(ar__,1)
+ SDSS,Tel,z,rmag_c,e_rmag_c,rmag_g,e_rmag_g,re,e_re,ell_g,pa_g,
+ T_g,fc,e_fc,u_fc,rMAG,Re_1
if (idig(ar__(24:29)).EQ.0) z = rNULL__
if (idig(ar__(31:35)).EQ.0) rmag_c = rNULL__
if(ar__(37:40) .EQ. '') e_rmag_c = rNULL__
if (idig(ar__(83:88)).EQ.0) rMAG = rNULL__
if (idig(ar__(90:93)).EQ.0) Re_1 = rNULL__
c ..............Just test output...........
write(6,1)
+ SDSS,Tel,z,rmag_c,e_rmag_c,rmag_g,e_rmag_g,re,e_re,ell_g,pa_g,
+ T_g,fc,e_fc,u_fc,rMAG,Re_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Main properties objects discussed in Heidt & Nilsson
* (2011, J/A+A/529/A162) and the present work
C Format for file interpretation
2 format(
+ A18,1X,A1,F6.4,1X,A1,1X,A1,F5.2,1X,A1,F4.2,1X,A1,F5.2,1X,F4.2,
+ 1X,F5.3,1X,A1,F5.3,1X,A1,F4.2,1X,F4.2,1X,F5.2,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,195
read(1,'(A90)')ar__1
read(ar__1,2)
+ SDSS_1,l_z,z_1,u_z,l_aox,aox,l_aro,aro,l_Pol,Pol,e_Pol,VL,
+ f_Vamp,Vamp,l_fc,fc_1,e_fc_1,lognuP,f_lognuP
if (idig(ar__1(21:26)).EQ.0) z_1 = rNULL__
if (idig(ar__1(38:41)).EQ.0) aro = rNULL__
if (idig(ar__1(44:48)).EQ.0) Pol = rNULL__
if (idig(ar__1(50:53)).EQ.0) e_Pol = rNULL__
if (idig(ar__1(74:77)).EQ.0) e_fc_1 = rNULL__
if (idig(ar__1(79:83)).EQ.0) lognuP = rNULL__
c ..............Just test output...........
write(6,2)
+ SDSS_1,l_z,z_1,u_z,l_aox,aox,l_aro,aro,l_Pol,Pol,e_Pol,VL,
+ f_Vamp,Vamp,l_fc,fc_1,e_fc_1,lognuP,f_lognuP
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end