Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/392/733 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/MNRAS/392/733 The NEWPS catalog from WMAP 5-yr survey (Massardi+, 2009)
*================================================================================
*Blind and non-blind source detection in WMAP 5-yr maps.
* Massardi M., Lopez-Caniego M., Gonzalez-Nuevo J., Herranz D., De Zotti G.,
* Sanz J.L.
* <Mon. Not. R. Astron. Soc., 392, 733-742 (2009)>
* =2009MNRAS.392..733M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! The NEWPS 5yr 5s Catalogue with a selection for
noise<1.5 times the median noise in the map
integer*4 nr__
parameter (nr__=516) ! Number of records
character*139 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 Seq ! Sequential number (1)
real*8 RAdeg ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg ! (deg) Declination in decimal degrees (J2000)
real*8 GLON ! (deg) Galactic longitude
real*8 GLAT ! (deg) Galactic latitude
integer*4 S23G ! (mJy) ?=0 Flux density at 23GHz
integer*4 S33G ! (mJy) ?=0 Flux density at 33GHz
integer*4 S41G ! (mJy) ?=0 Flux density at 41GHz
integer*4 S61G ! (mJy) ?=0 Flux density at 61GHz
integer*4 e_S23G ! (mJy) ?=0 rms uncertainty on S23GHz
integer*4 e_S33G ! (mJy) ?=0 rms uncertainty on S33GHz
integer*4 e_S41G ! (mJy) ?=0 rms uncertainty on S41GHz
integer*4 e_S61G ! (mJy) ?=0 rms uncertainty on S61GHz
real*4 N23 ! ?=- Ratio between the value of the noise map in
* the position of the source and the median
* noise in the map at 23GHz
real*4 N33 ! ?=- Ratio between the value of the noise map in
* the position of the source and the median
* noise in the map at 33GHz
real*4 N41 ! ?=- Ratio between the value of the noise map in
* the position of the source and the median
* noise in the map at 41GHz
real*4 N61 ! ?=- Ratio between the value of the noise map in
* the position of the source and the median
* noise in the map at 61GHz
character*1 f_S23G ! [BN.] Blind/non-blind detection at 23GHz
character*1 f_S33G ! [BN.] Blind/non-blind detection at 33GHz
character*1 f_S41G ! [BN.] Blind/non-blind detection at 41GHz
character*1 f_S61G ! [BN.] Blind/non-blind detection at 61GHz
character*3 fl1 ! [GLW.] Location flag (2)
character*1 fl2 ! [egu] Type flag (2)
real*8 z ! ?=- Redshift obtained from the NED database
* ('---' means that no value is available)
character*20 OName ! Identification in well known catalogues
* (blank means that no clear identification is
* available)
*Note (1): Sources identified as NEWPS55 DDD.dddd+DD.dddd in Simbad
*Note (2): Flags as follows:
* G = sources with 5{deg}<|b|<10{deg}
* L = sources within 5 degrees from the Large Magellanic Cloud
* (RA=5.3929 DE=-69.7561)
* W = objects within a beam unit (at 23GHz 21.35arcmin) from a
* WMAP_5yr source object identification
* e = extragalactic
* g = galactic
* u = not clearly identified objects
C=============================================================================
C Loading file 'table3.dat' ! The NEWPS 5yr 5s Catalogue with a selection for
* noise<1.5 times the median noise in the map
C Format for file interpretation
1 format(
+ I3,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,I5,1X,I5,1X,I5,1X,I5,1X,
+ I4,1X,I3,1X,I3,1X,I3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1,A1,
+ A1,A1,A3,A1,1X,F8.6,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,516
read(1,'(A139)')ar__
read(ar__,1)
+ Seq,RAdeg,DEdeg,GLON,GLAT,S23G,S33G,S41G,S61G,e_S23G,e_S33G,
+ e_S41G,e_S61G,N23,N33,N41,N61,f_S23G,f_S33G,f_S41G,f_S61G,fl1,
+ fl2,z,OName
if (idig(ar__(82:85)).EQ.0) N23 = rNULL__
if (idig(ar__(87:90)).EQ.0) N33 = rNULL__
if (idig(ar__(92:95)).EQ.0) N41 = rNULL__
if (idig(ar__(97:100)).EQ.0) N61 = rNULL__
if (idig(ar__(111:118)).EQ.0) z = rNULL__
c ..............Just test output...........
write(6,1)
+ Seq,RAdeg,DEdeg,GLON,GLAT,S23G,S33G,S41G,S61G,e_S23G,e_S33G,
+ e_S41G,e_S61G,N23,N33,N41,N61,f_S23G,f_S33G,f_S41G,f_S61G,fl1,
+ fl2,z,OName
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