Conversion of standardized ReadMe file for
catalog 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.5, on 2013-Jun-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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/MNRAS/416/559 PACO faint sources catalog (Bonavera+, 2011)
*================================================================================
*The Planck-ATCA coeval observations project: the faint sample.
* Bonavera L., Massardi M., Bonaldi A., Gonzalez-Nuevo J., De Zotti G.,
* Ekers R.D.
* <Mon. Not. R. Astron. Soc., 416, 559-566 (2011)>
* =2011MNRAS.416..559B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! The PACO faint catalogue
integer*4 nr__
parameter (nr__=674) ! Number of records
character*386 ar__ ! Full-size record
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 f1 ! [s-] Flag "s" to identify epochs within
* 10-days from Planck observations
character*1 f2 ! [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 ! ? Best fit parameter {nu}_0_ (1)
real*4 S0 ! (mJy) ? Best fit parameter S_0_ (1)
real*4 chi2 ! ? Best fit parameter {chi}^2^
*Note (1): Best fit parameters in the double power law:
* S({nu})=S_0_/[({nu}/{nu}_0_)^-a^ + ({nu}/{nu}_0_)^-n^], where
* {nu} is the frequency and S is the flux density in Jy.
C=============================================================================
C Loading file 'catalog.dat' ! The PACO faint 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,file='catalog.dat', status='old')
write(6,*) '....Loading file: catalog.dat'
do i__=1,674
read(1,'(A386)')ar__
read(ar__,1)
+ AT20G,RAhour,DEdeg,Date,f1,f2,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
if(ar__(60:64) .EQ. '') S4732 = iNULL__
if(ar__(66:70) .EQ. '') S5244 = iNULL__
if(ar__(72:76) .EQ. '') S5756 = iNULL__
if(ar__(78:82) .EQ. '') S6268 = iNULL__
if(ar__(84:88) .EQ. '') S8232 = iNULL__
if(ar__(90:94) .EQ. '') S8744 = iNULL__
if(ar__(96:100) .EQ. '') S9256 = iNULL__
if(ar__(102:106) .EQ. '') S9768 = iNULL__
if(ar__(108:112) .EQ. '') S17232 = iNULL__
if(ar__(114:118) .EQ. '') S17744 = iNULL__
if(ar__(120:124) .EQ. '') S18256 = iNULL__
if(ar__(126:130) .EQ. '') S18768 = iNULL__
if(ar__(132:136) .EQ. '') S23232 = iNULL__
if(ar__(138:142) .EQ. '') S23744 = iNULL__
if(ar__(144:148) .EQ. '') S24256 = iNULL__
if(ar__(150:154) .EQ. '') S24768 = iNULL__
if(ar__(156:160) .EQ. '') S32232 = iNULL__
if(ar__(162:166) .EQ. '') S32744 = iNULL__
if(ar__(168:172) .EQ. '') S33256 = iNULL__
if(ar__(174:178) .EQ. '') S33768 = iNULL__
if(ar__(180:184) .EQ. '') S38232 = iNULL__
if(ar__(186:190) .EQ. '') S38744 = iNULL__
if(ar__(192:196) .EQ. '') S39256 = iNULL__
if(ar__(198:202) .EQ. '') S39768 = iNULL__
if(ar__(204:208) .EQ. '') e_S4732 = rNULL__
if(ar__(210:214) .EQ. '') e_S5244 = rNULL__
if(ar__(216:220) .EQ. '') e_S5756 = rNULL__
if(ar__(222:226) .EQ. '') e_S6268 = rNULL__
if(ar__(228:232) .EQ. '') e_S8232 = rNULL__
if(ar__(234:238) .EQ. '') e_S8744 = rNULL__
if(ar__(240:244) .EQ. '') e_S9256 = rNULL__
if(ar__(246:250) .EQ. '') e_S9768 = rNULL__
if(ar__(252:256) .EQ. '') e_S17232 = rNULL__
if(ar__(258:262) .EQ. '') e_S17744 = rNULL__
if(ar__(264:268) .EQ. '') e_S18256 = rNULL__
if(ar__(270:274) .EQ. '') e_S18768 = rNULL__
if(ar__(276:280) .EQ. '') e_S23232 = rNULL__
if(ar__(282:286) .EQ. '') e_S23744 = rNULL__
if(ar__(288:292) .EQ. '') e_S24256 = rNULL__
if(ar__(294:298) .EQ. '') e_S24768 = rNULL__
if(ar__(300:304) .EQ. '') e_S32232 = rNULL__
if(ar__(306:310) .EQ. '') e_S32744 = rNULL__
if(ar__(312:316) .EQ. '') e_S33256 = rNULL__
if(ar__(318:322) .EQ. '') e_S33768 = rNULL__
if(ar__(324:328) .EQ. '') e_S38232 = rNULL__
if(ar__(330:334) .EQ. '') e_S38744 = rNULL__
if(ar__(336:340) .EQ. '') e_S39256 = rNULL__
if(ar__(342:346) .EQ. '') e_S39768 = rNULL__
if(ar__(349:354) .EQ. '') a = rNULL__
if(ar__(357:362) .EQ. '') b = rNULL__
if(ar__(365:370) .EQ. '') v0 = rNULL__
if(ar__(373:378) .EQ. '') S0 = rNULL__
if(ar__(382:386) .EQ. '') chi2 = rNULL__
c ..............Just test output...........
write(6,1)
+ AT20G,RAhour,DEdeg,Date,f1,f2,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