Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/435/521 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-12
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/435/521 Radio Sources in Low-Luminosity AGNs. IV (Nagar+, 2005)
*================================================================================
*Radio Sources in Low-Luminosity Active Galactic Nuclei.
*IV. Radio Luminosity Function, Importance of Jet Power,
*and Radio Properties of the Complete Palomar Sample
* Nagar N.M., Falcke H., Wilson A.S.
* <Astron. Astrophys. 435, 521 (2005)>
* =2005A&A...435..521N
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! *High-resolution 15GHz VLA imaging of Palomar
AGN & LLAGN (low-luminosity AGN)
integer*4 nr__
parameter (nr__=197) ! Number of records
character*103 ar__ ! Full-size record
C J2000.0 position composed of: RAObsh RAObsm RAObss DEObs- DEObsd DEObsm DEObss
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000.0
real*8 DEdeg (nr__) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*7 Name (nr__) ! Name from NGC or IC catalog
character*9 AType (nr__) ! Activity Type (G1)
real*4 TType (nr__) ! Morphological type
real*4 Dist (nr__) ! (Mpc) Distance
integer*4 RAObsh (nr__) ! (h) ? Observed right ascension (J2000.0)
integer*4 RAObsm (nr__) ! (min) ? Observed right ascension (J2000.0)
real*4 RAObss (nr__) ! (s) ? Observed right ascension (J2000.0)
character*1 DEObs_ (nr__) ! Observed declination sign (J2000.0)
integer*4 DEObsd (nr__) ! (deg) ? Observed declination (J2000.0)
integer*4 DEObsm (nr__) ! (arcmin) ? Observed declination (J2000.0)
real*4 DEObss (nr__) ! (arcsec) ? Observed declination (J2000.0)
real*4 Delta (nr__) ! (arcsec) ? Angular separation between nuclear radio and
* optical positions
real*4 Sp15GHz (nr__) ! (mJy) ? Peak VLA 15GHz flux
character*1 l_St15GHz (nr__) ! Limit sign for total flux
real*4 St15GHz (nr__) ! (mJy) ? Total VLA 15GHz flux
character*1 l_Lp15GHz (nr__) ! Limit sign for luminosity
real*4 Lp15GHz (nr__) ! (W/Hz) ? Corresponding peak 15GHz luminosity
integer*4 Note (nr__) ! ? Comments, in notes.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *VLBI observations of Palomar LLAGN & AGN
integer*4 nr__1
parameter (nr__1=44) ! Number of records
character*159 ar__1 ! Full-size record
C J2000.0 position composed of: RAObsh RAObsm RAObss DEObs- DEObsd DEObsm DEObss
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000.0
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000.0
C ---------------------------------- ! (position vector(s) in degrees)
character*7 Name_1 (nr__1) ! Name from NGC catalog
character*9 AType_1 (nr__1) ! Activity Type (from Ho et al., 1997,
* Cat. <J/ApJS/112/315>) (G1)
real*4 TType_1 (nr__1) ! Morphological type (from Ho et al., 1997,
* Cat. <J/ApJS/112/315>)
real*4 Dist_1 (nr__1) ! (Mpc) Distance (from Ho et al., 1997,
* Cat. <J/ApJS/112/315>)
real*4 LHalpha (nr__1) ! (10-7W) Halpha luminosity (from Ho et al., 1997,
* Cat. <J/ApJS/112/315>)
character*1 n_LHalpha (nr__1) ! [Lbc] Limit/quality flag for LHalpha (from Ho
* et al., 1997, Cat. <J/ApJS/112/315>) (1)
integer*4 RAObsh_1 (nr__1) ! (h) ? Observed right ascension (J2000.0)
integer*4 RAObsm_1 (nr__1) ! (min) ? Observed right ascension (J2000.0)
real*8 RAObss_1 (nr__1) ! (s) ? Observed right ascension (J2000.0)
character*1 DEObs__1 (nr__1) ! Observed declination sign (J2000.0)
integer*4 DEObsd_1 (nr__1) ! (deg) ? Observed declination (J2000.0)
integer*4 DEObsm_1 (nr__1) ! (arcmin) ? Observed declination (J2000.0)
real*4 DEObss_1 (nr__1) ! (arcsec) ? Observed declination (J2000.0)
real*4 Sp5GHz (nr__1) ! (mJy) ? Peak VLBA/I (5GHz) flux of the nuclear source
character*1 l_St5GHz (nr__1) ! Limit sign for Total VLBA/I flux of nuclear
* source
real*4 St5GHz (nr__1) ! (mJy) ? Total VLBA/I (5GHz) flux of the nuclear source
real*4 Stot (nr__1) ! (mJy) ? Total mas-scale flux
character*1 l_Lp15GHz_1(nr__1) ! Limit sign for luminosity
real*4 Lp15GHz_1 (nr__1) ! (W/Hz) 15GHz luminosity (derived from the Peak VLBA/I
* flux)
real*4 Size (nr__1) ! (mas) ? Angular size of VLBA/I mas-scale extended
* emission
integer*4 PA (nr__1) ! (deg) [-100/290]? Position angle of VLBA/I mas-scale
* extended emission
real*4 majB (nr__1) ! (mas) ? Beam size (major axis) in milliarcsec
real*4 minB (nr__1) ! (mas) ? Beam size (minor axis) in milliarcsec
integer*4 PAB (nr__1) ! (deg) ? Beam position angle (North to East)
character*2 l_logT (nr__1) ! [<>= ] limit sign to implied brightness
* temperature
real*4 logT (nr__1) ! ([K]) ? Implied brightness temperature
real*4 Scale (nr__1) ! (pc/mas) Linear scale at the distance of the source
character*5 Note_1 (nr__1) ! Comments, in notes.dat file
*Note (1): Limit/quality flag for LHalpha as follows:
* L = 3{sigma} lower limit
* b = Values with uncertainties of +/-30%-50%
* c = highly uncertain values with probable errors of +/-100%
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'notes.dat' ! Notes
integer*4 nr__2
parameter (nr__2=35) ! Number of records
character*80 ar__2 ! Full-size record
integer*4 Note_2 (nr__2) ! Note number
character*77 Text (nr__2) ! Text of the note
C=============================================================================
C Loading file 'table1.dat' ! *High-resolution 15GHz VLA imaging of Palomar
* AGN & LLAGN (low-luminosity AGN)
C Format for file interpretation
1 format(
+ A7,2X,A9,1X,F4.1,3X,F4.1,3X,I2,1X,I2,1X,F6.3,3X,A1,I2,1X,I2,
+ 1X,F5.2,2X,F4.1,7X,F6.1,2X,A1,F6.1,4X,A1,1X,F5.2,2X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,197
read(1,'(A103)')ar__
read(ar__,1)
+ Name(i__),AType(i__),TType(i__),Dist(i__),RAObsh(i__),
+ RAObsm(i__),RAObss(i__),DEObs_(i__),DEObsd(i__),DEObsm(i__),
+ DEObss(i__),Delta(i__),Sp15GHz(i__),l_St15GHz(i__),
+ St15GHz(i__),l_Lp15GHz(i__),Lp15GHz(i__),Note(i__)
if(ar__(34:35) .EQ. '') RAObsh(i__) = iNULL__
if(ar__(37:38) .EQ. '') RAObsm(i__) = iNULL__
if(ar__(40:45) .EQ. '') RAObss(i__) = rNULL__
if(ar__(50:51) .EQ. '') DEObsd(i__) = iNULL__
if(ar__(53:54) .EQ. '') DEObsm(i__) = iNULL__
if(ar__(56:60) .EQ. '') DEObss(i__) = rNULL__
if(ar__(63:66) .EQ. '') Delta(i__) = rNULL__
if(ar__(74:79) .EQ. '') Sp15GHz(i__) = rNULL__
if(ar__(83:88) .EQ. '') St15GHz(i__) = rNULL__
if(ar__(95:99) .EQ. '') Lp15GHz(i__) = rNULL__
if(ar__(102:103) .EQ. '') Note(i__) = iNULL__
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(RAObsh(i__) .GT. -180) RAdeg(i__)=RAObsh(i__)*15.
if(RAObsm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAObsm(i__)/4.
if(RAObss(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAObss(i__)/240.
if(DEObsd(i__) .GE. 0) DEdeg(i__)=DEObsd(i__)
if(DEObsm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEObsm(i__)/60.
if(DEObss(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEObss(i__)/3600.
if(DEObs_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),AType(i__),TType(i__),Dist(i__),RAObsh(i__),
+ RAObsm(i__),RAObss(i__),DEObs_(i__),DEObsd(i__),DEObsm(i__),
+ DEObss(i__),Delta(i__),Sp15GHz(i__),l_St15GHz(i__),
+ St15GHz(i__),l_Lp15GHz(i__),Lp15GHz(i__),Note(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 'table2.dat' ! *VLBI observations of Palomar LLAGN & AGN
C Format for file interpretation
2 format(
+ A7,2X,A9,1X,F4.1,3X,F4.1,2X,F4.1,1X,A1,2X,I2,1X,I2,1X,F7.4,1X,
+ A1,I2,1X,I2,1X,F6.3,2X,F6.1,3X,A1,1X,F4.1,3X,F5.1,2X,A1,1X,
+ F5.2,3X,F4.1,3X,I4,3X,F3.1,2X,F3.1,3X,I3,3X,A2,1X,F4.1,3X,
+ F4.2,2X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,44
read(1,'(A159)')ar__1
read(ar__1,2)
+ Name_1(i__),AType_1(i__),TType_1(i__),Dist_1(i__),
+ LHalpha(i__),n_LHalpha(i__),RAObsh_1(i__),RAObsm_1(i__),
+ RAObss_1(i__),DEObs__1(i__),DEObsd_1(i__),DEObsm_1(i__),
+ DEObss_1(i__),Sp5GHz(i__),l_St5GHz(i__),St5GHz(i__),Stot(i__),
+ l_Lp15GHz_1(i__),Lp15GHz_1(i__),Size(i__),PA(i__),majB(i__),
+ minB(i__),PAB(i__),l_logT(i__),logT(i__),Scale(i__),
+ Note_1(i__)
if(ar__1(41:42) .EQ. '') RAObsh_1(i__) = iNULL__
if(ar__1(44:45) .EQ. '') RAObsm_1(i__) = iNULL__
if(ar__1(47:53) .EQ. '') RAObss_1(i__) = rNULL__
if(ar__1(56:57) .EQ. '') DEObsd_1(i__) = iNULL__
if(ar__1(59:60) .EQ. '') DEObsm_1(i__) = iNULL__
if(ar__1(62:67) .EQ. '') DEObss_1(i__) = rNULL__
if(ar__1(70:75) .EQ. '') Sp5GHz(i__) = rNULL__
if(ar__1(81:84) .EQ. '') St5GHz(i__) = rNULL__
if(ar__1(88:92) .EQ. '') Stot(i__) = rNULL__
if(ar__1(105:108) .EQ. '') Size(i__) = rNULL__
if(ar__1(112:115) .EQ. '') PA(i__) = iNULL__
if(ar__1(119:121) .EQ. '') majB(i__) = rNULL__
if(ar__1(124:126) .EQ. '') minB(i__) = rNULL__
if(ar__1(130:132) .EQ. '') PAB(i__) = iNULL__
if(ar__1(139:142) .EQ. '') logT(i__) = rNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAObsh_1(i__) .GT. -180) RAdeg_1(i__)=RAObsh_1(i__)*15.
if(RAObsm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAObsm_1(i__)/4.
if(RAObss_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAObss_1(i__)/240.
if(DEObsd_1(i__) .GE. 0) DEdeg_1(i__)=DEObsd_1(i__)
if(DEObsm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEObsm_1(i__)/60.
if(DEObss_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEObss_1(i__)/3600.
if(DEObs__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),AType_1(i__),TType_1(i__),Dist_1(i__),
+ LHalpha(i__),n_LHalpha(i__),RAObsh_1(i__),RAObsm_1(i__),
+ RAObss_1(i__),DEObs__1(i__),DEObsd_1(i__),DEObsm_1(i__),
+ DEObss_1(i__),Sp5GHz(i__),l_St5GHz(i__),St5GHz(i__),Stot(i__),
+ l_Lp15GHz_1(i__),Lp15GHz_1(i__),Size(i__),PA(i__),majB(i__),
+ minB(i__),PAB(i__),l_logT(i__),logT(i__),Scale(i__),
+ Note_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'notes.dat' ! Notes
C Format for file interpretation
3 format(I2,1X,A77)
C Effective file loading
open(unit=1,status='old',file=
+'notes.dat')
write(6,*) '....Loading file: notes.dat'
do i__=1,35
read(1,'(A80)')ar__2
read(ar__2,3)Note_2(i__),Text(i__)
c ..............Just test output...........
write(6,3)Note_2(i__),Text(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end