Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/104/704 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-20
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/AJ/104/704 A VLA Search for Young Galactic Supernova Remnants (Sramek+, 1992)
*================================================================================
*A VLA Search for Young Galactic Supernova Remnants
* Sramek R.A., Cowen J.J., Roberts D.A., Goss W.M., Ekers R.D.
* <Astron. J. 104, 704 (1992)>
* =1992AJ....104..704S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! High resolution galactic plane radio SNR survey
integer*4 nr__
parameter (nr__=291) ! Number of records
character*52 ar__ ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension 1950
real*8 DEdeg ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
real*4 GLON ! (deg) Galactic Longitude (l2)
real*4 GLAT ! (deg) Galactic Latitude (b2)
integer*4 RAh ! (h) *Right Ascension hours (1950)
integer*4 RAm ! (min) *Right Ascension minutes (1950)
real*4 RAs ! (s) *Right Ascension seconds (1950)
character*1 DE_ ! *[+- ]Declination sign (1950)
integer*4 DEd ! (deg) *Declination degrees (1950)
integer*4 DEm ! (arcmin) *Declination minutes (1950)
integer*4 DEs ! (arcsec) *Declination seconds (1950)
character*1 l_s20cm ! *[ <]Limit flag for 20cm flux
real*4 s20cm ! (Jy) *Peak 20cm flux density in 3'x3' box
character*1 struct ! *[PSDCE]Radio structure of dominant source
character*2 obs ! *Epoch of observation
character*1 ref ! *Reference of original survey
character*2 newRAs ! (s) *?New Right Ascension (seconds only)
real*4 newDEm ! (arcmin) *?New Declination (minutes only)
character*1 note ! *?Notes
*Note on RAh,RAm,RAs,DE-,DEd,DEm,DEs:
* Right Ascension and Declination (1950.0) as listed in the original survey
* catalog (see Ref column for survey reference)
*Note on l_s20cm:
* < indicates s20cm is approximate 5 sigma upper limit for fields in which the
* source is resolved out.
*Note on s20cm:
* Peak 20cm flux density on VLA image in a 3.0x3.0 arcmin box centered on the
* catalogued position. For double and complex sources the peak flux density
* refers only to the brightest component. No attempt was made to integrate
* flux densities. See the notes for related information.
*Note on struct:
* The radio structure of the dominant high brightness compact source. The
* symbols' meanings are as follows
* P = point sources or slightly resolved
* S = single extended source
* D = double source
* C = complex structure with several knots and/or filaments
* E = empty field - the catalog source is probably resolved out
*Note on obs:
* "A" observations were toward the galactic anticenter;
* "B" observations were toward the galactic center
* A1 = 29 Nov. 1984
* A2 = 30 Nov. 1984
* B1 = 7 Dec. 1984
* B2 = 24 Dec. 1984
* B3 = 28 Dec. 1984
*Note on ref:
* 1 = Clark & Crawford 1974
* 2 = Haynes+ (1979)
* 3 = Altenhoff+ (1978)
*Note on newRAs:
* Improved right ascension seconds of time, if the peak of the compact
* structure is more than 0.5 arcmin from the catalogued position
*Note on newDEm:
* Improved declination arcminutes, if the peak of the compact structure
* is more than 0.5 arcmin from the catalogued position
*Note on note:
* A = improved positions and 6cm flux densities are given in Table 2
* B = low resolution, wide field images were made at 20cm; see Table 3
* 1 = diameter about 45 arcsec
* 2 = other components about 30 arcsec to SE; 7 and 5 mJy peak
* 3 = other components about 30 arcsec to S; 5 mJy
* 4 = listed source is same as listed at l=21.05
* 5 = extended source 4 mJy; peak 45 arcsec to N; another 90 arcsec to NE
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! 6cm observations of apparent shell structures
integer*4 nr__1
parameter (nr__1=14) ! Number of records
character*65 ar__1 ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension 1950
real*8 DEdeg_1 ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
real*4 GLON_1 ! (deg) *Galactic Longitude (l2)
real*4 GLAT_1 ! (deg) *Galactic Latitude (b2)
integer*4 RAh_1 ! (h) *?Right Ascension hours (1950)
integer*4 RAm_1 ! (min) *?Right Ascension minutes (1950)
real*4 RAs_1 ! (s) *?Right Ascension seconds (1950)
character*1 DE__1 ! *[+- ]Declination sign (1950)
integer*4 DEd_1 ! (deg) *?Declination degrees (1950)
integer*4 DEm_1 ! (arcmin) *?Declination minutes (1950)
integer*4 DEs_1 ! (arcsec) *?Declination seconds (1950)
real*4 s6cm ! (Jy) Peak 6cm flux density
real*4 SpIndex ! *Spectral index between 6 and 20cm
character*22 MType ! *?Description of morphology of source
*Note on GLON,GLAT:
* Galactic longitude and galactic latitude as derived from the original survey
* catalog
*Note on RAh,RAm,RAs,DE-,DEd,DEm,DEs:
* Right Ascension and Declination (1950.0) of the peak brightness of the source
* as determined by the 6cm observations
*Note on SpIndex:
* Spectral index between 6 and 20 cm wavelength using peak flux densities
* given in Table 1 and Table 2.
*Note on MType:
* Description of the source morphology, e.g. Shell, Point Source
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Low resolution galactic plane radio SNR survey
integer*4 nr__2
parameter (nr__2=62) ! Number of records
character*41 ar__2 ! Full-size record
C 1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 ! (deg) Right Ascension 1950
real*8 DEdeg_2 ! (deg) Declination 1950
C ---------------------------------- ! (position vector(s) in degrees)
real*4 GLON_2 ! (deg) *Galactic Longitude (l2)
real*4 GLAT_2 ! (deg) *Galactic Latitude (b2)
integer*4 RAh_2 ! (h) *?Right Ascension hours (1950)
integer*4 RAm_2 ! (min) *?Right Ascension minutes (1950)
real*4 RAs_2 ! (s) *?Right Ascension seconds (1950)
character*1 DE__2 ! *[+- ]Declination sign (1950)
integer*4 DEd_2 ! (deg) *?Declination degrees (1950)
integer*4 DEm_2 ! (arcmin) *?Declination minutes (1950)
integer*4 DEs_2 ! (arcsec) *?Declination seconds (1950)
character*1 l_s20cm_1 ! *[ <]Limit flag for 20cm flux
real*4 s20cm_1 ! (Jy) *Peak 20cm flux density in 3'x3' box
character*1 struct_1 ! *[PSDCE]Radio structure of dominant source
character*1 note1 ! *?Notes on general structure
integer*4 note2 ! *?Notes on individual sources; "listed
* source described in the table; other
* comments refer to additional sources in
* the image. Positions and flux densities
* are for the brightness maxima.
*Note on GLON,GLAT:
* Galactic longitude and galactic latitude as derived from the original survey
* catalog
*Note on RAh,RAm,RAs,DE-,DEd,DEm,DEs:
* Right Ascension and Declination (1950.0) of the peak brightness of the source
* as determined by the low resolution VLA 20cm observations
*Note on l_s20cm:
* < indicates s20cm is approximate 5 sigma upper limit for fields that remain
* completely resolved.
*Note on s20cm:
* Peak 20cm flux density. For double and complex sources the peak flux density
* refers only to the brightest component. No attempt was made to integrate
* flux densities. See the notes for related information.
*Note on struct:
* Classification of the radio structure of the dominant high brightness
* compact source. The symbols' meanings are as follows:
* P = point sources or slightly resolved
* S = single extended source
* D = double source
* C = complex structure with several knots and/or filaments
* E = empty field - the catalog source is probably resolved out
*Note on note1:
* A = Extended low brightness emission in addition to the compact
* structure indicated by the classification.
* B = Possible shell structure
*Note on note2:
* 1 = compact source, 1.40 Jy peak, at RA=18h00m37s, DE=24d22.9'
* 2 = double source, 0.25 Jy peak, at RA=18h03m19s, DE=-21d38'
* 3 = single, 0.162 Jy, at RA=18h06m46s, DE=-19d27.0'
* 4 = other weak knots nearby
* 5 = point source, 0.58 Jy, at RA=18h31m41s, DE=-07d57.1'
* 6 = point source, 0.15 Jy, at RA=18h37m57s
* 7 = another point source. 0.12 Jy, at RA=18h40m19s, DE=-04d17'
* 8 = listed double 30" separation, well resolved components.
* 9 = correction applied, FWZI about 1.1 arcmin.
* 10 = G29.7-0.3 is 15 arcmin to NW
* 11 = point source, 1.17 Jy, at RA=18h47m23s, DE=-01d36.4'
* 12 = point source, 0.43 Jy, at RA=18h44m59s, DE=-01d16.0'
* 13 = double source, 0.17 Jy, at RA=18h49m34s, DE=-00d04.5'
* 14 = point source, 0.53 Jy, at RA=18h49m14s, DE=-00d31.9'
* 15 = extended emission region 3' FWZI, offset 4'S of point source
* 16 = point source, 0.32 Jy, at RA=18h57m47s, DE=-03d59.0'
* tailed (shell?) source, .26 Jy peak, at RA 18h58m33s, DE=04d07.7'
* 17 = extended emission region 4' FWZI, offset 8' N of catalogued position
* 18 = listed double source has 0.5' separation
* 19 = listed double source has 0.6' separation
* 20 = listed double source has 1.5' separation
* 21 = point source, 0.23 Jy, at RA=07h18m27s, DE=-15d13.9'
* listed double source has 20" separation
* 22 = listed source has extended, 2.3' FWZI
* 23 = listed source has diameter 45" FWZI
* 24 = listed double source has 30" separation
* 25 = listed double source has 45" separation
* 26 = listed double source has 20" separation
* 27 = extended, partial shell source, 1.04 Jy peak, at RA=17h16m04s,
* DE=-37d07.8'
* 28 = extended, complex source, 1.1 Jy peak, centered on RA=17h17m30s,
* DE=-35d43'
* 29 = extended, complex source, 2.4 Jy peak, centered on RA=17h21m27s,
* DE=-34d08'
* 30 = broad structure with many knots and filaments seen on image;
* heavily resolved
* 31 = partial shell source, 0.15 Jy peak, at RA=17h28m11s, DE=-33d31.5',
* 2.5' diameter FWZI
* 32 = very extended diffuse emission over 8' diameter area
* 33 = heavily confused by SGR A * about 20' to north
C=============================================================================
C Loading file 'table1.dat' ! High resolution galactic plane radio SNR survey
C Format for file interpretation
1 format(
+ F6.2,F5.2,1X,I2,I2,F4.1,A1,I2,I2,I2,1X,A1,F5.3,1X,A1,1X,A2,1X,
+ A1,2X,A2,1X,F4.1,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,291
read(1,'(A52)')ar__
read(ar__,1)
+ GLON,GLAT,RAh,RAm,RAs,DE_,DEd,DEm,DEs,l_s20cm,s20cm,struct,
+ obs,ref,newRAs,newDEm,note
if(ar__(47:50) .EQ. '') newDEm = rNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ GLON,GLAT,RAh,RAm,RAs,DE_,DEd,DEm,DEs,l_s20cm,s20cm,struct,
+ obs,ref,newRAs,newDEm,note
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! 6cm observations of apparent shell structures
C Format for file interpretation
2 format(
+ F6.2,F5.2,1X,I2,I2,F4.1,A1,I2,I2,I2,2X,F5.3,2X,F5.2,2X,A22)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,14
read(1,'(A65)')ar__1
read(ar__1,2)
+ GLON_1,GLAT_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,s6cm,
+ SpIndex,MType
if(ar__1(13:14) .EQ. '') RAh_1 = iNULL__
if(ar__1(15:16) .EQ. '') RAm_1 = iNULL__
if(ar__1(17:20) .EQ. '') RAs_1 = rNULL__
if(ar__1(22:23) .EQ. '') DEd_1 = iNULL__
if(ar__1(24:25) .EQ. '') DEm_1 = iNULL__
if(ar__1(26:27) .EQ. '') DEs_1 = iNULL__
RAdeg_1 = rNULL__
DEdeg_1 = rNULL__
c Derive coordinates RAdeg_1 and DEdeg_1 from input data
c (RAdeg_1 and DEdeg_1 are set to rNULL__ when unknown)
if(RAh_1 .GT. -180) RAdeg_1=RAh_1*15.
if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4.
if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240.
if(DEd_1 .GE. 0) DEdeg_1=DEd_1
if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60.
if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1
c ..............Just test output...........
write(6,2)
+ GLON_1,GLAT_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,s6cm,
+ SpIndex,MType
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Low resolution galactic plane radio SNR survey
C Format for file interpretation
3 format(
+ F6.2,F5.2,1X,I2,I2,F4.1,A1,I2,I2,I2,1X,A1,F5.3,1X,A1,1X,A1,1X,
+ I2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,62
read(1,'(A41)')ar__2
read(ar__2,3)
+ GLON_2,GLAT_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,
+ l_s20cm_1,s20cm_1,struct_1,note1,note2
if(ar__2(13:14) .EQ. '') RAh_2 = iNULL__
if(ar__2(15:16) .EQ. '') RAm_2 = iNULL__
if(ar__2(17:20) .EQ. '') RAs_2 = rNULL__
if(ar__2(22:23) .EQ. '') DEd_2 = iNULL__
if(ar__2(24:25) .EQ. '') DEm_2 = iNULL__
if(ar__2(26:27) .EQ. '') DEs_2 = iNULL__
if(ar__2(40:41) .EQ. '') note2 = iNULL__
RAdeg_2 = rNULL__
DEdeg_2 = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2 .GT. -180) RAdeg_2=RAh_2*15.
if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4.
if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240.
if(DEd_2 .GE. 0) DEdeg_2=DEd_2
if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60.
if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600.
if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c ..............Just test output...........
write(6,3)
+ GLON_2,GLAT_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,
+ l_s20cm_1,s20cm_1,struct_1,note1,note2
write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end