Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/248/25 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 2026-Oct-12
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/ApJS/248/25 The SFRS: 3-33GHz VLA obs. of star-forming regions (Linden+, 2020)
*================================================================================
*The Star Formation in Radio Survey: 3-33GHz imaging of nearby galaxy nuclei and
*extranuclear star-forming regions.
* Linden S.T., Murphy E.J., Dong D., Momjian E., Kennicutt R.C., Meier D.S.,
* Schinnerer E., Turner J.L.
* <Astrophys. J. Suppl. Ser., 248, 25 (2020)>
* =2020ApJS..248...25L
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Galaxy properties
integer*4 nr__
parameter (nr__=50) ! Number of records
character*50 ar__ ! Full-size record
character*10 Name ! Galaxy name
character*6 Type ! Morphological type from RC3 (VII/155)
real*4 Dist ! (Mpc) [3.05/30.6] Redshift-independent distance (1)
character*6 Nuc ! Nuclear type (2)
character*1 f_Nuc ! [*] Flag on Nuc (3)
character*9 D25 ! (arcmin) Diameter from RC3 (VII/155)
integer*4 inc ! (deg) [21/83] Inclination
integer*4 PA ! (deg) [10/173] Position angle from RC3 (VII/155)
character*1 f_PA ! [d] d: PA from Jarrett+ (2003AJ....125..525J)
*Note (1): Redshift-independent distance taken from the list compiled by
* Kennicutt+ (2011PASP..123.1347K), except for the two non-KINGFISH
* galaxies NGC 5194 (Ciardullo+ 2002, J/ApJ/577/31) and
* NGC 2403 (Freedman+ 2001ApJ...553...47F).
*Note (2): Nuclear type based on optical spectroscopy. Code as follows:
* SF = star-forming;
* AGN = nonthermal emission as given in Table 5 of Moustakas+ J/ApJS/190/233
*Note (3):
* * = Nuclear type based on Table 4 of
* Lonsdale Persson & Helou (1987ApJ...314..513L).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! *Nuclear source positions and imaging characteristics
integer*4 nr__1
parameter (nr__1=50) ! Number of records
character*76 ar__1 ! 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*10 Name_1 ! Galaxy name
integer*4 RAh ! (h) Hour of right ascension (J2000)
integer*4 RAm ! (min) Minute of right ascension (J2000)
real*4 RAs ! (s) Second of right ascension (J2000)
character*1 DE_ ! Sign of declination (J2000)
integer*4 DEd ! (deg) Degree of declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs ! (arcsec) Arcsecond of declination (J2000)
character*9 Beam3GHz ! (arcsec) FWHM of the synthesized beam at 3GHz
real*4 Sig3GHz ! (uJy/beam) [12.8/45] 3GHZ point source sensitivity
* {sigma}
real*8 Tb3GHz ! (mK) [239/1431] 3GHz brightness temperature
* sensitivity {sigma}_Tb_
character*9 Beam15GHz ! (arcsec) FWHM of the synthesized beam at 15GHz
real*4 Sig15GHz ! (uJy/beam) [7.8/51] 15GHZ point source sensitivity
* {sigma}
real*4 Tb15GHz ! (mK) [10/133] 15GHz brightness temperature
* sensitivity {sigma}_Tb_
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'enuc.dat' ! Extranuclear source positions and imaging
characteristics (Table 3)
integer*4 nr__2
parameter (nr__2=62) ! Number of records
character*80 ar__2 ! 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*14 Name_2 ! Extranuclear source name
integer*4 RAh_1 ! (h) Hour of right ascension (J2000)
integer*4 RAm_1 ! (min) Minute of right ascension (J2000)
real*4 RAs_1 ! (s) Second of right ascension (J2000)
character*1 DE__1 ! Sign of declination (J2000)
integer*4 DEd_1 ! (deg) Degree of declination (J2000)
integer*4 DEm_1 ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs_1 ! (arcsec) Arcsecond of declination (J2000)
character*9 Beam3GHz_1 ! (arcsec) FWHM of the synthesized beam at 3GHz
real*4 Sig3GHz_1 ! (uJy/beam) [12.8/45] 3GHZ point source sensitivity
* {sigma}
real*8 Tb3GHz_1 ! (mK) [239/1431] 3GHz brightness temperature
* sensitivity {sigma}_Tb_
character*9 Beam15GHz_1 ! (arcsec) FWHM of the synthesized beam at 15GHz
real*4 Sig15GHz_1 ! (uJy/beam) [7.8/51] 15GHZ point source sensitivity
* {sigma}
real*4 Tb15GHz_1 ! (mK) [9.8/133] 15GHz brightness temperature
* sensitivity {sigma}_Tb_
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Region photometry and derived parameters
integer*4 nr__3
parameter (nr__3=377) ! Number of records
character*128 ar__3 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_2 ! (deg) Right Ascension J2000
real*8 DEdeg_2 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 OT ! Source type (G1)
character*14 Name_3 ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name ! Source identifier (blank if only one source
* in the field)
character*1 f_Name ! [a] Flag on Name (G2)
integer*4 RAh_2 ! (h) Hour of right ascension (J2000)
integer*4 RAm_2 ! (min) Minute of right ascension (J2000)
real*4 RAs_2 ! (s) Second of right ascension (J2000)
character*1 DE__2 ! Sign of declination (J2000)
integer*4 DEd_2 ! (deg) Degree of declination (J2000)
integer*4 DEm_2 ! (arcmin) Arcminute of declination (J2000)
real*4 DEs_2 ! (arcsec) Arcsecond of declination (J2000)
character*1 l_S3GHz ! Limit flag on S3GHz
real*4 S3GHz ! (mJy) [0.03/248.5] VLA 3GHz flux density
real*4 e_S3GHz ! (mJy) [0.01/1.2]? S3GHz uncertainty
character*1 l_S15GHz ! Limit flag on S15GHz
real*4 S15GHz ! (mJy) [0.02/45.4] VLA 15GHz flux density
real*4 e_S15GHz ! (mJy) [0.01/1.4]? S15GHz uncertainty
character*1 l_S33GHz ! Limit flag on S33GHz
real*4 S33GHz ! (mJy) [0.02/61.8] VLA 33GHz flux density
real*4 e_S33GHz ! (mJy) [0.01/1.1]? S33GHz uncertainty
real*4 Size ! (pc) [44/1556] Aperture size d_ap_
real*4 rG ! (kpc) [0/24] Deprojected galactocentric radius
real*4 alpha ! [-1.4/0.9]? Observed slope {alpha}
* from 3 to 33GHz
real*4 e_alpha ! [0.001/0.5]? alpha uncertainty
real*4 tfrac ! [0.05/1]? Thermal fraction at 33GHz
real*4 e_tfrac ! [0.002/0.9]? tfrac uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Region photometry and derived parameters at 7"
angular resolution
integer*4 nr__4
parameter (nr__4=180) ! Number of records
character*128 ar__4 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_3 ! (deg) Right Ascension J2000
real*8 DEdeg_3 ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 OT_1 ! Source type (G1)
character*14 Name_4 ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name_1 ! Source identifier (blank if only one source
* in the field)
character*1 f_Name_1 ! [a] Flag on Name (G2)
integer*4 RAh_3 ! (h) Hour of right ascension (J2000)
integer*4 RAm_3 ! (min) Minute of right ascension (J2000)
real*4 RAs_3 ! (s) Second of right ascension (J2000)
character*1 DE__3 ! Sign of declination (J2000)
integer*4 DEd_3 ! (deg) Degree of declination (J2000)
integer*4 DEm_3 ! (arcmin) Arcminute of declination (J2000)
real*4 DEs_3 ! (arcsec) Arcsecond of declination (J2000)
character*1 l_S3GHz_1 ! Limit flag on S3GHz
real*4 S3GHz_1 ! (mJy) [0.03/248.5] VLA 3GHz flux density
real*4 e_S3GHz_1 ! (mJy) [0.01/1.2]? S3GHz uncertainty
character*1 l_S15GHz_1 ! Limit flag on S15GHz
real*4 S15GHz_1 ! (mJy) [0.02/45.4] VLA 15GHz flux density
real*4 e_S15GHz_1 ! (mJy) [0.01/1.4]? S15GHz uncertainty
character*1 l_S33GHz_1 ! Limit flag on S33GHz
real*4 S33GHz_1 ! (mJy) [0.02/61.8] VLA 33GHz flux density
real*4 e_S33GHz_1 ! (mJy) [0.01/1.1]? S33GHz uncertainty
real*4 Size_1 ! (pc) [44/1556] Aperture size d_ap_
real*4 rG_1 ! (kpc) [0/24] Deprojected galactocentric radius
real*4 alpha_1 ! [-1.4/0.9]? Observed slope {alpha}
* from 3 to 33GHz
real*4 e_alpha_1 ! [0.001/0.5]? alpha uncertainty
real*4 tfrac_1 ! [0.05/1]? Thermal fraction at 33GHz
real*4 e_tfrac_1 ! [0.002/0.9]? tfrac uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Ancillary source photometry at 7" angular resolution
integer*4 nr__5
parameter (nr__5=180) ! Number of records
character*84 ar__5 ! Full-size record
character*3 OTa ! Source type, ancillary (G1)
character*14 Name_5 ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name_2 ! Source identifier (blank if only one
* source in the field)
character*1 l_F1528 ! Limit flag on F1528
real*8 F1528 ! (uJy) [0.2/1733]? GALEX FUV (1528{AA}) flux
* (f_{nu}_)
real*4 e_F1528 ! (uJy) [0.4/260]? F1528 uncertainty
character*1 l_F2271 ! Limit flag on F2271
real*8 F2271 ! (uJy) [0.4/2033]? GALEX NUV (2271{AA}) flux
* (f_{nu}_)
real*4 e_F2271 ! (uJy) [1/305]? F2271 uncertainty
character*1 l_FHa ! Limit flag on FHa
real*4 FHa ! (10-16W/m2) [0.01/8.2]? H{alpha} flux
* in 10^-13^erg/s/cm^2^ units (1)
real*4 e_FHa ! (10-16W/m2) [0.003/1.7]? FHa uncertainty
real*8 F24 ! (mJy) [0.3/2076]? Spitzer/MIPS 24um flux (f_{nu}_)
real*4 e_F24 ! (mJy) [0.02/104]? F24 uncertainty
character*1 f_F24 ! [a] Flag on F24 (2)
*Note (1): The H{alpha} data used in this analysis are taken from the
* compilation presented in Leroy+ (2012AJ....144....3L) as well as the
* SINGS archive. See section 3.1.
*Note (2):
* a = Photometry unreliable and not reported owing to saturation by the
* nucleus in the 24um image.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Local background measurements
integer*4 nr__6
parameter (nr__6=377) ! Number of records
character*61 ar__6 ! Full-size record
character*3 OT_2 ! Source type (G1)
character*14 Name_6 ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name_3 ! Source identifier (blank if only one
* source in the field)
character*1 f_Name_2 ! [a] Flag on Name (G2)
real*4 bg3GHz ! (mJy) [-0.14/16.8]? 3GHz local background value
real*4 bg15GHz ! (mJy) [-0.06/1.5]? 15GHz local background value
real*4 bg33GHz ! (mJy) [-0.2/1.4]? 33GHz local background value
real*4 f3GHzbg ! [0/0.37]? 3GHz fractional contribution (1)
real*4 f15GHzbg ! [0/0.33]? 15GHz fractional contribution (1)
real*4 f33GHzbg ! [0/0.26]? 33GHz fractional contribution (1)
*Note (1): The fractional contribution of the absolute value of the local
* background measurement to the aperture flux density.
C=============================================================================
C Loading file 'table1.dat' ! Galaxy properties
C Format for file interpretation
1 format(A10,1X,A6,1X,F5.2,1X,A6,A1,1X,A9,1X,I2,1X,I3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,50
read(1,'(A50)')ar__
read(ar__,1)Name,Type,Dist,Nuc,f_Nuc,D25,inc,PA,f_PA
c ..............Just test output...........
write(6,1)Name,Type,Dist,Nuc,f_Nuc,D25,inc,PA,f_PA
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! *Nuclear source positions and imaging characteristics
C Format for file interpretation
2 format(
+ A10,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A9,1X,F4.1,1X,
+ F7.2,1X,A9,1X,F4.1,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,50
read(1,'(A76)')ar__1
read(ar__1,2)
+ Name_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Beam3GHz,Sig3GHz,Tb3GHz,
+ Beam15GHz,Sig15GHz,Tb15GHz
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,2)
+ Name_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Beam3GHz,Sig3GHz,Tb3GHz,
+ Beam15GHz,Sig15GHz,Tb15GHz
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'enuc.dat' ! Extranuclear source positions and imaging
* characteristics (Table 3)
C Format for file interpretation
3 format(
+ A14,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A9,1X,F4.1,1X,
+ F7.2,1X,A9,1X,F4.1,1X,F6.2)
C Effective file loading
open(unit=1,status='old',file=
+'enuc.dat')
write(6,*) '....Loading file: enuc.dat'
do i__=1,62
read(1,'(A80)')ar__2
read(ar__2,3)
+ Name_2,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Beam3GHz_1,
+ Sig3GHz_1,Tb3GHz_1,Beam15GHz_1,Sig15GHz_1,Tb15GHz_1
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,3)
+ Name_2,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Beam3GHz_1,
+ Sig3GHz_1,Tb3GHz_1,Beam15GHz_1,Sig15GHz_1,Tb15GHz_1
write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Region photometry and derived parameters
C Format for file interpretation
4 format(
+ A3,1X,A14,1X,A1,1X,A1,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,
+ F6.3,1X,A1,1X,F6.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,
+ F5.2,1X,F4.2,1X,F6.1,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,377
read(1,'(A128)')ar__3
read(ar__3,4)
+ OT,Name_3,m_Name,f_Name,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,
+ DEs_2,l_S3GHz,S3GHz,e_S3GHz,l_S15GHz,S15GHz,e_S15GHz,l_S33GHz,
+ S33GHz,e_S33GHz,Size,rG,alpha,e_alpha,tfrac,e_tfrac
if(ar__3(60:63) .EQ. '') e_S3GHz = rNULL__
if(ar__3(73:76) .EQ. '') e_S15GHz = rNULL__
if(ar__3(86:89) .EQ. '') e_S33GHz = rNULL__
if(ar__3(105:110) .EQ. '') alpha = rNULL__
if(ar__3(112:116) .EQ. '') e_alpha = rNULL__
if(ar__3(118:122) .EQ. '') tfrac = rNULL__
if(ar__3(124:128) .EQ. '') e_tfrac = rNULL__
RAdeg_2 = rNULL__
DEdeg_2 = rNULL__
c Derive coordinates RAdeg_2 and DEdeg_2 from input data
c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown)
if(RAh_2 .GT. -180) RAdeg_2=RAh_2*15.
if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4.
if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240.
if(DEd_2 .GE. 0) DEdeg_2=DEd_2
if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60.
if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600.
if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2
c ..............Just test output...........
write(6,4)
+ OT,Name_3,m_Name,f_Name,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,
+ DEs_2,l_S3GHz,S3GHz,e_S3GHz,l_S15GHz,S15GHz,e_S15GHz,l_S33GHz,
+ S33GHz,e_S33GHz,Size,rG,alpha,e_alpha,tfrac,e_tfrac
write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Region photometry and derived parameters at 7"
* angular resolution
C Format for file interpretation
5 format(
+ A3,1X,A14,1X,A1,1X,A1,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,
+ F6.3,1X,A1,1X,F6.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,A1,1X,
+ F5.2,1X,F4.2,1X,F6.1,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,180
read(1,'(A128)')ar__4
read(ar__4,5)
+ OT_1,Name_4,m_Name_1,f_Name_1,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,
+ DEm_3,DEs_3,l_S3GHz_1,S3GHz_1,e_S3GHz_1,l_S15GHz_1,S15GHz_1,
+ e_S15GHz_1,l_S33GHz_1,S33GHz_1,e_S33GHz_1,Size_1,rG_1,alpha_1,
+ e_alpha_1,tfrac_1,e_tfrac_1
if(ar__4(60:63) .EQ. '') e_S3GHz_1 = rNULL__
if(ar__4(73:76) .EQ. '') e_S15GHz_1 = rNULL__
if(ar__4(86:89) .EQ. '') e_S33GHz_1 = rNULL__
if(ar__4(105:110) .EQ. '') alpha_1 = rNULL__
if(ar__4(112:116) .EQ. '') e_alpha_1 = rNULL__
if(ar__4(118:122) .EQ. '') tfrac_1 = rNULL__
if(ar__4(124:128) .EQ. '') e_tfrac_1 = rNULL__
RAdeg_3 = rNULL__
DEdeg_3 = rNULL__
c Derive coordinates RAdeg_3 and DEdeg_3 from input data
c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown)
if(RAh_3 .GT. -180) RAdeg_3=RAh_3*15.
if(RAm_3 .GT. -180) RAdeg_3=RAdeg_3+RAm_3/4.
if(RAs_3 .GT. -180) RAdeg_3=RAdeg_3+RAs_3/240.
if(DEd_3 .GE. 0) DEdeg_3=DEd_3
if(DEm_3 .GE. 0) DEdeg_3=DEdeg_3+DEm_3/60.
if(DEs_3 .GE. 0) DEdeg_3=DEdeg_3+DEs_3/3600.
if(DE__3.EQ.'-'.AND.DEdeg_3.GE.0) DEdeg_3=-DEdeg_3
c ..............Just test output...........
write(6,5)
+ OT_1,Name_4,m_Name_1,f_Name_1,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,
+ DEm_3,DEs_3,l_S3GHz_1,S3GHz_1,e_S3GHz_1,l_S15GHz_1,S15GHz_1,
+ e_S15GHz_1,l_S33GHz_1,S33GHz_1,e_S33GHz_1,Size_1,rG_1,alpha_1,
+ e_alpha_1,tfrac_1,e_tfrac_1
write(6,'(6H Pos: 2F8.4)') RAdeg_3,DEdeg_3
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea1.dat' ! Ancillary source photometry at 7" angular resolution
C Format for file interpretation
6 format(
+ A3,1X,A14,1X,A1,1X,A1,1X,F7.2,1X,F6.2,1X,A1,1X,F7.2,1X,F6.2,
+ 1X,A1,1X,F5.3,1X,F5.3,1X,F7.2,1X,F6.2,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,180
read(1,'(A84)')ar__5
read(ar__5,6)
+ OTa,Name_5,m_Name_2,l_F1528,F1528,e_F1528,l_F2271,F2271,
+ e_F2271,l_FHa,FHa,e_FHa,F24,e_F24,f_F24
if(ar__5(24:30) .EQ. '') F1528 = rNULL__
if(ar__5(32:37) .EQ. '') e_F1528 = rNULL__
if(ar__5(41:47) .EQ. '') F2271 = rNULL__
if(ar__5(49:54) .EQ. '') e_F2271 = rNULL__
if(ar__5(58:62) .EQ. '') FHa = rNULL__
if(ar__5(64:68) .EQ. '') e_FHa = rNULL__
if(ar__5(70:76) .EQ. '') F24 = rNULL__
if(ar__5(78:83) .EQ. '') e_F24 = rNULL__
c ..............Just test output...........
write(6,6)
+ OTa,Name_5,m_Name_2,l_F1528,F1528,e_F1528,l_F2271,F2271,
+ e_F2271,l_FHa,FHa,e_FHa,F24,e_F24,f_F24
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Local background measurements
C Format for file interpretation
7 format(
+ A3,1X,A14,1X,A1,1X,A1,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,377
read(1,'(A61)')ar__6
read(ar__6,7)
+ OT_2,Name_6,m_Name_3,f_Name_2,bg3GHz,bg15GHz,bg33GHz,f3GHzbg,
+ f15GHzbg,f33GHzbg
if(ar__6(24:29) .EQ. '') bg3GHz = rNULL__
if(ar__6(31:36) .EQ. '') bg15GHz = rNULL__
if(ar__6(38:43) .EQ. '') bg33GHz = rNULL__
if(ar__6(45:49) .EQ. '') f3GHzbg = rNULL__
if(ar__6(51:55) .EQ. '') f15GHzbg = rNULL__
if(ar__6(57:61) .EQ. '') f33GHzbg = rNULL__
c ..............Just test output...........
write(6,7)
+ OT_2,Name_6,m_Name_3,f_Name_2,bg3GHz,bg15GHz,bg33GHz,f3GHzbg,
+ f15GHzbg,f33GHzbg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end