Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/120/247 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 2026-Jun-08
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/ApJS/120/247 Supernova remnants in M33 (Gordon+, 1999)
*================================================================================
*A new sample of radio-selected and optically confirmed supernova remnants in M33
* Gordon S.M., Duric N., Kirshner R.P., Goss W.M., Viallefond F.
* <Astrophys. J. Suppl. Ser. 120, 247 (1999)>
* =1999ApJS..120..247G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Radio sources in M 33
integer*4 nr__
parameter (nr__=186) ! Number of records
character*68 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 v_GDK99_ (nr__) ! Identification number ([GDK99] NNN in Simbad)
integer*4 RAh (nr__) ! (h) Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Right ascension (J2000)
real*4 RAs (nr__) ! (s) Right ascension (J2000)
character*1 DE_ (nr__) ! Declination sign (J2000)
integer*4 DEd (nr__) ! (deg) Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Declination (J2000)
real*4 S20cm (nr__) ! (mJy) Total flux at 20cm (1.5GHz)
real*4 e_S20cm (nr__) ! (mJy) rms uncertainty on S20cm
character*1 l_S6cm (nr__) ! Limit flag on S6cm
real*4 S6cm (nr__) ! (mJy) Total flux at 6cm (5GHz)
real*4 e_S6cm (nr__) ! (mJy) ? rms uncertainty on S6cm
character*1 l_Sp_Index (nr__) ! Limit flag on Sp+Index
real*4 Sp_Index (nr__) ! Spectral index of the source
real*4 e_Sp_Index (nr__) ! ? rms uncertainty on Sp+Index
character*9 OptID (nr__) ! Optical identification (1)
*Note (1): the symbols have the following meaning:
* OI denotes the radiosurces that where outside the ield of view of
* the optical image
* SNR denotes the confirmed SNR
* SNRc denotes the SNR candidate
* H II denotes the H II region
* BGND-NOE denotes the background, no optical emission
* BGND-FOE denotes the background, faint optical emission
* BGND-STR denotes the background, stellar-like source
* GC denotes the galactic center
* * represent the H II regions with nonthermal spectral indices
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Radio-selected sample of supernova remnants
integer*4 nr__1
parameter (nr__1=53) ! Number of records
character*67 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 v_GDK99__1 (nr__1) ! Identification number ([GDK99] NNN in Simbad)
integer*4 M33SNR (nr__1) ! Gordon et al. (1998ApJS..117...89G)
* identification number (M33SNR NN in Simbad)
integer*4 RAh_1 (nr__1) ! (h) right ascension (J2000)
integer*4 RAm_1 (nr__1) ! (min) right ascension (J2000)
real*4 RAs_1 (nr__1) ! (s) right ascension (J2000)
character*1 DE__1 (nr__1) ! Declination sign (J2000)
integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000)
integer*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000)
integer*4 Diam (nr__1) ! (pc) Diameter (1)
character*1 n_Diam (nr__1) ! [*] Note on Diam (1)
character*1 l_S20cm (nr__1) ! Limit flag on S20cm
real*4 S20cm_1 (nr__1) ! (mJy) Total flux at 20cm (1.5GHz)
real*4 e_S20cm_1 (nr__1) ! (mJy) ? rms uncertainty on S20cm
character*1 l_S6cm_1 (nr__1) ! Limit flag on S6cm
real*4 S6cm_1 (nr__1) ! (mJy) Total flux at 6cm (5GHz)
real*4 e_S6cm_1 (nr__1) ! (mJy) ? rms uncertainty on S6cm
character*1 l_Sp_Index_1(nr__1) ! Limit flag on Sp+Index
real*4 Sp_Index_1 (nr__1) ! Spectral index
real*4 e_Sp_Index_1(nr__1) ! ? rms uncertainty on Sp+Index
character*1 Tc (nr__1) ! [x] Thermal correction (2)
character*4 SpConf (nr__1) ! Spectral confirmation (3)
*Note (1): Diameters as measured from optical images.
* Geometric mean for noncircular objects. * indicates that measurement
* is highly uncertain due to optical confusion. At 840 kpc, 1"=4.08pc
*Note (2): x indicates that the radio fluxes are corrected for thermal
* contamination
*Note (3): Sources for optical spectroscopy:
* G98: Gordon et al., 1998ApJS..117...89G (<M33SNR NN>)
* S93: Smith et al., 1993ApJ...407..564S
* BK85: Blair & Kirshner, 1985ApJ...289..582B
C=============================================================================
C Loading file 'table1.dat' ! Radio sources in M 33
C Format for file interpretation
1 format(
+ I3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.1,1X,F3.1,
+ 1X,A1,F4.1,1X,F3.1,1X,A1,F4.1,1X,F3.1,1X,A9)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,186
read(1,'(A68)')ar__
read(ar__,1)
+ v_GDK99_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),S20cm(i__),e_S20cm(i__),l_S6cm(i__),
+ S6cm(i__),e_S6cm(i__),l_Sp_Index(i__),Sp_Index(i__),
+ e_Sp_Index(i__),OptID(i__)
if(ar__(46:48) .EQ. '') e_S6cm(i__) = rNULL__
if(ar__(56:58) .EQ. '') e_Sp_Index(i__) = rNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ v_GDK99_(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),S20cm(i__),e_S20cm(i__),l_S6cm(i__),
+ S6cm(i__),e_S6cm(i__),l_Sp_Index(i__),Sp_Index(i__),
+ e_Sp_Index(i__),OptID(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Radio-selected sample of supernova remnants
C Format for file interpretation
2 format(
+ I3,1X,I2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,I3,A1,1X,
+ A1,F3.1,1X,F3.1,1X,A1,F3.1,1X,F3.1,1X,A1,F4.1,1X,F3.1,1X,A1,
+ 1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,53
read(1,'(A67)')ar__1
read(ar__1,2)
+ v_GDK99__1(i__),M33SNR(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Diam(i__),
+ n_Diam(i__),l_S20cm(i__),S20cm_1(i__),e_S20cm_1(i__),
+ l_S6cm_1(i__),S6cm_1(i__),e_S6cm_1(i__),l_Sp_Index_1(i__),
+ Sp_Index_1(i__),e_Sp_Index_1(i__),Tc(i__),SpConf(i__)
if(ar__1(39:41) .EQ. '') e_S20cm_1(i__) = rNULL__
if(ar__1(48:50) .EQ. '') e_S6cm_1(i__) = rNULL__
if(ar__1(58:60) .EQ. '') e_Sp_Index_1(i__) = rNULL__
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,2)
+ v_GDK99__1(i__),M33SNR(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Diam(i__),
+ n_Diam(i__),l_S20cm(i__),S20cm_1(i__),e_S20cm_1(i__),
+ l_S6cm_1(i__),S6cm_1(i__),e_S6cm_1(i__),l_Sp_Index_1(i__),
+ Sp_Index_1(i__),e_Sp_Index_1(i__),Tc(i__),SpConf(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end