Conversion of standardized ReadMe file for
file /./ftp/cats/VII/184 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-19 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. VII/184 Survey of 108 E-S0 galaxies (Michard+, 1993-94) *================================================================================ *Survey of 108 E-S0 RSA galaxies * Michard R., Marchal J. * <Astron. Astrophys. Suppl. Ser. 98, 29 (1993)> * <Astron. Astrophys. 283, 779 (1994)> * <Astron. Astrophys. Suppl. Ser. 105, 481 (1994)> * =1993A&AS...98...29M * =1994A&A...283..779M * =1994A&AS..105..481M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsafram' ! Observed galaxies and characteristics of the frame with available data integer*4 nr__ parameter (nr__=213) ! Number of records character*87 ar__ ! Full-size record integer*4 NGC (nr__) ! NGC number character*1 Color (nr__) ! Color of frame character*5 Code (nr__) ! Code number of frame character*8 Type (nr__) ! Hubble type from the Revised Shapley-Ames * Catalogue (RSA) character*3 Inst_ (nr__) ! Instrumentation (1) real*4 FWHM (nr__) ! (arcsec) FWHM of Point Spread Function in arc sec real*4 SkyBr (nr__) ! (mag/arcsec2) Sky brightness integer*4 Maj (nr__) ! (arcsec) Major axis of external fully analysed * isophote real*4 IsoBr (nr__) ! (mag/arcsec2) Brightness of the same isophote, character*30 Notes (nr__) ! Individual notes (2) *Note (1): hpc 120 cm in Haute-Provence with CCD * hpp id. with IIaO plates * t2r 2 m at Pic du Midi, f/10 or f/8 reducer with CCD * t2f id. , f/2.5 or f/3.5 reducer with CCD *Note (2): (1) Frames from the Toulouse team, kindly provided by P. Poulain. * (2) Photometric calibration from Michard and Simien, 1988. * (3) Photometric calibration from Michard, 1985. * (4) Galaxian center saturated. * (5) Frames kindly obtained for us by P. Prugniel. * (6) Frames kindly obtained for us by F. Simien. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsamorp' ! Main results of the morphological analysis of the studied galaxies. First list. (table2 of 1994A&AS..105..481M) integer*4 nr__1 parameter (nr__1=108) ! Number of records character*88 ar__1 ! Full-size record integer*4 NGC_1 (nr__1) ! NGC Number character*8 RSA (nr__1) ! Revised Shapley-Ames Catalogue Hubble type integer*4 T (nr__1) ! RC2 index (de Vaucouleurs Revised * Morphological Type) character*2 Ph (nr__1) ! Photometric evidence for bulk-disk * st = strong, ft = faint, * cl = clear, no = not seen real*4 qmin (nr__1) ! Minimum axis ratio real*4 log_aq (nr__1) ! ([arcsec]) Major axis of the isophote where qmin was * taken real*4 e_4 (nr__1) ! Carter coefficient of e4 real*4 e_6 (nr__1) ! Carter coefficient of e6 measured at the same * place as e*4 real*4 q25 (nr__1) ! Value of the axis ratio for the isophote * mu=25 in V light integer*4 DeltaPA (nr__1) ! Amplitude of the twist of the isophotes character*3 Bar (nr__1) ! Presence of a bar (bar or no) character*3 Dust (nr__1) ! Presence of dust * duR = dust ring, duL = dust lane * du- if of uncertain geometry, * ??? if uncertainty, no = no detection character*4 Disk (nr__1) ! Nature of disk * emDi = embedded disk, exDi = extended disk * miDi = mixed disk, no = no detection * ?Di if of dubful extension character*3 R_L (nr__1) ! Ring or lens presence * rl- = inner ring, -oL = outer lens, * roL = both seen, no = no detection character*3 Spi (nr__1) ! Spiral pattern presence * spP = possible spiral pattern, * no = no detection character*3 Env (nr__1) ! Envelope classification * spH = spheroidal halo, thD = thick disk * exD = extended disk, pec = peculiar envelope * ? if not classified character*4 Asy (nr__1) ! Miscellaneous structure * asyC = asymmetries with color effects * asyG = asymmetries without color effects * asy? = color class unknown * no = no significant detected structure * ? no suitable observations character*6 type_1 (nr__1) ! Revised Hubble type, completed for E objects. * diE = disky E, boE = boxyE, * unE if of undeterminate subclass * p added if the object has a peculiar envelope * SAO, SAB0, SB0 as in the RC2, except that the * r, s varieties are not noted c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsaphot' ! Main results of the morphological analysis of the studied galaxies. Second list. (table3 of 1994A&AS..105..481M) integer*4 nr__2 parameter (nr__2=108) ! Number of records character*75 ar__2 ! Full-size record integer*4 NGC_2 (nr__2) ! NGC number integer*4 T_1 (nr__2) ! RC2 index (de Vaucouleurs Revised * Morphological Type) character*5 Type_2 (nr__2) ! Revised and complete morphological * classification of rsamorp real*4 VT (nr__2) ! (mag) Asymptotic magnitude either measured from * our V frames, or sometimes converted form B * frames through the RC2 B-V colors real*4 log_r_e (nr__2) ! ([arcsec]) de Vaucouleurs radius real*4 log_l_e (nr__2) ! ([arcsec]) Corresponding semi-major axis real*4 mu_e_V (nr__2) ! (mag/arcmin2) Mean surface brightness of the effective * isophote in V light real*4 mu0 (nr__2) ! (mag) Adopted distance modulus real*4 BT0 (nr__2) ! (mag) Asymptotic magnitude converted to B light real*4 MT_B (nr__2) ! (mag) Absolute magnitude real*4 log_R_e_1 (nr__2) ! ([kpc]) Mean radius of the effective isophote real*4 mu_e_B (nr__2) ! (mag/arcmin2) Mean surface brightness of the effective * isophote in B light real*4 logD25 (nr__2) ! ([0.1arcmin]) RC2-like parameter real*4 logR25 (nr__2) ! RC2-like parameter c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsamuls1' ! Initial parameters for photometric measurements integer*4 nr__3 parameter (nr__3=214) ! Number of records character*111 ar__3 ! Full-size record integer*4 NGC_3 (nr__3) ! NGC number character*1 Color_1 (nr__3) ! Color of the observation character*6 Frame (nr__3) ! Serial number for the frame character*6 OtherName (nr__3) ! Other name real*4 size (nr__3) ! (arcsec) pixel size real*8 FWHM_1 (nr__3) ! (pix) FWHM of the PSF integer*4 Nharm (nr__3) ! number of harmonics in Carter's isophotal * representation (always 10) integer*4 Niso (nr__3) ! number of isophotes in the table real*8 p1 (nr__3) ! Parameter used to describe the sky level (1) real*4 p2 (nr__3) ! Parameter used to describe the sky level (1) real*8 cal (nr__3) ! (mag) calibration constant real*8 pSuBr (nr__3) ! (mag) peak surface brightness of the galaxy real*8 SuBr (nr__3) ! (mag) surface brightness of the faintest measured * isophote *Note (1): Those parameters are in frame units and of "house keeping" interest c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsamuls2' ! Average magnitudes on axis integer*4 nr__4 parameter (nr__4=15479) ! Number of records character*46 ar__4 ! Full-size record integer*4 NGC_4 (nr__4) ! NGC number character*1 Color_2 (nr__4) ! Color of the observation character*6 Frame_1 (nr__4) ! Serial number for the frame real*4 a (nr__4) ! (arcsec) major axis of contour real*4 SuBra (nr__4) ! (mag) surface brightness corresponding to major axis * of contour real*4 b (nr__4) ! (arcsec) minor axis of contour real*4 SuBrb (nr__4) ! (mag) surface brightness corresponding to minor axis * of contour c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsamuls3' ! Differences between sectors integer*4 nr__5 parameter (nr__5=15396) ! Number of records character*49 ar__5 ! Full-size record integer*4 NGC_5 (nr__5) ! NGC number character*1 Color_3 (nr__5) ! Color of the observation character*6 Frame_2 (nr__5) ! Serial number for the frame real*4 Dmaj (nr__5) ! (arcsec) Distance from center along major axis real*4 Dsubrmaj (nr__5) ! (mag) Difference of surface brightness at opposite * major axis positions of the given central * distance real*4 Dmin (nr__5) ! (arcsec) Distance from center along minor axis real*4 Dsubrmin (nr__5) ! (mag) Difference of surface brightness at opposite * minor axis positions of the given central * distance c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsamuls4' ! Gradients on the two axis integer*4 nr__6 parameter (nr__6=14967) ! Number of records character*50 ar__6 ! Full-size record integer*4 NGC_6 (nr__6) ! NGC number character*1 Color_4 (nr__6) ! Color of the observation character*6 Frame_3 (nr__6) ! Serial number for the frame real*4 a_1 (nr__6) ! (arcsec) major axis of contour real*8 GradSuBra (nr__6) ! (mag/arcsec) Gradient of surface brightness along major * axis real*4 b_1 (nr__6) ! (arcsec) minor axis of contour real*8 GradSuBrb (nr__6) ! (mag/arcsec) Gradient of surface brightness along minor * axis c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsaresu1' ! Initial parameters for isophotal analysis integer*4 nr__7 parameter (nr__7=214) ! Number of records character*103 ar__7 ! Full-size record integer*4 NGC_7 (nr__7) ! NGC number character*1 Color_5 (nr__7) ! Color of the observation character*6 Frame_4 (nr__7) ! Serial number for the frame real*4 size_1 (nr__7) ! (arcsec) Pixel size real*8 FWHM_2 (nr__7) ! (pix) FWHM of the PSF integer*4 Nharm_1 (nr__7) ! Number of harmonics in Carter's isophotal * representation (always 10) integer*4 Niso_1 (nr__7) ! Number of isophotes in the table real*8 p1_1 (nr__7) ! Parameter used to describe the sky level (1) integer*4 p2_1 (nr__7) ! Parameter used to describe the sky level (1) real*4 p3 (nr__7) ! Parameter used to describe the sky level (1) real*8 Cal_1 (nr__7) ! (mag) Calibration constant real*8 Pint (nr__7) ! Peak intensity of the galaxy (2) real*4 isof (nr__7) ! Level of faintest measures isophote (2) real*8 PSuBr_1 (nr__7) ! (mag) Peak surface brightness of the galaxy real*8 SuBrf (nr__7) ! (mag) []? Surface brightness of the faintest * measured isophote *Note (1): Those parameters are in frame units and of "house keeping" interest *Note (2): Those parameters are in frame units c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rsaresu2' ! Isophote parameters integer*4 nr__8 parameter (nr__8=15182) ! Number of records character*93 ar__8 ! Full-size record integer*4 NGC_8 (nr__8) ! NGC number character*1 Color_6 (nr__8) ! Color of the observation character*6 Frame_5 (nr__8) ! Serial number for the frame integer*4 Npix (nr__8) ! Number of pixels used to describe the isophote real*4 MaRE (nr__8) ! (arcsec) Major axis of the Reference Ellipse (RE) integer*4 rRE (nr__8) ! (10-3) Axis ratio of the RE real*4 PArRE (nr__8) ! (deg) Position angle of the major axis of the RE integer*4 C4 (nr__8) ! (10-4) Coefficient of the 4th cosine harmonic integer*4 C6 (nr__8) ! (10-4) Coefficient of the 6th cosine harmonic integer*4 C8 (nr__8) ! (10-4) Coefficient of the 8th cosine harmonic integer*4 C10 (nr__8) ! (10-4) Coefficient of the 10th cosine harmonic integer*4 C3 (nr__8) ! (10-4) Coefficient of the 3rd cosine harmonic integer*4 S3 (nr__8) ! (10-4) Coefficient of the 3rd sine harmonic integer*4 S2 (nr__8) ! (10-4) Coefficient of the 2nd sine harmonic. * This coefficient is expected to be zero if * the isophote is reasonnably well defined. * Significant values occur at low surface * brightness and indicate a poor isophote * representation. integer*4 S4 (nr__8) ! (10-4) Coefficient of the 4th sine harmonic integer*4 S6 (nr__8) ! (10-4) Coefficient of the 6th sine harmonic real*4 SuBr_1 (nr__8) ! (mag/arcsec2) Surface brightness of the isophote C============================================================================= C Loading file 'rsafram' ! Observed galaxies and characteristics of the * frame with available data C Format for file interpretation 1 format( + 4X,I4,2X,A1,2X,A5,1X,A8,3X,A3,2X,F4.2,1X,F5.2,1X,I3,1X,F5.2, + 2X,A30) C Effective file loading open(unit=1,status='old',file= +'rsafram') write(6,*) '....Loading file: rsafram' do i__=1,213 read(1,'(A87)')ar__ read(ar__,1) + NGC(i__),Color(i__),Code(i__),Type(i__),Inst_(i__),FWHM(i__), + SkyBr(i__),Maj(i__),IsoBr(i__),Notes(i__) c ..............Just test output........... write(6,1) + NGC(i__),Color(i__),Code(i__),Type(i__),Inst_(i__),FWHM(i__), + SkyBr(i__),Maj(i__),IsoBr(i__),Notes(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsamorp' ! Main results of the morphological analysis of * the studied galaxies. First list. * (table2 of 1994A&AS..105..481M) C Format for file interpretation 2 format( + 1X,I4,1X,A8,1X,I2,1X,A2,1X,F4.2,1X,F4.2,1X,F4.1,1X,F4.1,1X, + F4.2,1X,I2,2X,A3,1X,A3,1X,A4,1X,A3,1X,A3,2X,A3,1X,A4,2X,A6) C Effective file loading open(unit=1,status='old',file= +'rsamorp') write(6,*) '....Loading file: rsamorp' do i__=1,108 read(1,'(A88)')ar__1 read(ar__1,2) + NGC_1(i__),RSA(i__),T(i__),Ph(i__),qmin(i__),log_aq(i__), + e_4(i__),e_6(i__),q25(i__),DeltaPA(i__),Bar(i__),Dust(i__), + Disk(i__),R_L(i__),Spi(i__),Env(i__),Asy(i__),type_1(i__) c ..............Just test output........... write(6,2) + NGC_1(i__),RSA(i__),T(i__),Ph(i__),qmin(i__),log_aq(i__), + e_4(i__),e_6(i__),q25(i__),DeltaPA(i__),Bar(i__),Dust(i__), + Disk(i__),R_L(i__),Spi(i__),Env(i__),Asy(i__),type_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsaphot' ! Main results of the morphological analysis of * the studied galaxies. Second list. * (table3 of 1994A&AS..105..481M) C Format for file interpretation 3 format( + 1X,I4,2X,I2,1X,A5,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X, + F5.2,F6.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F3.2) C Effective file loading open(unit=1,status='old',file= +'rsaphot') write(6,*) '....Loading file: rsaphot' do i__=1,108 read(1,'(A75)')ar__2 read(ar__2,3) + NGC_2(i__),T_1(i__),Type_2(i__),VT(i__),log_r_e(i__), + log_l_e(i__),mu_e_V(i__),mu0(i__),BT0(i__),MT_B(i__), + log_R_e_1(i__),mu_e_B(i__),logD25(i__),logR25(i__) c ..............Just test output........... write(6,3) + NGC_2(i__),T_1(i__),Type_2(i__),VT(i__),log_r_e(i__), + log_l_e(i__),mu_e_V(i__),mu0(i__),BT0(i__),MT_B(i__), + log_R_e_1(i__),mu_e_B(i__),logD25(i__),logR25(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsamuls1' ! Initial parameters for photometric measurements C Format for file interpretation 4 format( + 1X,I4,A1,1X,A6,4X,A6,7X,F6.4,5X,F8.5,1X,I2,1X,I2,1X,F7.2,3X, + E13.3,1X,F7.4,1X,F7.4,4X,F8.5) C Effective file loading open(unit=1,status='old',file= +'rsamuls1') write(6,*) '....Loading file: rsamuls1' do i__=1,214 read(1,'(A111)')ar__3 read(ar__3,4) + NGC_3(i__),Color_1(i__),Frame(i__),OtherName(i__),size(i__), + FWHM_1(i__),Nharm(i__),Niso(i__),p1(i__),p2(i__),cal(i__), + pSuBr(i__),SuBr(i__) c ..............Just test output........... write(6,4) + NGC_3(i__),Color_1(i__),Frame(i__),OtherName(i__),size(i__), + FWHM_1(i__),Nharm(i__),Niso(i__),p1(i__),p2(i__),cal(i__), + pSuBr(i__),SuBr(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsamuls2' ! Average magnitudes on axis C Format for file interpretation 5 format(1X,I4,A1,1X,A6,1X,F6.2,2X,F6.3,3X,F6.2,2X,F6.3) C Effective file loading open(unit=1,status='old',file= +'rsamuls2') write(6,*) '....Loading file: rsamuls2' do i__=1,15479 read(1,'(A46)')ar__4 read(ar__4,5) + NGC_4(i__),Color_2(i__),Frame_1(i__),a(i__),SuBra(i__),b(i__), + SuBrb(i__) c ..............Just test output........... write(6,5) + NGC_4(i__),Color_2(i__),Frame_1(i__),a(i__),SuBra(i__),b(i__), + SuBrb(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsamuls3' ! Differences between sectors C Format for file interpretation 6 format(1X,I4,A1,1X,A6,5X,F6.2,2X,F6.3,3X,F6.2,2X,F6.3) C Effective file loading open(unit=1,status='old',file= +'rsamuls3') write(6,*) '....Loading file: rsamuls3' do i__=1,15396 read(1,'(A49)')ar__5 read(ar__5,6) + NGC_5(i__),Color_3(i__),Frame_2(i__),Dmaj(i__),Dsubrmaj(i__), + Dmin(i__),Dsubrmin(i__) c ..............Just test output........... write(6,6) + NGC_5(i__),Color_3(i__),Frame_2(i__),Dmaj(i__),Dsubrmaj(i__), + Dmin(i__),Dsubrmin(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsamuls4' ! Gradients on the two axis C Format for file interpretation 7 format(1X,I4,A1,1X,A6,7X,F6.2,1X,F7.4,2X,F6.2,1X,F7.4) C Effective file loading open(unit=1,status='old',file= +'rsamuls4') write(6,*) '....Loading file: rsamuls4' do i__=1,14967 read(1,'(A50)')ar__6 read(ar__6,7) + NGC_6(i__),Color_4(i__),Frame_3(i__),a_1(i__),GradSuBra(i__), + b_1(i__),GradSuBrb(i__) c ..............Just test output........... write(6,7) + NGC_6(i__),Color_4(i__),Frame_3(i__),a_1(i__),GradSuBra(i__), + b_1(i__),GradSuBrb(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsaresu1' ! Initial parameters for isophotal analysis C Format for file interpretation 8 format( + 1X,I4,A1,1X,A6,2X,F6.4,1X,F8.5,1X,I2,1X,I2,1X,F7.2,1X,I2,1X, + E13.5,1X,F7.4,1X,F8.2,1X,F5.2,1X,F7.4,1X,F10.7) C Effective file loading open(unit=1,status='old',file= +'rsaresu1') write(6,*) '....Loading file: rsaresu1' do i__=1,214 read(1,'(A103)')ar__7 read(ar__7,8) + NGC_7(i__),Color_5(i__),Frame_4(i__),size_1(i__),FWHM_2(i__), + Nharm_1(i__),Niso_1(i__),p1_1(i__),p2_1(i__),p3(i__), + Cal_1(i__),Pint(i__),isof(i__),PSuBr_1(i__),SuBrf(i__) if(ar__7(94:103) .EQ. '') SuBrf(i__) = rNULL__ c ..............Just test output........... write(6,8) + NGC_7(i__),Color_5(i__),Frame_4(i__),size_1(i__),FWHM_2(i__), + Nharm_1(i__),Niso_1(i__),p1_1(i__),p2_1(i__),p3(i__), + Cal_1(i__),Pint(i__),isof(i__),PSuBr_1(i__),SuBrf(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rsaresu2' ! Isophote parameters C Format for file interpretation 9 format( + 1X,I4,A1,1X,A6,4X,I4,1X,F6.2,1X,I3,2X,F6.2,1X,I5,1X,I4,1X,I4, + 1X,I4,1X,I4,1X,I4,1X,I4,1X,I4,1X,I4,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'rsaresu2') write(6,*) '....Loading file: rsaresu2' do i__=1,15182 read(1,'(A93)')ar__8 read(ar__8,9) + NGC_8(i__),Color_6(i__),Frame_5(i__),Npix(i__),MaRE(i__), + rRE(i__),PArRE(i__),C4(i__),C6(i__),C8(i__),C10(i__),C3(i__), + S3(i__),S2(i__),S4(i__),S6(i__),SuBr_1(i__) c ..............Just test output........... write(6,9) + NGC_8(i__),Color_6(i__),Frame_5(i__),Npix(i__),MaRE(i__), + rRE(i__),PArRE(i__),C4(i__),C6(i__),C8(i__),C10(i__),C3(i__), + S3(i__),S2(i__),S4(i__),S6(i__),SuBr_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end