Conversion of standardized ReadMe file for
catalog 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-May-23
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. VIII/19 Catalogues from a deep 327 MHz Westerbork Survey (Wieringa 1991,1993)
*================================================================================
*Catalogues from a deep 327 MHz Westerbork Survey
* Wieringa M.H.: 1991, 1993
* <Bull. Inf. CDS 43, 17 (1993), and Thesis, Leiden Univ. (1991)>
* =1993BICDS..43...17W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '64w92cm.dat' ! Survey 64W2, 327 MHz (in Lynx)
integer*4 nr__
parameter (nr__=706) ! Number of records
character*110 ar__ ! Full-size record
real*8 RA_ (nr__) ! (deg) Right Ascension B1950
real*8 Dec (nr__) ! (deg) Declination B1950
character*8 NAME (nr__) ! Westerbork survey name + source number (1)
integer*4 RAh (nr__) ! (h) Right Ascension B1950 (hours)
integer*4 RAm (nr__) ! (min) Right Ascension B1950 (minutes)
real*4 RAs (nr__) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs (nr__) ! (s) Mean error on RA
character*1 DE_ (nr__) ! Declination B1950 (sign)
integer*4 DEd (nr__) ! (deg) Declination B1950 (degrees)
integer*4 DEm (nr__) ! (arcmin) Declination B1950 (minutes)
real*4 DEs (nr__) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs (nr__) ! (arcsec) Mean error on DE
real*8 S (nr__) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S (nr__) ! (mJy) Mean error on S
character*1 Res (nr__) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis (nr__) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis (nr__) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis (nr__) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA (nr__) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA (nr__) ! (deg) Mean error in PA (2)
real*4 MinAxis (nr__) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R (nr__) ! Ratio of total flux to fitted peak flux
real*8 SN (nr__) ! Peak signal to noise ratio of source in map
real*4 Atten (nr__) ! Primary beam attenuation (>1)
real*4 W (nr__) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '69w50cm.dat' ! Survey 69W, 608MHz (in Draco)
integer*4 nr__1
parameter (nr__1=794) ! Number of records
character*110 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*8 NAME_1 (nr__1) ! Westerbork survey name + source number (1)
integer*4 RAh_1 (nr__1) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_1 (nr__1) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_1 (nr__1) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_1 (nr__1) ! (s) Mean error on RA
character*1 DE__1 (nr__1) ! Declination B1950 (sign)
integer*4 DEd_1 (nr__1) ! (deg) Declination B1950 (degrees)
integer*4 DEm_1 (nr__1) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_1 (nr__1) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_1 (nr__1) ! (arcsec) Mean error on DE
real*8 S_1 (nr__1) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_1 (nr__1) ! (mJy) Mean error on S
character*1 Res_1 (nr__1) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_1(nr__1) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_1 (nr__1) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_1(nr__1) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_1 (nr__1) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_1 (nr__1) ! (deg) Mean error in PA (2)
real*4 MinAxis_1 (nr__1) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_1 (nr__1) ! Ratio of total flux to fitted peak flux
real*8 SN_1 (nr__1) ! Peak signal to noise ratio of source in map
real*4 Atten_1 (nr__1) ! Primary beam attenuation (>1)
real*4 W_1 (nr__1) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '69w92cm.dat' ! Survey 69W, 327MHz (in Draco)
integer*4 nr__2
parameter (nr__2=1199) ! Number of records
character*110 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*8 NAME_2 (nr__2) ! Westerbork survey name + source number (1)
integer*4 RAh_2 (nr__2) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_2 (nr__2) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_2 (nr__2) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_2 (nr__2) ! (s) Mean error on RA
character*1 DE__2 (nr__2) ! Declination B1950 (sign)
integer*4 DEd_2 (nr__2) ! (deg) Declination B1950 (degrees)
integer*4 DEm_2 (nr__2) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_2 (nr__2) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_2 (nr__2) ! (arcsec) Mean error on DE
real*8 S_2 (nr__2) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_2 (nr__2) ! (mJy) Mean error on S
character*1 Res_2 (nr__2) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_2(nr__2) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_2 (nr__2) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_2(nr__2) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_2 (nr__2) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_2 (nr__2) ! (deg) Mean error in PA (2)
real*4 MinAxis_2 (nr__2) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_2 (nr__2) ! Ratio of total flux to fitted peak flux
real*8 SN_2 (nr__2) ! Peak signal to noise ratio of source in map
real*4 Atten_2 (nr__2) ! Primary beam attenuation (>1)
real*4 W_2 (nr__2) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '70w92cm.dat' ! Survey 70W, 327MHz (in UMi)
integer*4 nr__3
parameter (nr__3=595) ! Number of records
character*110 ar__3 ! Full-size record
real*8 RA__3 (nr__3) ! (deg) Right Ascension B1950
real*8 Dec_3 (nr__3) ! (deg) Declination B1950
character*8 NAME_3 (nr__3) ! Westerbork survey name + source number (1)
integer*4 RAh_3 (nr__3) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_3 (nr__3) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_3 (nr__3) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_3 (nr__3) ! (s) Mean error on RA
character*1 DE__3 (nr__3) ! Declination B1950 (sign)
integer*4 DEd_3 (nr__3) ! (deg) Declination B1950 (degrees)
integer*4 DEm_3 (nr__3) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_3 (nr__3) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_3 (nr__3) ! (arcsec) Mean error on DE
real*8 S_3 (nr__3) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_3 (nr__3) ! (mJy) Mean error on S
character*1 Res_3 (nr__3) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_3(nr__3) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_3 (nr__3) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_3(nr__3) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_3 (nr__3) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_3 (nr__3) ! (deg) Mean error in PA (2)
real*4 MinAxis_3 (nr__3) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_3 (nr__3) ! Ratio of total flux to fitted peak flux
real*8 SN_3 (nr__3) ! Peak signal to noise ratio of source in map
real*4 Atten_3 (nr__3) ! Primary beam attenuation (>1)
real*4 W_3 (nr__3) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '75w50cm.dat' ! Survey 75W, 608MHz (OH 471)
integer*4 nr__4
parameter (nr__4=414) ! Number of records
character*110 ar__4 ! Full-size record
real*8 RA__4 (nr__4) ! (deg) Right Ascension B1950
real*8 Dec_4 (nr__4) ! (deg) Declination B1950
character*8 NAME_4 (nr__4) ! Westerbork survey name + source number (1)
integer*4 RAh_4 (nr__4) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_4 (nr__4) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_4 (nr__4) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_4 (nr__4) ! (s) Mean error on RA
character*1 DE__4 (nr__4) ! Declination B1950 (sign)
integer*4 DEd_4 (nr__4) ! (deg) Declination B1950 (degrees)
integer*4 DEm_4 (nr__4) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_4 (nr__4) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_4 (nr__4) ! (arcsec) Mean error on DE
real*8 S_4 (nr__4) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_4 (nr__4) ! (mJy) Mean error on S
character*1 Res_4 (nr__4) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_4(nr__4) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_4 (nr__4) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_4(nr__4) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_4 (nr__4) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_4 (nr__4) ! (deg) Mean error in PA (2)
real*4 MinAxis_4 (nr__4) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_4 (nr__4) ! Ratio of total flux to fitted peak flux
real*8 SN_4 (nr__4) ! Peak signal to noise ratio of source in map
real*4 Atten_4 (nr__4) ! Primary beam attenuation (>1)
real*4 W_4 (nr__4) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '75w92cm.dat' ! Survey 75W, 327MHz (OH 471)
integer*4 nr__5
parameter (nr__5=806) ! Number of records
character*110 ar__5 ! Full-size record
real*8 RA__5 (nr__5) ! (deg) Right Ascension B1950
real*8 Dec_5 (nr__5) ! (deg) Declination B1950
character*8 NAME_5 (nr__5) ! Westerbork survey name + source number (1)
integer*4 RAh_5 (nr__5) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_5 (nr__5) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_5 (nr__5) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_5 (nr__5) ! (s) Mean error on RA
character*1 DE__5 (nr__5) ! Declination B1950 (sign)
integer*4 DEd_5 (nr__5) ! (deg) Declination B1950 (degrees)
integer*4 DEm_5 (nr__5) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_5 (nr__5) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_5 (nr__5) ! (arcsec) Mean error on DE
real*8 S_5 (nr__5) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_5 (nr__5) ! (mJy) Mean error on S
character*1 Res_5 (nr__5) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_5(nr__5) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_5 (nr__5) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_5(nr__5) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_5 (nr__5) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_5 (nr__5) ! (deg) Mean error in PA (2)
real*4 MinAxis_5 (nr__5) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_5 (nr__5) ! Ratio of total flux to fitted peak flux
real*8 SN_5 (nr__5) ! Peak signal to noise ratio of source in map
real*4 Atten_5 (nr__5) ! Primary beam attenuation (>1)
real*4 W_5 (nr__5) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '76w50cm.dat' ! Survey 76W, 608MHz (in Cam)
integer*4 nr__6
parameter (nr__6=485) ! Number of records
character*110 ar__6 ! Full-size record
real*8 RA__6 (nr__6) ! (deg) Right Ascension B1950
real*8 Dec_6 (nr__6) ! (deg) Declination B1950
character*8 NAME_6 (nr__6) ! Westerbork survey name + source number (1)
integer*4 RAh_6 (nr__6) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_6 (nr__6) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_6 (nr__6) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_6 (nr__6) ! (s) Mean error on RA
character*1 DE__6 (nr__6) ! Declination B1950 (sign)
integer*4 DEd_6 (nr__6) ! (deg) Declination B1950 (degrees)
integer*4 DEm_6 (nr__6) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_6 (nr__6) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_6 (nr__6) ! (arcsec) Mean error on DE
real*8 S_6 (nr__6) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_6 (nr__6) ! (mJy) Mean error on S
character*1 Res_6 (nr__6) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_6(nr__6) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_6 (nr__6) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_6(nr__6) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_6 (nr__6) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_6 (nr__6) ! (deg) Mean error in PA (2)
real*4 MinAxis_6 (nr__6) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_6 (nr__6) ! Ratio of total flux to fitted peak flux
real*8 SN_6 (nr__6) ! Peak signal to noise ratio of source in map
real*4 Atten_6 (nr__6) ! Primary beam attenuation (>1)
real*4 W_6 (nr__6) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for '76w92cm.dat' ! Survey 76W, 327MHz (in Cam)
integer*4 nr__7
parameter (nr__7=851) ! Number of records
character*110 ar__7 ! Full-size record
real*8 RA__7 (nr__7) ! (deg) Right Ascension B1950
real*8 Dec_7 (nr__7) ! (deg) Declination B1950
character*8 NAME_7 (nr__7) ! Westerbork survey name + source number (1)
integer*4 RAh_7 (nr__7) ! (h) Right Ascension B1950 (hours)
integer*4 RAm_7 (nr__7) ! (min) Right Ascension B1950 (minutes)
real*4 RAs_7 (nr__7) ! (s) Right Ascension B1950 (seconds)
real*4 e_RAs_7 (nr__7) ! (s) Mean error on RA
character*1 DE__7 (nr__7) ! Declination B1950 (sign)
integer*4 DEd_7 (nr__7) ! (deg) Declination B1950 (degrees)
integer*4 DEm_7 (nr__7) ! (arcmin) Declination B1950 (minutes)
real*4 DEs_7 (nr__7) ! (arcsec) [0/60] Declination B1950 (seconds)
real*4 e_DEs_7 (nr__7) ! (arcsec) Mean error on DE
real*8 S_7 (nr__7) ! (mJy) Flux density at: 327 MHz for *92cm.dat files
* 608 MHz for *50cm.dat files
real*8 e_S_7 (nr__7) ! (mJy) Mean error on S
character*1 Res_7 (nr__7) ! [URE] indicates whether the source is
* unresolved (2)
* resolved (R, fitted gaussian parameters
* follow),
* extended (E, multiple components found and
* given on subsequent lines)
character*1 l_MajAxis_7(nr__7) ! '<' indicates that an estimated upper limit
* to the size follows (for Unresolved sources)
real*4 MajAxis_7 (nr__7) ! (arcsec) The major axis (FWHM) of the fitted gaussian
* (for resolved sources) or maximum component
* separation (for extended sources) (2)
real*4 e_MajAxis_7(nr__7) ! (arcsec) *]0/30]?=-1.0 Mean error on MajAxis
real*4 PA_7 (nr__7) ! (deg) [0/180] Position angle (North through East)(2)
real*4 e_PA_7 (nr__7) ! (deg) Mean error in PA (2)
real*4 MinAxis_7 (nr__7) ! (arcsec) The minor axis (FWHM, arcsec) of the fitted
* gaussian (2)
real*4 R_7 (nr__7) ! Ratio of total flux to fitted peak flux
real*8 SN_7 (nr__7) ! Peak signal to noise ratio of source in map
real*4 Atten_7 (nr__7) ! Primary beam attenuation (>1)
real*4 W_7 (nr__7) ! Statistical weight for use in source counts
*Note on e_MajAxis: -1.0 indicates an undefined mean error
*Note (1): The NAME has the form nnWmmmmM, where
* nn is Westerbork survey number,
* mmmm the source sequence number and
* M is blank for single sources,
* '*' for multiples,
* 'A', 'B', 'C',... for components
*Note (2): The 'Res' flag takes the values
* 'U' when the source is unresolved; in this case, the values of
* MajAxis, PA and MinAxis (as well as their mean errors)
* should be ignored (see the *.tex files)
* 'R' when the source is resolved; in this case, the fitted gaussian
* parameters follow
* 'E' for extended source; multiple components are found and given
* on subsequent lines.
C=============================================================================
C Loading file '64w92cm.dat' ! Survey 64W2, 327 MHz (in Lynx)
C Format for file interpretation
1 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='64w92cm.dat', status='old')
write(6,*) '....Loading file: 64w92cm.dat'
do i__=1,706
read(1,'(A110)')ar__
read(ar__,1)
+ NAME(i__),RAh(i__),RAm(i__),RAs(i__),e_RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),e_DEs(i__),S(i__),e_S(i__),
+ Res(i__),l_MajAxis(i__),MajAxis(i__),e_MajAxis(i__),PA(i__),
+ e_PA(i__),MinAxis(i__),R(i__),SN(i__),Atten(i__),W(i__)
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__),e_RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),e_DEs(i__),S(i__),e_S(i__),
+ Res(i__),l_MajAxis(i__),MajAxis(i__),e_MajAxis(i__),PA(i__),
+ e_PA(i__),MinAxis(i__),R(i__),SN(i__),Atten(i__),W(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 '69w50cm.dat' ! Survey 69W, 608MHz (in Draco)
C Format for file interpretation
2 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='69w50cm.dat', status='old')
write(6,*) '....Loading file: 69w50cm.dat'
do i__=1,794
read(1,'(A110)')ar__1
read(ar__1,2)
+ NAME_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),e_RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),e_DEs_1(i__),
+ S_1(i__),e_S_1(i__),Res_1(i__),l_MajAxis_1(i__),
+ MajAxis_1(i__),e_MajAxis_1(i__),PA_1(i__),e_PA_1(i__),
+ MinAxis_1(i__),R_1(i__),SN_1(i__),Atten_1(i__),W_1(i__)
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__),e_RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),e_DEs_1(i__),
+ S_1(i__),e_S_1(i__),Res_1(i__),l_MajAxis_1(i__),
+ MajAxis_1(i__),e_MajAxis_1(i__),PA_1(i__),e_PA_1(i__),
+ MinAxis_1(i__),R_1(i__),SN_1(i__),Atten_1(i__),W_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 '69w92cm.dat' ! Survey 69W, 327MHz (in Draco)
C Format for file interpretation
3 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='69w92cm.dat', status='old')
write(6,*) '....Loading file: 69w92cm.dat'
do i__=1,1199
read(1,'(A110)')ar__2
read(ar__2,3)
+ NAME_2(i__),RAh_2(i__),RAm_2(i__),RAs_2(i__),e_RAs_2(i__),
+ DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),e_DEs_2(i__),
+ S_2(i__),e_S_2(i__),Res_2(i__),l_MajAxis_2(i__),
+ MajAxis_2(i__),e_MajAxis_2(i__),PA_2(i__),e_PA_2(i__),
+ MinAxis_2(i__),R_2(i__),SN_2(i__),Atten_2(i__),W_2(i__)
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__),e_RAs_2(i__),
+ DE__2(i__),DEd_2(i__),DEm_2(i__),DEs_2(i__),e_DEs_2(i__),
+ S_2(i__),e_S_2(i__),Res_2(i__),l_MajAxis_2(i__),
+ MajAxis_2(i__),e_MajAxis_2(i__),PA_2(i__),e_PA_2(i__),
+ MinAxis_2(i__),R_2(i__),SN_2(i__),Atten_2(i__),W_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=============================================================================
C Loading file '70w92cm.dat' ! Survey 70W, 327MHz (in UMi)
C Format for file interpretation
4 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='70w92cm.dat', status='old')
write(6,*) '....Loading file: 70w92cm.dat'
do i__=1,595
read(1,'(A110)')ar__3
read(ar__3,4)
+ NAME_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),e_RAs_3(i__),
+ DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),e_DEs_3(i__),
+ S_3(i__),e_S_3(i__),Res_3(i__),l_MajAxis_3(i__),
+ MajAxis_3(i__),e_MajAxis_3(i__),PA_3(i__),e_PA_3(i__),
+ MinAxis_3(i__),R_3(i__),SN_3(i__),Atten_3(i__),W_3(i__)
c Derive coordinates RA__3 and Dec_3 from input data
c (RA__3 and Dec_3 are set to rNULL_ when unknown)
RA__3(i__) = RAh_3(i__)
if(RAh_3(i__) .GE. 0) RA__3(i__)=RAh_3(i__)*15.
if(RAm_3(i__) .GE. 0) RA__3(i__)=RA__3(i__)+RAm_3(i__)/4.
if(RAs_3(i__) .GE. 0) RA__3(i__)=RA__3(i__)+RAs_3(i__)/240.
Dec_3(i__) = DEd_3(i__)
if(DEm_3(i__) .GE. 0) Dec_3(i__)=Dec_3(i__)+DEm_3(i__)/60.
if(DEs_3(i__) .GE. 0) Dec_3(i__)=Dec_3(i__)+DEs_3(i__)/3600.
if(DE__3(i__).EQ.'-'.AND.Dec_3(i__).GE.0) Dec_3(i__)=-Dec_3(i__)
c ..............Just test output...........
write(6,4)
+ NAME_3(i__),RAh_3(i__),RAm_3(i__),RAs_3(i__),e_RAs_3(i__),
+ DE__3(i__),DEd_3(i__),DEm_3(i__),DEs_3(i__),e_DEs_3(i__),
+ S_3(i__),e_S_3(i__),Res_3(i__),l_MajAxis_3(i__),
+ MajAxis_3(i__),e_MajAxis_3(i__),PA_3(i__),e_PA_3(i__),
+ MinAxis_3(i__),R_3(i__),SN_3(i__),Atten_3(i__),W_3(i__)
write(6,'(6H Pos: 2F8.4)') RA__3(i__),Dec_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '75w50cm.dat' ! Survey 75W, 608MHz (OH 471)
C Format for file interpretation
5 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='75w50cm.dat', status='old')
write(6,*) '....Loading file: 75w50cm.dat'
do i__=1,414
read(1,'(A110)')ar__4
read(ar__4,5)
+ NAME_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),e_RAs_4(i__),
+ DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),e_DEs_4(i__),
+ S_4(i__),e_S_4(i__),Res_4(i__),l_MajAxis_4(i__),
+ MajAxis_4(i__),e_MajAxis_4(i__),PA_4(i__),e_PA_4(i__),
+ MinAxis_4(i__),R_4(i__),SN_4(i__),Atten_4(i__),W_4(i__)
c Derive coordinates RA__4 and Dec_4 from input data
c (RA__4 and Dec_4 are set to rNULL_ when unknown)
RA__4(i__) = RAh_4(i__)
if(RAh_4(i__) .GE. 0) RA__4(i__)=RAh_4(i__)*15.
if(RAm_4(i__) .GE. 0) RA__4(i__)=RA__4(i__)+RAm_4(i__)/4.
if(RAs_4(i__) .GE. 0) RA__4(i__)=RA__4(i__)+RAs_4(i__)/240.
Dec_4(i__) = DEd_4(i__)
if(DEm_4(i__) .GE. 0) Dec_4(i__)=Dec_4(i__)+DEm_4(i__)/60.
if(DEs_4(i__) .GE. 0) Dec_4(i__)=Dec_4(i__)+DEs_4(i__)/3600.
if(DE__4(i__).EQ.'-'.AND.Dec_4(i__).GE.0) Dec_4(i__)=-Dec_4(i__)
c ..............Just test output...........
write(6,5)
+ NAME_4(i__),RAh_4(i__),RAm_4(i__),RAs_4(i__),e_RAs_4(i__),
+ DE__4(i__),DEd_4(i__),DEm_4(i__),DEs_4(i__),e_DEs_4(i__),
+ S_4(i__),e_S_4(i__),Res_4(i__),l_MajAxis_4(i__),
+ MajAxis_4(i__),e_MajAxis_4(i__),PA_4(i__),e_PA_4(i__),
+ MinAxis_4(i__),R_4(i__),SN_4(i__),Atten_4(i__),W_4(i__)
write(6,'(6H Pos: 2F8.4)') RA__4(i__),Dec_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '75w92cm.dat' ! Survey 75W, 327MHz (OH 471)
C Format for file interpretation
6 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='75w92cm.dat', status='old')
write(6,*) '....Loading file: 75w92cm.dat'
do i__=1,806
read(1,'(A110)')ar__5
read(ar__5,6)
+ NAME_5(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),e_RAs_5(i__),
+ DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),e_DEs_5(i__),
+ S_5(i__),e_S_5(i__),Res_5(i__),l_MajAxis_5(i__),
+ MajAxis_5(i__),e_MajAxis_5(i__),PA_5(i__),e_PA_5(i__),
+ MinAxis_5(i__),R_5(i__),SN_5(i__),Atten_5(i__),W_5(i__)
c Derive coordinates RA__5 and Dec_5 from input data
c (RA__5 and Dec_5 are set to rNULL_ when unknown)
RA__5(i__) = RAh_5(i__)
if(RAh_5(i__) .GE. 0) RA__5(i__)=RAh_5(i__)*15.
if(RAm_5(i__) .GE. 0) RA__5(i__)=RA__5(i__)+RAm_5(i__)/4.
if(RAs_5(i__) .GE. 0) RA__5(i__)=RA__5(i__)+RAs_5(i__)/240.
Dec_5(i__) = DEd_5(i__)
if(DEm_5(i__) .GE. 0) Dec_5(i__)=Dec_5(i__)+DEm_5(i__)/60.
if(DEs_5(i__) .GE. 0) Dec_5(i__)=Dec_5(i__)+DEs_5(i__)/3600.
if(DE__5(i__).EQ.'-'.AND.Dec_5(i__).GE.0) Dec_5(i__)=-Dec_5(i__)
c ..............Just test output...........
write(6,6)
+ NAME_5(i__),RAh_5(i__),RAm_5(i__),RAs_5(i__),e_RAs_5(i__),
+ DE__5(i__),DEd_5(i__),DEm_5(i__),DEs_5(i__),e_DEs_5(i__),
+ S_5(i__),e_S_5(i__),Res_5(i__),l_MajAxis_5(i__),
+ MajAxis_5(i__),e_MajAxis_5(i__),PA_5(i__),e_PA_5(i__),
+ MinAxis_5(i__),R_5(i__),SN_5(i__),Atten_5(i__),W_5(i__)
write(6,'(6H Pos: 2F8.4)') RA__5(i__),Dec_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '76w50cm.dat' ! Survey 76W, 608MHz (in Cam)
C Format for file interpretation
7 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='76w50cm.dat', status='old')
write(6,*) '....Loading file: 76w50cm.dat'
do i__=1,485
read(1,'(A110)')ar__6
read(ar__6,7)
+ NAME_6(i__),RAh_6(i__),RAm_6(i__),RAs_6(i__),e_RAs_6(i__),
+ DE__6(i__),DEd_6(i__),DEm_6(i__),DEs_6(i__),e_DEs_6(i__),
+ S_6(i__),e_S_6(i__),Res_6(i__),l_MajAxis_6(i__),
+ MajAxis_6(i__),e_MajAxis_6(i__),PA_6(i__),e_PA_6(i__),
+ MinAxis_6(i__),R_6(i__),SN_6(i__),Atten_6(i__),W_6(i__)
c Derive coordinates RA__6 and Dec_6 from input data
c (RA__6 and Dec_6 are set to rNULL_ when unknown)
RA__6(i__) = RAh_6(i__)
if(RAh_6(i__) .GE. 0) RA__6(i__)=RAh_6(i__)*15.
if(RAm_6(i__) .GE. 0) RA__6(i__)=RA__6(i__)+RAm_6(i__)/4.
if(RAs_6(i__) .GE. 0) RA__6(i__)=RA__6(i__)+RAs_6(i__)/240.
Dec_6(i__) = DEd_6(i__)
if(DEm_6(i__) .GE. 0) Dec_6(i__)=Dec_6(i__)+DEm_6(i__)/60.
if(DEs_6(i__) .GE. 0) Dec_6(i__)=Dec_6(i__)+DEs_6(i__)/3600.
if(DE__6(i__).EQ.'-'.AND.Dec_6(i__).GE.0) Dec_6(i__)=-Dec_6(i__)
c ..............Just test output...........
write(6,7)
+ NAME_6(i__),RAh_6(i__),RAm_6(i__),RAs_6(i__),e_RAs_6(i__),
+ DE__6(i__),DEd_6(i__),DEm_6(i__),DEs_6(i__),e_DEs_6(i__),
+ S_6(i__),e_S_6(i__),Res_6(i__),l_MajAxis_6(i__),
+ MajAxis_6(i__),e_MajAxis_6(i__),PA_6(i__),e_PA_6(i__),
+ MinAxis_6(i__),R_6(i__),SN_6(i__),Atten_6(i__),W_6(i__)
write(6,'(6H Pos: 2F8.4)') RA__6(i__),Dec_6(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file '76w92cm.dat' ! Survey 76W, 327MHz (in Cam)
C Format for file interpretation
8 format(
+ A8,1X,I2,I2,F5.2,1X,F4.2,1X,A1,I2,I2,F4.1,1X,F4.1,1X,F8.2,1X,
+ F7.2,1X,A1,A1,1X,F5.1,1X,F4.1,1X,F4.0,1X,F3.0,1X,F5.1,1X,F5.2,
+ F8.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,file='76w92cm.dat', status='old')
write(6,*) '....Loading file: 76w92cm.dat'
do i__=1,851
read(1,'(A110)')ar__7
read(ar__7,8)
+ NAME_7(i__),RAh_7(i__),RAm_7(i__),RAs_7(i__),e_RAs_7(i__),
+ DE__7(i__),DEd_7(i__),DEm_7(i__),DEs_7(i__),e_DEs_7(i__),
+ S_7(i__),e_S_7(i__),Res_7(i__),l_MajAxis_7(i__),
+ MajAxis_7(i__),e_MajAxis_7(i__),PA_7(i__),e_PA_7(i__),
+ MinAxis_7(i__),R_7(i__),SN_7(i__),Atten_7(i__),W_7(i__)
c Derive coordinates RA__7 and Dec_7 from input data
c (RA__7 and Dec_7 are set to rNULL_ when unknown)
RA__7(i__) = RAh_7(i__)
if(RAh_7(i__) .GE. 0) RA__7(i__)=RAh_7(i__)*15.
if(RAm_7(i__) .GE. 0) RA__7(i__)=RA__7(i__)+RAm_7(i__)/4.
if(RAs_7(i__) .GE. 0) RA__7(i__)=RA__7(i__)+RAs_7(i__)/240.
Dec_7(i__) = DEd_7(i__)
if(DEm_7(i__) .GE. 0) Dec_7(i__)=Dec_7(i__)+DEm_7(i__)/60.
if(DEs_7(i__) .GE. 0) Dec_7(i__)=Dec_7(i__)+DEs_7(i__)/3600.
if(DE__7(i__).EQ.'-'.AND.Dec_7(i__).GE.0) Dec_7(i__)=-Dec_7(i__)
c ..............Just test output...........
write(6,8)
+ NAME_7(i__),RAh_7(i__),RAm_7(i__),RAs_7(i__),e_RAs_7(i__),
+ DE__7(i__),DEd_7(i__),DEm_7(i__),DEs_7(i__),e_DEs_7(i__),
+ S_7(i__),e_S_7(i__),Res_7(i__),l_MajAxis_7(i__),
+ MajAxis_7(i__),e_MajAxis_7(i__),PA_7(i__),e_PA_7(i__),
+ MinAxis_7(i__),R_7(i__),SN_7(i__),Atten_7(i__),W_7(i__)
write(6,'(6H Pos: 2F8.4)') RA__7(i__),Dec_7(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end