Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/248/25 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.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 (nr__) ! Galaxy name
character*6 Type (nr__) ! Morphological type from RC3 (VII/155)
real*4 Dist (nr__) ! (Mpc) [3.05/30.6] Redshift-independent distance (1)
character*6 Nuc (nr__) ! Nuclear type (2)
character*1 f_Nuc (nr__) ! [*] Flag on Nuc (3)
character*9 D25 (nr__) ! (arcmin) Diameter from RC3 (VII/155)
integer*4 inc (nr__) ! (deg) [21/83] Inclination
integer*4 PA (nr__) ! (deg) [10/173] Position angle from RC3 (VII/155)
character*1 f_PA (nr__) ! [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 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*10 Name_1 (nr__1) ! Galaxy name
integer*4 RAh (nr__1) ! (h) Hour of right ascension (J2000)
integer*4 RAm (nr__1) ! (min) Minute of right ascension (J2000)
real*4 RAs (nr__1) ! (s) Second of right ascension (J2000)
character*1 DE_ (nr__1) ! Sign of declination (J2000)
integer*4 DEd (nr__1) ! (deg) Degree of declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs (nr__1) ! (arcsec) Arcsecond of declination (J2000)
character*9 Beam3GHz (nr__1) ! (arcsec) FWHM of the synthesized beam at 3GHz
real*4 Sig3GHz (nr__1) ! (uJy/beam) [12.8/45] 3GHZ point source sensitivity
* {sigma}
real*8 Tb3GHz (nr__1) ! (mK) [239/1431] 3GHz brightness temperature
* sensitivity {sigma}_Tb_
character*9 Beam15GHz (nr__1) ! (arcsec) FWHM of the synthesized beam at 15GHz
real*4 Sig15GHz (nr__1) ! (uJy/beam) [7.8/51] 15GHZ point source sensitivity
* {sigma}
real*4 Tb15GHz (nr__1) ! (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 (nr__2) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__2) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*14 Name_2 (nr__2) ! Extranuclear source name
integer*4 RAh_1 (nr__2) ! (h) Hour of right ascension (J2000)
integer*4 RAm_1 (nr__2) ! (min) Minute of right ascension (J2000)
real*4 RAs_1 (nr__2) ! (s) Second of right ascension (J2000)
character*1 DE__1 (nr__2) ! Sign of declination (J2000)
integer*4 DEd_1 (nr__2) ! (deg) Degree of declination (J2000)
integer*4 DEm_1 (nr__2) ! (arcmin) Arcminute of declination (J2000)
integer*4 DEs_1 (nr__2) ! (arcsec) Arcsecond of declination (J2000)
character*9 Beam3GHz_1 (nr__2) ! (arcsec) FWHM of the synthesized beam at 3GHz
real*4 Sig3GHz_1 (nr__2) ! (uJy/beam) [12.8/45] 3GHZ point source sensitivity
* {sigma}
real*8 Tb3GHz_1 (nr__2) ! (mK) [239/1431] 3GHz brightness temperature
* sensitivity {sigma}_Tb_
character*9 Beam15GHz_1(nr__2) ! (arcsec) FWHM of the synthesized beam at 15GHz
real*4 Sig15GHz_1 (nr__2) ! (uJy/beam) [7.8/51] 15GHZ point source sensitivity
* {sigma}
real*4 Tb15GHz_1 (nr__2) ! (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 (nr__3) ! (deg) Right Ascension J2000
real*8 DEdeg_2 (nr__3) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 OT (nr__3) ! Source type (G1)
character*14 Name_3 (nr__3) ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name (nr__3) ! Source identifier (blank if only one source
* in the field)
character*1 f_Name (nr__3) ! [a] Flag on Name (G2)
integer*4 RAh_2 (nr__3) ! (h) Hour of right ascension (J2000)
integer*4 RAm_2 (nr__3) ! (min) Minute of right ascension (J2000)
real*4 RAs_2 (nr__3) ! (s) Second of right ascension (J2000)
character*1 DE__2 (nr__3) ! Sign of declination (J2000)
integer*4 DEd_2 (nr__3) ! (deg) Degree of declination (J2000)
integer*4 DEm_2 (nr__3) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs_2 (nr__3) ! (arcsec) Arcsecond of declination (J2000)
character*1 l_S3GHz (nr__3) ! Limit flag on S3GHz
real*4 S3GHz (nr__3) ! (mJy) [0.03/248.5] VLA 3GHz flux density
real*4 e_S3GHz (nr__3) ! (mJy) [0.01/1.2]? S3GHz uncertainty
character*1 l_S15GHz (nr__3) ! Limit flag on S15GHz
real*4 S15GHz (nr__3) ! (mJy) [0.02/45.4] VLA 15GHz flux density
real*4 e_S15GHz (nr__3) ! (mJy) [0.01/1.4]? S15GHz uncertainty
character*1 l_S33GHz (nr__3) ! Limit flag on S33GHz
real*4 S33GHz (nr__3) ! (mJy) [0.02/61.8] VLA 33GHz flux density
real*4 e_S33GHz (nr__3) ! (mJy) [0.01/1.1]? S33GHz uncertainty
real*4 Size (nr__3) ! (pc) [44/1556] Aperture size d_ap_
real*4 rG (nr__3) ! (kpc) [0/24] Deprojected galactocentric radius
real*4 alpha (nr__3) ! [-1.4/0.9]? Observed slope {alpha}
* from 3 to 33GHz
real*4 e_alpha (nr__3) ! [0.001/0.5]? alpha uncertainty
real*4 tfrac (nr__3) ! [0.05/1]? Thermal fraction at 33GHz
real*4 e_tfrac (nr__3) ! [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 (nr__4) ! (deg) Right Ascension J2000
real*8 DEdeg_3 (nr__4) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*3 OT_1 (nr__4) ! Source type (G1)
character*14 Name_4 (nr__4) ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name_1 (nr__4) ! Source identifier (blank if only one source
* in the field)
character*1 f_Name_1 (nr__4) ! [a] Flag on Name (G2)
integer*4 RAh_3 (nr__4) ! (h) Hour of right ascension (J2000)
integer*4 RAm_3 (nr__4) ! (min) Minute of right ascension (J2000)
real*4 RAs_3 (nr__4) ! (s) Second of right ascension (J2000)
character*1 DE__3 (nr__4) ! Sign of declination (J2000)
integer*4 DEd_3 (nr__4) ! (deg) Degree of declination (J2000)
integer*4 DEm_3 (nr__4) ! (arcmin) Arcminute of declination (J2000)
real*4 DEs_3 (nr__4) ! (arcsec) Arcsecond of declination (J2000)
character*1 l_S3GHz_1 (nr__4) ! Limit flag on S3GHz
real*4 S3GHz_1 (nr__4) ! (mJy) [0.03/248.5] VLA 3GHz flux density
real*4 e_S3GHz_1 (nr__4) ! (mJy) [0.01/1.2]? S3GHz uncertainty
character*1 l_S15GHz_1 (nr__4) ! Limit flag on S15GHz
real*4 S15GHz_1 (nr__4) ! (mJy) [0.02/45.4] VLA 15GHz flux density
real*4 e_S15GHz_1 (nr__4) ! (mJy) [0.01/1.4]? S15GHz uncertainty
character*1 l_S33GHz_1 (nr__4) ! Limit flag on S33GHz
real*4 S33GHz_1 (nr__4) ! (mJy) [0.02/61.8] VLA 33GHz flux density
real*4 e_S33GHz_1 (nr__4) ! (mJy) [0.01/1.1]? S33GHz uncertainty
real*4 Size_1 (nr__4) ! (pc) [44/1556] Aperture size d_ap_
real*4 rG_1 (nr__4) ! (kpc) [0/24] Deprojected galactocentric radius
real*4 alpha_1 (nr__4) ! [-1.4/0.9]? Observed slope {alpha}
* from 3 to 33GHz
real*4 e_alpha_1 (nr__4) ! [0.001/0.5]? alpha uncertainty
real*4 tfrac_1 (nr__4) ! [0.05/1]? Thermal fraction at 33GHz
real*4 e_tfrac_1 (nr__4) ! [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 (nr__5) ! Source type, ancillary (G1)
character*14 Name_5 (nr__5) ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name_2 (nr__5) ! Source identifier (blank if only one
* source in the field)
character*1 l_F1528 (nr__5) ! Limit flag on F1528
real*8 F1528 (nr__5) ! (uJy) [0.2/1733]? GALEX FUV (1528{AA}) flux
* (f_{nu}_)
real*4 e_F1528 (nr__5) ! (uJy) [0.4/260]? F1528 uncertainty
character*1 l_F2271 (nr__5) ! Limit flag on F2271
real*8 F2271 (nr__5) ! (uJy) [0.4/2033]? GALEX NUV (2271{AA}) flux
* (f_{nu}_)
real*4 e_F2271 (nr__5) ! (uJy) [1/305]? F2271 uncertainty
character*1 l_FHa (nr__5) ! Limit flag on FHa
real*4 FHa (nr__5) ! (10-16W/m2) [0.01/8.2]? H{alpha} flux
* in 10^-13^erg/s/cm^2^ units (1)
real*4 e_FHa (nr__5) ! (10-16W/m2) [0.003/1.7]? FHa uncertainty
real*8 F24 (nr__5) ! (mJy) [0.3/2076]? Spitzer/MIPS 24um flux (f_{nu}_)
real*4 e_F24 (nr__5) ! (mJy) [0.02/104]? F24 uncertainty
character*1 f_F24 (nr__5) ! [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 (nr__6) ! Source type (G1)
character*14 Name_6 (nr__6) ! Field identifier (galaxy nuclei or
* extranuclear HII regions)
character*1 m_Name_3 (nr__6) ! Source identifier (blank if only one
* source in the field)
character*1 f_Name_2 (nr__6) ! [a] Flag on Name (G2)
real*4 bg3GHz (nr__6) ! (mJy) [-0.14/16.8]? 3GHz local background value
real*4 bg15GHz (nr__6) ! (mJy) [-0.06/1.5]? 15GHz local background value
real*4 bg33GHz (nr__6) ! (mJy) [-0.2/1.4]? 33GHz local background value
real*4 f3GHzbg (nr__6) ! [0/0.37]? 3GHz fractional contribution (1)
real*4 f15GHzbg (nr__6) ! [0/0.33]? 15GHz fractional contribution (1)
real*4 f33GHzbg (nr__6) ! [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(i__),Type(i__),Dist(i__),Nuc(i__),f_Nuc(i__),D25(i__),
+ inc(i__),PA(i__),f_PA(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),Type(i__),Dist(i__),Nuc(i__),f_Nuc(i__),D25(i__),
+ inc(i__),PA(i__),f_PA(i__)
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(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Beam3GHz(i__),Sig3GHz(i__),Tb3GHz(i__),
+ Beam15GHz(i__),Sig15GHz(i__),Tb15GHz(i__)
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),Beam3GHz(i__),Sig3GHz(i__),Tb3GHz(i__),
+ Beam15GHz(i__),Sig15GHz(i__),Tb15GHz(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),Beam3GHz_1(i__),
+ Sig3GHz_1(i__),Tb3GHz_1(i__),Beam15GHz_1(i__),Sig15GHz_1(i__),
+ Tb15GHz_1(i__)
RAdeg_1(i__) = rNULL__
DEdeg_1(i__) = 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(i__) .GT. -180) RAdeg_1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GT. -180) RAdeg_1(i__)=RAdeg_1(i__)+RAs_1(i__)/240.
if(DEd_1(i__) .GE. 0) DEdeg_1(i__)=DEd_1(i__)
if(DEm_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) DEdeg_1(i__)=DEdeg_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.DEdeg_1(i__).GE.0) DEdeg_1(i__)=-DEdeg_1(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),Beam3GHz_1(i__),
+ Sig3GHz_1(i__),Tb3GHz_1(i__),Beam15GHz_1(i__),Sig15GHz_1(i__),
+ Tb15GHz_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_1(i__),DEdeg_1(i__)
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(i__),Name_3(i__),m_Name(i__),f_Name(i__),RAh_2(i__),
+ RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),
+ DEs_2(i__),l_S3GHz(i__),S3GHz(i__),e_S3GHz(i__),l_S15GHz(i__),
+ S15GHz(i__),e_S15GHz(i__),l_S33GHz(i__),S33GHz(i__),
+ e_S33GHz(i__),Size(i__),rG(i__),alpha(i__),e_alpha(i__),
+ tfrac(i__),e_tfrac(i__)
if(ar__3(60:63) .EQ. '') e_S3GHz(i__) = rNULL__
if(ar__3(73:76) .EQ. '') e_S15GHz(i__) = rNULL__
if(ar__3(86:89) .EQ. '') e_S33GHz(i__) = rNULL__
if(ar__3(105:110) .EQ. '') alpha(i__) = rNULL__
if(ar__3(112:116) .EQ. '') e_alpha(i__) = rNULL__
if(ar__3(118:122) .EQ. '') tfrac(i__) = rNULL__
if(ar__3(124:128) .EQ. '') e_tfrac(i__) = rNULL__
RAdeg_2(i__) = rNULL__
DEdeg_2(i__) = 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(i__) .GT. -180) RAdeg_2(i__)=RAh_2(i__)*15.
if(RAm_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAm_2(i__)/4.
if(RAs_2(i__) .GT. -180) RAdeg_2(i__)=RAdeg_2(i__)+RAs_2(i__)/240.
if(DEd_2(i__) .GE. 0) DEdeg_2(i__)=DEd_2(i__)
if(DEm_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEm_2(i__)/60.
if(DEs_2(i__) .GE. 0) DEdeg_2(i__)=DEdeg_2(i__)+DEs_2(i__)/3600.
if(DE__2(i__).EQ.'-'.AND.DEdeg_2(i__).GE.0) DEdeg_2(i__)=-DEdeg_2(i__)
c ..............Just test output...........
write(6,4)
+ OT(i__),Name_3(i__),m_Name(i__),f_Name(i__),RAh_2(i__),
+ RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__),DEm_2(i__),
+ DEs_2(i__),l_S3GHz(i__),S3GHz(i__),e_S3GHz(i__),l_S15GHz(i__),
+ S15GHz(i__),e_S15GHz(i__),l_S33GHz(i__),S33GHz(i__),
+ e_S33GHz(i__),Size(i__),rG(i__),alpha(i__),e_alpha(i__),
+ tfrac(i__),e_tfrac(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__)
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(i__),Name_4(i__),m_Name_1(i__),f_Name_1(i__),RAh_3(i__),
+ RAm_3(i__),RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),
+ DEs_3(i__),l_S3GHz_1(i__),S3GHz_1(i__),e_S3GHz_1(i__),
+ l_S15GHz_1(i__),S15GHz_1(i__),e_S15GHz_1(i__),l_S33GHz_1(i__),
+ S33GHz_1(i__),e_S33GHz_1(i__),Size_1(i__),rG_1(i__),
+ alpha_1(i__),e_alpha_1(i__),tfrac_1(i__),e_tfrac_1(i__)
if(ar__4(60:63) .EQ. '') e_S3GHz_1(i__) = rNULL__
if(ar__4(73:76) .EQ. '') e_S15GHz_1(i__) = rNULL__
if(ar__4(86:89) .EQ. '') e_S33GHz_1(i__) = rNULL__
if(ar__4(105:110) .EQ. '') alpha_1(i__) = rNULL__
if(ar__4(112:116) .EQ. '') e_alpha_1(i__) = rNULL__
if(ar__4(118:122) .EQ. '') tfrac_1(i__) = rNULL__
if(ar__4(124:128) .EQ. '') e_tfrac_1(i__) = rNULL__
RAdeg_3(i__) = rNULL__
DEdeg_3(i__) = 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(i__) .GT. -180) RAdeg_3(i__)=RAh_3(i__)*15.
if(RAm_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAm_3(i__)/4.
if(RAs_3(i__) .GT. -180) RAdeg_3(i__)=RAdeg_3(i__)+RAs_3(i__)/240.
if(DEd_3(i__) .GE. 0) DEdeg_3(i__)=DEd_3(i__)
if(DEm_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEm_3(i__)/60.
if(DEs_3(i__) .GE. 0) DEdeg_3(i__)=DEdeg_3(i__)+DEs_3(i__)/3600.
if(DE__3(i__).EQ.'-'.AND.DEdeg_3(i__).GE.0) DEdeg_3(i__)=-DEdeg_3(i__)
c ..............Just test output...........
write(6,5)
+ OT_1(i__),Name_4(i__),m_Name_1(i__),f_Name_1(i__),RAh_3(i__),
+ RAm_3(i__),RAs_3(i__),DE__3(i__),DEd_3(i__),DEm_3(i__),
+ DEs_3(i__),l_S3GHz_1(i__),S3GHz_1(i__),e_S3GHz_1(i__),
+ l_S15GHz_1(i__),S15GHz_1(i__),e_S15GHz_1(i__),l_S33GHz_1(i__),
+ S33GHz_1(i__),e_S33GHz_1(i__),Size_1(i__),rG_1(i__),
+ alpha_1(i__),e_alpha_1(i__),tfrac_1(i__),e_tfrac_1(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg_3(i__),DEdeg_3(i__)
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(i__),Name_5(i__),m_Name_2(i__),l_F1528(i__),F1528(i__),
+ e_F1528(i__),l_F2271(i__),F2271(i__),e_F2271(i__),l_FHa(i__),
+ FHa(i__),e_FHa(i__),F24(i__),e_F24(i__),f_F24(i__)
if(ar__5(24:30) .EQ. '') F1528(i__) = rNULL__
if(ar__5(32:37) .EQ. '') e_F1528(i__) = rNULL__
if(ar__5(41:47) .EQ. '') F2271(i__) = rNULL__
if(ar__5(49:54) .EQ. '') e_F2271(i__) = rNULL__
if(ar__5(58:62) .EQ. '') FHa(i__) = rNULL__
if(ar__5(64:68) .EQ. '') e_FHa(i__) = rNULL__
if(ar__5(70:76) .EQ. '') F24(i__) = rNULL__
if(ar__5(78:83) .EQ. '') e_F24(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ OTa(i__),Name_5(i__),m_Name_2(i__),l_F1528(i__),F1528(i__),
+ e_F1528(i__),l_F2271(i__),F2271(i__),e_F2271(i__),l_FHa(i__),
+ FHa(i__),e_FHa(i__),F24(i__),e_F24(i__),f_F24(i__)
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(i__),Name_6(i__),m_Name_3(i__),f_Name_2(i__),bg3GHz(i__),
+ bg15GHz(i__),bg33GHz(i__),f3GHzbg(i__),f15GHzbg(i__),
+ f33GHzbg(i__)
if(ar__6(24:29) .EQ. '') bg3GHz(i__) = rNULL__
if(ar__6(31:36) .EQ. '') bg15GHz(i__) = rNULL__
if(ar__6(38:43) .EQ. '') bg33GHz(i__) = rNULL__
if(ar__6(45:49) .EQ. '') f3GHzbg(i__) = rNULL__
if(ar__6(51:55) .EQ. '') f15GHzbg(i__) = rNULL__
if(ar__6(57:61) .EQ. '') f33GHzbg(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ OT_2(i__),Name_6(i__),m_Name_3(i__),f_Name_2(i__),bg3GHz(i__),
+ bg15GHz(i__),bg33GHz(i__),f3GHzbg(i__),f15GHzbg(i__),
+ f33GHzbg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end