Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/505/1049 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/A+A/505/1049 Parkes HI observations behind Milky Way II (Schroeder+, 2009)
*================================================================================
*Parkes HI observations of galaxies behind the southern Milky Way.
*II. The Crux and Great Attractor regions (l~289{deg} to 338{deg}).
* Schroeder A.C., Kraan-Korteweg R.C., Henning P.A.
* <Astron. Astrophys. 505, 1049 (2009)>
* =2009A&A...505.1049S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! HI-detections in the Crux and Great Attractor region
integer*4 nr__
parameter (nr__=162) ! Number of records
character*154 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*11 Name ! Source name (WKKNNNN, ICNNNN or YNNN-N)
character*3 f_Name ! Uncertain identification flag (1)
character*13 OName ! Other name
character*2 IR ! [MID ] Identification in the infrared (IR) and
* near-infrared (NIR) (2)
integer*4 RAh ! (h) ? Right ascension (J2000)
integer*4 RAm ! (min) ? Right ascension (J2000)
real*4 RAs ! (s) ? Right ascension (J2000)
character*1 DE_ ! Declination sign (J2000)
integer*4 DEd ! (deg) ? Declination (J2000)
integer*4 DEm ! (arcmin) ? Declination (J2000)
integer*4 DEs ! (arcsec) [0/60]? Declination (J2000)
character*1 u_DEs ! Uncertainty flag on Position
real*4 GLON ! (deg) ? Galactic longitude
real*4 GLAT ! (deg) ? Galactic latitude
character*5 MType ! Morphological type
integer*4 D ! (arcsec) ? Major axis diameter
integer*4 d_1 ! (arcsec) ? Minor axis diameter
real*4 Bj ! (mag) ? Apparent B_J_ magnitude
character*1 n_Bj ! [NR] Unexplained meaning
real*4 E_B_V ! (mag) Galactic reddening at the position
* of the galaxy
character*1 l_HRV ! Limit flag on HRV
integer*4 HRV ! (km/s) Heliocentric HI radial velocity (3)
character*1 u_HRV ! Uncertainty flag on HRV
integer*4 DV50 ! (km/s) ? Velocity width of the HI profile measured
* at the 50% level of the peak intensity
character*1 u_DV50 ! Uncertainty flag on DV50
integer*4 DV20 ! (km/s) ? Velocity width of the HI profile measured
* at the 20% level of the peak intensity
character*1 u_DV20 ! Uncertainty flag on DV20
real*4 Sint ! (Jy.km/s) ? Integrated HI flux density uncorrected
* for finite beam size
character*1 u_Sint ! Uncertainty flag on integrated HI flux
real*4 rms ! ? rms noise level
character*1 hann ! [n] n for no Hanning-smoothing (4)
character*1 n_Name ! [*] * indicates a note (5)
real*4 Dist ! (arcmin) ? Angular distance of the detected galaxy from
* the centre of the beam in arcminutes
real*4 Sc ! (Jy.km/s) ? Corrected flux density for off-centre
* detection
character*14 RFI ! Excised RFI (radio frequency interference) on
* or near the detected HI profile
*Note (1): Note as follows:
* ? = uncertain identification of the HI signal
* + = more than one signal was found in the pointing or in the
* associated OFF-observation
*Note (2): Flag as follows:
* I = indicates an entry in the IRAS Point Source Catalog (II/125)
* according to the precepts explained in Woudt & Kraan-Korteweg
* (2001, Cat. J/A+A/380/441)
* M = an entry in the 2MASS Extended Objects Catalogue (2MASX, 2003,
* Cat. VII/233)
* D = indicates an entry in the DENIS catalogue by Vauglin et al.
* (2002, Cat. J/A+A/387/1
*Note (3): Heliocentric HI radial velocity taken at the midpoint of the
* HI profile at the 20% level.
*Note (4): Most spectra have been Hanning-smoothed, except when the line
* width was smaller than 100km/s or for other reasons (see below),
* which is indicated with an 'n' for 'no Hanning-smoothing'.
*Note (5): Notes on individual sources as follows:
* WKK0969 = see plot of WKK1117
* WKK2576 = the signal is blended with WKK2595 and WKK2597 (cf. text),
* the width and velocity of the narrow peak are given here
* WKK2597 = see plot of WKK2595
* WKK2863 = see plot of WKK2844
* WKK4016 = the signal is blended with WKK4022, the velocity
* measurements refer to the narrow peak at v~4640km/s only
* WKK5299 = the parameters are uncertain due to an RFI at the edge of
* the profile at v~4500km/s
* WKK5616 = see plot of WKK5659
* WKK5636 = Hanning-smoothed once despite the narrow line, due to
* ringing in the spectrum
* WKK5829 = not Hanning-smoothed due to an RFI close to the profile
* at v=4550km/s
* WKK5993 = the profile of WKK5999 (see separate pointing) has been
* excluded from the profile measurement
* WKK6680 = once Hanning-smoothed for the plot
* SKK209-5 = the RFI is next to the profile and possibly affects the line width
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb2.dat' ! HI-detections in the Vela region
integer*4 nr__1
parameter (nr__1=14) ! Number of records
character*154 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 ! (deg) Right Ascension J2000
real*8 DEdeg_1 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*11 Name_1 ! Source name (WKKNNNN, ICNNNN or YNNN-N)
character*3 f_Name_1 ! Uncertain identification flag (1)
character*13 OName_1 ! Other name
character*2 IR_1 ! [MID ] Identification in the infrared (IR) and
* near-infrared (NIR) (2)
integer*4 RAh_1 ! (h) ? Right ascension (J2000)
integer*4 RAm_1 ! (min) ? Right ascension (J2000)
real*4 RAs_1 ! (s) ? Right ascension (J2000)
character*1 DE__1 ! Declination sign (J2000)
integer*4 DEd_1 ! (deg) ? Declination (J2000)
integer*4 DEm_1 ! (arcmin) ? Declination (J2000)
integer*4 DEs_1 ! (arcsec) [0/60]? Declination (J2000)
character*1 u_DEs_1 ! Uncertainty flag on Position
real*4 GLON_1 ! (deg) ? Galactic longitude
real*4 GLAT_1 ! (deg) ? Galactic latitude
character*5 MType_1 ! Morphological type
integer*4 D_2 ! (arcsec) ? Major axis diameter
integer*4 d_3 ! (arcsec) ? Minor axis diameter
real*4 Bj_1 ! (mag) ? Apparent B_J_ magnitude
character*1 n_Bj_1 ! [NR] Unexplained meaning
real*4 E_B_V_1 ! (mag) Galactic reddening at the position
* of the galaxy
character*1 l_HRV_1 ! Limit flag on HRV
integer*4 HRV_1 ! (km/s) Heliocentric HI radial velocity (3)
character*1 u_HRV_1 ! Uncertainty flag on HRV
integer*4 DV50_1 ! (km/s) ? Velocity width of the HI profile measured
* at the 50% level of the peak intensity
character*1 u_DV50_1 ! Uncertainty flag on DV50
integer*4 DV20_1 ! (km/s) ? Velocity width of the HI profile measured
* at the 20% level of the peak intensity
character*1 u_DV20_1 ! Uncertainty flag on DV20
real*4 Sint_1 ! (Jy.km/s) ? Integrated HI flux density uncorrected
* for finite beam size
character*1 u_Sint_1 ! Uncertainty flag on integrated HI flux
real*4 rms_1 ! ? rms noise level
character*1 hann_1 ! [n] n for no Hanning-smoothing (4)
character*1 n_Name_1 ! [*] * indicates a note (5)
real*4 Dist_1 ! (arcmin) ? Angular distance of the detected galaxy from
* the centre of the beam in arcminutes
real*4 Sc_1 ! (Jy.km/s) ? Corrected flux density for off-centre
* detection
character*14 RFI_1 ! Excised RFI (radio frequency interference) on
* or near the detected HI profile
*Note (1): Note as follows:
* ? = uncertain identification of the HI signal
* + = more than one signal was found in the pointing or in the
* associated OFF-observation
*Note (2): Flag as follows:
* I = indicates an entry in the IRAS Point Source Catalog (II/125)
* according to the precepts explained in Woudt & Kraan-Korteweg
* (2001, Cat. J/A+A/380/441)
* M = an entry in the 2MASS Extended Objects Catalogue (2MASX, 2003,
* Cat. VII/233)
* D = indicates an entry in the DENIS catalogue by Vauglin et al.
* (2002, Cat. J/A+A/387/1
*Note (3): Heliocentric HI radial velocity taken at the midpoint of the
* HI profile at the 20% level.
*Note (4): Most spectra have been Hanning-smoothed, except when the line
* width was smaller than 100km/s or for other reasons (see below),
* which is indicated with an 'n' for 'no Hanning-smoothing'.
*Note (5): Notes on individual sources as follows:
* WKK0969 = see plot of WKK1117
* WKK2576 = the signal is blended with WKK2595 and WKK2597 (cf. text),
* the width and velocity of the narrow peak are given here
* WKK2597 = see plot of WKK2595
* WKK2863 = see plot of WKK2844
* WKK4016 = the signal is blended with WKK4022, the velocity
* measurements refer to the narrow peak at v~4640km/s only
* WKK5299 = the parameters are uncertain due to an RFI at the edge of
* the profile at v~4500km/s
* WKK5616 = see plot of WKK5659
* WKK5636 = Hanning-smoothed once despite the narrow line, due to
* ringing in the spectrum
* WKK5829 = not Hanning-smoothed due to an RFI close to the profile
* at v=4550km/s
* WKK5993 = the profile of WKK5999 (see separate pointing) has been
* excluded from the profile measurement
* WKK6680 = once Hanning-smoothed for the plot
* SKK209-5 = the RFI is next to the profile and possibly affects the line width
C=============================================================================
C Loading file 'table1.dat' ! HI-detections in the Crux and Great Attractor region
C Format for file interpretation
1 format(
+ A11,A3,1X,A13,1X,A2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,
+ A1,1X,F6.2,1X,F6.2,1X,A5,1X,I3,1X,I3,1X,F4.1,A1,1X,F4.2,1X,A1,
+ I5,A1,1X,I4,A1,1X,I3,A1,1X,F6.2,A1,1X,F4.1,1X,A1,1X,A1,1X,
+ F4.1,1X,F5.1,1X,A14)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,162
read(1,'(A154)')ar__
read(ar__,1)
+ Name,f_Name,OName,IR,RAh,RAm,RAs,DE_,DEd,DEm,DEs,u_DEs,GLON,
+ GLAT,MType,D,d_1,Bj,n_Bj,E_B_V,l_HRV,HRV,u_HRV,DV50,u_DV50,
+ DV20,u_DV20,Sint,u_Sint,rms,hann,n_Name,Dist,Sc,RFI
if(ar__(33:34) .EQ. '') RAh = iNULL__
if(ar__(36:37) .EQ. '') RAm = iNULL__
if(ar__(39:42) .EQ. '') RAs = rNULL__
if(ar__(45:46) .EQ. '') DEd = iNULL__
if(ar__(48:49) .EQ. '') DEm = iNULL__
if(ar__(51:52) .EQ. '') DEs = iNULL__
if(ar__(55:60) .EQ. '') GLON = rNULL__
if(ar__(62:67) .EQ. '') GLAT = rNULL__
if(ar__(75:77) .EQ. '') D = iNULL__
if(ar__(79:81) .EQ. '') d_1 = iNULL__
if(ar__(83:86) .EQ. '') Bj = rNULL__
if(ar__(102:105) .EQ. '') DV50 = iNULL__
if(ar__(108:110) .EQ. '') DV20 = iNULL__
if(ar__(113:118) .EQ. '') Sint = rNULL__
if(ar__(121:124) .EQ. '') rms = rNULL__
if(ar__(130:133) .EQ. '') Dist = rNULL__
if(ar__(135:139) .EQ. '') Sc = 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)
+ Name,f_Name,OName,IR,RAh,RAm,RAs,DE_,DEd,DEm,DEs,u_DEs,GLON,
+ GLAT,MType,D,d_1,Bj,n_Bj,E_B_V,l_HRV,HRV,u_HRV,DV50,u_DV50,
+ DV20,u_DV20,Sint,u_Sint,rms,hann,n_Name,Dist,Sc,RFI
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb2.dat' ! HI-detections in the Vela region
C Format for file interpretation
2 format(
+ A11,A3,1X,A13,1X,A2,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,
+ A1,1X,F6.2,1X,F6.2,1X,A5,1X,I3,1X,I3,1X,F4.1,A1,1X,F4.2,1X,A1,
+ I5,A1,1X,I4,A1,1X,I3,A1,1X,F6.2,A1,1X,F4.1,1X,A1,1X,A1,1X,
+ F4.1,1X,F5.1,1X,A14)
C Effective file loading
open(unit=1,status='old',file=
+'tableb2.dat')
write(6,*) '....Loading file: tableb2.dat'
do i__=1,14
read(1,'(A154)')ar__1
read(ar__1,2)
+ Name_1,f_Name_1,OName_1,IR_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,
+ DEm_1,DEs_1,u_DEs_1,GLON_1,GLAT_1,MType_1,D_2,d_3,Bj_1,n_Bj_1,
+ E_B_V_1,l_HRV_1,HRV_1,u_HRV_1,DV50_1,u_DV50_1,DV20_1,u_DV20_1,
+ Sint_1,u_Sint_1,rms_1,hann_1,n_Name_1,Dist_1,Sc_1,RFI_1
if(ar__1(33:34) .EQ. '') RAh_1 = iNULL__
if(ar__1(36:37) .EQ. '') RAm_1 = iNULL__
if(ar__1(39:42) .EQ. '') RAs_1 = rNULL__
if(ar__1(45:46) .EQ. '') DEd_1 = iNULL__
if(ar__1(48:49) .EQ. '') DEm_1 = iNULL__
if(ar__1(51:52) .EQ. '') DEs_1 = iNULL__
if(ar__1(55:60) .EQ. '') GLON_1 = rNULL__
if(ar__1(62:67) .EQ. '') GLAT_1 = rNULL__
if(ar__1(75:77) .EQ. '') D_2 = iNULL__
if(ar__1(79:81) .EQ. '') d_3 = iNULL__
if(ar__1(83:86) .EQ. '') Bj_1 = rNULL__
if(ar__1(102:105) .EQ. '') DV50_1 = iNULL__
if(ar__1(108:110) .EQ. '') DV20_1 = iNULL__
if(ar__1(113:118) .EQ. '') Sint_1 = rNULL__
if(ar__1(121:124) .EQ. '') rms_1 = rNULL__
if(ar__1(130:133) .EQ. '') Dist_1 = rNULL__
if(ar__1(135:139) .EQ. '') Sc_1 = rNULL__
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)
+ Name_1,f_Name_1,OName_1,IR_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,
+ DEm_1,DEs_1,u_DEs_1,GLON_1,GLAT_1,MType_1,D_2,d_3,Bj_1,n_Bj_1,
+ E_B_V_1,l_HRV_1,HRV_1,u_HRV_1,DV50_1,u_DV50_1,DV20_1,u_DV20_1,
+ Sint_1,u_Sint_1,rms_1,hann_1,n_Name_1,Dist_1,Sc_1,RFI_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end