Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/76/317 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-May-24
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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/A+AS/76/317 Possible new planetary nebulae in IRAS PSC (Preite-Martinez, 1988)
*================================================================================
*Possible new planetary nebulae in the IRAS point source catalogue.
* Preite-Martinez A.
* <Astron. Astrophys. Suppl. Ser., 76, 317-330 (1988)>
* =1988A&AS...76..317P
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Possible new planetary nebulae in the
IRAS Point Source Catalogue
integer*4 nr__
parameter (nr__=340) ! Number of records
character*108 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension B1950
real*8 Dec ! (deg) Declination B1950
integer*4 L ! [1,2] List: 1=Possible new PN (table4),
* 2=Suspected non-PN (table5)
integer*4 Seq ! Sequential number (2)
character*11 IRAS ! IRAS name (HHMMm+DDMMA, B1950)
integer*4 RAh ! (h) Right ascension (B1950)
integer*4 RAm ! (min) Right ascension (B1950)
real*4 RAs ! (s) Right ascension (B1950)
character*1 DE_ ! Declination sign (B1950)
integer*4 DEd ! (deg) Declination (B1950)
integer*4 DEm ! (arcmin) Declination (B1950)
integer*4 DEs ! (arcsec) Declination (B1950)
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
real*4 F12 ! (Jy) IRAS/12{mu}m band flux as in IRAS PSC
real*4 F25 ! (Jy) IRAS/25{mu}m band flux as in IRAS PSC
real*4 F60 ! (Jy) IRAS/60{mu}m band flux as in IRAS PSC
real*4 F100 ! (Jy) IRAS/100{mu}m band flux as in IRAS PSC
character*4 Qual ! [123] Quality flags (3=good, 1=upper limit),
* 1 per band (as in IRAS PSC Cats)
integer*4 Tdust ! (K) Dust temperature (assuming the dust is
* the unique source of the observed FIR)
integer*4 FIR ! (10-14W/m2) Far-IR flux
real*4 logFHb ! ([mW/m2]) ? Log of estimated H-beta flux (1)
integer*4 S5GHz ! (mJy) ? Estimated radio continuum flux at 5GHz (1)
real*4 Dist ! (kpc) ? Estimated heliocentric distance (1)
*Note (1): Only in table5.
*Note (2): <PN PM 1-NNN> and <PN PM 2-NN> in Simbad for new possible PN
* and suspected non-PN respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Suspected non-PN in the selected sample of
possible PN
integer*4 nr__1
parameter (nr__1=48) ! Number of records
character*108 ar__1 ! Full-size record
real*8 RA__1 ! (deg) Right Ascension B1950
real*8 Dec_1 ! (deg) Declination B1950
integer*4 L_1 ! [1,2] List: 1=Possible new PN (table4),
* 2=Suspected non-PN (table5)
integer*4 Seq_1 ! Sequential number (2)
character*11 IRAS_1 ! IRAS name (HHMMm+DDMMA, B1950)
integer*4 RAh_1 ! (h) Right ascension (B1950)
integer*4 RAm_1 ! (min) Right ascension (B1950)
real*4 RAs_1 ! (s) Right ascension (B1950)
character*1 DE__1 ! Declination sign (B1950)
integer*4 DEd_1 ! (deg) Declination (B1950)
integer*4 DEm_1 ! (arcmin) Declination (B1950)
integer*4 DEs_1 ! (arcsec) Declination (B1950)
real*4 GLON_1 ! (deg) Galactic longitude
real*4 GLAT_1 ! (deg) Galactic latitude
real*4 F12_1 ! (Jy) IRAS/12{mu}m band flux as in IRAS PSC
real*4 F25_1 ! (Jy) IRAS/25{mu}m band flux as in IRAS PSC
real*4 F60_1 ! (Jy) IRAS/60{mu}m band flux as in IRAS PSC
real*4 F100_1 ! (Jy) IRAS/100{mu}m band flux as in IRAS PSC
character*4 Qual_1 ! [123] Quality flags (3=good, 1=upper limit),
* 1 per band (as in IRAS PSC Cats)
integer*4 Tdust_1 ! (K) Dust temperature (assuming the dust is
* the unique source of the observed FIR)
integer*4 FIR_1 ! (10-14W/m2) Far-IR flux
real*4 logFHb_1 ! ([mW/m2]) ? Log of estimated H-beta flux (1)
integer*4 S5GHz_1 ! (mJy) ? Estimated radio continuum flux at 5GHz (1)
real*4 Dist_1 ! (kpc) ? Estimated heliocentric distance (1)
*Note (1): Only in table5.
*Note (2): <PN PM 1-NNN> and <PN PM 2-NN> in Simbad for new possible PN
* and suspected non-PN respectively.
C=============================================================================
C Loading file 'table4.dat' ! Possible new planetary nebulae in the
* IRAS Point Source Catalogue
C Format for file interpretation
1 format(
+ 3X,I1,1X,I3,1X,A11,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,
+ 1X,F5.1,1X,F5.1,1X,F5.2,1X,F6.2,1X,F6.2,1X,F5.1,1X,A4,1X,I3,
+ 1X,I4,1X,F6.2,1X,I3,1X,F3.1)
C Effective file loading
open(unit=1,file='table4.dat', status='old')
write(6,*) '....Loading file: table4.dat'
do i__=1,340
read(1,'(A108)')ar__
read(ar__,1)
+ L,Seq,IRAS,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,F12,F25,F60,
+ F100,Qual,Tdust,FIR,logFHb,S5GHz,Dist
if(ar__(95:100) .EQ. '') logFHb = rNULL__
if(ar__(102:104) .EQ. '') S5GHz = iNULL__
if(ar__(106:108) .EQ. '') Dist = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ L,Seq,IRAS,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,F12,F25,F60,
+ F100,Qual,Tdust,FIR,logFHb,S5GHz,Dist
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Suspected non-PN in the selected sample of
* possible PN
C Format for file interpretation
2 format(
+ 3X,I1,1X,I3,1X,A11,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,
+ 1X,F5.1,1X,F5.1,1X,F5.2,1X,F6.2,1X,F6.2,1X,F5.1,1X,A4,1X,I3,
+ 1X,I4,1X,F6.2,1X,I3,1X,F3.1)
C Effective file loading
open(unit=1,file='table5.dat', status='old')
write(6,*) '....Loading file: table5.dat'
do i__=1,48
read(1,'(A108)')ar__1
read(ar__1,2)
+ L_1,Seq_1,IRAS_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ GLON_1,GLAT_1,F12_1,F25_1,F60_1,F100_1,Qual_1,Tdust_1,FIR_1,
+ logFHb_1,S5GHz_1,Dist_1
if(ar__1(95:100) .EQ. '') logFHb_1 = rNULL__
if(ar__1(102:104) .EQ. '') S5GHz_1 = iNULL__
if(ar__1(106:108) .EQ. '') Dist_1 = rNULL__
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1 = RAh_1
if(RAh_1 .GE. 0) RA__1=RAh_1*15.
if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
Dec_1 = DEd_1
if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
if(DEs_1 .GE. 0) Dec_1=Dec_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c ..............Just test output...........
write(6,2)
+ L_1,Seq_1,IRAS_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ GLON_1,GLAT_1,F12_1,F25_1,F60_1,F100_1,Qual_1,Tdust_1,FIR_1,
+ logFHb_1,S5GHz_1,Dist_1
write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end