Conversion of standardized ReadMe file for
file /./ftp/cats/VIII/53 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.5, on 2013-Jun-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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. VIII/53 The Green Bank Survey (Maslowski+ 1972-87)
*================================================================================
*A Green Bank sky survey in search of radio sources at 1400 MHz. III.
*Positions and flux densities of the GB radio sources.
* Maslowski J.
* <Acta Astron. 22, 227 (1972)>
* =1972AcA....22..227M
*The Green Bank second survey of extragalactic radio sources at 1400 MHz.
*The catalogue of sources.
* Machalski J.
* <Acta Astron. 28, 367 (1978)>
* =1978AcA....28..367M
*The Green Bank third (GB3) survey of extragalactic radiosources at 1400 MHz.
* Rys S., Machalski J.
* <Acta Astron. 37, 163 (1987)>
* =1987AcA....37..163R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gb1.dat' ! Green Bank I Survey at 1400 MHz (Dec: 45d to 52d)
integer*4 nr__
parameter (nr__=1086) ! Number of records
character*59 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension B1950
real*8 Dec (nr__) ! (deg) Declination B1950
character*13 Name (nr__) ! Coordinate-designated B1950 name (1)
integer*4 RAh (nr__) ! (h) Right Ascension (B1950) (hour) (2)
integer*4 RAm (nr__) ! (min) Right Ascension (B1950) (min)
real*4 RAs (nr__) ! (s) Right Ascension (B1950) (sec)
character*1 DE_ (nr__) ! Declination sign (B1950) (2)
integer*4 DEd (nr__) ! (deg) Declination (B1950) (deg)
integer*4 DEm (nr__) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs (nr__) ! (arcsec) Declination (B1950) (arcsec)
character*1 pos_flg (nr__) ! [ :] Flagged ':' if position is less certain (3)
real*4 S1400 (nr__) ! (Jy) Integrated 1400 MHz flux density in Jy (4)
character*1 u_S1400 (nr__) ! [ :] ':' if flux is less certain (4)
real*4 S1400p (nr__) ! (Jy) ? Peak 1400 MHz flux density in Jy (4)
character*4 flags (nr__) ! Source flags (5)
real*4 Sp_Index (nr__) ! ? Spectral index between 1400 MHz and 408 MHz (6)
character*2 u_Sp_Index (nr__) ! [ :*] ':' if spectral index is less certain (6)
*Note (1): This column contains the name of the source, derived from its
* B1950 coordinates. Names have been designated with a 'GBx', where 'x'
* is 1 , 2, or 3, followed by the hours and minutes of RA, declination
* sign, then degrees and truncated minutes (e.g. 'GB2 0714+384'). For a
* few objects close to each other, a duplication of names was avoided by
* appending an "A" or "B" to the name.
*Note (2): These columns give the B1950 coordinates of the source, in
* degrees and in radians. Positional uncertainties for GB3 in RA and DEC
* are estimated to be:
*
* ra_err = +/-[(0.45)^2^ + (120/S)^2^]^1/2^ arcmin
*
* dec_err = +/-[(0.39)^2^ + (165/S)^2^]^1/2^ arcmin
*
* where flux density S is expressed in mJy.
*Note (3): pos_flg is flagged with a colon ':' if the position is less
* certain. In these cases the positional uncertainty is larger than that
* expressed above.
*Note (4): The integrated 1400 MHz flux is given in the S1400 column, in
* Jy. The u_S1400 column is marked with a colon ':' if the integrated flux
* is uncertain. The S1400p column gives the peak 1400 MHz flux density,
* when known.
*Note (5): This column gives flags on the source. An 'X' indicates that the
* source-fitting routine found an extension in RA. A 'Y' indicates a
* corresponding extension in declination. It is likely that in most
* cases such extension is caused by confusing faint sources. An
* evidently confused source is indicated by 'C'. Note that the X is
* always in the first column; the Y is always in the second column; and
* the C is always in the third column. Thus, a confused source has this
* column equal to ' C'. A question mark indicates uncertainty.
*Note (6): The 'Sp-Index' column gives the spectral index between the 1400 MHz
* flux density and the 408 MHz flux density from the B2 survey, GV
* survey (Grueff and Vigotti 1968ApL.....2..113G), or the WK survey
* (Windram and Kenderdine 1969MNRAS.146..265W).
*
* A colon in the 'u_Sp-Index' column indicates a less certain spectral
* index. An asterisk ('*') in the u_Sp-Index column corresponds to
* indices given in brackets in the published paper, i.e. very uncertain
* spectral indices. Note that no corrections were made to the original
* flux density data in the other catalogs.
*
* 'Sp-Index' and 'u_Sp-Index' appear in table 'gb2.dat' only.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gb2.dat' ! Green Bank II Survey at 1400 MHz (Dec: 32d to 40d)
integer*4 nr__1
parameter (nr__1=2022) ! Number of records
character*59 ar__1 ! Full-size record
real*8 RA__1 (nr__1) ! (deg) Right Ascension B1950
real*8 Dec_1 (nr__1) ! (deg) Declination B1950
character*13 Name_1 (nr__1) ! Coordinate-designated B1950 name (1)
integer*4 RAh_1 (nr__1) ! (h) Right Ascension (B1950) (hour) (2)
integer*4 RAm_1 (nr__1) ! (min) Right Ascension (B1950) (min)
real*4 RAs_1 (nr__1) ! (s) Right Ascension (B1950) (sec)
character*1 DE__1 (nr__1) ! Declination sign (B1950) (2)
integer*4 DEd_1 (nr__1) ! (deg) Declination (B1950) (deg)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs_1 (nr__1) ! (arcsec) Declination (B1950) (arcsec)
character*1 pos_flg_1 (nr__1) ! [ :] Flagged ':' if position is less certain (3)
real*4 S1400_1 (nr__1) ! (Jy) Integrated 1400 MHz flux density in Jy (4)
character*1 u_S1400_1 (nr__1) ! [ :] ':' if flux is less certain (4)
real*4 S1400p_1 (nr__1) ! (Jy) ? Peak 1400 MHz flux density in Jy (4)
character*4 flags_1 (nr__1) ! Source flags (5)
real*4 Sp_Index_1 (nr__1) ! ? Spectral index between 1400 MHz and 408 MHz (6)
character*2 u_Sp_Index_1(nr__1) ! [ :*] ':' if spectral index is less certain (6)
*Note (1): This column contains the name of the source, derived from its
* B1950 coordinates. Names have been designated with a 'GBx', where 'x'
* is 1 , 2, or 3, followed by the hours and minutes of RA, declination
* sign, then degrees and truncated minutes (e.g. 'GB2 0714+384'). For a
* few objects close to each other, a duplication of names was avoided by
* appending an "A" or "B" to the name.
*Note (2): These columns give the B1950 coordinates of the source, in
* degrees and in radians. Positional uncertainties for GB3 in RA and DEC
* are estimated to be:
*
* ra_err = +/-[(0.45)^2^ + (120/S)^2^]^1/2^ arcmin
*
* dec_err = +/-[(0.39)^2^ + (165/S)^2^]^1/2^ arcmin
*
* where flux density S is expressed in mJy.
*Note (3): pos_flg is flagged with a colon ':' if the position is less
* certain. In these cases the positional uncertainty is larger than that
* expressed above.
*Note (4): The integrated 1400 MHz flux is given in the S1400 column, in
* Jy. The u_S1400 column is marked with a colon ':' if the integrated flux
* is uncertain. The S1400p column gives the peak 1400 MHz flux density,
* when known.
*Note (5): This column gives flags on the source. An 'X' indicates that the
* source-fitting routine found an extension in RA. A 'Y' indicates a
* corresponding extension in declination. It is likely that in most
* cases such extension is caused by confusing faint sources. An
* evidently confused source is indicated by 'C'. Note that the X is
* always in the first column; the Y is always in the second column; and
* the C is always in the third column. Thus, a confused source has this
* column equal to ' C'. A question mark indicates uncertainty.
*Note (6): The 'Sp-Index' column gives the spectral index between the 1400 MHz
* flux density and the 408 MHz flux density from the B2 survey, GV
* survey (Grueff and Vigotti 1968ApL.....2..113G), or the WK survey
* (Windram and Kenderdine 1969MNRAS.146..265W).
*
* A colon in the 'u_Sp-Index' column indicates a less certain spectral
* index. An asterisk ('*') in the u_Sp-Index column corresponds to
* indices given in brackets in the published paper, i.e. very uncertain
* spectral indices. Note that no corrections were made to the original
* flux density data in the other catalogs.
*
* 'Sp-Index' and 'u_Sp-Index' appear in table 'gb2.dat' only.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'gb3.dat' ! Green Bank III Survey at 1400 MHz (Dec: 70d to 77d)
integer*4 nr__2
parameter (nr__2=676) ! Number of records
character*59 ar__2 ! Full-size record
real*8 RA__2 (nr__2) ! (deg) Right Ascension B1950
real*8 Dec_2 (nr__2) ! (deg) Declination B1950
character*13 Name_2 (nr__2) ! Coordinate-designated B1950 name (1)
integer*4 RAh_2 (nr__2) ! (h) Right Ascension (B1950) (hour) (2)
integer*4 RAm_2 (nr__2) ! (min) Right Ascension (B1950) (min)
real*4 RAs_2 (nr__2) ! (s) Right Ascension (B1950) (sec)
character*1 DE__2 (nr__2) ! Declination sign (B1950) (2)
integer*4 DEd_2 (nr__2) ! (deg) Declination (B1950) (deg)
integer*4 DEm_2 (nr__2) ! (arcmin) Declination (B1950) (arcmin)
integer*4 DEs_2 (nr__2) ! (arcsec) Declination (B1950) (arcsec)
character*1 pos_flg_2 (nr__2) ! [ :] Flagged ':' if position is less certain (3)
real*4 S1400_2 (nr__2) ! (Jy) Integrated 1400 MHz flux density in Jy (4)
character*1 u_S1400_2 (nr__2) ! [ :] ':' if flux is less certain (4)
real*4 S1400p_2 (nr__2) ! (Jy) ? Peak 1400 MHz flux density in Jy (4)
character*4 flags_2 (nr__2) ! Source flags (5)
real*4 Sp_Index_2 (nr__2) ! ? Spectral index between 1400 MHz and 408 MHz (6)
character*2 u_Sp_Index_2(nr__2) ! [ :*] ':' if spectral index is less certain (6)
*Note (1): This column contains the name of the source, derived from its
* B1950 coordinates. Names have been designated with a 'GBx', where 'x'
* is 1 , 2, or 3, followed by the hours and minutes of RA, declination
* sign, then degrees and truncated minutes (e.g. 'GB2 0714+384'). For a
* few objects close to each other, a duplication of names was avoided by
* appending an "A" or "B" to the name.
*Note (2): These columns give the B1950 coordinates of the source, in
* degrees and in radians. Positional uncertainties for GB3 in RA and DEC
* are estimated to be:
*
* ra_err = +/-[(0.45)^2^ + (120/S)^2^]^1/2^ arcmin
*
* dec_err = +/-[(0.39)^2^ + (165/S)^2^]^1/2^ arcmin
*
* where flux density S is expressed in mJy.
*Note (3): pos_flg is flagged with a colon ':' if the position is less
* certain. In these cases the positional uncertainty is larger than that
* expressed above.
*Note (4): The integrated 1400 MHz flux is given in the S1400 column, in
* Jy. The u_S1400 column is marked with a colon ':' if the integrated flux
* is uncertain. The S1400p column gives the peak 1400 MHz flux density,
* when known.
*Note (5): This column gives flags on the source. An 'X' indicates that the
* source-fitting routine found an extension in RA. A 'Y' indicates a
* corresponding extension in declination. It is likely that in most
* cases such extension is caused by confusing faint sources. An
* evidently confused source is indicated by 'C'. Note that the X is
* always in the first column; the Y is always in the second column; and
* the C is always in the third column. Thus, a confused source has this
* column equal to ' C'. A question mark indicates uncertainty.
*Note (6): The 'Sp-Index' column gives the spectral index between the 1400 MHz
* flux density and the 408 MHz flux density from the B2 survey, GV
* survey (Grueff and Vigotti 1968ApL.....2..113G), or the WK survey
* (Windram and Kenderdine 1969MNRAS.146..265W).
*
* A colon in the 'u_Sp-Index' column indicates a less certain spectral
* index. An asterisk ('*') in the u_Sp-Index column corresponds to
* indices given in brackets in the published paper, i.e. very uncertain
* spectral indices. Note that no corrections were made to the original
* flux density data in the other catalogs.
*
* 'Sp-Index' and 'u_Sp-Index' appear in table 'gb2.dat' only.
C=============================================================================
C Loading file 'gb1.dat' ! Green Bank I Survey at 1400 MHz (Dec: 45d to 52d)
C Format for file interpretation
1 format(
+ A13,2X,I2,I2,F4.1,A1,I2,I2,I2,A1,2X,F6.3,A1,1X,F4.2,2X,A4,2X,
+ F4.2,A2)
C Effective file loading
open(unit=1,file='gb1.dat', status='old')
write(6,*) '....Loading file: gb1.dat'
do i__=1,1086
read(1,'(A59)')ar__
read(ar__,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),pos_flg(i__),S1400(i__),u_S1400(i__),
+ S1400p(i__),flags(i__),Sp_Index(i__),u_Sp_Index(i__)
if(ar__(42:45) .EQ. '') S1400p(i__) = rNULL__
if(ar__(54:57) .EQ. '') Sp_Index(i__) = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),pos_flg(i__),S1400(i__),u_S1400(i__),
+ S1400p(i__),flags(i__),Sp_Index(i__),u_Sp_Index(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gb2.dat' ! Green Bank II Survey at 1400 MHz (Dec: 32d to 40d)
C Format for file interpretation
2 format(
+ A13,2X,I2,I2,F4.1,A1,I2,I2,I2,A1,2X,F6.3,A1,1X,F4.2,2X,A4,2X,
+ F4.2,A2)
C Effective file loading
open(unit=1,file='gb2.dat', status='old')
write(6,*) '....Loading file: gb2.dat'
do i__=1,2022
read(1,'(A59)')ar__1
read(ar__1,2)
+ Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),pos_flg_1(i__),S1400_1(i__),
+ u_S1400_1(i__),S1400p_1(i__),flags_1(i__),Sp_Index_1(i__),
+ u_Sp_Index_1(i__)
if(ar__1(42:45) .EQ. '') S1400p_1(i__) = rNULL__
if(ar__1(54:57) .EQ. '') Sp_Index_1(i__) = 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(i__) = RAh_1(i__)
if(RAh_1(i__) .GE. 0) RA__1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAs_1(i__)/240.
Dec_1(i__) = DEd_1(i__)
if(DEm_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),pos_flg_1(i__),S1400_1(i__),
+ u_S1400_1(i__),S1400p_1(i__),flags_1(i__),Sp_Index_1(i__),
+ u_Sp_Index_1(i__)
write(6,'(6H Pos: 2F8.4)') RA__1(i__),Dec_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'gb3.dat' ! Green Bank III Survey at 1400 MHz (Dec: 70d to 77d)
C Format for file interpretation
3 format(
+ A13,2X,I2,I2,F4.1,A1,I2,I2,I2,A1,2X,F6.3,A1,1X,F4.2,2X,A4,2X,
+ F4.2,A2)
C Effective file loading
open(unit=1,file='gb3.dat', status='old')
write(6,*) '....Loading file: gb3.dat'
do i__=1,676
read(1,'(A59)')ar__2
read(ar__2,3)
+ Name_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),pos_flg_2(i__),S1400_2(i__),
+ u_S1400_2(i__),S1400p_2(i__),flags_2(i__),Sp_Index_2(i__),
+ u_Sp_Index_2(i__)
if(ar__2(42:45) .EQ. '') S1400p_2(i__) = rNULL__
if(ar__2(54:57) .EQ. '') Sp_Index_2(i__) = rNULL__
c Derive coordinates RA__2 and Dec_2 from input data
c (RA__2 and Dec_2 are set to rNULL_ when unknown)
RA__2(i__) = RAh_2(i__)
if(RAh_2(i__) .GE. 0) RA__2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GE. 0) RA__2(i__)=RA__2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GE. 0) RA__2(i__)=RA__2(i__)+RAs_2(i__)/240.
Dec_2(i__) = DEd_2(i__)
if(DEm_2(i__) .GE. 0) Dec_2(i__)=Dec_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) Dec_2(i__)=Dec_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.Dec_2(i__).GE.0) Dec_2(i__)=-Dec_2(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),
+ DEd_2(i__),DEm_2(i__),DEs_2(i__),pos_flg_2(i__),S1400_2(i__),
+ u_S1400_2(i__),S1400p_2(i__),flags_2(i__),Sp_Index_2(i__),
+ u_Sp_Index_2(i__)
write(6,'(6H Pos: 2F8.4)') RA__2(i__),Dec_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end