Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/415/1597 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/MNRAS/415/1597 PACO catalogue of bright sources (Massardi+, 2011)
*================================================================================
*The Planck-ATCA Co-eval Observations project: the bright sample.
* Massardi M., Bonaldi A., Bonavera L., Lopez-Caniego M., De Zotti G.,
* Ekers R.D.
* <Mon. Not. R. Astron. Soc., 415, 1597-1610 (2011)>
* =2011MNRAS.415.1597M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tables1.dat' ! *The PACO bright catalogue
integer*4 nr__
parameter (nr__=1004) ! Number of records
character*386 ar__ ! Full-size record
C J2000 position composed of: RAhour DEdeg
character*14 AT20G ! AT20G name (JHHMMSS+DDMMSS)
real*8 RAhour ! (h) Right ascension in decimal hours (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
character*10 Date ! ("YYYY-MM-DD") Observation date
character*1 fs ! [s-] Flag "s" to identify epochs within
* 10-days from Planck observations
character*1 fe ! [e-] Flag "e" for extended sources
character*1 l_S4732 ! Limit flag on S4732
integer*4 S4732 ! (mJy) PACO flux at 4732MHz
character*1 l_S5244 ! Limit flag on S5244
integer*4 S5244 ! (mJy) PACO flux at 5244MHz
character*1 l_S5756 ! Limit flag on S5756
integer*4 S5756 ! (mJy) PACO flux at 5756MHz
character*1 l_S6268 ! Limit flag on S6268
integer*4 S6268 ! (mJy) PACO flux at 6268MHz
character*1 l_S8232 ! Limit flag on S8232
integer*4 S8232 ! (mJy) PACO flux at 8232MHz
character*1 l_S8744 ! Limit flag on S8744
integer*4 S8744 ! (mJy) PACO flux at 8744MHz
character*1 l_S9256 ! Limit flag on S9256
integer*4 S9256 ! (mJy) PACO flux at 9256MHz
character*1 l_S9768 ! Limit flag on S9768
integer*4 S9768 ! (mJy) PACO flux at 9768MHz
character*1 l_S17232 ! Limit flag on S17232
integer*4 S17232 ! (mJy) PACO flux at 17232MHz
character*1 l_S17744 ! Limit flag on S17744
integer*4 S17744 ! (mJy) PACO flux at 17744MHz
character*1 l_S18256 ! Limit flag on S18256
integer*4 S18256 ! (mJy) PACO flux at 18256MHz
character*1 l_S18768 ! Limit flag on S18768
integer*4 S18768 ! (mJy) PACO flux at 18768MHz
character*1 l_S23232 ! Limit flag on S23232
integer*4 S23232 ! (mJy) PACO flux at 23232MHz
character*1 l_S23744 ! Limit flag on S23744
integer*4 S23744 ! (mJy) PACO flux at 23744MHz
character*1 l_S24256 ! Limit flag on S24256
integer*4 S24256 ! (mJy) PACO flux at 24256MHz
character*1 l_S24768 ! Limit flag on S24768
integer*4 S24768 ! (mJy) PACO flux at 24768MHz
character*1 l_S32232 ! Limit flag on S32232
integer*4 S32232 ! (mJy) PACO flux at 32232MHz
character*1 l_S32744 ! Limit flag on S32744
integer*4 S32744 ! (mJy) PACO flux at 32744MHz
character*1 l_S33256 ! Limit flag on S33256
integer*4 S33256 ! (mJy) PACO flux at 33256MHz
character*1 l_S33768 ! Limit flag on S33768
integer*4 S33768 ! (mJy) PACO flux at 33768MHz
character*1 l_S38232 ! Limit flag on S38232
integer*4 S38232 ! (mJy) PACO flux at 38232MHz
character*1 l_S38744 ! Limit flag on S38744
integer*4 S38744 ! (mJy) PACO flux at 38744MHz
character*1 l_S39256 ! Limit flag on S39256
integer*4 S39256 ! (mJy) PACO flux at 39256MHz
character*1 l_S39768 ! Limit flag on S39768
integer*4 S39768 ! (mJy) PACO flux at 39768MHz
real*4 e_S4732 ! (mJy) rms uncertainty on S4732
real*4 e_S5244 ! (mJy) rms uncertainty on S5244
real*4 e_S5756 ! (mJy) rms uncertainty on S5756
real*4 e_S6268 ! (mJy) rms uncertainty on S6268
real*4 e_S8232 ! (mJy) rms uncertainty on S8232
real*4 e_S8744 ! (mJy) rms uncertainty on S8744
real*4 e_S9256 ! (mJy) rms uncertainty on S9256
real*4 e_S9768 ! (mJy) rms uncertainty on S9768
real*4 e_S17232 ! (mJy) rms uncertainty on S17232
real*4 e_S17744 ! (mJy) rms uncertainty on S17744
real*4 e_S18256 ! (mJy) rms uncertainty on S18256
real*4 e_S18768 ! (mJy) rms uncertainty on S18768
real*4 e_S23232 ! (mJy) rms uncertainty on S23232
real*4 e_S23744 ! (mJy) rms uncertainty on S23744
real*4 e_S24256 ! (mJy) rms uncertainty on S24256
real*4 e_S24768 ! (mJy) rms uncertainty on S24768
real*4 e_S32232 ! (mJy) rms uncertainty on S32232
real*4 e_S32744 ! (mJy) rms uncertainty on S32744
real*4 e_S33256 ! (mJy) rms uncertainty on S33256
real*4 e_S33768 ! (mJy) rms uncertainty on S33768
real*4 e_S38232 ! (mJy) rms uncertainty on S38232
real*4 e_S38744 ! (mJy) rms uncertainty on S38744
real*4 e_S39256 ! (mJy) rms uncertainty on S39256
real*4 e_S39768 ! (mJy) rms uncertainty on S39768
real*4 a ! Best fit parameter a (1)
real*4 b ! Best fit parameter b (1)
real*4 v0 ! (GHz) Best fit parameter v0 (1)
real*4 S0 ! (Jy) Best fit parameter S0 (1)
real*4 chi2 ! Best fit parameter chi^2^
*Note (1): Best fit parameters in the double power law:
* S(v)=S0/[(v/v0)^-a^+(v/v0)^-b^], where
* v is the frequency and S is the flux density in Jy.
C=============================================================================
C Loading file 'tables1.dat' ! *The PACO bright catalogue
C Format for file interpretation
1 format(
+ 5X,A14,1X,F10.7,1X,F11.7,1X,A10,1X,A1,1X,A1,1X,A1,I5,A1,I5,A1,
+ I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,
+ I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,I5,A1,
+ I5,A1,I5,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,
+ F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,
+ 1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,
+ F5.1,1X,F5.1,1X,F5.1,2X,F6.3,2X,F6.3,2X,F6.3,2X,F6.3,3X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'tables1.dat')
write(6,*) '....Loading file: tables1.dat'
do i__=1,1004
read(1,'(A386)')ar__
read(ar__,1)
+ AT20G,RAhour,DEdeg,Date,fs,fe,l_S4732,S4732,l_S5244,S5244,
+ l_S5756,S5756,l_S6268,S6268,l_S8232,S8232,l_S8744,S8744,
+ l_S9256,S9256,l_S9768,S9768,l_S17232,S17232,l_S17744,S17744,
+ l_S18256,S18256,l_S18768,S18768,l_S23232,S23232,l_S23744,
+ S23744,l_S24256,S24256,l_S24768,S24768,l_S32232,S32232,
+ l_S32744,S32744,l_S33256,S33256,l_S33768,S33768,l_S38232,
+ S38232,l_S38744,S38744,l_S39256,S39256,l_S39768,S39768,
+ e_S4732,e_S5244,e_S5756,e_S6268,e_S8232,e_S8744,e_S9256,
+ e_S9768,e_S17232,e_S17744,e_S18256,e_S18768,e_S23232,e_S23744,
+ e_S24256,e_S24768,e_S32232,e_S32744,e_S33256,e_S33768,
+ e_S38232,e_S38744,e_S39256,e_S39768,a,b,v0,S0,chi2
c ..............Just test output...........
write(6,1)
+ AT20G,RAhour,DEdeg,Date,fs,fe,l_S4732,S4732,l_S5244,S5244,
+ l_S5756,S5756,l_S6268,S6268,l_S8232,S8232,l_S8744,S8744,
+ l_S9256,S9256,l_S9768,S9768,l_S17232,S17232,l_S17744,S17744,
+ l_S18256,S18256,l_S18768,S18768,l_S23232,S23232,l_S23744,
+ S23744,l_S24256,S24256,l_S24768,S24768,l_S32232,S32232,
+ l_S32744,S32744,l_S33256,S33256,l_S33768,S33768,l_S38232,
+ S38232,l_S38744,S38744,l_S39256,S39256,l_S39768,S39768,
+ e_S4732,e_S5244,e_S5756,e_S6268,e_S8232,e_S8744,e_S9256,
+ e_S9768,e_S17232,e_S17744,e_S18256,e_S18768,e_S23232,e_S23744,
+ e_S24256,e_S24768,e_S32232,e_S32744,e_S33256,e_S33768,
+ e_S38232,e_S38744,e_S39256,e_S39768,a,b,v0,S0,chi2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end