Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/392/733 into FORTRAN code for loading all data files into arrays.
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-17 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 (nr__) ! Sequential number (1) real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000) real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000) real*8 GLON (nr__) ! (deg) Galactic longitude real*8 GLAT (nr__) ! (deg) Galactic latitude integer*4 S23G (nr__) ! (mJy) ?=0 Flux density at 23GHz integer*4 S33G (nr__) ! (mJy) ?=0 Flux density at 33GHz integer*4 S41G (nr__) ! (mJy) ?=0 Flux density at 41GHz integer*4 S61G (nr__) ! (mJy) ?=0 Flux density at 61GHz integer*4 e_S23G (nr__) ! (mJy) ?=0 rms uncertainty on S23GHz integer*4 e_S33G (nr__) ! (mJy) ?=0 rms uncertainty on S33GHz integer*4 e_S41G (nr__) ! (mJy) ?=0 rms uncertainty on S41GHz integer*4 e_S61G (nr__) ! (mJy) ?=0 rms uncertainty on S61GHz real*4 N23 (nr__) ! ?=- 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 (nr__) ! ?=- 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 (nr__) ! ?=- 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 (nr__) ! ?=- 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 (nr__) ! [BN.] Blind/non-blind detection at 23GHz character*1 f_S33G (nr__) ! [BN.] Blind/non-blind detection at 33GHz character*1 f_S41G (nr__) ! [BN.] Blind/non-blind detection at 41GHz character*1 f_S61G (nr__) ! [BN.] Blind/non-blind detection at 61GHz character*3 fl1 (nr__) ! [GLW.] Location flag (2) character*1 fl2 (nr__) ! [egu] Type flag (2) real*8 z (nr__) ! ?=- Redshift obtained from the NED database * ('---' means that no value is available) character*20 OName (nr__) ! 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(i__),RAdeg(i__),DEdeg(i__),GLON(i__),GLAT(i__),S23G(i__), + S33G(i__),S41G(i__),S61G(i__),e_S23G(i__),e_S33G(i__), + e_S41G(i__),e_S61G(i__),N23(i__),N33(i__),N41(i__),N61(i__), + f_S23G(i__),f_S33G(i__),f_S41G(i__),f_S61G(i__),fl1(i__), + fl2(i__),z(i__),OName(i__) if (idig(ar__(82:85)).EQ.0) N23(i__) = rNULL__ if (idig(ar__(87:90)).EQ.0) N33(i__) = rNULL__ if (idig(ar__(92:95)).EQ.0) N41(i__) = rNULL__ if (idig(ar__(97:100)).EQ.0) N61(i__) = rNULL__ if (idig(ar__(111:118)).EQ.0) z(i__) = rNULL__ c ..............Just test output........... write(6,1) + Seq(i__),RAdeg(i__),DEdeg(i__),GLON(i__),GLAT(i__),S23G(i__), + S33G(i__),S41G(i__),S61G(i__),e_S23G(i__),e_S33G(i__), + e_S41G(i__),e_S61G(i__),N23(i__),N33(i__),N41(i__),N61(i__), + f_S23G(i__),f_S33G(i__),f_S41G(i__),f_S61G(i__),fl1(i__), + fl2(i__),z(i__),OName(i__) 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