Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/132/211 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.5, on 2013-May-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. J/A+AS/132/211 Interstellar SiO sources (Harju+ 1998)
*================================================================================
*A survey of SiO emission towards interstellar masers.
*I. SiO line characteristics.
* Harju J., Lehtinen K., Booth R.S., Zinchenko I.
* <Astron. Astrophys. Suppl. Ser. 132, 211 (1998)>
* =1998A&AS..132..211H (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sest2-1.dat' ! Observed SiO J=2-1 line parameters, SEST
integer*4 nr__
parameter (nr__=184) ! Number of records
character*200 ar__ ! Full-size record
real*8 RA_ ! (deg) Right Ascension 1950
real*8 Dec ! (deg) Declination 1950
integer*4 Source ! Source number
character*15 Name ! Source name
integer*4 RAh ! (h) Right ascension (1950)
integer*4 RAm ! (min) Right ascension (1950)
real*4 RAs ! (s) Right ascension (1950)
character*1 DE_ ! Declination sign
integer*4 DEd ! (deg) Declination (1950)
integer*4 DEm ! (arcmin) Declination (1950)
integer*4 DEs ! (arcsec) Declination (1950)
real*4 GLON ! (deg) Galactic longitude
real*4 GLAT ! (deg) Galactic latitude
real*4 Dist ! (kpc) ? Distance
character*1 r_Dist ! [abcefghi] Distance reference (G1)
character*5 Line ! SiO line (2-->1 or 3-->2)
character*10 FigFile ! Name of the figure with the spectrum in
* subdirectory fig
real*4 TA_ ! (K) ? Peak antenna temperature of the SiO line
real*4 e_TA_ ! (K) ? rms uncertainty on TA*
real*4 Vpeak ! (km/s) ? Peak velocity
real*4 Vmin ! (km/s) ? Minimum velocity
real*4 Vmax ! (km/s) ? Maximum velocity
real*4 Vcen ! (km/s) ? Centroid velocity (=(V_min + V_max)/2)
real*4 Vmean ! (km/s) ? Intensity weighted average velocity
real*4 e_Vmean ! (km/s) ? rms uncertainty on Vmean
real*4 FW ! (km/s) ? Full width of the detected SiO line (2)
real*4 Width ! (km/s) ? SiO line width (variance or the second moment)
real*4 e_Width ! (km/s) ? rms uncertainty on Width
real*4 Area ! (K.km/s) ? Integrated intensity of the SiO line
real*4 e_Area ! (K.km/s) ? rms uncertainty on Area
real*4 P ! ? Asymmetry parameter of the SiO line (3)
character*40 Rem ! Remarks (G2)
*Note (2): FW = V_max - V_min
*Note (3): P = 2*(V_peak - V_cen)/FW
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sest3-2.dat' ! Observed SiO J=3-2 line parameters, SEST
integer*4 nr__1
parameter (nr__1=184) ! Number of records
character*200 ar__1 ! Full-size record
real*8 RA__1 ! (deg) Right Ascension 1950
real*8 Dec_1 ! (deg) Declination 1950
integer*4 Source_1 ! Source number
character*15 Name_1 ! Source name
integer*4 RAh_1 ! (h) Right ascension (1950)
integer*4 RAm_1 ! (min) Right ascension (1950)
real*4 RAs_1 ! (s) Right ascension (1950)
character*1 DE__1 ! Declination sign
integer*4 DEd_1 ! (deg) Declination (1950)
integer*4 DEm_1 ! (arcmin) Declination (1950)
integer*4 DEs_1 ! (arcsec) Declination (1950)
real*4 GLON_1 ! (deg) Galactic longitude
real*4 GLAT_1 ! (deg) Galactic latitude
real*4 Dist_1 ! (kpc) ? Distance
character*1 r_Dist_1 ! [abcefghi] Distance reference (G1)
character*5 Line_1 ! SiO line (2-->1 or 3-->2)
character*10 FigFile_1 ! Name of the figure with the spectrum in
* subdirectory fig
real*4 TA__1 ! (K) ? Peak antenna temperature of the SiO line
real*4 e_TA__1 ! (K) ? rms uncertainty on TA*
real*4 Vpeak_1 ! (km/s) ? Peak velocity
real*4 Vmin_1 ! (km/s) ? Minimum velocity
real*4 Vmax_1 ! (km/s) ? Maximum velocity
real*4 Vcen_1 ! (km/s) ? Centroid velocity (=(V_min + V_max)/2)
real*4 Vmean_1 ! (km/s) ? Intensity weighted average velocity
real*4 e_Vmean_1 ! (km/s) ? rms uncertainty on Vmean
real*4 FW_1 ! (km/s) ? Full width of the detected SiO line (2)
real*4 Width_1 ! (km/s) ? SiO line width (variance or the second moment)
real*4 e_Width_1 ! (km/s) ? rms uncertainty on Width
real*4 Area_1 ! (K.km/s) ? Integrated intensity of the SiO line
real*4 e_Area_1 ! (K.km/s) ? rms uncertainty on Area
real*4 P_1 ! ? Asymmetry parameter of the SiO line (3)
character*40 Rem_1 ! Remarks (G2)
*Note (2): FW = V_max - V_min
*Note (3): P = 2*(V_peak - V_cen)/FW
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'onsa2-1.dat' ! Observed SiO J=2-1 line parameters, Onsala
integer*4 nr__2
parameter (nr__2=188) ! Number of records
character*200 ar__2 ! Full-size record
real*8 RA__2 ! (deg) Right Ascension 1950
real*8 Dec_2 ! (deg) Declination 1950
integer*4 Source_2 ! Source number
character*15 Name_2 ! Source name
integer*4 RAh_2 ! (h) Right ascension (1950)
integer*4 RAm_2 ! (min) Right ascension (1950)
real*4 RAs_2 ! (s) Right ascension (1950)
character*1 DE__2 ! Declination sign
integer*4 DEd_2 ! (deg) Declination (1950)
integer*4 DEm_2 ! (arcmin) Declination (1950)
integer*4 DEs_2 ! (arcsec) Declination (1950)
real*4 GLON_2 ! (deg) Galactic longitude
real*4 GLAT_2 ! (deg) Galactic latitude
real*4 Dist_2 ! (kpc) ? Distance
character*1 r_Dist_2 ! [abcefghi] Distance reference (G1)
character*5 Line_2 ! SiO line (2-->1 or 3-->2)
character*10 FigFile_2 ! Name of the figure with the spectrum in
* subdirectory fig
real*4 TA__2 ! (K) ? Peak antenna temperature of the SiO line
real*4 e_TA__2 ! (K) ? rms uncertainty on TA*
real*4 Vpeak_2 ! (km/s) ? Peak velocity
real*4 Vmin_2 ! (km/s) ? Minimum velocity
real*4 Vmax_2 ! (km/s) ? Maximum velocity
real*4 Vcen_2 ! (km/s) ? Centroid velocity (=(V_min + V_max)/2)
real*4 Vmean_2 ! (km/s) ? Intensity weighted average velocity
real*4 e_Vmean_2 ! (km/s) ? rms uncertainty on Vmean
real*4 FW_2 ! (km/s) ? Full width of the detected SiO line (2)
real*4 Width_2 ! (km/s) ? SiO line width (variance or the second moment)
real*4 e_Width_2 ! (km/s) ? rms uncertainty on Width
real*4 Area_2 ! (K.km/s) ? Integrated intensity of the SiO line
real*4 e_Area_2 ! (K.km/s) ? rms uncertainty on Area
real*4 P_2 ! ? Asymmetry parameter of the SiO line (3)
character*40 Rem_2 ! Remarks (G2)
*Note (2): FW = V_max - V_min
*Note (3): P = 2*(V_peak - V_cen)/FW
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sestcat.dat' ! SEST source catalogue
integer*4 nr__3
parameter (nr__3=184) ! Number of records
character*171 ar__3 ! Full-size record
real*8 RA__3 ! (deg) Right Ascension 1950
real*8 Dec_3 ! (deg) Declination 1950
integer*4 Source_3 ! Source number
character*15 Name_3 ! Source name
integer*4 RAh_3 ! (h) Right ascension (1950)
integer*4 RAm_3 ! (min) Right ascension (1950)
real*4 RAs_3 ! (s) Right ascension (1950)
character*1 DE__3 ! Declination sign
integer*4 DEd_3 ! (deg) Declination (1950)
integer*4 DEm_3 ! (arcmin) Declination (1950)
integer*4 DEs_3 ! (arcsec) Declination (1950)
real*4 GLON_3 ! (deg) Galactic longitude
real*4 GLAT_3 ! (deg) Galactic latitude
real*4 Dist_3 ! (kpc) ? Distance
character*1 r_Dist_3 ! [abcefghi] Distance reference (G1)
character*5 Cat ! Catalogue designation
real*8 Fpeak ! (Jy) ? Maximun flux density of the H2O maser (2)
real*4 e_Fpeak ! (Jy) ? rms uncertainty on Fpeak
real*8 Vmaser ! (km/s) ? Velocity of the strongest maser component
real*4 Vmin_3 ! (km/s) ? Minimum velocity
real*4 Vmax_3 ! (km/s) ? Maximum velocity
real*4 Vcen_3 ! (km/s) ? Centroid velocity ((=(V_min + V_max)/2)
real*4 Vmean_3 ! (km/s) ? Intensity weighted average velocity
real*4 FW_3 ! (km/s) ? Full width of the maser emission feature (3)
real*8 S ! (Jy.km/s) ? Integrated flux density of the maser line
real*4 Vcloud ! (km/s) ? Velocity of the associated molecular cloud
integer*4 r_Vcloud ! ? [1/21] Cloud velocity reference (4)
character*33 Rem_3 ! Remarks (G2)
*Note (2): H_2_O maser or sometimes OH maser, in which case clearly indicated
*Note (3): FW = V_max - V_min
*Note (4): The cloud velocity refers usually to the velocity of the
* associated dense core as estimated from CS data. When the latter has
* not been available, a radial velocity determined from NH3 or CO
* observations has been used.
*
* If no molecular line measurements has been found in the literature we
* have used for the kinematic distance estimate the peak radial velocity
* of the SiO(J=2-1) line from the present survey or the peak velocity of
* the associated maser. In this case, there is no velocity cloud data.
*
* Cloud velocity references:
* 1: Bronfman et al. (1996, Cat. <J/A+AS/115/81>)
* 2: Plume et al. (1992ApJS...78..505P)
* 3: Molinari et al. (1996, cat. <J/A+A/308/573>)
* 4: Anglada et al. (1996ApJ...463..205A)
* 5: Zinchenko et al. (1995, Cat. <J/A+AS/111/95>)
* 6: Juvela (1996, Cat. <J/A+AS/118/191>)
* 7: Wouterloot & Brand (1989, Cat. <J/A+AS/80/149>)
* 8: Churchwell et al. (1990A&AS...83..119C)
* 9: Jenness et al. (1995MNRAS.276.1024J)
* 10: Wouterloot et al. (1988A&A...203..367W)
* 11: Brand & Blitz (1993, Cat. <J/A+A/275/67>)
* 12: Fich et al. (1989ApJ...342..272F)
* 13: Lada et al. (1991ApJ...368..432L)
* 14: Haschick et al. (ApJ...265..281H)
* 15: Gardner & Whiteoak (1978MNRAS.183..711G)
* 16: Codella & Felli (1995A&A...302..521C)
* 17: Cesaroni et al. (1991A&A...252..278C)
* 18: Joncas et al. (1992ApJ...387..591J)
* 19: Haschick et al. (1983ApJ...265..281H)
* 20: Clemens & Barvainis (1988ApJS...68..257C)
* 21: Zinchenko et al. (1994, Cat. <J/A+A/288/601>)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'onsacat.dat' ! Onsala source catalogue
integer*4 nr__4
parameter (nr__4=189) ! Number of records
character*171 ar__4 ! Full-size record
real*8 RA__4 ! (deg) Right Ascension 1950
real*8 Dec_4 ! (deg) Declination 1950
integer*4 Source_4 ! Source number
character*15 Name_4 ! Source name
integer*4 RAh_4 ! (h) Right ascension (1950)
integer*4 RAm_4 ! (min) Right ascension (1950)
real*4 RAs_4 ! (s) Right ascension (1950)
character*1 DE__4 ! Declination sign
integer*4 DEd_4 ! (deg) Declination (1950)
integer*4 DEm_4 ! (arcmin) Declination (1950)
integer*4 DEs_4 ! (arcsec) Declination (1950)
real*4 GLON_4 ! (deg) Galactic longitude
real*4 GLAT_4 ! (deg) Galactic latitude
real*4 Dist_4 ! (kpc) ? Distance
character*1 r_Dist_4 ! [abcefghi] Distance reference (G1)
character*5 Cat_1 ! Catalogue designation
real*8 Fpeak_1 ! (Jy) ? Maximun flux density of the H2O maser (2)
real*4 e_Fpeak_1 ! (Jy) ? rms uncertainty on Fpeak
real*8 Vmaser_1 ! (km/s) ? Velocity of the strongest maser component
real*4 Vmin_4 ! (km/s) ? Minimum velocity
real*4 Vmax_4 ! (km/s) ? Maximum velocity
real*4 Vcen_4 ! (km/s) ? Centroid velocity ((=(V_min + V_max)/2)
real*4 Vmean_4 ! (km/s) ? Intensity weighted average velocity
real*4 FW_4 ! (km/s) ? Full width of the maser emission feature (3)
real*8 S_1 ! (Jy.km/s) ? Integrated flux density of the maser line
real*4 Vcloud_1 ! (km/s) ? Velocity of the associated molecular cloud
integer*4 r_Vcloud_1 ! ? [1/21] Cloud velocity reference (4)
character*33 Rem_4 ! Remarks (G2)
*Note (2): H_2_O maser or sometimes OH maser, in which case clearly indicated
*Note (3): FW = V_max - V_min
*Note (4): The cloud velocity refers usually to the velocity of the
* associated dense core as estimated from CS data. When the latter has
* not been available, a radial velocity determined from NH3 or CO
* observations has been used.
*
* If no molecular line measurements has been found in the literature we
* have used for the kinematic distance estimate the peak radial velocity
* of the SiO(J=2-1) line from the present survey or the peak velocity of
* the associated maser. In this case, there is no velocity cloud data.
*
* Cloud velocity references:
* 1: Bronfman et al. (1996, Cat. <J/A+AS/115/81>)
* 2: Plume et al. (1992ApJS...78..505P)
* 3: Molinari et al. (1996, cat. <J/A+A/308/573>)
* 4: Anglada et al. (1996ApJ...463..205A)
* 5: Zinchenko et al. (1995, Cat. <J/A+AS/111/95>)
* 6: Juvela (1996, Cat. <J/A+AS/118/191>)
* 7: Wouterloot & Brand (1989, Cat. <J/A+AS/80/149>)
* 8: Churchwell et al. (1990A&AS...83..119C)
* 9: Jenness et al. (1995MNRAS.276.1024J)
* 10: Wouterloot et al. (1988A&A...203..367W)
* 11: Brand & Blitz (1993, Cat. <J/A+A/275/67>)
* 12: Fich et al. (1989ApJ...342..272F)
* 13: Lada et al. (1991ApJ...368..432L)
* 14: Haschick et al. (ApJ...265..281H)
* 15: Gardner & Whiteoak (1978MNRAS.183..711G)
* 16: Codella & Felli (1995A&A...302..521C)
* 17: Cesaroni et al. (1991A&A...252..278C)
* 18: Joncas et al. (1992ApJ...387..591J)
* 19: Haschick et al. (1983ApJ...265..281H)
* 20: Clemens & Barvainis (1988ApJS...68..257C)
* 21: Zinchenko et al. (1994, Cat. <J/A+A/288/601>)
C=============================================================================
C Loading file 'sest2-1.dat' ! Observed SiO J=2-1 line parameters, SEST
C Format for file interpretation
1 format(
+ I3,1X,A15,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,
+ F6.2,1X,F5.2,A1,A5,1X,A10,1X,F4.2,1X,F4.2,1X,F5.1,1X,F6.1,1X,
+ F5.1,1X,F5.1,1X,F6.2,1X,F5.2,1X,F5.1,1X,F5.2,1X,F4.2,1X,F6.2,
+ 1X,F4.2,1X,F4.1,1X,A40)
C Effective file loading
open(unit=1,file='sest2-1.dat', status='old')
write(6,*) '....Loading file: sest2-1.dat'
do i__=1,184
read(1,'(A200)')ar__
read(ar__,1)
+ Source,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,Dist,r_Dist,
+ Line,FigFile,TA_,e_TA_,Vpeak,Vmin,Vmax,Vcen,Vmean,e_Vmean,FW,
+ Width,e_Width,Area,e_Area,P,Rem
if(ar__(56:60) .EQ. '') Dist = rNULL__
if(ar__(79:82) .EQ. '') TA_ = rNULL__
if(ar__(84:87) .EQ. '') e_TA_ = rNULL__
if(ar__(89:93) .EQ. '') Vpeak = rNULL__
if(ar__(95:100) .EQ. '') Vmin = rNULL__
if(ar__(102:106) .EQ. '') Vmax = rNULL__
if(ar__(108:112) .EQ. '') Vcen = rNULL__
if(ar__(114:119) .EQ. '') Vmean = rNULL__
if(ar__(121:125) .EQ. '') e_Vmean = rNULL__
if(ar__(127:131) .EQ. '') FW = rNULL__
if(ar__(133:137) .EQ. '') Width = rNULL__
if(ar__(139:142) .EQ. '') e_Width = rNULL__
if(ar__(144:149) .EQ. '') Area = rNULL__
if(ar__(151:154) .EQ. '') e_Area = rNULL__
if(ar__(156:159) .EQ. '') P = rNULL__
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_ = RAh
if(RAh .GE. 0) RA_=RAh*15.
if(RAm .GE. 0) RA_=RA_+RAm/4.
if(RAs .GE. 0) RA_=RA_+RAs/240.
Dec = DEd
if(DEm .GE. 0) Dec=Dec+DEm/60.
if(DEs .GE. 0) Dec=Dec+DEs/3600.
if(DE_.EQ.'-'.AND.Dec.GE.0) Dec=-Dec
c ..............Just test output...........
write(6,1)
+ Source,Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,GLON,GLAT,Dist,r_Dist,
+ Line,FigFile,TA_,e_TA_,Vpeak,Vmin,Vmax,Vcen,Vmean,e_Vmean,FW,
+ Width,e_Width,Area,e_Area,P,Rem
write(6,'(6H Pos: 2F8.4)') RA_,Dec
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'sest3-2.dat' ! Observed SiO J=3-2 line parameters, SEST
C Format for file interpretation
2 format(
+ I3,1X,A15,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,
+ F6.2,1X,F5.2,A1,A5,1X,A10,1X,F4.2,1X,F4.2,1X,F5.1,1X,F6.1,1X,
+ F5.1,1X,F5.1,1X,F6.2,1X,F5.2,1X,F5.1,1X,F5.2,1X,F4.2,1X,F6.2,
+ 1X,F4.2,1X,F4.1,1X,A40)
C Effective file loading
open(unit=1,file='sest3-2.dat', status='old')
write(6,*) '....Loading file: sest3-2.dat'
do i__=1,184
read(1,'(A200)')ar__1
read(ar__1,2)
+ Source_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ GLON_1,GLAT_1,Dist_1,r_Dist_1,Line_1,FigFile_1,TA__1,e_TA__1,
+ Vpeak_1,Vmin_1,Vmax_1,Vcen_1,Vmean_1,e_Vmean_1,FW_1,Width_1,
+ e_Width_1,Area_1,e_Area_1,P_1,Rem_1
if(ar__1(56:60) .EQ. '') Dist_1 = rNULL__
if(ar__1(79:82) .EQ. '') TA__1 = rNULL__
if(ar__1(84:87) .EQ. '') e_TA__1 = rNULL__
if(ar__1(89:93) .EQ. '') Vpeak_1 = rNULL__
if(ar__1(95:100) .EQ. '') Vmin_1 = rNULL__
if(ar__1(102:106) .EQ. '') Vmax_1 = rNULL__
if(ar__1(108:112) .EQ. '') Vcen_1 = rNULL__
if(ar__1(114:119) .EQ. '') Vmean_1 = rNULL__
if(ar__1(121:125) .EQ. '') e_Vmean_1 = rNULL__
if(ar__1(127:131) .EQ. '') FW_1 = rNULL__
if(ar__1(133:137) .EQ. '') Width_1 = rNULL__
if(ar__1(139:142) .EQ. '') e_Width_1 = rNULL__
if(ar__1(144:149) .EQ. '') Area_1 = rNULL__
if(ar__1(151:154) .EQ. '') e_Area_1 = rNULL__
if(ar__1(156:159) .EQ. '') P_1 = 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 = RAh_1
if(RAh_1 .GE. 0) RA__1=RAh_1*15.
if(RAm_1 .GE. 0) RA__1=RA__1+RAm_1/4.
if(RAs_1 .GE. 0) RA__1=RA__1+RAs_1/240.
Dec_1 = DEd_1
if(DEm_1 .GE. 0) Dec_1=Dec_1+DEm_1/60.
if(DEs_1 .GE. 0) Dec_1=Dec_1+DEs_1/3600.
if(DE__1.EQ.'-'.AND.Dec_1.GE.0) Dec_1=-Dec_1
c ..............Just test output...........
write(6,2)
+ Source_1,Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,
+ GLON_1,GLAT_1,Dist_1,r_Dist_1,Line_1,FigFile_1,TA__1,e_TA__1,
+ Vpeak_1,Vmin_1,Vmax_1,Vcen_1,Vmean_1,e_Vmean_1,FW_1,Width_1,
+ e_Width_1,Area_1,e_Area_1,P_1,Rem_1
write(6,'(6H Pos: 2F8.4)') RA__1,Dec_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'onsa2-1.dat' ! Observed SiO J=2-1 line parameters, Onsala
C Format for file interpretation
3 format(
+ I3,1X,A15,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,
+ F6.2,1X,F5.2,A1,A5,1X,A10,1X,F4.2,1X,F4.2,1X,F5.1,1X,F6.1,1X,
+ F5.1,1X,F5.1,1X,F6.2,1X,F5.2,1X,F5.1,1X,F5.2,1X,F4.2,1X,F6.2,
+ 1X,F4.2,1X,F4.1,1X,A40)
C Effective file loading
open(unit=1,file='onsa2-1.dat', status='old')
write(6,*) '....Loading file: onsa2-1.dat'
do i__=1,188
read(1,'(A200)')ar__2
read(ar__2,3)
+ Source_2,Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,
+ GLON_2,GLAT_2,Dist_2,r_Dist_2,Line_2,FigFile_2,TA__2,e_TA__2,
+ Vpeak_2,Vmin_2,Vmax_2,Vcen_2,Vmean_2,e_Vmean_2,FW_2,Width_2,
+ e_Width_2,Area_2,e_Area_2,P_2,Rem_2
if(ar__2(56:60) .EQ. '') Dist_2 = rNULL__
if(ar__2(79:82) .EQ. '') TA__2 = rNULL__
if(ar__2(84:87) .EQ. '') e_TA__2 = rNULL__
if(ar__2(89:93) .EQ. '') Vpeak_2 = rNULL__
if(ar__2(95:100) .EQ. '') Vmin_2 = rNULL__
if(ar__2(102:106) .EQ. '') Vmax_2 = rNULL__
if(ar__2(108:112) .EQ. '') Vcen_2 = rNULL__
if(ar__2(114:119) .EQ. '') Vmean_2 = rNULL__
if(ar__2(121:125) .EQ. '') e_Vmean_2 = rNULL__
if(ar__2(127:131) .EQ. '') FW_2 = rNULL__
if(ar__2(133:137) .EQ. '') Width_2 = rNULL__
if(ar__2(139:142) .EQ. '') e_Width_2 = rNULL__
if(ar__2(144:149) .EQ. '') Area_2 = rNULL__
if(ar__2(151:154) .EQ. '') e_Area_2 = rNULL__
if(ar__2(156:159) .EQ. '') P_2 = 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 = RAh_2
if(RAh_2 .GE. 0) RA__2=RAh_2*15.
if(RAm_2 .GE. 0) RA__2=RA__2+RAm_2/4.
if(RAs_2 .GE. 0) RA__2=RA__2+RAs_2/240.
Dec_2 = DEd_2
if(DEm_2 .GE. 0) Dec_2=Dec_2+DEm_2/60.
if(DEs_2 .GE. 0) Dec_2=Dec_2+DEs_2/3600.
if(DE__2.EQ.'-'.AND.Dec_2.GE.0) Dec_2=-Dec_2
c ..............Just test output...........
write(6,3)
+ Source_2,Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,
+ GLON_2,GLAT_2,Dist_2,r_Dist_2,Line_2,FigFile_2,TA__2,e_TA__2,
+ Vpeak_2,Vmin_2,Vmax_2,Vcen_2,Vmean_2,e_Vmean_2,FW_2,Width_2,
+ e_Width_2,Area_2,e_Area_2,P_2,Rem_2
write(6,'(6H Pos: 2F8.4)') RA__2,Dec_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'sestcat.dat' ! SEST source catalogue
C Format for file interpretation
4 format(
+ I3,1X,A15,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,
+ F6.2,1X,F5.2,A1,A5,1X,F9.2,1X,F4.2,1X,F7.2,1X,F5.1,1X,F5.1,1X,
+ F5.1,1X,F4.1,1X,F5.1,1X,F9.2,1X,F5.1,1X,I2,1X,A33)
C Effective file loading
open(unit=1,file='sestcat.dat', status='old')
write(6,*) '....Loading file: sestcat.dat'
do i__=1,184
read(1,'(A171)')ar__3
read(ar__3,4)
+ Source_3,Name_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3,
+ GLON_3,GLAT_3,Dist_3,r_Dist_3,Cat,Fpeak,e_Fpeak,Vmaser,Vmin_3,
+ Vmax_3,Vcen_3,Vmean_3,FW_3,S,Vcloud,r_Vcloud,Rem_3
if(ar__3(56:60) .EQ. '') Dist_3 = rNULL__
if(ar__3(68:76) .EQ. '') Fpeak = rNULL__
if(ar__3(78:81) .EQ. '') e_Fpeak = rNULL__
if(ar__3(83:89) .EQ. '') Vmaser = rNULL__
if(ar__3(91:95) .EQ. '') Vmin_3 = rNULL__
if(ar__3(97:101) .EQ. '') Vmax_3 = rNULL__
if(ar__3(103:107) .EQ. '') Vcen_3 = rNULL__
if(ar__3(109:112) .EQ. '') Vmean_3 = rNULL__
if(ar__3(114:118) .EQ. '') FW_3 = rNULL__
if(ar__3(120:128) .EQ. '') S = rNULL__
if(ar__3(130:134) .EQ. '') Vcloud = rNULL__
if(ar__3(136:137) .EQ. '') r_Vcloud = iNULL__
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 = RAh_3
if(RAh_3 .GE. 0) RA__3=RAh_3*15.
if(RAm_3 .GE. 0) RA__3=RA__3+RAm_3/4.
if(RAs_3 .GE. 0) RA__3=RA__3+RAs_3/240.
Dec_3 = DEd_3
if(DEm_3 .GE. 0) Dec_3=Dec_3+DEm_3/60.
if(DEs_3 .GE. 0) Dec_3=Dec_3+DEs_3/3600.
if(DE__3.EQ.'-'.AND.Dec_3.GE.0) Dec_3=-Dec_3
c ..............Just test output...........
write(6,4)
+ Source_3,Name_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3,
+ GLON_3,GLAT_3,Dist_3,r_Dist_3,Cat,Fpeak,e_Fpeak,Vmaser,Vmin_3,
+ Vmax_3,Vcen_3,Vmean_3,FW_3,S,Vcloud,r_Vcloud,Rem_3
write(6,'(6H Pos: 2F8.4)') RA__3,Dec_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'onsacat.dat' ! Onsala source catalogue
C Format for file interpretation
5 format(
+ I3,1X,A15,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,
+ F6.2,1X,F5.2,A1,A5,1X,F9.2,1X,F4.2,1X,F7.2,1X,F5.1,1X,F5.1,1X,
+ F5.1,1X,F4.1,1X,F5.1,1X,F9.2,1X,F5.1,1X,I2,1X,A33)
C Effective file loading
open(unit=1,file='onsacat.dat', status='old')
write(6,*) '....Loading file: onsacat.dat'
do i__=1,189
read(1,'(A171)')ar__4
read(ar__4,5)
+ Source_4,Name_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4,
+ GLON_4,GLAT_4,Dist_4,r_Dist_4,Cat_1,Fpeak_1,e_Fpeak_1,
+ Vmaser_1,Vmin_4,Vmax_4,Vcen_4,Vmean_4,FW_4,S_1,Vcloud_1,
+ r_Vcloud_1,Rem_4
if(ar__4(56:60) .EQ. '') Dist_4 = rNULL__
if(ar__4(68:76) .EQ. '') Fpeak_1 = rNULL__
if(ar__4(78:81) .EQ. '') e_Fpeak_1 = rNULL__
if(ar__4(83:89) .EQ. '') Vmaser_1 = rNULL__
if(ar__4(91:95) .EQ. '') Vmin_4 = rNULL__
if(ar__4(97:101) .EQ. '') Vmax_4 = rNULL__
if(ar__4(103:107) .EQ. '') Vcen_4 = rNULL__
if(ar__4(109:112) .EQ. '') Vmean_4 = rNULL__
if(ar__4(114:118) .EQ. '') FW_4 = rNULL__
if(ar__4(120:128) .EQ. '') S_1 = rNULL__
if(ar__4(130:134) .EQ. '') Vcloud_1 = rNULL__
if(ar__4(136:137) .EQ. '') r_Vcloud_1 = iNULL__
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 = RAh_4
if(RAh_4 .GE. 0) RA__4=RAh_4*15.
if(RAm_4 .GE. 0) RA__4=RA__4+RAm_4/4.
if(RAs_4 .GE. 0) RA__4=RA__4+RAs_4/240.
Dec_4 = DEd_4
if(DEm_4 .GE. 0) Dec_4=Dec_4+DEm_4/60.
if(DEs_4 .GE. 0) Dec_4=Dec_4+DEs_4/3600.
if(DE__4.EQ.'-'.AND.Dec_4.GE.0) Dec_4=-Dec_4
c ..............Just test output...........
write(6,5)
+ Source_4,Name_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4,
+ GLON_4,GLAT_4,Dist_4,r_Dist_4,Cat_1,Fpeak_1,e_Fpeak_1,
+ Vmaser_1,Vmin_4,Vmax_4,Vcen_4,Vmean_4,FW_4,S_1,Vcloud_1,
+ r_Vcloud_1,Rem_4
write(6,'(6H Pos: 2F8.4)') RA__4,Dec_4
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end