Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/458/2221 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 2024-Mar-28 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/MNRAS/458/2221 ATLAS3D Project. XXXI (Nyland+, 2016) *================================================================================ *The ATLAS^3D^ Project. *XXXI. Nuclear radio emission in nearby early-type galaxies. * Nyland K., Young L.M., Wrobel J.M., Sarzi M., Morganti R., Alatalo K., * Blitz L., Bournaud F., Bureau M., Cappellari M., Crocker A.F., Davies R.L., * Davis T.A., De Zeeuw P.T., Duc P.-A., Emsellem E., Khochfar S., * Krajnovic D., Kuntschner H., McDermid R.M., Naab T., Oosterloo T., * Scott N., Serra P., Weijmans A.-M. * <Mon. Not. R. Astron. Soc., 458, 2221-2268 (2016)> * =2016MNRAS.458.2221N (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea1.dat' ! VLA 5GHz sample and flux density measurements integer*4 nr__ parameter (nr__=125) ! Number of records character*76 ar__ ! Full-size record character*1 f_Name ! [a] Flag on Name (1) character*9 Name ! Galaxy name character*1 n_Name ! [b] Note on Name (1) real*4 Dist ! (Mpc) Atlas3D distance (Paper I, * Cappellari et al. 2011, Cat. J/MNRAS/413/813) integer*4 Virgo ! [0/1] Virgo membership character*1 F_S ! kinematic class (Paper III, Emsellem et al., * 2011, Cat. J/MNRAS/414/888) (G1) real*4 logMJAM ! ([Msun]) ?=- Dynamically modelled stellar mass * (Paper XV, Cappellari et al., * 2013MNRAS.432.1709C) character*1 l_logMH2 ! Limit flag on logMH2 real*4 logMH2 ! ([Msun]) Molecular hydrogen mass (Paper IV, * Young et al., 2011, Cat. J/MNRAS/414/940) character*1 l_SW91 ! Limit flag on SW91 real*4 SW91 ! (mJy) ?=- Integrated flux density from the 5GHz * study of nearby ETGs carried-out by * Wrobel & Heeschen (1991AJ....101..148W) * at 5 arcsec resolution real*4 e_SW91 ! (mJy) ? rms uncertainty on SW91 integer*4 rms ! (uJy/beam) ?=- Average rms noise (3) character*1 l_Speak ! Limit flag on Speak (3) real*4 Speak ! (mJy/beam) ?=- Peak flux density (3) real*4 e_Speak ! (mJy/beam) ? rms uncertainty on Speak (3) real*4 Sint ! (mJy) ?=- Integrated flux density (3)(4) real*4 e_Sint ! (mJy) ? rms uncertainty on Sint (3) character*1 l_logL ! Limit flag on logL (3) real*4 logL ! ([W/Hz]) ?=- 5 GHz radio luminosity (3)(5) *Note (1): Notes as follows: * a = Extended source not well represented by a single two-dimensional Gaussian * model. The peak and integrated flux densities were calculated by drawing * an aperture at the 3xrms level around the source in the casa Viewer * and then using the IMSTAT task to determine the flux parameters. * b = Multicomponent source. The integrated flux density refers to the sum of * all components. See Table A3 for information on the properties of * individual components. *Note (3): refer to data from our new high-resolution 5GHz VLA observations. *Note (4): Note that measurements of the integrated flux density are only given * for sources that are formally resolved by JMFIT. *Note (5): When an integrated flux density is given, L is based on the integrated * flux density. If only a peak flux density is given (either a measurement or * an upper limit), then L is based on the peak flux density. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea2.dat' ! 5GHz spatial parameters of detections integer*4 nr__1 parameter (nr__1=53) ! Number of records character*108 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*1 f_Name_1 ! [a] Flag on Name (1) character*9 Name_1 ! galaxy name character*1 n_Name_1 ! [b] Note on Name (1) character*1 Morph ! [RU] Radio morphology based on the output of * the JMFIT task in aips (G2) integer*4 RAh ! (h) Right ascension (J2000) of the emission at * the peak flux density (G3) integer*4 RAm ! (min) Right ascension (J2000) of the emission * at the peak flux density (G3) real*4 RAs ! (s) Right ascension (J2000) of the emission * at the peak flux density (G3) character*1 DE_ ! Declination sign (J2000) of the central * position of the emission (G3) integer*4 DEd ! (deg) Declination of the central position of * the emission (G3) integer*4 DEm ! (arcmin) Declination of the central position of * the emission (G3) real*4 DEs ! (arcsec) Declination of the central position of * the emission (G3) real*4 Beam1 ! (arcsec) Angular dimension of the synthesized beam * (major axis) real*4 Beam2 ! (arcsec) Angular dimension of the synthesized beam * (minor axis) real*4 BPA ! (deg) Beam position angle, measured anti clockwise * from North character*1 l_thetaM ! Limit flag on thetaM real*4 thetaM ! (arcsec) Radius upper limit or major axis of the * emission (4) real*4 e_thetaM ! (arcsec) ? rms uncertainty on thetaM real*4 thetam_1 ! (arcsec) ? Minor axis of the emission (4) real*4 e_thetam_1 ! (arcsec) ? rms uncertainty on thetam real*4 PA ! (deg) ?=- Position angle of the emission from JMFIT real*4 e_PA ! (deg) ? rms uncertainty on PA character*1 l_Maj ! Limit flag on Maj real*8 Maj ! (pc) Linear radius upper limit or major axis of * the emission in physical units real*8 Min ! (pc) ? Linear minor axis of the emission * in physical units *Note (1): Notes as follows: * a = Extended source not well-represented by a two-dimensional Gaussian model. * The source dimensions were measured using the casa Viewer. * b = Multicomponent source. The spatial dimensions refer to the source closest * to the centre of the galaxy based on ground-based optical measurements. * See Table A4 for information on individual components. *Note (4): Angular dimensions of the emission (major x minor axis). If JMFIT * was only able to deconvolve the major axis of the source, then the minor * axis extent is given as 0.00. The errors are from JMFIT and are only given * if the emission was successfully deconvolved in at least one dimension and * categorized as resolved. For non-Gaussian sources, source dimensions were * determined using the casa Viewer and no error is reported. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea3.dat' ! 5GHz image properties of sources with multiple components integer*4 nr__2 parameter (nr__2=16) ! Number of records character*82 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*8 Name_2 ! Galaxy name character*19 Comp ! Radio component name. The brightest * component in each galaxy is listed first integer*4 RAh_1 ! (h) Right ascension (J2000) at the location of * the peak flux density integer*4 RAm_1 ! (min) Right ascension (J2000) at the location of * the peak flux density real*4 RAs_1 ! (s) Right ascension (J2000) at the location of * the peak flux density character*1 DE__1 ! Declination sign (J2000) at the location of * the peak flux density integer*4 DEd_1 ! (deg) Declination (J2000) at the location of the * peak flux density integer*4 DEm_1 ! (arcmin) Declination (J2000) at the location of the * peak flux density real*4 DEs_1 ! (arcsec) Declination (J2000) at the location of the * peak flux density real*4 Speak_1 ! (mJy/beam) Peak flux density real*4 e_Speak_1 ! (mJy/beam) rms uncertainty on Speak real*4 Sint_1 ! (mJy) ?=- Integrated flux density real*4 e_Sint_1 ! (mJy) ? rms uncertainty on Sint real*4 logL_1 ! ([W/Hz]) Radio luminosity (1) *Note (1): When an integrated flux density is given, L is based on the integrated * flux density. If only a peak flux density is given (either a measurement or * an upper limit), then L is based on the peak flux density. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea4.dat' ! 5GHz spatial properties of sources with multiple components integer*4 nr__3 parameter (nr__3=16) ! Number of records character*80 ar__3 ! Full-size record character*8 Name_3 ! Galaxy name character*19 Comp_1 ! Radio component name. The brightest component * in each galaxy is listed first character*1 Morph_1 ! [RU] Radio morphology as reported by JMFIT (G2) character*1 l_thetaM_1 ! Limit flag on thetaM real*4 thetaM_2 ! (arcsec) Radius upper limit or major axis of the * emission (G3) real*4 e_thetaM_2 ! (arcsec) ? rms uncertainty on thetaM real*4 thetam_3 ! (arcsec) ? Minor axis of the emission (G3) real*4 e_thetam_3 ! (arcsec) ? rms uncertainty on thetam real*4 PA_1 ! (deg) ?=- Position angle of the emission from JMFIT real*4 e_PA_1 ! (deg) ? rms uncertainty on PA character*1 l_Maj_1 ! Limit flag on Maj real*8 Maj_1 ! (pc) Linear radius upper limit or major axis of * the emission in physical units real*4 Min_1 ! (pc) ? Linear minor axis of the emission in * physical units c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea5.dat' ! *Archival high-resolution radio continuum data integer*4 nr__4 parameter (nr__4=29) ! Number of records character*76 ar__4 ! Full-size record character*8 Name_4 ! Galaxy name character*1 n_Name_2 ! [abc] Note on Name (1) real*4 Dist_1 ! (Mpc) Atlas3D distance (Paper I, * Cappellari et al. 2011, Cat. J/MNRAS/413/813)) integer*4 Virgo_1 ! [0/1] Virgo membership character*1 F_S_1 ! Kinematic class (Paper III, Emsellem et al., * 2011, Cat. J/MNRAS/414/888) (G1) real*4 logMJAM_1 ! ([Msun]) ?=- Dynamically modelled stellar mass * (Paper XV, Cappellari et al., * 2013MNRAS.432.1709C)) character*1 l_logMH2_1 ! Limit flag on logMH2 real*4 logMH2_1 ! ([Msun]) Molecular hydrogen mass (Paper IV, * Young et al., 2011, Cat. J/MNRAS/414/940) real*4 SW91_1 ! (mJy) ?=- Integrated flux density from the 5GHz study * of nearby ETGs carried-out by Wrobel & * Heeschen (1991AJ....101..148W) * at 5 arcsec resolution real*4 e_SW91_1 ! (mJy) ? rms uncertainty on SW91 real*4 Freq ! (GHz) Original observing frequency real*4 Res ! (arcsec) Approximate synthesized beam major axis character*1 l_S5GHz ! Limit flag on S5GHz real*8 S5GHz ! (mJy) ?=- Integrated nuclear flux density at 5GHz (2) real*4 e_S5GHz ! (mJy/beam) ? rms uncertainty on S5GHz character*1 l_logL_1 ! Limit flag on logL real*4 logL_2 ! ([W/Hz]) ?=- Nuclear radio luminosity at 5GHz character*1 Ref ! References (3) *Note (1): Notes as follows: * a = For NGC 4261, a milliarcsecond-scale spatial resolution VLBA radio flux * density is reported. * b = Although NGC 4697 and NGC 4477 were included in our new 5GHz VLA * observations, the sensitivity and quality of the final images were not * sufficient to detect them. We include a literature detection of NGC 4697 * (Wrobel et al., 2008ApJ...686..838W) as well as a detection of NGC 4477 * from project 12B-191 based on our own independent analysis. * c = For NGC 5846, the flux density has been extrapolated from to 5GHz using * a spectral index of -0.03. This spectral index was calculated based on * high-resolution measurements at 8.4 and 15GHz in Filho et al. * (2004A&A...418..429F) and Nagar et al. (2005, Cat. J/A+A/435/521), * respectively. *Note (2): For measurements originally at frequencies other than 5GHz, the flux * densities listed here have been scaled to 5GHz using the source spectral * index reported in the literature. If no such spectral index information is * available, the flux density is scaled to 5GHz assuming a flat spectral index * of {alpha}=-0.1, where S~{nu}^{alpha}. Uncertainties are reported as given in * the literature when possible. If errors were not reported in the literature, * we estimate the flux density uncertainty as described in Section 5. *Note (3): References as follows: * 1 = Filho, Barthel & Ho (2006A&A...451...71F) * 2 = Filho et al. (2004A&A...418..429F) * 3 = Kharb et al. (2012AJ....143...78K) * 4 = Nagar et al. (2001ApJ...559L..87N) * 5 = Nagar et al. (2005, Cat. J/A+A/435/521) * 6 = Project 12B-191 * 7 = Wrobel et al. (2008ApJ...686..838W) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablea6.dat' ! Additional galaxy properties integer*4 nr__5 parameter (nr__5=148) ! Number of records character*59 ar__5 ! Full-size record character*9 Name_5 ! Galaxy name real*4 logMBH ! ([Msun]) Black hole mass integer*4 r_logMBH ! ? Black hole mass reference (1) character*1 l_logLX ! Limit flag on logLX real*4 logLX ! ([10-7W]) ?=- 2-10keV X-ray luminosity from * the literature integer*4 r_logLX ! ?=- X-ray luminosity reference (2) character*1 l_log_LBol_LEdd ! Limit flag on log(LBol/LEdd) real*4 log_LBol_LEdd ! ?=- Eddington ratio (3) character*1 l_logRX ! Limit flag on logRX real*4 logRX ! ([-]) ?=- Radio-X-ray ratio as defined in * Terashima & Wilson * (2003ApJ...583..145T) real*4 logL_OIII_ ! ([10-7J]) ?=- Nuclear [OIII] luminosity real*4 EW_OIII_ ! (0.1nm) ?=- Nuclear [OIII] equivalent width character*5 Class ! Nuclear emission line classification *Note (1): for galaxies with dynamical mass estimates available in the * literature. If no reference is given, the black hole mass is estimated from * the MBH-{sigma}* relation (McConnell & Ma, 2013ApJ...764..184M). * References as follows: * 1 = Krajnovic et al. (2002A&A...390..423K) * 2 = Schulze & Gebhardt (2011ApJ...729...21S) * 3 = Bower et al. (2001ApJ...550...75B) * 4 = Barth et al. (2001ApJ...555..685B) * 5 = van den Bosch & de Zeeuw (2010MNRAS.401.1770V) * 6 = Nowak et al. (2010MNRAS.403..646N) * 7 = Gultekin et al. (2009ApJ...695.1577G) * 8 = Walsh et al. (2012ApJ...753...79W) * 9 = Ferrarese, Ford & Jaffe (1996ApJ...470..444F) * 10 = Walsh, Barth & Sarzi (2010ApJ...721..762W) * 11 = Sarzi et al. (2001ApJ...550...65S) * 12 = Rusli et al. (2013AJ....146..160R) * 13 = Gebhardt et al. (2011ApJ...729..119G) * 14 = Davis et al. (2014MNRAS.444.3427D) *Note (2): Details on the computation of this parameter are provided in * Section 6.3.2. *Note (3): References as follows: * 15 = Pellegrini (2010ApJ...717..640P) * 16 = Miller et al. (2012ApJ...747...57M) * 17 = Kharb et al. (2012AJ....143...78K) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'tablec1.dat' ! Relative contour levels in the 5GHz continuum maps integer*4 nr__6 parameter (nr__6=53) ! Number of records character*55 ar__6 ! Full-size record character*9 Name_6 ! Galaxy name integer*4 rms_1 ! (uJy/beam) rms noise character*42 RC ! Relative contours C============================================================================= C Loading file 'tablea1.dat' ! VLA 5GHz sample and flux density measurements C Format for file interpretation 1 format( + A1,A9,A1,1X,F4.1,1X,I1,1X,A1,1X,F5.2,1X,A1,F4.2,1X,A1,F4.1,1X, + F3.1,1X,I2,1X,A1,F5.2,1X,F4.2,2X,F5.2,1X,F4.2,1X,A1,F5.2) C Effective file loading open(unit=1,status='old',file= +'tablea1.dat') write(6,*) '....Loading file: tablea1.dat' do i__=1,125 read(1,'(A76)')ar__ read(ar__,1) + f_Name,Name,n_Name,Dist,Virgo,F_S,logMJAM,l_logMH2,logMH2, + l_SW91,SW91,e_SW91,rms,l_Speak,Speak,e_Speak,Sint,e_Sint, + l_logL,logL if (idig(ar__(22:26)).EQ.0) logMJAM = rNULL__ if (idig(ar__(35:38)).EQ.0) SW91 = rNULL__ if(ar__(40:42) .EQ. '') e_SW91 = rNULL__ if (rms .EQ. 45) rms = iNULL__ if (idig(ar__(48:52)).EQ.0) Speak = rNULL__ if(ar__(54:57) .EQ. '') e_Speak = rNULL__ if (idig(ar__(60:64)).EQ.0) Sint = rNULL__ if(ar__(66:69) .EQ. '') e_Sint = rNULL__ if (idig(ar__(72:76)).EQ.0) logL = rNULL__ c ..............Just test output........... write(6,1) + f_Name,Name,n_Name,Dist,Virgo,F_S,logMJAM,l_logMH2,logMH2, + l_SW91,SW91,e_SW91,rms,l_Speak,Speak,e_Speak,Sint,e_Sint, + l_logL,logL c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablea2.dat' ! 5GHz spatial parameters of detections C Format for file interpretation 2 format( + A1,A9,A1,1X,A1,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X, + F4.2,1X,F4.2,1X,F6.2,1X,A1,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X, + F6.2,1X,F5.2,1X,A1,F7.2,1X,F7.2) C Effective file loading open(unit=1,status='old',file= +'tablea2.dat') write(6,*) '....Loading file: tablea2.dat' do i__=1,53 read(1,'(A108)')ar__1 read(ar__1,2) + f_Name_1,Name_1,n_Name_1,Morph,RAh,RAm,RAs,DE_,DEd,DEm,DEs, + Beam1,Beam2,BPA,l_thetaM,thetaM,e_thetaM,thetam_1,e_thetam_1, + PA,e_PA,l_Maj,Maj,Min if(ar__1(65:68) .EQ. '') e_thetaM = rNULL__ if(ar__1(70:73) .EQ. '') thetam_1 = rNULL__ if(ar__1(75:78) .EQ. '') e_thetam_1 = rNULL__ if (idig(ar__1(80:85)).EQ.0) PA = rNULL__ if(ar__1(87:91) .EQ. '') e_PA = rNULL__ if(ar__1(102:108) .EQ. '') Min = rNULL__ 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) + f_Name_1,Name_1,n_Name_1,Morph,RAh,RAm,RAs,DE_,DEd,DEm,DEs, + Beam1,Beam2,BPA,l_thetaM,thetaM,e_thetaM,thetam_1,e_thetam_1, + PA,e_PA,l_Maj,Maj,Min write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablea3.dat' ! 5GHz image properties of sources with multiple * components C Format for file interpretation 3 format( + A8,1X,A19,1X,I2,1X,I2,1X,F6.3,1X,A1,I2,1X,I2,1X,F5.2,1X,F5.2, + 1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'tablea3.dat') write(6,*) '....Loading file: tablea3.dat' do i__=1,16 read(1,'(A82)')ar__2 read(ar__2,3) + Name_2,Comp,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Speak_1, + e_Speak_1,Sint_1,e_Sint_1,logL_1 if (idig(ar__2(67:71)).EQ.0) Sint_1 = rNULL__ if(ar__2(73:76) .EQ. '') e_Sint_1 = rNULL__ 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,Comp,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Speak_1, + e_Speak_1,Sint_1,e_Sint_1,logL_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 'tablea4.dat' ! 5GHz spatial properties of sources with * multiple components C Format for file interpretation 4 format( + A8,1X,A19,1X,A1,1X,A1,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.2,1X, + F5.2,1X,A1,F7.2,1X,F6.2) C Effective file loading open(unit=1,status='old',file= +'tablea4.dat') write(6,*) '....Loading file: tablea4.dat' do i__=1,16 read(1,'(A80)')ar__3 read(ar__3,4) + Name_3,Comp_1,Morph_1,l_thetaM_1,thetaM_2,e_thetaM_2,thetam_3, + e_thetam_3,PA_1,e_PA_1,l_Maj_1,Maj_1,Min_1 if(ar__3(38:41) .EQ. '') e_thetaM_2 = rNULL__ if(ar__3(43:46) .EQ. '') thetam_3 = rNULL__ if(ar__3(48:51) .EQ. '') e_thetam_3 = rNULL__ if (idig(ar__3(53:58)).EQ.0) PA_1 = rNULL__ if(ar__3(60:64) .EQ. '') e_PA_1 = rNULL__ if(ar__3(75:80) .EQ. '') Min_1 = rNULL__ c ..............Just test output........... write(6,4) + Name_3,Comp_1,Morph_1,l_thetaM_1,thetaM_2,e_thetaM_2,thetam_3, + e_thetam_3,PA_1,e_PA_1,l_Maj_1,Maj_1,Min_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablea5.dat' ! *Archival high-resolution radio continuum data C Format for file interpretation 5 format( + A8,A1,1X,F4.1,1X,I1,1X,A1,1X,F5.2,1X,A1,F4.2,1X,F5.2,1X,F4.2, + 1X,F4.1,1X,F5.3,1X,A1,F7.2,1X,F5.2,1X,A1,F5.2,1X,A1) C Effective file loading open(unit=1,status='old',file= +'tablea5.dat') write(6,*) '....Loading file: tablea5.dat' do i__=1,29 read(1,'(A76)')ar__4 read(ar__4,5) + Name_4,n_Name_2,Dist_1,Virgo_1,F_S_1,logMJAM_1,l_logMH2_1, + logMH2_1,SW91_1,e_SW91_1,Freq,Res,l_S5GHz,S5GHz,e_S5GHz, + l_logL_1,logL_2,Ref if (idig(ar__4(20:24)).EQ.0) logMJAM_1 = rNULL__ if (idig(ar__4(32:36)).EQ.0) SW91_1 = rNULL__ if(ar__4(38:41) .EQ. '') e_SW91_1 = rNULL__ if (idig(ar__4(55:61)).EQ.0) S5GHz = rNULL__ if(ar__4(63:67) .EQ. '') e_S5GHz = rNULL__ if (idig(ar__4(70:74)).EQ.0) logL_2 = rNULL__ c ..............Just test output........... write(6,5) + Name_4,n_Name_2,Dist_1,Virgo_1,F_S_1,logMJAM_1,l_logMH2_1, + logMH2_1,SW91_1,e_SW91_1,Freq,Res,l_S5GHz,S5GHz,e_S5GHz, + l_logL_1,logL_2,Ref c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablea6.dat' ! Additional galaxy properties C Format for file interpretation 6 format( + A9,1X,F4.2,1X,I2,1X,A1,F5.2,1X,I2,1X,A1,F5.2,1X,A1,F5.2,1X, + F5.2,1X,F5.2,1X,A5) C Effective file loading open(unit=1,status='old',file= +'tablea6.dat') write(6,*) '....Loading file: tablea6.dat' do i__=1,148 read(1,'(A59)')ar__5 read(ar__5,6) + Name_5,logMBH,r_logMBH,l_logLX,logLX,r_logLX,l_log_LBol_LEdd, + log_LBol_LEdd,l_logRX,logRX,logL_OIII_,EW_OIII_,Class if(ar__5(16:17) .EQ. '') r_logMBH = iNULL__ if (idig(ar__5(20:24)).EQ.0) logLX = rNULL__ if (r_logLX .EQ. 45) r_logLX = iNULL__ if (idig(ar__5(30:34)).EQ.0) log_LBol_LEdd = rNULL__ if (idig(ar__5(37:41)).EQ.0) logRX = rNULL__ if (idig(ar__5(43:47)).EQ.0) logL_OIII_ = rNULL__ if (idig(ar__5(49:53)).EQ.0) EW_OIII_ = rNULL__ c ..............Just test output........... write(6,6) + Name_5,logMBH,r_logMBH,l_logLX,logLX,r_logLX,l_log_LBol_LEdd, + log_LBol_LEdd,l_logRX,logRX,logL_OIII_,EW_OIII_,Class c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'tablec1.dat' ! Relative contour levels in the 5GHz continuum maps C Format for file interpretation 7 format(A9,1X,I2,1X,A42) C Effective file loading open(unit=1,status='old',file= +'tablec1.dat') write(6,*) '....Loading file: tablec1.dat' do i__=1,53 read(1,'(A55)')ar__6 read(ar__6,7)Name_6,rms_1,RC c ..............Just test output........... write(6,7)Name_6,rms_1,RC c .......End.of.Just test output........... end do close(1) C============================================================================= stop end C Locate position of first digit in string; or return 0 integer function idig(c) character*(*) c character*1 c1 integer lc,i lc=len(c) idig=0 do i=1,lc if(c(i:i).ne.' ') go to 1 end do 1 if(i.gt.lc) return c1=c(i:i) if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1 if(i.gt.lc) return c1=c(i:i) if(c1.ge.'0'.and.c1.le.'9') idig=i return end