Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/439/487 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 2024-Apr-16
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/439/487 Nuker law fits of radio galaxies (De Ruiter+, 2005)
*================================================================================
*Are radio galaxies and quiescent galaxies different?
*Results from the analysis of HST brightness profiles.
* De Ruiter H.R., Parma P., Capetti A., Fanti R., Morganti R., Santantonio L.
* <Astron. Astrophys., 439, 487-496 (2005)>
* =2005A&A...439..487D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Nuker law parameters of B2 radio galaxies
with central dust
integer*4 nr__
parameter (nr__=17) ! Number of records
character*102 ar__ ! Full-size record
character*10 Name ! Name of radio source
real*4 alpha ! Nuker parameter {alpha} (1)
real*4 e_alpha ! 1-{sigme} uncertainty of alpha
real*4 beta ! Nuker parameter {beta} (1)
real*4 e_beta ! 1-{sigme} uncertainty of beta
real*4 gamma ! Nuker parameter {gamma} (1)
real*4 e_gamma ! 1-{sigme} uncertainty of gamma
real*4 rb ! (arcsec) Nuker parameter r_b_ (1)
real*4 e_rb ! (arcsec) 1-{sigme} uncertainty of rb
real*4 mub ! (mag/arcsec2) Nuker parameter {mu}_b_ (1)
real*4 e_mub ! (mag/arcsec2) 1-{sigme} uncertainty of mub
real*4 z ! Redshift
real*4 log_Pt ! ([W/Hz]) logarithm of total radio power (at 1.4GHz)
real*4 VMAG ! (mag) Absolute V magnitude of galaxy
real*4 b4 ! (%) ?=- Deviation from pure ellipse in %
real*4 e_b4 ! (%) ?=- 1-{sigme} uncertainty of b4 in %
character*8 OName ! Alternative NGC or 3C name
*Note (1): The "Nuker profile" is an empirical law introduced by Lauer et al.
* (1995AJ....110.2622L). In terms of magnitudes the Nuker law reads:
* {mu}(r) = {mu}_b_ - [({beta}-{gamma})/{alpha}]*2.5log2
* + 2.5*{gamma}*log(r/r_b_)
* + 2.5*[({beta}-{gamma})/{alpha}]*log(1+(r/r_b_)^{alpha}^)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Nuker law parameters of B2 radio galaxies
without central dust
integer*4 nr__1
parameter (nr__1=22) ! Number of records
character*102 ar__1 ! Full-size record
character*10 Name_1 ! Name of radio source
real*4 alpha_1 ! Nuker parameter {alpha} (1)
real*4 e_alpha_1 ! 1-{sigme} uncertainty of alpha
real*4 beta_1 ! Nuker parameter {beta} (1)
real*4 e_beta_1 ! 1-{sigme} uncertainty of beta
real*4 gamma_1 ! Nuker parameter {gamma} (1)
real*4 e_gamma_1 ! 1-{sigme} uncertainty of gamma
real*4 rb_1 ! (arcsec) Nuker parameter r_b_ (1)
real*4 e_rb_1 ! (arcsec) 1-{sigme} uncertainty of rb
real*4 mub_1 ! (mag/arcsec2) Nuker parameter {mu}_b_ (1)
real*4 e_mub_1 ! (mag/arcsec2) 1-{sigme} uncertainty of mub
real*4 z_1 ! Redshift
real*4 log_Pt_1 ! ([W/Hz]) logarithm of total radio power (at 1.4GHz)
real*4 VMAG_1 ! (mag) Absolute V magnitude of galaxy
real*4 b4_1 ! (%) ?=- Deviation from pure ellipse in %
real*4 e_b4_1 ! (%) ?=- 1-{sigme} uncertainty of b4 in %
character*8 OName_1 ! Alternative NGC or 3C name
*Note (1): The "Nuker profile" is an empirical law introduced by Lauer et al.
* (1995AJ....110.2622L). In terms of magnitudes the Nuker law reads:
* {mu}(r) = {mu}_b_ - [({beta}-{gamma})/{alpha}]*2.5log2
* + 2.5*{gamma}*log(r/r_b_)
* + 2.5*[({beta}-{gamma})/{alpha}]*log(1+(r/r_b_)^{alpha}^)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Nuker law parameters of 3CR radio galaxies
of FR type II
integer*4 nr__2
parameter (nr__2=9) ! Number of records
character*102 ar__2 ! Full-size record
character*10 Name_2 ! Name of radio source
real*4 alpha_2 ! Nuker parameter {alpha} (1)
real*4 e_alpha_2 ! 1-{sigme} uncertainty of alpha
real*4 beta_2 ! Nuker parameter {beta} (1)
real*4 e_beta_2 ! 1-{sigme} uncertainty of beta
real*4 gamma_2 ! Nuker parameter {gamma} (1)
real*4 e_gamma_2 ! 1-{sigme} uncertainty of gamma
real*4 rb_2 ! (arcsec) Nuker parameter r_b_ (1)
real*4 e_rb_2 ! (arcsec) 1-{sigme} uncertainty of rb
real*4 mub_2 ! (mag/arcsec2) Nuker parameter {mu}_b_ (1)
real*4 e_mub_2 ! (mag/arcsec2) 1-{sigme} uncertainty of mub
real*4 z_2 ! Redshift
real*4 log_Pt_1 ! ([W/Hz]) logarithm of total radio power (at 1.4GHz)
real*4 VMAG_2 ! (mag) Absolute V magnitude of galaxy
real*4 b4_2 ! (%) ?=- Deviation from pure ellipse in %
real*4 e_b4_2 ! (%) ?=- 1-{sigme} uncertainty of b4 in %
character*8 OName_2 ! Alternative NGC or 3C name
*Note (1): The "Nuker profile" is an empirical law introduced by Lauer et al.
* (1995AJ....110.2622L). In terms of magnitudes the Nuker law reads:
* {mu}(r) = {mu}_b_ - [({beta}-{gamma})/{alpha}]*2.5log2
* + 2.5*{gamma}*log(r/r_b_)
* + 2.5*[({beta}-{gamma})/{alpha}]*log(1+(r/r_b_)^{alpha}^)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Faber et al. (1997AJ.....114.1771F) and
Laine et al. (2003AJ.....125..478L) basic data
integer*4 nr__3
parameter (nr__3=33) ! Number of records
character*46 ar__3 ! Full-size record
character*9 Galaxy ! Name of galaxy (1)
real*4 z_3 ! redshift
real*4 VMAG_3 ! (mag) Absolute magnitude of galaxy
character*1 l_log_Pt ! Upper limit flag (<) of radio power (2)
real*4 log_Pt_1 ! ([W/Hz]) ?=- logarithm of total radio power (at 1.4GHz)
character*1 Sample ! F: Faber+ (1997); L: Laine+ (2003) (1)
character*13 Info ! Information on radio source (3)
*Note (1): Galaxy names taken from the original papers by Faber et al.
* (1997AJ.....114.1771F) and Laine et al. (2003AJ.....125..478L).
*Note (2): 3-{sigme} detection limit is 1.35mJy (NVSS).
*Note (3): Some objects could not be searched for radio emission (e.g. because
* too far to the south for NVSS, or because a precise position is
* lacking). Galaxies with VMAG>-17 were excluded from the analysis.
* Abell 168 was excluded because of a large discrepancy between the fit
* parameter {gamma} and the slope directly measured at 0.05 arcsec.
C=============================================================================
C Loading file 'table1.dat' ! Nuker law parameters of B2 radio galaxies
* with central dust
C Format for file interpretation
1 format(
+ A10,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.3,
+ 1X,F5.3,1X,F5.2,1X,F4.2,1X,F6.4,1X,F5.2,1X,F6.2,1X,F4.1,1X,
+ F3.1,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,17
read(1,'(A102)')ar__
read(ar__,1)
+ Name,alpha,e_alpha,beta,e_beta,gamma,e_gamma,rb,e_rb,mub,
+ e_mub,z,log_Pt,VMAG,b4,e_b4,OName
if (idig(ar__(86:89)).EQ.0) b4 = rNULL__
if (idig(ar__(91:93)).EQ.0) e_b4 = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,alpha,e_alpha,beta,e_beta,gamma,e_gamma,rb,e_rb,mub,
+ e_mub,z,log_Pt,VMAG,b4,e_b4,OName
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Nuker law parameters of B2 radio galaxies
* without central dust
C Format for file interpretation
2 format(
+ A10,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.3,
+ 1X,F5.3,1X,F5.2,1X,F4.2,1X,F6.4,1X,F5.2,1X,F6.2,1X,F4.1,1X,
+ F3.1,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,22
read(1,'(A102)')ar__1
read(ar__1,2)
+ Name_1,alpha_1,e_alpha_1,beta_1,e_beta_1,gamma_1,e_gamma_1,
+ rb_1,e_rb_1,mub_1,e_mub_1,z_1,log_Pt_1,VMAG_1,b4_1,e_b4_1,
+ OName_1
if (idig(ar__1(86:89)).EQ.0) b4_1 = rNULL__
if (idig(ar__1(91:93)).EQ.0) e_b4_1 = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,alpha_1,e_alpha_1,beta_1,e_beta_1,gamma_1,e_gamma_1,
+ rb_1,e_rb_1,mub_1,e_mub_1,z_1,log_Pt_1,VMAG_1,b4_1,e_b4_1,
+ OName_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Nuker law parameters of 3CR radio galaxies
* of FR type II
C Format for file interpretation
3 format(
+ A10,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.3,
+ 1X,F5.3,1X,F5.2,1X,F4.2,1X,F6.4,1X,F5.2,1X,F6.2,1X,F4.1,1X,
+ F3.1,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,9
read(1,'(A102)')ar__2
read(ar__2,3)
+ Name_2,alpha_2,e_alpha_2,beta_2,e_beta_2,gamma_2,e_gamma_2,
+ rb_2,e_rb_2,mub_2,e_mub_2,z_2,log_Pt_1,VMAG_2,b4_2,e_b4_2,
+ OName_2
if (idig(ar__2(86:89)).EQ.0) b4_2 = rNULL__
if (idig(ar__2(91:93)).EQ.0) e_b4_2 = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,alpha_2,e_alpha_2,beta_2,e_beta_2,gamma_2,e_gamma_2,
+ rb_2,e_rb_2,mub_2,e_mub_2,z_2,log_Pt_1,VMAG_2,b4_2,e_b4_2,
+ OName_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Faber et al. (1997AJ.....114.1771F) and
* Laine et al. (2003AJ.....125..478L) basic data
C Format for file interpretation
4 format(A9,1X,F6.4,1X,F6.2,1X,A1,F5.2,1X,A1,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,33
read(1,'(A46)')ar__3
read(ar__3,4)Galaxy,z_3,VMAG_3,l_log_Pt,log_Pt_1,Sample,Info
if (idig(ar__3(26:30)).EQ.0) log_Pt_1 = rNULL__
c ..............Just test output...........
write(6,4)Galaxy,z_3,VMAG_3,l_log_Pt,log_Pt_1,Sample,Info
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