Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/269/3 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-13
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/269/3 Siena Galaxy Atlas 2020 (SGA-2020) (Moustakas+, 2023)
*================================================================================
*Siena Galaxy Atlas 2020.
* Moustakas J., Lang D., Dey A., Juneau S., Meisner A., Myers A.D.,
* Schlafly E.F., Schlegel D.J., Valdes F., Weaver B.A., Zhou R.
* <Astrophys. J. Suppl. Ser., 269, 3 (2023)>
* =2023ApJS..269....3M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'sga2020.dat' ! Ellipse-fitting results
integer*4 nr__
parameter (nr__=383620) ! Number of records
character*1136 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID (nr__) ! [2/5005247] SGA unique
* identification number (SGA_ID)
character*16 SGA2020 (nr__) ! SGA-2020 name (SGA-2020 @{ID})
* (SGA_GALAXY)
character*29 Name (nr__) ! Unique galaxy name (GALAXY)
integer*4 PGC (nr__) ! [1/6740716]? Unique identifier from
* the Principal Catalogue of Galaxies
* (Paturel+ 1989, VII/119 ;
* obsoleted by VII/237) (PGC)
real*8 RALdeg (nr__) ! (deg) Right ascension (J2000) from the
* reference indicated in Ref (RA_LEDA)
real*8 DELdeg (nr__) ! (deg) Declination (J2000) from the
* reference indicated in Ref
* (DEC_LEDA)
character*9 MType (nr__) ! Visual morphological type from
* HyperLeda (if available) (MORPHTYPE)
real*4 PA_LEDA (nr__) ! (deg) Galaxy position angle (measured
* east of north); taken from the
* reference indicated in Ref (PA_LEDA)
real*8 D25_LEDA (nr__) ! (arcmin) [0.16/63] Major-axis diameter at
* the 25mag/arcsec^2^ (optical)
* surface-brightness isophote; taken
* from the reference indicated in Ref
* (D25_LEDA)
real*4 b_a_LEDA (nr__) ! [0.03/1] Ratio of the semiminor
* axis to the semimajor axis; taken
* from the reference indicated in Ref
* (BA_LEDA)
real*8 z_LEDA (nr__) ! [-0.003/4.4]? Heliocentric redshift
* from HyperLeda; a NULL value does
* not imply that no redshift exists
* (Z_LEDA)
real*8 SBD25_LEDA (nr__) ! (mag/arcsec2) [14.5/31.7] Mean surface brightness
* based on D25-LEDA and mag-LEDA
* (SB_D25_LEDA)
real*8 mag_LEDA (nr__) ! (mag) [6.3/25] Apparent b_t_-band Vega
* magnitude (MAG_LEDA) (1)
integer*4 ByHand (nr__) ! [0/1] Boolean flag indicating
* whether one or more of RALdeg,
* DELdeg, D25-LEDA, PA-LEDA,
* b/a-LEDA, or mag-LEDA were changed
* from their published values
* (usually via visual inspection)
* while building the parent sample
* (BYHAND)
character*13 Ref (nr__) ! Reference indicating the origin of
* the object (REF) (2)
integer*4 GroupID (nr__) ! [0/531676] Unique group
* identification number (GROUP_ID)
character*29 Group (nr__) ! Group name, constructed from the
* name of its largest member (see
* Section 2.2). For isolated
* galaxies, identical to Name
* (GROUP_NAME)
integer*4 Nm (nr__) ! [1/13] Number of group members
* (i.e., group multiplicity)
* (GROUP_MULT)
integer*4 Prim (nr__) ! [0/1] Boolean flag indicating the
* primary (i.e., largest) group
* member (GROUP_PRIMARY)
real*8 RAGdeg (nr__) ! (deg) Mean R.A. of the group weighted by
* D25-LEDA (GROUP_RA)
real*8 DEGdeg (nr__) ! (deg) Mean decl. of the group weighted by
* D25-LEDA (GROUP_DEC)
real*8 GrDiam (nr__) ! (arcmin) [0.16/63] Approximate group
* diameter (see Section 2.2)
* (GROUP_DIAMETER)
character*8 BrickName (nr__) ! Name of custom Tractor "brick,"
* encoding the sky position (e.g.,
* "1126p222" is centered on
* RA=112.6,DE=+22.2) (BRICKNAME)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000) based on
* The Tractor model fit (RA)
real*8 DEdeg (nr__) ! (deg) Declination (J2000) based on The
* Tractor model fit (DEC)
real*8 D26 (nr__) ! (arcmin) [0.06/78] Major axis diameter at
* the {mu}=26mag/arcsec^2^ r-band
* isophote (D26)
character*4 r_D26 (nr__) ! Reference indicating the origin of
* the DIAM measurement: SB26, SB25,
* or LEDA (D26_REF) (3)
real*4 PA (nr__) ! (deg) Galaxy position angle (measured
* east of north), as measured from
* the ellipse moments (see Section
* 3.2.3); equivalent to PA-LEDA if
* the ellipse moments could not be
* measured (PA)
real*4 b_a (nr__) ! [0.06/1] Minor-to-major axis ratio,
* as measured from the ellipse
* moments; equivalent to b/a-LEDA if
* the ellipse moments could not be
* measured (BA)
real*8 RAMdeg (nr__) ! (deg) Light-weighted R.A. (J2000), as
* measured from the ellipse moments
* (RA_MOMENT) (4)
real*8 DEMdeg (nr__) ! (deg) Light-weighted decl. (J2000), as
* measured from the ellipse moments
* (DEC_MOMENT) (4)
real*8 SMA_Moment (nr__) ! (arcsec) [0.3/1084]? Second moment of the
* light distribution along the major
* axis based on the measured ellipse
* moments (SMA_MOMENT) (5)
real*4 gSMA50 (nr__) ! (arcsec) [0/93493]? Half-light semimajor
* axis length based on Equation-2 in
* g-band (G_SMA50)
real*4 rSMA50 (nr__) ! (arcsec) [0/62331]? Half-light semimajor
* axis length based on Equation-2 in
* r-band (R_SMA50)
real*4 zSMA50 (nr__) ! (arcsec) [0/7e+9]? Half-light semimajor
* axis length based on Equation-2 in
* z-band (Z_SMA50)
real*8 SMASB22 (nr__) ! (arcsec) [0.03/1901]? Semimajor axis length
* at the r-band 22mag/arcsec^2^
* isophote (SMA_SB22)
real*8 SMASB22_5 (nr__) ! (arcsec) [0.003/104140]? Semimajor axis
* length at the r-band
* 22.5mag/arcsec^2^ isophote
* (SMA_SB22.5)
real*8 SMASB23 (nr__) ! (arcsec) [0.005/147363]? Semimajor axis
* length at the r-band
* 23mag/arcsec^2^ isophote (SMA_SB23)
real*8 SMASB23_5 (nr__) ! (arcsec) [0.003/245391]? Semimajor axis
* length at the r-band
* 23.5mag/arcsec^2^ isophote
* (SMA_SB23.5)
real*8 SMASB24 (nr__) ! (arcsec) [0.1/733]? Semimajor axis length at
* the r-band 24mag/arcsec^2^ isophote
* (SMA_SB24)
real*8 SMASB24_5 (nr__) ! (arcsec) [0.2/867]? Semimajor axis length at
* the r-band 24.5mag/arcsec^2^
* isophote (SMA_SB24.5)
real*8 SMASB25 (nr__) ! (arcsec) [0.19/702]? Semimajor axis length
* at the r-band 25mag/arcsec^2^
* isophote (SMA_SB25)
real*8 SMASB25_5 (nr__) ! (arcsec) [0.9/774]? Semimajor axis length at
* the r-band 25.5mag/arcsec^2^
* isophote (SMA_SB25.5)
real*8 SMASB26 (nr__) ! (arcsec) [1.9/828]? Semimajor axis length at
* the r-band 26mag/arcsec^2^ isophote
* (SMA_SB26)
real*8 gmagSB22 (nr__) ! (mag) [8.8/51]? Cumulative g-band
* brightness measured within SMASB22
* (G_MAG_SB22)
real*8 rmagSB22 (nr__) ! (mag) [8/28.4]? Cumulative r-band
* brightness measured within SMASB22
* (R_MAG_SB22)
real*8 zmagSB22 (nr__) ! (mag) [7.6/43.4]? Cumulative z-band
* brightness measured within SMASB22
* (Z_MAG_SB22)
real*8 gmagSB22_5 (nr__) ! (mag) [8.8/49.1]? Cumulative g-band
* brightness measured within
* SMASB22.5 (G_MAG_SB22.5)
real*8 rmagSB22_5 (nr__) ! (mag) [8/36]? Cumulative r-band
* brightness measured within
* SMASB22.5 (R_MAG_SB22.5)
real*8 zmagSB22_5 (nr__) ! (mag) [7.5/46.3]? Cumulative z-band
* brightness measured within
* SMASB22.5 (Z_MAG_SB22.5)
real*8 gmagSB23 (nr__) ! (mag) [8.8/47]? Cumulative g-band
* brightness measured within SMASB23
* (G_MAG_SB23)
real*8 rmagSB23 (nr__) ! (mag) [8/34.7]? Cumulative r-band
* brightness measured within SMASB23
* (R_MAG_SB23)
real*8 zmagSB23 (nr__) ! (mag) [7.5/51.6]? Cumulative z-band
* brightness measured within SMASB23
* (Z_MAG_SB23)
real*8 gmagSB23_5 (nr__) ! (mag) [8.7/48.1]? Cumulative g-band
* brightness measured within
* SMASB23.5 (G_MAG_SB23.5)
real*8 rmagSB23_5 (nr__) ! (mag) [8/35]? Cumulative r-band
* brightness measured within
* SMASB23.5 (R_MAG_SB23.5)
real*8 zmagSB23_5 (nr__) ! (mag) [7.5/43.4]? Cumulative z-band
* brightness measured within
* SMASB23.5 (Z_MAG_SB23.5)
real*8 gmagSB24 (nr__) ! (mag) [8.7/46]? Cumulative g-band
* brightness measured within SMASB24
* (G_MAG_SB24)
real*8 rmagSB24 (nr__) ! (mag) [8/28.5]? Cumulative r-band
* brightness measured within SMASB24
* (R_MAG_SB24)
real*8 zmagSB24 (nr__) ! (mag) [7.4/43.4]? Cumulative z-band
* brightness measured within SMASB24
* (Z_MAG_SB24)
real*8 gmagSB24_5 (nr__) ! (mag) [8.7/45.3]? Cumulative g-band
* brightness measured within
* SMASB24.5 (G_MAG_SB24.5)
real*8 rmagSB24_5 (nr__) ! (mag) [8/27]? Cumulative r-band
* brightness measured within
* SMASB24.5 (R_MAG_SB24.5)
real*8 zmagSB24_5 (nr__) ! (mag) [7.4/44.3]? Cumulative z-band
* brightness measured within
* SMASB24.5 (Z_MAG_SB24.5)
real*8 gmagSB25 (nr__) ! (mag) [8.7/51.2]? Cumulative g-band
* brightness measured within SMASB25
* (G_MAG_SB25)
real*8 rmagSB25 (nr__) ! (mag) [8/27.3]? Cumulative r-band
* brightness measured within SMASB25
* (R_MAG_SB25)
real*8 zmagSB25 (nr__) ! (mag) [7.6/49]? Cumulative z-band
* brightness measured within SMASB25
* (Z_MAG_SB25)
real*8 gmagSB25_5 (nr__) ! (mag) [8.7/46]? Cumulative g-band
* brightness measured within
* SMASB25.5 (G_MAG_SB25.5)
real*8 rmagSB25_5 (nr__) ! (mag) [8/23.1]? Cumulative r-band
* brightness measured within
* SMASB25.5 (R_MAG_SB25.5)
real*8 zmagSB25_5 (nr__) ! (mag) [7.6/42.3]? Cumulative z-band
* brightness measured within
* SMASB25.5 (Z_MAG_SB25.5)
real*8 gmagSB26 (nr__) ! (mag) [8.6/47.1]? Cumulative g-band
* brightness measured within SMASB26
* (G_MAG_SB26)
real*8 rmagSB26 (nr__) ! (mag) [8/24.5]? Cumulative r-band
* brightness measured within SMASB26
* (R_MAG_SB26)
real*8 zmagSB26 (nr__) ! (mag) [7.6/23]? Cumulative z-band
* brightness measured within SMASB26
* (Z_MAG_SB26)
real*8 e_SMASB22 (nr__) ! (arcsec) [0.001/3763]? SMASB22 1{sigma}
* uncertainty (SMA_SB22_ERR)
real*8 e_SMASB22_5(nr__) ! (arcsec) [0.002/349384]? SMASB22.5 1{sigma}
* uncertainty (SMA_SB22.5_ERR)
real*8 e_SMASB23 (nr__) ! (arcsec) [0.002/318385]? SMASB23 1{sigma}
* uncertainty (SMA_SB23_ERR)
real*8 e_SMASB23_5(nr__) ! (arcsec) [0.003/380471]? SMASB23.5 1{sigma}
* uncertainty (SMA_SB23.5_ERR)
real*8 e_SMASB24 (nr__) ! (arcsec) [0.004/6]? SMASB24 1{sigma}
* uncertainty (SMA_SB24_ERR)
real*8 e_SMASB24_5(nr__) ! (arcsec) [0.006/8]? SMASB24.5 1{sigma}
* uncertainty (SMA_SB24.5_ERR)
real*8 e_SMASB25 (nr__) ! (arcsec) [0.009/7]? SMASB25 1{sigma}
* uncertainty (SMA_SB25_ERR)
real*8 e_SMASB25_5(nr__) ! (arcsec) [0.01/13]? SMASB25.5 1{sigma}
* uncertainty (SMA_SB25.5_ERR)
real*8 e_SMASB26 (nr__) ! (arcsec) [0.01/12]? SMASB26 1{sigma}
* uncertainty (SMA_SB26_ERR)
real*8 e_gmagSB22 (nr__) ! (mag) [0.0002/27]? gmagSB22 1{sigma}
* uncertainty (G_MAG_SB22_ERR)
real*8 e_rmagSB22 (nr__) ! (mag) [0.0002/3]? rmagSB22 1{sigma}
* uncertainty (R_MAG_SB22_ERR)
real*8 e_zmagSB22 (nr__) ! (mag) [0/4]? zmagSB22 1{sigma}
* uncertainty (Z_MAG_SB22_ERR)
real*8 e_gmagSB22_5(nr__) ! (mag) [0.0002/1]? gmagSB22.5 1{sigma}
* uncertainty (G_MAG_SB22.5_ERR)
real*8 e_rmagSB22_5(nr__) ! (mag) [0.0002/12]? rmagSB22.5 1{sigma}
* uncertainty (R_MAG_SB22.5_ERR)
real*8 e_zmagSB22_5(nr__) ! (mag) [0.0002/8]? zmagSB22.5 1{sigma}
* uncertainty (Z_MAG_SB22.5_ERR)
real*8 e_gmagSB23 (nr__) ! (mag) [0.0002/10]? gmagSB23 1{sigma}
* uncertainty (G_MAG_SB23_ERR)
real*8 e_rmagSB23 (nr__) ! (mag) [0.0002/7]? rmagSB23 1{sigma}
* uncertainty (R_MAG_SB23_ERR)
real*8 e_zmagSB23 (nr__) ! (mag) [0.0002/64]? zmagSB23 1{sigma}
* uncertainty (Z_MAG_SB23_ERR)
real*8 e_gmagSB23_5(nr__) ! (mag) [0.0002/7]? gmagSB23.5 1{sigma}
* uncertainty (G_MAG_SB23.5_ERR)
real*8 e_rmagSB23_5(nr__) ! (mag) [0.0001/9]? rmagSB23.5 1{sigma}
* uncertainty (R_MAG_SB23.5_ERR)
real*8 e_zmagSB23_5(nr__) ! (mag) [0.0002/9]? zmagSB23.5 1{sigma}
* uncertainty (Z_MAG_SB23.5_ERR)
real*8 e_gmagSB24 (nr__) ! (mag) [0.0002/2]? gmagSB24 1{sigma}
* uncertainty (G_MAG_SB24_ERR)
real*8 e_rmagSB24 (nr__) ! (mag) [0.0001/3]? rmagSB24 1{sigma}
* uncertainty (R_MAG_SB24_ERR)
real*8 e_zmagSB24 (nr__) ! (mag) [0.0002/74]? zmagSB24 1{sigma}
* uncertainty (Z_MAG_SB24_ERR)
real*8 e_gmagSB24_5(nr__) ! (mag) [0.0002/18]? gmagSB24.5 1{sigma}
* uncertainty (G_MAG_SB24.5_ERR)
real*8 e_rmagSB24_5(nr__) ! (mag) [0.0001/3]? rmagSB24.5 1{sigma}
* uncertainty (R_MAG_SB24.5_ERR)
real*8 e_zmagSB24_5(nr__) ! (mag) [0.0001/54]? zmagSB24.5 1{sigma}
* uncertainty (Z_MAG_SB24.5_ERR)
real*8 e_gmagSB25 (nr__) ! (mag) [0.0002/16]? gmagSB25 1{sigma}
* uncertainty (G_MAG_SB25_ERR)
real*8 e_rmagSB25 (nr__) ! (mag) [0.0001/5]? rmagSB25 1{sigma}
* uncertainty (R_MAG_SB25_ERR)
real*8 e_zmagSB25 (nr__) ! (mag) [0.0001/9]? zmagSB25 1{sigma}
* uncertainty (Z_MAG_SB25_ERR)
real*8 e_gmagSB25_5(nr__) ! (mag) [0.0002/23]? gmagSB25.5 1{sigma}
* uncertainty (G_MAG_SB25.5_ERR)
real*8 e_rmagSB25_5(nr__) ! (mag) [0.0001/1]? rmagSB25.5 1{sigma}
* uncertainty (R_MAG_SB25.5_ERR)
real*8 e_zmagSB25_5(nr__) ! (mag) [0.0001/12]? zmagSB25.5 1{sigma}
* uncertainty (Z_MAG_SB25.5_ERR)
real*8 e_gmagSB26 (nr__) ! (mag) [0.0002/28]? gmagSB26 1{sigma}
* uncertainty (G_MAG_SB26_ERR)
real*8 e_rmagSB26 (nr__) ! (mag) [0.0001/7]? rmagSB26 1{sigma}
* uncertainty (R_MAG_SB26_ERR)
real*8 e_zmagSB26 (nr__) ! (mag) [0.0001/133]? zmagSB26 1{sigma}
* uncertainty (Z_MAG_SB26_ERR)
real*8 gCOG_m1 (nr__) ! (mag) [7.6/45.1]? Best-fitting
* curve-of-growth parameter m_tot_
* from Equation-1 in g-band
* (G_COG_PARAMS_MTOT) (6)
real*4 gCOG_m0 (nr__) ! (mag) [1.5e-07/4.3e+11]? Best-fitting
* curve-of-growth parameter m0 from
* Equation-1 in g-band
* (G_COG_PARAMS_M0) (6)
real*4 gCOG_alp1 (nr__) ! [3.5e-10/6.2e+15]? Best-fitting
* curve-of-growth parameter {alpha}1
* from Equation-1 in g-band
* (G_COG_PARAMS_ALPHA1) (6)
real*8 gCOG_alp2 (nr__) ! [0/169]? Best-fitting
* curve-of-growth parameter {alpha}2
* from Equation-1 in g-band
* (G_COG_PARAMS_ALPHA2) (6)
real*8 gCOG_Chi2 (nr__) ! [0/463626]? {Chi}^2^ of the fit to
* the curve of growth in g-band
* (see Section 3.2.3)
* (G_COG_PARAMS_CHI2)
real*8 rCOG_m1 (nr__) ! (mag) [7.5/23.5]? Best-fitting
* curve-of-growth parameter m_tot_
* from Equation-1 in r-band
* (R_COG_PARAMS_MTOT) (6)
real*4 rCOG_m0 (nr__) ! (mag) [5.6e-7/2.1e+13]? Best-fitting
* curve-of-growth parameter m0 from
* Equation-1 in r-band
* (R_COG_PARAMS_M0) (6)
real*4 rCOG_alp1 (nr__) ! [7.5e-9/1.5e+16]? Best-fitting
* curve-of-growth parameter {alpha}1
* from Equation-1 in r-band
* (R_COG_PARAMS_ALPHA1) (6)
real*8 rCOG_alp2 (nr__) ! [0/183]? Best-fitting
* curve-of-growth parameter {alpha}2
* from Equation-1 in r-band
* (R_COG_PARAMS_ALPHA2) (6)
real*8 rCOG_Chi2 (nr__) ! [0/572291]? {Chi}^2^ of the fit to
* the curve of growth in r-band
* (see Section 3.2.3)
* (R_COG_PARAMS_CHI2)
real*8 zCOG_m1 (nr__) ! (mag) [6.9/40.4]? Best-fitting
* curve-of-growth parameter m_tot_
* from Equation-1 in z-band
* (Z_COG_PARAMS_MTOT) (6)
real*4 zCOG_m0 (nr__) ! (mag) [7.6e-8/4e+12]? Best-fitting
* curve-of-growth parameter m0 from
* Equation-1 in z-band
* (Z_COG_PARAMS_M0) (6)
real*4 zCOG_alp1 (nr__) ! [2.7e-10/6e+15]? Best-fitting
* curve-of-growth parameter {alpha}1
* from Equation-1 in z-band
* (Z_COG_PARAMS_ALPHA1) (6)
real*8 zCOG_alp2 (nr__) ! [0/197]? Best-fitting
* curve-of-growth parameter {alpha}2
* from Equation-1 in z-band
* (Z_COG_PARAMS_ALPHA2) (6)
real*8 zCOG_Chi2 (nr__) ! [0/405272]? {Chi}^2^ of the fit to
* the curve of growth in z-band
* (see Section 3.2.3)
* (Z_COG_PARAMS_CHI2)
integer*4 Fit (nr__) ! [0/32] Ellipse-fitting Bitmasks
* described in Table 3 (ELLIPSEBIT)
* (7)
*Note (1): This magnitude estimate is heterogeneous in both bandpass and
* aperture but, for most galaxies, it is measured in the B band; use
* with care.
*Note (2): References as follows:
* LEDA-20181114 = the HyperLeda extragalactic database;
* http://leda.univ-lyon1.fr/ (380532 occurrences)
* DR8 = the Legacy Surveys Data Release 8;
* http://legacysurvey.org/dr8 (2889 occurrences)
* OpenNGC = OpenNGC catalog; http://github.com/mattiaverga/OpenNGC
* (142 occurrences)
* LGDWARFS = Local Group dwarf galaxies from McConnachie (2012, J/AJ/144/4);
* (29 occurrences)
* RC3 = The third Reference Cat. of Bright Galaxies; see VII/155
* (28 occurrences)
*
* See Section 2.1.
*Note (3): By default, we infer D26 from SMASB26. However, if the r-band
* surface-brightness profile could not be measured at this level, we
* estimate D26 as 2.5xSMASB25 or 1.5xD25-LEDA, in that order of priority.
*Note (4): Light-weighted R.A. and decl. (J2000), as measured from the ellipse
* moments. Equivalent to RA_X0 in Table 5 but set to RALdeg if ellipse
* fitting failed or was not carried out.
*Note (5): Second moment of the light distribution along the major axis based
* on the measured ellipse moments. Equivalent to MAJORAXIS in Table 5
* but converted to arcseconds. If ellipse-fitting or curve-of-growth
* modeling failed or was not attempted, then the column value is NULL.
*Note (6): We fit the curve of growth, m(r), using the following empirical model:
*
* m(r)=m_tot_+m_0_log_e_[1+{alpha}1(r/r_0_)^-{alpha}2^]; Equation-1
*
* where m_tot_, m_0_, {alpha}1, {alpha}2, and r_0_ are constant
* parameters of the model, and r is the semimajor axis in arcseconds.
* See Section 3.2.3.
*Note (7): Ellipse-fitting Bitmasks as in Table 3:
* ---------------------------------------------------------------------------
* Bit Bit name Definition
* Number
* ---------------------------------------------------------------------------
* 0 --- Not used; ignore
* 1 REX_TOOSMALL Ellipse fit skipped; galaxy classified as
* too-small type REX
* 2 NOTREX_TOOSMALL Ellipse fit skipped; galaxy classified as
* too-small type EXP, DEV, or SER
* 3 FAILED Ellipse fitting was attempted but failed
* 4 NOTFIT Ellipse fitting was not attempted
* 5 REJECTED Ellipse-fitting results were rejected based
* on visual inspection
* ---------------------------------------------------------------------------
* Bitmask encoding various reasons why ellipse fitting failed or was not
* attempted; see Section 3.2.3 for additional details.
* See also: http://portal.nersc.gov/project/cosmo/temp/ioannis/sga-dr9/sga/
* sga2020/#bitmasks
C=============================================================================
C Loading file 'sga2020.dat' ! Ellipse-fitting results
C Format for file interpretation
1 format(
+ I7,1X,A16,1X,A29,1X,I7,1X,F11.7,1X,F11.7,1X,A9,1X,F6.2,1X,
+ F7.4,1X,F6.4,1X,F9.5,1X,F7.4,1X,F7.4,1X,I1,1X,A13,1X,I6,1X,
+ A29,1X,I2,1X,I1,1X,F11.7,1X,F11.7,1X,F7.4,1X,A8,1X,F11.7,1X,
+ F11.7,1X,F7.4,1X,A4,1X,F6.2,1X,F6.4,1X,F11.7,1X,F11.7,1X,F9.4,
+ 1X,E11.6,1X,E11.6,1X,E11.6,1X,F9.4,1X,F11.4,1X,F11.4,1X,F11.4,
+ 1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F7.4,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F9.4,1X,F11.4,1X,F11.4,1X,F11.4,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
+ F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
+ F8.4,1X,F7.4,1X,E11.6,1X,E11.6,1X,F8.4,1X,F11.4,1X,F7.4,1X,
+ E11.6,1X,E11.6,1X,F8.4,1X,F11.4,1X,F7.4,1X,E11.6,1X,E11.6,1X,
+ F8.4,1X,F11.4,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'sga2020.dat')
write(6,*) '....Loading file: sga2020.dat'
do i__=1,383620
read(1,'(A1136)')ar__
read(ar__,1)
+ ID(i__),SGA2020(i__),Name(i__),PGC(i__),RALdeg(i__),
+ DELdeg(i__),MType(i__),PA_LEDA(i__),D25_LEDA(i__),
+ b_a_LEDA(i__),z_LEDA(i__),SBD25_LEDA(i__),mag_LEDA(i__),
+ ByHand(i__),Ref(i__),GroupID(i__),Group(i__),Nm(i__),
+ Prim(i__),RAGdeg(i__),DEGdeg(i__),GrDiam(i__),BrickName(i__),
+ RAdeg(i__),DEdeg(i__),D26(i__),r_D26(i__),PA(i__),b_a(i__),
+ RAMdeg(i__),DEMdeg(i__),SMA_Moment(i__),gSMA50(i__),
+ rSMA50(i__),zSMA50(i__),SMASB22(i__),SMASB22_5(i__),
+ SMASB23(i__),SMASB23_5(i__),SMASB24(i__),SMASB24_5(i__),
+ SMASB25(i__),SMASB25_5(i__),SMASB26(i__),gmagSB22(i__),
+ rmagSB22(i__),zmagSB22(i__),gmagSB22_5(i__),rmagSB22_5(i__),
+ zmagSB22_5(i__),gmagSB23(i__),rmagSB23(i__),zmagSB23(i__),
+ gmagSB23_5(i__),rmagSB23_5(i__),zmagSB23_5(i__),gmagSB24(i__),
+ rmagSB24(i__),zmagSB24(i__),gmagSB24_5(i__),rmagSB24_5(i__),
+ zmagSB24_5(i__),gmagSB25(i__),rmagSB25(i__),zmagSB25(i__),
+ gmagSB25_5(i__),rmagSB25_5(i__),zmagSB25_5(i__),gmagSB26(i__),
+ rmagSB26(i__),zmagSB26(i__),e_SMASB22(i__),e_SMASB22_5(i__),
+ e_SMASB23(i__),e_SMASB23_5(i__),e_SMASB24(i__),
+ e_SMASB24_5(i__),e_SMASB25(i__),e_SMASB25_5(i__),
+ e_SMASB26(i__),e_gmagSB22(i__),e_rmagSB22(i__),
+ e_zmagSB22(i__),e_gmagSB22_5(i__),e_rmagSB22_5(i__),
+ e_zmagSB22_5(i__),e_gmagSB23(i__),e_rmagSB23(i__),
+ e_zmagSB23(i__),e_gmagSB23_5(i__),e_rmagSB23_5(i__),
+ e_zmagSB23_5(i__),e_gmagSB24(i__),e_rmagSB24(i__),
+ e_zmagSB24(i__),e_gmagSB24_5(i__),e_rmagSB24_5(i__),
+ e_zmagSB24_5(i__),e_gmagSB25(i__),e_rmagSB25(i__),
+ e_zmagSB25(i__),e_gmagSB25_5(i__),e_rmagSB25_5(i__),
+ e_zmagSB25_5(i__),e_gmagSB26(i__),e_rmagSB26(i__),
+ e_zmagSB26(i__),gCOG_m1(i__),gCOG_m0(i__),gCOG_alp1(i__),
+ gCOG_alp2(i__),gCOG_Chi2(i__),rCOG_m1(i__),rCOG_m0(i__),
+ rCOG_alp1(i__),rCOG_alp2(i__),rCOG_Chi2(i__),zCOG_m1(i__),
+ zCOG_m0(i__),zCOG_alp1(i__),zCOG_alp2(i__),zCOG_Chi2(i__),
+ Fit(i__)
if(ar__(56:62) .EQ. '') PGC(i__) = iNULL__
if(ar__(120:128) .EQ. '') z_LEDA(i__) = rNULL__
if(ar__(320:328) .EQ. '') SMA_Moment(i__) = rNULL__
if(ar__(330:340) .EQ. '') gSMA50(i__) = rNULL__
if(ar__(342:352) .EQ. '') rSMA50(i__) = rNULL__
if(ar__(354:364) .EQ. '') zSMA50(i__) = rNULL__
if(ar__(366:374) .EQ. '') SMASB22(i__) = rNULL__
if(ar__(376:386) .EQ. '') SMASB22_5(i__) = rNULL__
if(ar__(388:398) .EQ. '') SMASB23(i__) = rNULL__
if(ar__(400:410) .EQ. '') SMASB23_5(i__) = rNULL__
if(ar__(412:419) .EQ. '') SMASB24(i__) = rNULL__
if(ar__(421:428) .EQ. '') SMASB24_5(i__) = rNULL__
if(ar__(430:437) .EQ. '') SMASB25(i__) = rNULL__
if(ar__(439:446) .EQ. '') SMASB25_5(i__) = rNULL__
if(ar__(448:455) .EQ. '') SMASB26(i__) = rNULL__
if(ar__(457:463) .EQ. '') gmagSB22(i__) = rNULL__
if(ar__(465:471) .EQ. '') rmagSB22(i__) = rNULL__
if(ar__(473:479) .EQ. '') zmagSB22(i__) = rNULL__
if(ar__(481:487) .EQ. '') gmagSB22_5(i__) = rNULL__
if(ar__(489:495) .EQ. '') rmagSB22_5(i__) = rNULL__
if(ar__(497:503) .EQ. '') zmagSB22_5(i__) = rNULL__
if(ar__(505:511) .EQ. '') gmagSB23(i__) = rNULL__
if(ar__(513:519) .EQ. '') rmagSB23(i__) = rNULL__
if(ar__(521:527) .EQ. '') zmagSB23(i__) = rNULL__
if(ar__(529:535) .EQ. '') gmagSB23_5(i__) = rNULL__
if(ar__(537:543) .EQ. '') rmagSB23_5(i__) = rNULL__
if(ar__(545:551) .EQ. '') zmagSB23_5(i__) = rNULL__
if(ar__(553:559) .EQ. '') gmagSB24(i__) = rNULL__
if(ar__(561:567) .EQ. '') rmagSB24(i__) = rNULL__
if(ar__(569:575) .EQ. '') zmagSB24(i__) = rNULL__
if(ar__(577:583) .EQ. '') gmagSB24_5(i__) = rNULL__
if(ar__(585:591) .EQ. '') rmagSB24_5(i__) = rNULL__
if(ar__(593:599) .EQ. '') zmagSB24_5(i__) = rNULL__
if(ar__(601:607) .EQ. '') gmagSB25(i__) = rNULL__
if(ar__(609:615) .EQ. '') rmagSB25(i__) = rNULL__
if(ar__(617:623) .EQ. '') zmagSB25(i__) = rNULL__
if(ar__(625:631) .EQ. '') gmagSB25_5(i__) = rNULL__
if(ar__(633:639) .EQ. '') rmagSB25_5(i__) = rNULL__
if(ar__(641:647) .EQ. '') zmagSB25_5(i__) = rNULL__
if(ar__(649:655) .EQ. '') gmagSB26(i__) = rNULL__
if(ar__(657:663) .EQ. '') rmagSB26(i__) = rNULL__
if(ar__(665:671) .EQ. '') zmagSB26(i__) = rNULL__
if(ar__(673:681) .EQ. '') e_SMASB22(i__) = rNULL__
if(ar__(683:693) .EQ. '') e_SMASB22_5(i__) = rNULL__
if(ar__(695:705) .EQ. '') e_SMASB23(i__) = rNULL__
if(ar__(707:717) .EQ. '') e_SMASB23_5(i__) = rNULL__
if(ar__(719:725) .EQ. '') e_SMASB24(i__) = rNULL__
if(ar__(727:733) .EQ. '') e_SMASB24_5(i__) = rNULL__
if(ar__(735:741) .EQ. '') e_SMASB25(i__) = rNULL__
if(ar__(743:749) .EQ. '') e_SMASB25_5(i__) = rNULL__
if(ar__(751:757) .EQ. '') e_SMASB26(i__) = rNULL__
if(ar__(759:765) .EQ. '') e_gmagSB22(i__) = rNULL__
if(ar__(767:773) .EQ. '') e_rmagSB22(i__) = rNULL__
if(ar__(775:781) .EQ. '') e_zmagSB22(i__) = rNULL__
if(ar__(783:789) .EQ. '') e_gmagSB22_5(i__) = rNULL__
if(ar__(791:797) .EQ. '') e_rmagSB22_5(i__) = rNULL__
if(ar__(799:805) .EQ. '') e_zmagSB22_5(i__) = rNULL__
if(ar__(807:813) .EQ. '') e_gmagSB23(i__) = rNULL__
if(ar__(815:821) .EQ. '') e_rmagSB23(i__) = rNULL__
if(ar__(823:829) .EQ. '') e_zmagSB23(i__) = rNULL__
if(ar__(831:837) .EQ. '') e_gmagSB23_5(i__) = rNULL__
if(ar__(839:845) .EQ. '') e_rmagSB23_5(i__) = rNULL__
if(ar__(847:853) .EQ. '') e_zmagSB23_5(i__) = rNULL__
if(ar__(855:861) .EQ. '') e_gmagSB24(i__) = rNULL__
if(ar__(863:869) .EQ. '') e_rmagSB24(i__) = rNULL__
if(ar__(871:877) .EQ. '') e_zmagSB24(i__) = rNULL__
if(ar__(879:885) .EQ. '') e_gmagSB24_5(i__) = rNULL__
if(ar__(887:893) .EQ. '') e_rmagSB24_5(i__) = rNULL__
if(ar__(895:901) .EQ. '') e_zmagSB24_5(i__) = rNULL__
if(ar__(903:909) .EQ. '') e_gmagSB25(i__) = rNULL__
if(ar__(911:917) .EQ. '') e_rmagSB25(i__) = rNULL__
if(ar__(919:925) .EQ. '') e_zmagSB25(i__) = rNULL__
if(ar__(927:933) .EQ. '') e_gmagSB25_5(i__) = rNULL__
if(ar__(935:941) .EQ. '') e_rmagSB25_5(i__) = rNULL__
if(ar__(943:949) .EQ. '') e_zmagSB25_5(i__) = rNULL__
if(ar__(951:957) .EQ. '') e_gmagSB26(i__) = rNULL__
if(ar__(959:965) .EQ. '') e_rmagSB26(i__) = rNULL__
if(ar__(967:974) .EQ. '') e_zmagSB26(i__) = rNULL__
if(ar__(976:982) .EQ. '') gCOG_m1(i__) = rNULL__
if(ar__(984:994) .EQ. '') gCOG_m0(i__) = rNULL__
if(ar__(996:1006) .EQ. '') gCOG_alp1(i__) = rNULL__
if(ar__(1008:1015) .EQ. '') gCOG_alp2(i__) = rNULL__
if(ar__(1017:1027) .EQ. '') gCOG_Chi2(i__) = rNULL__
if(ar__(1029:1035) .EQ. '') rCOG_m1(i__) = rNULL__
if(ar__(1037:1047) .EQ. '') rCOG_m0(i__) = rNULL__
if(ar__(1049:1059) .EQ. '') rCOG_alp1(i__) = rNULL__
if(ar__(1061:1068) .EQ. '') rCOG_alp2(i__) = rNULL__
if(ar__(1070:1080) .EQ. '') rCOG_Chi2(i__) = rNULL__
if(ar__(1082:1088) .EQ. '') zCOG_m1(i__) = rNULL__
if(ar__(1090:1100) .EQ. '') zCOG_m0(i__) = rNULL__
if(ar__(1102:1112) .EQ. '') zCOG_alp1(i__) = rNULL__
if(ar__(1114:1121) .EQ. '') zCOG_alp2(i__) = rNULL__
if(ar__(1123:1133) .EQ. '') zCOG_Chi2(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),SGA2020(i__),Name(i__),PGC(i__),RALdeg(i__),
+ DELdeg(i__),MType(i__),PA_LEDA(i__),D25_LEDA(i__),
+ b_a_LEDA(i__),z_LEDA(i__),SBD25_LEDA(i__),mag_LEDA(i__),
+ ByHand(i__),Ref(i__),GroupID(i__),Group(i__),Nm(i__),
+ Prim(i__),RAGdeg(i__),DEGdeg(i__),GrDiam(i__),BrickName(i__),
+ RAdeg(i__),DEdeg(i__),D26(i__),r_D26(i__),PA(i__),b_a(i__),
+ RAMdeg(i__),DEMdeg(i__),SMA_Moment(i__),gSMA50(i__),
+ rSMA50(i__),zSMA50(i__),SMASB22(i__),SMASB22_5(i__),
+ SMASB23(i__),SMASB23_5(i__),SMASB24(i__),SMASB24_5(i__),
+ SMASB25(i__),SMASB25_5(i__),SMASB26(i__),gmagSB22(i__),
+ rmagSB22(i__),zmagSB22(i__),gmagSB22_5(i__),rmagSB22_5(i__),
+ zmagSB22_5(i__),gmagSB23(i__),rmagSB23(i__),zmagSB23(i__),
+ gmagSB23_5(i__),rmagSB23_5(i__),zmagSB23_5(i__),gmagSB24(i__),
+ rmagSB24(i__),zmagSB24(i__),gmagSB24_5(i__),rmagSB24_5(i__),
+ zmagSB24_5(i__),gmagSB25(i__),rmagSB25(i__),zmagSB25(i__),
+ gmagSB25_5(i__),rmagSB25_5(i__),zmagSB25_5(i__),gmagSB26(i__),
+ rmagSB26(i__),zmagSB26(i__),e_SMASB22(i__),e_SMASB22_5(i__),
+ e_SMASB23(i__),e_SMASB23_5(i__),e_SMASB24(i__),
+ e_SMASB24_5(i__),e_SMASB25(i__),e_SMASB25_5(i__),
+ e_SMASB26(i__),e_gmagSB22(i__),e_rmagSB22(i__),
+ e_zmagSB22(i__),e_gmagSB22_5(i__),e_rmagSB22_5(i__),
+ e_zmagSB22_5(i__),e_gmagSB23(i__),e_rmagSB23(i__),
+ e_zmagSB23(i__),e_gmagSB23_5(i__),e_rmagSB23_5(i__),
+ e_zmagSB23_5(i__),e_gmagSB24(i__),e_rmagSB24(i__),
+ e_zmagSB24(i__),e_gmagSB24_5(i__),e_rmagSB24_5(i__),
+ e_zmagSB24_5(i__),e_gmagSB25(i__),e_rmagSB25(i__),
+ e_zmagSB25(i__),e_gmagSB25_5(i__),e_rmagSB25_5(i__),
+ e_zmagSB25_5(i__),e_gmagSB26(i__),e_rmagSB26(i__),
+ e_zmagSB26(i__),gCOG_m1(i__),gCOG_m0(i__),gCOG_alp1(i__),
+ gCOG_alp2(i__),gCOG_Chi2(i__),rCOG_m1(i__),rCOG_m0(i__),
+ rCOG_alp1(i__),rCOG_alp2(i__),rCOG_Chi2(i__),zCOG_m1(i__),
+ zCOG_m0(i__),zCOG_alp1(i__),zCOG_alp2(i__),zCOG_Chi2(i__),
+ Fit(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end