Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/141/113 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 2024-Apr-20 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/A+AS/141/113 Coma extensive H photometry (Andreon+, 2000) *================================================================================ *Extensive near-infrared (H-band) photometry in Coma * Andreon S., Pello R., Davoust E., Dominguez R., Poulain P. * <Astron. Astrophys. Suppl. Ser. 141, 113 (2000)> * =2000A&AS..141..113A (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! *List of detected sources in the region "Patch1" integer*4 nr__ parameter (nr__=77) ! Number of records character*143 ar__ ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 RAh (nr__) ! (h) Right ascension (J2000) integer*4 RAm (nr__) ! (min) Right ascension (J2000) real*4 RAs (nr__) ! (s) Right ascension (J2000) character*1 DE_ (nr__) ! Declination sign (J2000) integer*4 DEd (nr__) ! (deg) Declination (J2000) integer*4 DEm (nr__) ! (arcmin) Declination (J2000) integer*4 DEs (nr__) ! (arcsec) Declination (J2000) real*4 Hmag5 (nr__) ! (mag) H magnitude within the aperture 5" real*4 e_Hmag5 (nr__) ! (mag) rms uncertainty on Hmag5 (1) real*4 Hmag10 (nr__) ! (mag) H magnitude within the aperture 10" real*4 e_Hmag10 (nr__) ! (mag) rms uncertainty on Hmag10 (1) real*4 Hmag13 (nr__) ! (mag) ?=99.99 H magnitude within the aperture 13" (2) real*4 e_Hmag13 (nr__) ! (mag) ?=9.99 rms uncertainty on Hmag13 (1) real*4 Hmag14_5 (nr__) ! (mag) ?=99.99 H magnitude within the aperture 14.5" real*4 e_Hmag14_5 (nr__) ! (mag) ?=9.99 rms uncertainty on Hmag14.5 (1) real*4 Hmag18_5 (nr__) ! (mag) ?=99.99 H magnitude within the aperture 18.5" real*4 e_Hmag18_5 (nr__) ! (mag) ?=9.99 rms uncertainty on Hmag18.5 (1) real*4 Isomag (nr__) ! (mag) Isophotal magnitude at {mu}=22mag/arcsec^2^ real*4 e_Isomag (nr__) ! (mag) rms uncertainty on Isomag real*4 Kronmag (nr__) ! (mag) Kron magnitude (3) real*4 e_Kronmag (nr__) ! (mag) rms uncertainty on Kronmag integer*4 Area (nr__) ! (arcsec) Isophotal area (at {mu}=22mag/arcsec^2^) real*4 ell (nr__) ! Ellipticity integer*4 PA (nr__) ! (deg) Position angle (from North to East) real*4 Radius (nr__) ! (arcsec) Kron radius, or 1.8 arcsec if the * latter is larger real*4 Stellar (nr__) ! Stellarity index as defined by Sextractor (4) character*19 Name (nr__) ! Name (5) *Note (1): Errors computed by Sextractor. Listed errors do not include * photometric zero-point errors (see Sect. 3.7), those arising from the * compositing, and implicitly assume a gaussian or poissonian error * distribution. *Note (2): "99.99" means that no measure is available (usually because the * object is too small or too faint for its flux to be measured through a * large aperture). *Note (3): Kron magnitudes are computed by integrating the flux in an area * 2.5 times larger than the Kron area, with a minimum aperture radius of * 1.8 arcsec. *Note (4): When Index <~0.3, objects are identified as galaxies, whereas * Index >~0.7 (for 0404 and 0504) and Index >~0.8 (for patch1) * correspond to stars. The difference between fields is due to seeing * conditions. *Note (5): the acronyms are: * 5C: Wilson, 1970MNRAS.151....1W * CSO: Sanduleak & Pesch, 1990ApJS...72..291S * GMP: Godwin, Metcalfe & Peach (1983, Cat. <VII/42>) * NGC: Dreyer, See Cat. <VII/118> * NGP9: Odewahn & Aldering, 1995, Cat. <VII/204> * USNO: Monet et al., 1999, The USNO-A1 guide star catalog, Washington DC, * US Naval Observatory, see Cat. <I/243>, USNO-A1 * and Cat. <I/252>, USNO-A2.0 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! *List of detected sources in the region "0404" integer*4 nr__1 parameter (nr__1=124) ! Number of records character*143 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 RAh_1 (nr__1) ! (h) Right ascension (J2000) integer*4 RAm_1 (nr__1) ! (min) Right ascension (J2000) real*4 RAs_1 (nr__1) ! (s) Right ascension (J2000) character*1 DE__1 (nr__1) ! Declination sign (J2000) integer*4 DEd_1 (nr__1) ! (deg) Declination (J2000) integer*4 DEm_1 (nr__1) ! (arcmin) Declination (J2000) integer*4 DEs_1 (nr__1) ! (arcsec) Declination (J2000) real*4 Hmag5_1 (nr__1) ! (mag) H magnitude within the aperture 5" real*4 e_Hmag5_1 (nr__1) ! (mag) rms uncertainty on Hmag5 (1) real*4 Hmag10_1 (nr__1) ! (mag) H magnitude within the aperture 10" real*4 e_Hmag10_1 (nr__1) ! (mag) rms uncertainty on Hmag10 (1) real*4 Hmag13_1 (nr__1) ! (mag) ?=99.99 H magnitude within the aperture 13" (2) real*4 e_Hmag13_1 (nr__1) ! (mag) ?=9.99 rms uncertainty on Hmag13 (1) real*4 Hmag14_5_1 (nr__1) ! (mag) ?=99.99 H magnitude within the aperture 14.5" real*4 e_Hmag14_5_1(nr__1) ! (mag) ?=9.99 rms uncertainty on Hmag14.5 (1) real*4 Hmag18_5_1 (nr__1) ! (mag) ?=99.99 H magnitude within the aperture 18.5" real*4 e_Hmag18_5_1(nr__1) ! (mag) ?=9.99 rms uncertainty on Hmag18.5 (1) real*4 Isomag_1 (nr__1) ! (mag) Isophotal magnitude at {mu}=22mag/arcsec^2^ real*4 e_Isomag_1 (nr__1) ! (mag) rms uncertainty on Isomag real*4 Kronmag_1 (nr__1) ! (mag) Kron magnitude (3) real*4 e_Kronmag_1(nr__1) ! (mag) rms uncertainty on Kronmag integer*4 Area_1 (nr__1) ! (arcsec) Isophotal area (at {mu}=22mag/arcsec^2^) real*4 ell_1 (nr__1) ! Ellipticity integer*4 PA_1 (nr__1) ! (deg) Position angle (from North to East) real*4 Radius_1 (nr__1) ! (arcsec) Kron radius, or 1.8 arcsec if the * latter is larger real*4 Stellar_1 (nr__1) ! Stellarity index as defined by Sextractor (4) character*19 Name_1 (nr__1) ! Name (5) *Note (1): Errors computed by Sextractor. Listed errors do not include * photometric zero-point errors (see Sect. 3.7), those arising from the * compositing, and implicitly assume a gaussian or poissonian error * distribution. *Note (2): "99.99" means that no measure is available (usually because the * object is too small or too faint for its flux to be measured through a * large aperture). *Note (3): Kron magnitudes are computed by integrating the flux in an area * 2.5 times larger than the Kron area, with a minimum aperture radius of * 1.8 arcsec. *Note (4): When Index <~0.3, objects are identified as galaxies, whereas * Index >~0.7 (for 0404 and 0504) and Index >~0.8 (for patch1) * correspond to stars. The difference between fields is due to seeing * conditions. *Note (5): the acronyms are: * 5C: Wilson, 1970MNRAS.151....1W * CSO: Sanduleak & Pesch, 1990ApJS...72..291S * GMP: Godwin, Metcalfe & Peach (1983, Cat. <VII/42>) * NGC: Dreyer, See Cat. <VII/118> * NGP9: Odewahn & Aldering, 1995, Cat. <VII/204> * USNO: Monet et al., 1999, The USNO-A1 guide star catalog, Washington DC, * US Naval Observatory, see Cat. <I/243>, USNO-A1 * and Cat. <I/252>, USNO-A2.0 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table5.dat' ! *List of detected sources in the region "0504" integer*4 nr__2 parameter (nr__2=88) ! Number of records character*143 ar__2 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension J2000 real*8 DEdeg_2 (nr__2) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) integer*4 RAh_2 (nr__2) ! (h) Right ascension (J2000) integer*4 RAm_2 (nr__2) ! (min) Right ascension (J2000) real*4 RAs_2 (nr__2) ! (s) Right ascension (J2000) character*1 DE__2 (nr__2) ! Declination sign (J2000) integer*4 DEd_2 (nr__2) ! (deg) Declination (J2000) integer*4 DEm_2 (nr__2) ! (arcmin) Declination (J2000) integer*4 DEs_2 (nr__2) ! (arcsec) Declination (J2000) real*4 Hmag5_2 (nr__2) ! (mag) H magnitude within the aperture 5" real*4 e_Hmag5_2 (nr__2) ! (mag) rms uncertainty on Hmag5 (1) real*4 Hmag10_2 (nr__2) ! (mag) H magnitude within the aperture 10" real*4 e_Hmag10_2 (nr__2) ! (mag) rms uncertainty on Hmag10 (1) real*4 Hmag13_2 (nr__2) ! (mag) ?=99.99 H magnitude within the aperture 13" (2) real*4 e_Hmag13_2 (nr__2) ! (mag) ?=9.99 rms uncertainty on Hmag13 (1) real*4 Hmag14_5_2 (nr__2) ! (mag) ?=99.99 H magnitude within the aperture 14.5" real*4 e_Hmag14_5_2(nr__2) ! (mag) ?=9.99 rms uncertainty on Hmag14.5 (1) real*4 Hmag18_5_2 (nr__2) ! (mag) ?=99.99 H magnitude within the aperture 18.5" real*4 e_Hmag18_5_2(nr__2) ! (mag) ?=9.99 rms uncertainty on Hmag18.5 (1) real*4 Isomag_2 (nr__2) ! (mag) Isophotal magnitude at {mu}=22mag/arcsec^2^ real*4 e_Isomag_2 (nr__2) ! (mag) rms uncertainty on Isomag real*4 Kronmag_2 (nr__2) ! (mag) Kron magnitude (3) real*4 e_Kronmag_2(nr__2) ! (mag) rms uncertainty on Kronmag integer*4 Area_2 (nr__2) ! (arcsec) Isophotal area (at {mu}=22mag/arcsec^2^) real*4 ell_2 (nr__2) ! Ellipticity integer*4 PA_2 (nr__2) ! (deg) Position angle (from North to East) real*4 Radius_2 (nr__2) ! (arcsec) Kron radius, or 1.8 arcsec if the * latter is larger real*4 Stellar_2 (nr__2) ! Stellarity index as defined by Sextractor (4) character*19 Name_2 (nr__2) ! Name (5) *Note (1): Errors computed by Sextractor. Listed errors do not include * photometric zero-point errors (see Sect. 3.7), those arising from the * compositing, and implicitly assume a gaussian or poissonian error * distribution. *Note (2): "99.99" means that no measure is available (usually because the * object is too small or too faint for its flux to be measured through a * large aperture). *Note (3): Kron magnitudes are computed by integrating the flux in an area * 2.5 times larger than the Kron area, with a minimum aperture radius of * 1.8 arcsec. *Note (4): When Index <~0.3, objects are identified as galaxies, whereas * Index >~0.7 (for 0404 and 0504) and Index >~0.8 (for patch1) * correspond to stars. The difference between fields is due to seeing * conditions. *Note (5): the acronyms are: * 5C: Wilson, 1970MNRAS.151....1W * CSO: Sanduleak & Pesch, 1990ApJS...72..291S * GMP: Godwin, Metcalfe & Peach (1983, Cat. <VII/42>) * NGC: Dreyer, See Cat. <VII/118> * NGP9: Odewahn & Aldering, 1995, Cat. <VII/204> * USNO: Monet et al., 1999, The USNO-A1 guide star catalog, Washington DC, * US Naval Observatory, see Cat. <I/243>, USNO-A1 * and Cat. <I/252>, USNO-A2.0 C============================================================================= C Loading file 'table3.dat' ! *List of detected sources in the region "Patch1" C Format for file interpretation 1 format( + 1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F5.2,1X,F4.2,1X, + F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2, + 1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,2X,I3,1X,F4.2,2X,I3,1X,F3.1, + 2X,F4.2,1X,A19) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,77 read(1,'(A143)')ar__ read(ar__,1) + RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__), + DEs(i__),Hmag5(i__),e_Hmag5(i__),Hmag10(i__),e_Hmag10(i__), + Hmag13(i__),e_Hmag13(i__),Hmag14_5(i__),e_Hmag14_5(i__), + Hmag18_5(i__),e_Hmag18_5(i__),Isomag(i__),e_Isomag(i__), + Kronmag(i__),e_Kronmag(i__),Area(i__),ell(i__),PA(i__), + Radius(i__),Stellar(i__),Name(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,1) + RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),DEm(i__), + DEs(i__),Hmag5(i__),e_Hmag5(i__),Hmag10(i__),e_Hmag10(i__), + Hmag13(i__),e_Hmag13(i__),Hmag14_5(i__),e_Hmag14_5(i__), + Hmag18_5(i__),e_Hmag18_5(i__),Isomag(i__),e_Isomag(i__), + Kronmag(i__),e_Kronmag(i__),Area(i__),ell(i__),PA(i__), + Radius(i__),Stellar(i__),Name(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 'table4.dat' ! *List of detected sources in the region "0404" C Format for file interpretation 2 format( + 1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F5.2,1X,F4.2,1X, + F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2, + 1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,2X,I3,1X,F4.2,2X,I3,1X,F3.1, + 2X,F4.2,1X,A19) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,124 read(1,'(A143)')ar__1 read(ar__1,2) + RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__), + DEm_1(i__),DEs_1(i__),Hmag5_1(i__),e_Hmag5_1(i__), + Hmag10_1(i__),e_Hmag10_1(i__),Hmag13_1(i__),e_Hmag13_1(i__), + Hmag14_5_1(i__),e_Hmag14_5_1(i__),Hmag18_5_1(i__), + e_Hmag18_5_1(i__),Isomag_1(i__),e_Isomag_1(i__), + Kronmag_1(i__),e_Kronmag_1(i__),Area_1(i__),ell_1(i__), + PA_1(i__),Radius_1(i__),Stellar_1(i__),Name_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,2) + RAh_1(i__),RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__), + DEm_1(i__),DEs_1(i__),Hmag5_1(i__),e_Hmag5_1(i__), + Hmag10_1(i__),e_Hmag10_1(i__),Hmag13_1(i__),e_Hmag13_1(i__), + Hmag14_5_1(i__),e_Hmag14_5_1(i__),Hmag18_5_1(i__), + e_Hmag18_5_1(i__),Isomag_1(i__),e_Isomag_1(i__), + Kronmag_1(i__),e_Kronmag_1(i__),Area_1(i__),ell_1(i__), + PA_1(i__),Radius_1(i__),Stellar_1(i__),Name_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 'table5.dat' ! *List of detected sources in the region "0504" C Format for file interpretation 3 format( + 1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F5.2,1X,F4.2,1X, + F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2, + 1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,2X,I3,1X,F4.2,2X,I3,1X,F3.1, + 2X,F4.2,1X,A19) C Effective file loading open(unit=1,status='old',file= +'table5.dat') write(6,*) '....Loading file: table5.dat' do i__=1,88 read(1,'(A143)')ar__2 read(ar__2,3) + RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__), + DEm_2(i__),DEs_2(i__),Hmag5_2(i__),e_Hmag5_2(i__), + Hmag10_2(i__),e_Hmag10_2(i__),Hmag13_2(i__),e_Hmag13_2(i__), + Hmag14_5_2(i__),e_Hmag14_5_2(i__),Hmag18_5_2(i__), + e_Hmag18_5_2(i__),Isomag_2(i__),e_Isomag_2(i__), + Kronmag_2(i__),e_Kronmag_2(i__),Area_2(i__),ell_2(i__), + PA_2(i__),Radius_2(i__),Stellar_2(i__),Name_2(i__) 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,3) + RAh_2(i__),RAm_2(i__),RAs_2(i__),DE__2(i__),DEd_2(i__), + DEm_2(i__),DEs_2(i__),Hmag5_2(i__),e_Hmag5_2(i__), + Hmag10_2(i__),e_Hmag10_2(i__),Hmag13_2(i__),e_Hmag13_2(i__), + Hmag14_5_2(i__),e_Hmag14_5_2(i__),Hmag18_5_2(i__), + e_Hmag18_5_2(i__),Isomag_2(i__),e_Isomag_2(i__), + Kronmag_2(i__),e_Kronmag_2(i__),Area_2(i__),ell_2(i__), + PA_2(i__),Radius_2(i__),Stellar_2(i__),Name_2(i__) write(6,'(6H Pos: 2F8.4)') RAdeg_2(i__),DEdeg_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end