Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/144/125 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/AJ/144/125 Photometry and spectroscopy of Markarian 266 (Mazzarella+, 2012)
*================================================================================
*Investigation of dual active nuclei, outflows, shock-heated gas, and young star
*clusters in Markarian 266.
* Mazzarella J.M., Iwasawa K., Vavilkin T., Armus L., Kim D.-C., Bothun G.,
* Evans A.S., Spoon H.W.W., Haan S., Howell J.H., Lord S., Marshall J.A.,
* Ishida C.M., Xu C.K., Petric A., Sanders D.B., Surace J.A., Appleton P.,
* Chan B.H.P., Frayer D.T., Inami H., Khachikian E.Y., Madore B.F.,
* Privon G.C., Sturm E., U V., Veilleux S.
* <Astron. J., 144, 125 (2012)>
* =2012AJ....144..125M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Measurements from low-resolution spectral mapping
integer*4 nr__
parameter (nr__=38) ! Number of records
character*57 ar__ ! Full-size record
character*5 Region (nr__) ! Mrk 266 region (G1)
character*1 f_Region (nr__) ! [ad] Flag on Region (1)
real*4 lambda (nr__) ! (um) [6/35] Spectral line wavelength ({mu}m)
character*16 Line (nr__) ! Spectral line identification
real*8 Flux (nr__) ! (10-16W/m2) Line flux (in 10^-20^W/cm^2^)
real*4 e_Flux (nr__) ! (10-16W/m2) Uncertainty in Flux
character*1 f_Flux (nr__) ! [bc] Flag on Flux (2)
real*4 EW (nr__) ! (um) ? Equivalent width; in {mu}m
real*4 e_EW (nr__) ! (um) ? Uncertainty in EW; in {mu}m
character*1 f_EW (nr__) ! [bc] Flag on EW (2)
*Note (1): Flag on Mrk 266 component as follows:
* a = Spectral measurements for the integrated emission from both galaxies in
* Mrk 266. Features in the 5.2-14.5{mu}m (Short-Low, SL) and 14.5-38{mu}m
* (Long-Low, LL) regions were measured from the 1D spectrum plotted in
* Figure 11; this region is 4*4 LL pixels in size, where each pixel is
* 5.1"*5.1", as outlined in brown in Figure 10(a).
* d = Spectral features measured from the Short-Low (SL) spectra of Mrk 266 SW
* and NE as plotted in Figures 12(c) and (b). These 1D spectra were extracted
* from the blue and red regions of the spectral map shown in Figure 10, which
* have dimensions of 5*5 SL pixels, where each pixel is 1.8"*1.8".
*Note (2): Flag as follows:
* b = These line fluxes were measured using PAHFIT and include an extinction
* correction based on a fit with a {tau}_9.7{mu}m_=0.30 dust screen, which
* corresponds to A_V_{approx}5.7mag.
* c = Equivalent width measurements for the PAH features were measured using a
* spline-fitting technique for estimating the local continuum. Spline fitting
* was also used to estimate fluxes of features in the Short-Low (SL) spectra
* of the individual galaxies which could not be measured with PAHFIT.
* See Section 2.7 for more details.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Elliptical apertures for regions of interest
integer*4 nr__1
parameter (nr__1=5) ! Number of records
character*54 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*16 Region_1 (nr__1) ! Region used to perform aperture photometry (G1)
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 the Declination
integer*4 DEd (nr__1) ! (deg) Degree of Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
real*4 a (nr__1) ! (arcsec) First parameter for elliptical aperture (1)
real*4 b (nr__1) ! (arcsec) Second parameter for elliptical aperture (1)
integer*4 PA (nr__1) ! (deg) ? Position angle
*Note (1): Used for photometric measurements, as illustrated in Figure 31.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Photometric data for Mrk 266 - Total system
integer*4 nr__2
parameter (nr__2=35) ! Number of records
character*70 ar__2 ! Full-size record
character*18 Region_2 (nr__2) ! Mrk 266 region (G1)
character*20 Band (nr__2) ! Bandpass
character*1 f_Band (nr__2) ! [a] The 70 and 160{mu}m flux densities for the
* South-West (SW) and North-East (NE) galaxies
* were estimated using the 8 and 24{mu}m
* measurements as described in Section 2.6.2
real*4 Freq (nr__2) ! (Hz) [1.5e8/1.5e18] Frequency ({nu})
real*4 Flux_1 (nr__2) ! (Jy) Flux density at Freq (f_{nu}_)
real*4 e_Flux_1 (nr__2) ! (Jy) The 1{sigma} uncertainty in Flux {sigma}_f__{nu}_
integer*4 Ref (nr__2) ! Reference; in refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! Photometric data for Mrk 266 - Components
integer*4 nr__3
parameter (nr__3=55) ! Number of records
character*70 ar__3 ! Full-size record
character*18 Region_3 (nr__3) ! Mrk 266 region (G1)
character*20 Band_1 (nr__3) ! Bandpass
character*1 f_Band_1 (nr__3) ! [a] The 70 and 160{mu}m flux densities for the
* South-West (SW) and North-East (NE) galaxies
* were estimated using the 8 and 24{mu}m
* measurements as described in Section 2.6.2
real*4 Freq_1 (nr__3) ! (Hz) [1.5e8/1.5e18] Frequency ({nu})
real*4 Flux_2 (nr__3) ! (Jy) Flux density at Freq (f_{nu}_)
real*4 e_Flux_2 (nr__3) ! (Jy) The 1{sigma} uncertainty in Flux {sigma}_f__{nu}_
integer*4 Ref_1 (nr__3) ! Reference; in refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table14.dat' ! Regions covered in published IRS measurements
integer*4 nr__4
parameter (nr__4=42) ! Number of records
character*86 ar__4 ! Full-size record
character*27 Region_4 (nr__4) ! Mrk 266 region (G1)
real*4 lambda_1 (nr__4) ! (um) [6/35] Spectral line wavelength ({mu}m)
character*16 Line_1 (nr__4) ! Spectral line identification
character*2 IRS (nr__4) ! Infrared Spectrograph (IRS) module used (1)
character*1 f_Line (nr__4) ! [ab] Flag on Line (2)
character*1 l_Flux (nr__4) ! Upper limit on Flux
real*4 Flux_3 (nr__4) ! (W/m2) ? Flux
real*4 e_Flux_3 (nr__4) ! (W/m2) ? Uncertainty in Flux
real*4 EW_1 (nr__4) ! (um) ? Equivalent width; in {mu}m
character*3 Ref_2 (nr__4) ! Reference; in refs.dat file
*Note (1): IRS module as follows:
* SL = Short-Low (5.2-14.5{mu}m, 3.7"*57", R~60-127);
* LL = Long-Low (14.0-38.0{mu}m, 10.5"*168.3", R~57-126);
* SH = 9.9-19.6{mu}m, 4.7"*11.3", R~600;
* LH = 18.7-37.3{mu}m, 11.1"*22.3", R~600.
*Note (2): Flag on spectral line identification as follows:
* a = These features are visible in the low-res spectrum published by
* Brandl et al. (2006ApJ...653.1129B);
* b = This is consistent with an [NeV] 24.32{mu}m flux of
* 1.19(+/-0.06)*10^-20^W.cm-2 measured from the same LH aperture by
* Dudik et al. (2007ApJ...664...71D; Table 2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table15.dat' ! Star cluster photometry
integer*4 nr__5
parameter (nr__5=119) ! Number of records
character*67 ar__5 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__5) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__5) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq (nr__5) ! Star cluster (candidate) number
* <[MIV2012] Mrk 266 SCNNN> in Simbad
integer*4 RAh_1 (nr__5) ! (h) Hour of right ascension (J2000) (1)
integer*4 RAm_1 (nr__5) ! (min) Minute of right ascension (J2000) (1)
real*4 RAs_1 (nr__5) ! (s) Second of right ascension (J2000) (1)
character*1 DE__1 (nr__5) ! Sign of the declination (J2000)
integer*4 DEd_1 (nr__5) ! (deg) Degree of declination (J2000) (1)
integer*4 DEm_1 (nr__5) ! (arcmin) Arcminute of declination (J2000) (1)
real*4 DEs_1 (nr__5) ! (arcsec) Arcsecond of declination (J2000) (1)
real*4 Bmag (nr__5) ! (mag) Apparent mag in the HST ACS F435W (B) filter
real*4 e_Bmag (nr__5) ! (mag) 1{sigma} uncertainty in Bmag
real*4 BMag_1 (nr__5) ! (mag) Absolute B magnitude (2)
real*4 Imag (nr__5) ! (mag) Apparent mag in the HST ACS F814W (I) filter
real*4 e_Imag (nr__5) ! (mag) 1{sigma} uncertainty in Imag
real*4 B_I (nr__5) ! (mag) B-I color index (2)
real*4 e_B_I (nr__5) ! (mag) 1{sigma} uncertainty in B-I
*Note (1): The 1{sigma} uncertainty in the coordinates is +/-0.3".
*Note (2): Corrected for foreground Galactic extinction using values from
* Schlegel et al. (1998ApJ...500..525S) as provided by NED:
* A_B_=0.056mag; A_t_=0.025mag. These data have not been corrected for
* extinction within Mrk 266 (see Section 3.6.2).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__6
parameter (nr__6=14) ! Number of records
character*77 ar__6 ! Full-size record
character*3 Ref_3 (nr__6) ! Reference code
character*19 BibCode (nr__6) ! Bibliographic code
character*20 Aut (nr__6) ! Author's name
character*32 Com (nr__6) ! Comments
C=============================================================================
C Loading file 'table3.dat' ! Measurements from low-resolution spectral mapping
C Format for file interpretation
1 format(
+ A5,1X,A1,1X,F5.2,1X,A16,F7.3,1X,F5.3,1X,A1,1X,F4.2,1X,F4.2,1X,
+ A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,38
read(1,'(A57)')ar__
read(ar__,1)
+ Region(i__),f_Region(i__),lambda(i__),Line(i__),Flux(i__),
+ e_Flux(i__),f_Flux(i__),EW(i__),e_EW(i__),f_EW(i__)
if(ar__(47:50) .EQ. '') EW(i__) = rNULL__
if(ar__(52:55) .EQ. '') e_EW(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Region(i__),f_Region(i__),lambda(i__),Line(i__),Flux(i__),
+ e_Flux(i__),f_Flux(i__),EW(i__),e_EW(i__),f_EW(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Elliptical apertures for regions of interest
C Format for file interpretation
2 format(
+ A16,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F4.1,1X,
+ F4.1,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,5
read(1,'(A54)')ar__1
read(ar__1,2)
+ Region_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),a(i__),b(i__),PA(i__)
if(ar__1(52:54) .EQ. '') PA(i__) = iNULL__
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)
+ Region_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),a(i__),b(i__),PA(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 'table12.dat' ! Photometric data for Mrk 266 - Total system
C Format for file interpretation
3 format(A18,1X,A20,1X,A1,1X,E8.3,1X,E8.3,1X,E7.2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,35
read(1,'(A70)')ar__2
read(ar__2,3)
+ Region_2(i__),Band(i__),f_Band(i__),Freq(i__),Flux_1(i__),
+ e_Flux_1(i__),Ref(i__)
c ..............Just test output...........
write(6,3)
+ Region_2(i__),Band(i__),f_Band(i__),Freq(i__),Flux_1(i__),
+ e_Flux_1(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! Photometric data for Mrk 266 - Components
C Format for file interpretation
4 format(A18,1X,A20,1X,A1,1X,E8.3,1X,E8.3,1X,E7.2,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,55
read(1,'(A70)')ar__3
read(ar__3,4)
+ Region_3(i__),Band_1(i__),f_Band_1(i__),Freq_1(i__),
+ Flux_2(i__),e_Flux_2(i__),Ref_1(i__)
c ..............Just test output...........
write(6,4)
+ Region_3(i__),Band_1(i__),f_Band_1(i__),Freq_1(i__),
+ Flux_2(i__),e_Flux_2(i__),Ref_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table14.dat' ! Regions covered in published IRS measurements
C Format for file interpretation
5 format(
+ A27,1X,F5.2,1X,A16,1X,A2,1X,A1,1X,A1,1X,E9.4,1X,E8.3,1X,F5.3,
+ 1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table14.dat')
write(6,*) '....Loading file: table14.dat'
do i__=1,42
read(1,'(A86)')ar__4
read(ar__4,5)
+ Region_4(i__),lambda_1(i__),Line_1(i__),IRS(i__),f_Line(i__),
+ l_Flux(i__),Flux_3(i__),e_Flux_3(i__),EW_1(i__),Ref_2(i__)
if(ar__4(59:67) .EQ. '') Flux_3(i__) = rNULL__
if(ar__4(69:76) .EQ. '') e_Flux_3(i__) = rNULL__
if(ar__4(78:82) .EQ. '') EW_1(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Region_4(i__),lambda_1(i__),Line_1(i__),IRS(i__),f_Line(i__),
+ l_Flux(i__),Flux_3(i__),e_Flux_3(i__),EW_1(i__),Ref_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table15.dat' ! Star cluster photometry
C Format for file interpretation
6 format(
+ I3,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,F4.2,
+ 1X,F6.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table15.dat')
write(6,*) '....Loading file: table15.dat'
do i__=1,119
read(1,'(A67)')ar__5
read(ar__5,6)
+ Seq(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),Bmag(i__),e_Bmag(i__),
+ BMag_1(i__),Imag(i__),e_Imag(i__),B_I(i__),e_B_I(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,6)
+ Seq(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),
+ DEd_1(i__),DEm_1(i__),DEs_1(i__),Bmag(i__),e_Bmag(i__),
+ BMag_1(i__),Imag(i__),e_Imag(i__),B_I(i__),e_B_I(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 'refs.dat' ! References
C Format for file interpretation
7 format(A3,1X,A19,1X,A20,1X,A32)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,14
read(1,'(A77)')ar__6
read(ar__6,7)Ref_3(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,7)Ref_3(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end