Conversion of standardized ReadMe file for
file /./ftp/cats/J/AZh/86/35 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-19
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/AZh/86/35 Radio source 111MHz interplan. scintillation. II (Tyul'bashev, 2009)
*================================================================================
*Study of Compact Radio Sources Using Interplanetary Scintillations at 111 MHz.
*The Pearson-Readhead Sample.
* Tyul'bashev S.A.
* <Astron. Zh. 86, 35 (2009)>
* =2009AZh....86...35T
* =2009ARep...53...30T
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Estimated flux densities and compactness of
the studied sources
integer*4 nr__
parameter (nr__=32) ! Number of records
character*117 ar__ ! Full-size record
character*9 Source ! Source name (BHHMM+DDd)
character*11 ONames ! Other name(s) of the source
character*2 Type ! [QG BL] source type (Quasar, Galaxy or BL Lac)
character*2 l_DS ! [<=] Limit flag on DS
real*4 DS ! (Jy) ? Estimated 111MHz flux-density fluctuation
character*2 l_DS2 ! [<=] Limit flag on DS2
real*4 DS2 ! (Jy) ? Second value of 111MHz flux-density
* fluctuation when interval
real*4 eps ! (deg) ?=- Mean solar elongation of the source during
* the observation
character*2 l_Ssct ! [<=] Limit flag on Ssct
real*4 Ssct ! (Jy) ? Estimated 111MHz flux density of the
* scintillating component
real*4 Ssct2 ! (Jy) ? Second value of 111MHz flux density of the
* scintillating component (1)
character*2 l_Sint ! [<=] Limit flag on Sint
real*4 Sint ! (Jy) ? Estimated 111MHz peak flux density
real*4 e_Sint ! (Jy) ? rms uncertainty on Sint
character*1 u_Sint ! [?] Uncertainty flag on Sint
character*2 l_Sint2 ! [<=] Limit flag on Sint2
real*4 Sint2 ! (Jy) ? Second value of 111MHz peak flux density (1)
real*4 Sint2b ! (Jy) ? Upper value of 111MHz peak flux density
* when interval (1)
character*1 u_Sint2 ! [?] Uncertainty flag on Sint2
character*7 Date ! ("MM/YYYY") Observation date (2)
character*7 Date2 ! ("MM/YYYY") Second observation date when interval
character*2 l_R111 ! [<=] Limit flag on R111
real*4 R111 ! ? Compactness (R=Scomp/Sint) at 111MHz
character*2 l_R5000 ! [<=] Limit flag on R5000
real*4 R5000 ! ? Compactness (R=Scomp/Sint) at 5GHz
integer*4 r_R111 ! ? Compactness parameter reference, in
* refs.dat file
*Note (1): Estimates for sources with low flux densities obtained by
* extrapolating the NED high-frequency data
*Note (2): The 1996 observations where carried out at 102.5MHz.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Estimated flux density and compactness of the
previously studied sources
integer*4 nr__1
parameter (nr__1=27) ! Number of records
character*67 ar__1 ! Full-size record
character*9 Name ! Source name (BHHMM+DDd)
character*8 OName ! Other name(s) of the source
character*2 Type_1 ! [QG BL] source type (Quasar, Galaxy or BL Lac)
character*2 l_Ssct_1 ! [<= ] Limit flag on Ssct
real*4 Ssct_1 ! (Jy) ? Estimated 111MHz flux density of the
* scintillating component
character*1 u_Ssct ! [?] Uncertainty flag on Ssct
character*2 r_Ssct ! References for Ssct, in refs.dat file
character*2 l_Sint_1 ! [<= ] Limit flag on Sint
real*4 Sint_1 ! (Jy) Estimated 111MHz peak flux density
character*1 u_Sint_1 ! [?] Uncertainty flag on Sint
real*4 e_Sint_1 ! (Jy) ? rms uncertainty on Sint
character*1 u_e_Sint ! [?] Uncertainty flag on e_Sint
character*2 l_R111_1 ! [<= ] Limit flag on R111
real*4 R111_1 ! ?=- Compactness (R=Scomp/Sint ) at 111MHz
character*2 l_R5000_1 ! [>= ] Limit flag on R5000
real*4 R5000_1 ! ?=- Compactness (R=Scomp/Sint ) at 5GHz
integer*4 r_R5000 ! ? Reference for R5000, in refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__2
parameter (nr__2=22) ! Number of records
character*71 ar__2 ! Full-size record
integer*4 Ref ! Reference code
character*19 BibCode ! BibCode
character*23 Aut ! Auhtor's name
character*20 Com ! Comments
C=============================================================================
C Loading file 'table1.dat' ! Estimated flux densities and compactness of
* the studied sources
C Format for file interpretation
1 format(
+ A9,2X,A11,2X,A2,1X,A2,F4.2,2X,A2,F3.1,1X,F5.1,1X,A2,F5.2,1X,
+ F3.1,A2,F5.1,1X,F5.2,A1,A2,F4.2,1X,F3.1,A1,2X,A7,1X,A7,1X,A2,
+ F4.2,1X,A2,F4.2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,32
read(1,'(A117)')ar__
read(ar__,1)
+ Source,ONames,Type,l_DS,DS,l_DS2,DS2,eps,l_Ssct,Ssct,Ssct2,
+ l_Sint,Sint,e_Sint,u_Sint,l_Sint2,Sint2,Sint2b,u_Sint2,Date,
+ Date2,l_R111,R111,l_R5000,R5000,r_R111
if(ar__(30:33) .EQ. '') DS = rNULL__
if(ar__(38:40) .EQ. '') DS2 = rNULL__
if (idig(ar__(42:46)).EQ.0) eps = rNULL__
if(ar__(50:54) .EQ. '') Ssct = rNULL__
if(ar__(56:58) .EQ. '') Ssct2 = rNULL__
if(ar__(61:65) .EQ. '') Sint = rNULL__
if(ar__(67:71) .EQ. '') e_Sint = rNULL__
if(ar__(75:78) .EQ. '') Sint2 = rNULL__
if(ar__(80:82) .EQ. '') Sint2b = rNULL__
if(ar__(104:107) .EQ. '') R111 = rNULL__
if(ar__(111:114) .EQ. '') R5000 = rNULL__
if(ar__(116:117) .EQ. '') r_R111 = iNULL__
c ..............Just test output...........
write(6,1)
+ Source,ONames,Type,l_DS,DS,l_DS2,DS2,eps,l_Ssct,Ssct,Ssct2,
+ l_Sint,Sint,e_Sint,u_Sint,l_Sint2,Sint2,Sint2b,u_Sint2,Date,
+ Date2,l_R111,R111,l_R5000,R5000,r_R111
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Estimated flux density and compactness of the
* previously studied sources
C Format for file interpretation
2 format(
+ A9,1X,A8,1X,A2,1X,A2,F5.2,A1,1X,A2,1X,A2,F5.1,A1,1X,F4.2,A1,
+ 1X,A2,F5.3,1X,A2,F5.3,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,27
read(1,'(A67)')ar__1
read(ar__1,2)
+ Name,OName,Type_1,l_Ssct_1,Ssct_1,u_Ssct,r_Ssct,l_Sint_1,
+ Sint_1,u_Sint_1,e_Sint_1,u_e_Sint,l_R111_1,R111_1,l_R5000_1,
+ R5000_1,r_R5000
if(ar__1(25:29) .EQ. '') Ssct_1 = rNULL__
if(ar__1(44:47) .EQ. '') e_Sint_1 = rNULL__
if (idig(ar__1(52:56)).EQ.0) R111_1 = rNULL__
if (idig(ar__1(60:64)).EQ.0) R5000_1 = rNULL__
if(ar__1(66:67) .EQ. '') r_R5000 = iNULL__
c ..............Just test output...........
write(6,2)
+ Name,OName,Type_1,l_Ssct_1,Ssct_1,u_Ssct,r_Ssct,l_Sint_1,
+ Sint_1,u_Sint_1,e_Sint_1,u_e_Sint,l_R111_1,R111_1,l_R5000_1,
+ R5000_1,r_R5000
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
3 format(I2,1X,A19,1X,A23,5X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,22
read(1,'(A71)')ar__2
read(ar__2,3)Ref,BibCode,Aut,Com
c ..............Just test output...........
write(6,3)Ref,BibCode,Aut,Com
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