Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_AS/116/429 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 2026-Aug-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/A+AS/116/429 Classification of Coma early galaxies (Andreon+ 1996)
*================================================================================
*Morphological classification and structural parameters for early-type galaxies
*in the Coma cluster.
* Andreon S., Davoust E., Michard R., Nieto J.-L., Poulain P.
* <Astron. Astrophys. Suppl. Ser. 116, 429 (1996)>
* =1996A&AS..116..429A (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1' ! Observing log
integer*4 nr__
parameter (nr__=151) ! Number of records
character*43 ar__ ! Full-size record
character*5 Name ! Abridged GMP or NGC number
integer*4 Exp ! (s) ? Exposure time
character*1 Filt ! Filter
integer*4 Run ! Observing run
real*4 FWHM ! (arcsec) ? Seeing
character*1 u_FWHM ! [?] Uncertainty flag on FWHM
real*4 Sky ! (mag/arcsec+2) ? Sky brightness
real*4 Isolim ! (mag/arcsec+2) ? Level of last measured isophote
character*4 Notes ! Notes (1)
*Note (1): (1) presumably absorbed weather;
* (2) not analysed;
* (3) spiral;
* (4) calibration from mean sky brightness;
* (5) elongated stellar images;
* (6) calibration uncertain, mean of previous night;
* (7) near edge of image;
* (8) very doubtful identity;
* (9) blended binary;
* (10) inadequate field;
* (11) Superimposed on NGC 4889, measured after subtraction of a model
* image of the cD;
* (12) underexposed image with bright sky;
* (13) companion to the W;
* (14) halo of a very bright star;
* (15) suspiciously high seeing;
* (16) very poor seeing;
* (17) calibrator of KP2225.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2a' ! CCD Telescope Bernard Lyot (TBL) subsample
photometric parameters and classification
integer*4 nr__1
parameter (nr__1=105) ! Number of records
character*133 ar__1 ! Full-size record
integer*4 GMPa ! Number in the abridged version of GMP
* catalogue <VII/42>
integer*4 GMP ! Number in the unabridged version of GMP
* catalogue <VII/42>
integer*4 Dress ! ? Number in Dressler catalogue <VII/174>
character*7 Name_1 ! NGC, Ic or RB (Rood & Baum 1967) number.
character*1 m_Name ! Multiplicity index on Name
integer*4 Xpos ! (arcsec) X differential position (1)
integer*4 Ypos ! (arcsec) Y differential position (1)
real*4 Rmag ! (mag) ? Asymptotic magnitude in R_c_
real*4 log_Reff ! ([arcsec]) ? Effective radius (10)
real*4 log_a ! ([arcsec]) ? Semi-major axis of the effective
* isophote (10)
real*4 SuBr ! (mag/arcmin+2) ? Mean surface brightness inside the
* effective isophote, in R magnitude
character*2 Diskp ! Photometric evidence for a disk (2)
real*4 b_a ! ? Typical axis ratio (9)
real*4 e4 ! (%) ? Typical e4 parameter (9)
character*2 n_e4 ! Location where the e4 parameter was
* estimated (3)
real*4 b_ae ! ? Axis ratio at the envelope, i.e. at the
* isophote {mu}_r_=24mag/arcsec^2^
integer*4 Ampli ! (deg) ? Amplitude of isophotal twist in the range
* of reliable measurements
character*4 Bar ! Detection of bar (4)
character*4 Disk ! Detection and classification of a disk (5)
character*5 Spiral ! Detection of a spiral pattern (6)
character*4 Envelope ! Classification of an envelope (7)
character*8 MType ! Morphological classification (8)
character*1 Note ! [*] Individual notes (see notes file)
*Note (1): X and Y positions in arc seconds relative to the approximate
* center of the cluster, increasing westwards and northwards
* respectively. The origin of these coordinates is at
* RA=12h57.3 Dec=+28d14.4'(1950.0).
*Note (2): st: strong;
* cl: clear;
* ft: faint;
* no: none
*Note (3): ex: at its extremum,
* re: at the effective isophote,
* co: if the value is the same at both locations
*Note (4): bar: bar seen,
* bar?: bar suspected,
* no: no bar seen
*Note (5): emDi: embedded disk,
* miDi: mixed disk,
* exDi: extended disk,
* -?Di: detected but unclassified disk,
* -no- (no disk seen)
*Note (6): spiP: spiral pattern seen,
* spiP?: spiral pattern suspected,
* -no-: no spiral pattern seen
*Note (7): spH: spheroidal halo,
* thD: thick disk,
* exD: extended disk,
* pec: peculiar envelope,
* -?-: unclassified envelope
*Note (8): boE: boxy E,
* unE: undetermined E,
* diE: disky E,
* SA0, SAB0, SB0, Sa, etc.,
* S...: spiral of unknown stage)
*Note (9): Typical value, either its minimum value, if clearly defined,
* or its value at the effective isophote otherwise
*Note (10): The units are erroneously indicated as 0.1arcmin in the
* printed paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2b' ! CCD CFH subsample photometric parameters and
classification
integer*4 nr__2
parameter (nr__2=17) ! Number of records
character*133 ar__2 ! Full-size record
integer*4 GMPa_1 ! Number in the abridged version of GMP
* catalogue <VII/42>
integer*4 GMP_1 ! Number in the unabridged version of GMP
* catalogue <VII/42>
integer*4 Dress_1 ! ? Number in Dressler catalogue <VII/174>
character*7 Name_2 ! NGC, Ic or RB (Rood & Baum 1967) number.
character*1 m_Name_1 ! Multiplicity index on Name
integer*4 Xpos_1 ! (arcsec) X differential position (1)
integer*4 Ypos_1 ! (arcsec) Y differential position (1)
real*4 Rmag_1 ! (mag) ? Asymptotic magnitude in R_c_
real*4 log_Reff_1 ! ([arcsec]) ? Effective radius (10)
real*4 log_a_1 ! ([arcsec]) ? Semi-major axis of the effective
* isophote (10)
real*4 SuBr_1 ! (mag/arcmin+2) ? Mean surface brightness inside the
* effective isophote, in R magnitude
character*2 Diskp_1 ! Photometric evidence for a disk (2)
real*4 b_a_1 ! ? Typical axis ratio (9)
real*4 e4_1 ! (%) ? Typical e4 parameter (9)
character*2 n_e4_1 ! Location where the e4 parameter was
* estimated (3)
real*4 b_ae_1 ! ? Axis ratio at the envelope, i.e. at the
* isophote {mu}_r_=24mag/arcsec^2^
integer*4 Ampli_1 ! (deg) ? Amplitude of isophotal twist in the range
* of reliable measurements
character*4 Bar_1 ! Detection of bar (4)
character*4 Disk_1 ! Detection and classification of a disk (5)
character*5 Spiral_1 ! Detection of a spiral pattern (6)
character*4 Envelope_1 ! Classification of an envelope (7)
character*8 MType_1 ! Morphological classification (8)
character*1 Note_1 ! [*] Individual notes (see notes file)
*Note (1): X and Y positions in arc seconds relative to the approximate
* center of the cluster, increasing westwards and northwards
* respectively. The origin of these coordinates is at
* RA=12h57.3 Dec=+28d14.4'(1950.0).
*Note (2): st: strong;
* cl: clear;
* ft: faint;
* no: none
*Note (3): ex: at its extremum,
* re: at the effective isophote,
* co: if the value is the same at both locations
*Note (4): bar: bar seen,
* bar?: bar suspected,
* no: no bar seen
*Note (5): emDi: embedded disk,
* miDi: mixed disk,
* exDi: extended disk,
* -?Di: detected but unclassified disk,
* -no- (no disk seen)
*Note (6): spiP: spiral pattern seen,
* spiP?: spiral pattern suspected,
* -no-: no spiral pattern seen
*Note (7): spH: spheroidal halo,
* thD: thick disk,
* exD: extended disk,
* pec: peculiar envelope,
* -?-: unclassified envelope
*Note (8): boE: boxy E,
* unE: undetermined E,
* diE: disky E,
* SA0, SAB0, SB0, Sa, etc.,
* S...: spiral of unknown stage)
*Note (9): Typical value, either its minimum value, if clearly defined,
* or its value at the effective isophote otherwise
*Note (10): The units are erroneously indicated as 0.1arcmin in the
* printed paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3' ! KP1608 sample photometric parameters and
classification
integer*4 nr__3
parameter (nr__3=113) ! Number of records
character*133 ar__3 ! Full-size record
integer*4 GMPa_2 ! Number in the abridged version of GMP
* catalogue <VII/42>
integer*4 GMP_2 ! Number in the unabridged version of GMP
* catalogue <VII/42>
integer*4 Dress_2 ! ? Number in Dressler catalogue <VII/174>
character*7 Name_3 ! NGC, Ic or RB (Rood & Baum 1967) number.
character*1 m_Name_2 ! Multiplicity index on Name
integer*4 Xpos_2 ! (arcsec) X differential position (1)
integer*4 Ypos_2 ! (arcsec) Y differential position (1)
real*4 Rmag_2 ! (mag) ? Asymptotic magnitude in R_c_
real*4 log_Reff_1 ! ([arcsec]) ? Effective radius (10)
real*4 log_a_1 ! ([arcsec]) ? Semi-major axis of the effective
* isophote (10)
real*4 SuBr_2 ! (mag/arcmin+2) ? Mean surface brightness inside the
* effective isophote, in R magnitude
character*2 Diskp_2 ! Photometric evidence for a disk (2)
real*4 b_a_2 ! ? Typical axis ratio (9)
real*4 e4_2 ! (%) ? Typical e4 parameter (9)
character*2 n_e4_2 ! Location where the e4 parameter was
* estimated (3)
real*4 b_ae_2 ! ? Axis ratio at the envelope, i.e. at the
* isophote {mu}_r_=24mag/arcsec^2^
integer*4 Ampli_2 ! (deg) ? Amplitude of isophotal twist in the range
* of reliable measurements
character*4 Bar_2 ! Detection of bar (4)
character*4 Disk_2 ! Detection and classification of a disk (5)
character*5 Spiral_2 ! Detection of a spiral pattern (6)
character*4 Envelope_2 ! Classification of an envelope (7)
character*8 MType_2 ! Morphological classification (8)
character*1 Note_2 ! [*] Individual notes (see notes file)
*Note (1): X and Y positions in arc seconds relative to the approximate
* center of the cluster, increasing westwards and northwards
* respectively. The origin of these coordinates is at
* RA=12h57.3 Dec=+28d14.4'(1950.0).
*Note (2): st: strong;
* cl: clear;
* ft: faint;
* no: none
*Note (3): ex: at its extremum,
* re: at the effective isophote,
* co: if the value is the same at both locations
*Note (4): bar: bar seen,
* bar?: bar suspected,
* no: no bar seen
*Note (5): emDi: embedded disk,
* miDi: mixed disk,
* exDi: extended disk,
* -?Di: detected but unclassified disk,
* -no- (no disk seen)
*Note (6): spiP: spiral pattern seen,
* spiP?: spiral pattern suspected,
* -no-: no spiral pattern seen
*Note (7): spH: spheroidal halo,
* thD: thick disk,
* exD: extended disk,
* pec: peculiar envelope,
* -?-: unclassified envelope
*Note (8): boE: boxy E,
* unE: undetermined E,
* diE: disky E,
* SA0, SAB0, SB0, Sa, etc.,
* S...: spiral of unknown stage)
*Note (9): Typical value, either its minimum value, if clearly defined,
* or its value at the effective isophote otherwise
*Note (10): The units are erroneously indicated as 0.1arcmin in the
* printed paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4' ! FW sample photometric parameters and
classification
integer*4 nr__4
parameter (nr__4=47) ! Number of records
character*121 ar__4 ! Full-size record
C 1950.0 position composed of: RAh RAm RAs DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension 1950.0
real*8 DEdeg ! (deg) Declination 1950.0
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 RAh ! (h) Right ascension (1950.0)
integer*4 RAm ! (min) Right ascension (1950.0)
integer*4 RAs ! (s) Right ascension (1950.0)
integer*4 DEd ! (deg) Declination (1950.0)
integer*4 DEm ! (arcmin) Declination (1950.0)
integer*4 DEs ! (arcsec) Declination (1950.0)
character*7 Name_4 ! NGC, IC or RB (Rood & Baum 1967) number.
* FW and D from Strom & Strom (1978),
* anon = anonymous (m_R_<=16.05),
* supll = anonymous (m_R_> 16.05)
character*1 m_Name_3 ! Multiplicity index on Name
real*4 Rmag_3 ! (mag) ? Asymptotic magnitude in R_c_
real*4 log_Reff_1 ! ([arcsec]) ? Effective radius (10)
real*4 log_a_1 ! ([arcsec]) ? Semi-major axis of the effective
* isophote (10)
real*4 SuBr_3 ! (mag/arcmin+2) ? Mean surface brightness inside the
* effective isophote, in R magnitude
character*2 Diskp_3 ! Photometric evidence for a disk (1)
real*4 b_a_3 ! ? Typical axis ratio (8)
integer*4 e4_3 ! (%) ? Typical e4 parameter (8)
character*2 n_e4_3 ! Location where the e4 parameter was
* estimated (2)
real*4 b_ae_3 ! ? Axis ratio at the envelope, i.e. at the
* isophote {mu}_r_=24mag/arcsec^2^
integer*4 Ampli_3 ! (deg) ? Amplitude of isophotal twist in the
* range of reliable measurements
character*4 Bar_3 ! Detection of bar (3)
character*4 Disk_3 ! Detection and classification of a disk (4)
character*5 Spiral_3 ! Detection of a spiral pattern (5)
character*4 Envelope_3 ! Classification of an envelope (6)
character*8 MType_3 ! Morphological classification (7)
character*1 Note_3 ! [*] Individual notes (see notes file)
*Note (1): st (strong); cl (clear); ft (faint); no (none)
*Note (2): ex (at its extremum), re (at the effective isophote),
* co (if the value is the same at both locations)
*Note (3): bar (bar seen), bar? (bar suspected), no (no bar seen)
*Note (4): emDi (embedded disk), miDi (mixed disk), exDi (extended disk),
* -?Di (detected but unclassified disk), -no- (no disk seen)
*Note (5): spiP (spiral pattern seen), spiP? (spiral pattern suspected),
* -no- (no spiral pattern seen)
*Note (6): spH (spheroidal halo), thD (thick disk), exD (extended disk),
* pec (peculiar envelope), -?- (unclassified envelope)
*Note (7): boE (boxy E), unE (undetermined E), diE (disky E), SA0, SAB0, SB0,
* Sa, etc., S... (spiral of unknown stage)
*Note (8): Typical value, either its minimum value, if clearly defined, or its
* value at the effective isophote otherwise
*Note (10): The units are erroneously indicated as 0.1arcmin in the
* printed paper.
C=============================================================================
C Loading file 'table1' ! Observing log
C Format for file interpretation
1 format(A5,2X,I2,2X,A1,1X,I1,2X,F4.2,A1,2X,F5.2,2X,F5.2,1X,A4)
C Effective file loading
open(unit=1,status='old',file=
+'table1')
write(6,*) '....Loading file: table1'
do i__=1,151
read(1,'(A43)')ar__
read(ar__,1)Name,Exp,Filt,Run,FWHM,u_FWHM,Sky,Isolim,Notes
if(ar__(8:9) .EQ. '') Exp = iNULL__
if(ar__(17:20) .EQ. '') FWHM = rNULL__
if(ar__(24:28) .EQ. '') Sky = rNULL__
if(ar__(31:35) .EQ. '') Isolim = rNULL__
c ..............Just test output...........
write(6,1)Name,Exp,Filt,Run,FWHM,u_FWHM,Sky,Isolim,Notes
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2a' ! CCD Telescope Bernard Lyot (TBL) subsample
* photometric parameters and classification
C Format for file interpretation
2 format(
+ I4,2X,I4,2X,I3,2X,A7,A1,2X,I5,2X,I5,2X,F5.2,2X,F4.2,2X,F4.2,
+ 2X,F5.2,2X,A2,2X,F4.2,2X,F4.1,2X,A2,2X,F4.2,2X,I2,2X,A4,2X,A4,
+ 2X,A5,2X,A4,2X,A8,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2a')
write(6,*) '....Loading file: table2a'
do i__=1,105
read(1,'(A133)')ar__1
read(ar__1,2)
+ GMPa,GMP,Dress,Name_1,m_Name,Xpos,Ypos,Rmag,log_Reff,log_a,
+ SuBr,Diskp,b_a,e4,n_e4,b_ae,Ampli,Bar,Disk,Spiral,Envelope,
+ MType,Note
if(ar__1(13:15) .EQ. '') Dress = iNULL__
if(ar__1(42:46) .EQ. '') Rmag = rNULL__
if(ar__1(49:52) .EQ. '') log_Reff = rNULL__
if(ar__1(55:58) .EQ. '') log_a = rNULL__
if(ar__1(61:65) .EQ. '') SuBr = rNULL__
if(ar__1(72:75) .EQ. '') b_a = rNULL__
if(ar__1(78:81) .EQ. '') e4 = rNULL__
if(ar__1(88:91) .EQ. '') b_ae = rNULL__
if(ar__1(94:95) .EQ. '') Ampli = iNULL__
c ..............Just test output...........
write(6,2)
+ GMPa,GMP,Dress,Name_1,m_Name,Xpos,Ypos,Rmag,log_Reff,log_a,
+ SuBr,Diskp,b_a,e4,n_e4,b_ae,Ampli,Bar,Disk,Spiral,Envelope,
+ MType,Note
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2b' ! CCD CFH subsample photometric parameters and
* classification
C Format for file interpretation
3 format(
+ I4,2X,I4,2X,I3,2X,A7,A1,2X,I5,2X,I5,2X,F5.2,2X,F4.2,2X,F4.2,
+ 2X,F5.2,2X,A2,2X,F4.2,2X,F4.1,2X,A2,2X,F4.2,2X,I2,2X,A4,2X,A4,
+ 2X,A5,2X,A4,2X,A8,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2b')
write(6,*) '....Loading file: table2b'
do i__=1,17
read(1,'(A133)')ar__2
read(ar__2,3)
+ GMPa_1,GMP_1,Dress_1,Name_2,m_Name_1,Xpos_1,Ypos_1,Rmag_1,
+ log_Reff_1,log_a_1,SuBr_1,Diskp_1,b_a_1,e4_1,n_e4_1,b_ae_1,
+ Ampli_1,Bar_1,Disk_1,Spiral_1,Envelope_1,MType_1,Note_1
if(ar__2(13:15) .EQ. '') Dress_1 = iNULL__
if(ar__2(42:46) .EQ. '') Rmag_1 = rNULL__
if(ar__2(49:52) .EQ. '') log_Reff_1 = rNULL__
if(ar__2(55:58) .EQ. '') log_a_1 = rNULL__
if(ar__2(61:65) .EQ. '') SuBr_1 = rNULL__
if(ar__2(72:75) .EQ. '') b_a_1 = rNULL__
if(ar__2(78:81) .EQ. '') e4_1 = rNULL__
if(ar__2(88:91) .EQ. '') b_ae_1 = rNULL__
if(ar__2(94:95) .EQ. '') Ampli_1 = iNULL__
c ..............Just test output...........
write(6,3)
+ GMPa_1,GMP_1,Dress_1,Name_2,m_Name_1,Xpos_1,Ypos_1,Rmag_1,
+ log_Reff_1,log_a_1,SuBr_1,Diskp_1,b_a_1,e4_1,n_e4_1,b_ae_1,
+ Ampli_1,Bar_1,Disk_1,Spiral_1,Envelope_1,MType_1,Note_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3' ! KP1608 sample photometric parameters and
* classification
C Format for file interpretation
4 format(
+ I4,2X,I4,2X,I3,2X,A7,A1,2X,I5,2X,I5,2X,F5.2,2X,F4.2,2X,F4.2,
+ 2X,F5.2,2X,A2,2X,F4.2,2X,F4.1,2X,A2,2X,F4.2,2X,I2,2X,A4,2X,A4,
+ 2X,A5,2X,A4,2X,A8,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3')
write(6,*) '....Loading file: table3'
do i__=1,113
read(1,'(A133)')ar__3
read(ar__3,4)
+ GMPa_2,GMP_2,Dress_2,Name_3,m_Name_2,Xpos_2,Ypos_2,Rmag_2,
+ log_Reff_1,log_a_1,SuBr_2,Diskp_2,b_a_2,e4_2,n_e4_2,b_ae_2,
+ Ampli_2,Bar_2,Disk_2,Spiral_2,Envelope_2,MType_2,Note_2
if(ar__3(13:15) .EQ. '') Dress_2 = iNULL__
if(ar__3(42:46) .EQ. '') Rmag_2 = rNULL__
if(ar__3(49:52) .EQ. '') log_Reff_1 = rNULL__
if(ar__3(55:58) .EQ. '') log_a_1 = rNULL__
if(ar__3(61:65) .EQ. '') SuBr_2 = rNULL__
if(ar__3(72:75) .EQ. '') b_a_2 = rNULL__
if(ar__3(78:81) .EQ. '') e4_2 = rNULL__
if(ar__3(88:91) .EQ. '') b_ae_2 = rNULL__
if(ar__3(94:95) .EQ. '') Ampli_2 = iNULL__
c ..............Just test output...........
write(6,4)
+ GMPa_2,GMP_2,Dress_2,Name_3,m_Name_2,Xpos_2,Ypos_2,Rmag_2,
+ log_Reff_1,log_a_1,SuBr_2,Diskp_2,b_a_2,e4_2,n_e4_2,b_ae_2,
+ Ampli_2,Bar_2,Disk_2,Spiral_2,Envelope_2,MType_2,Note_2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4' ! FW sample photometric parameters and
* classification
C Format for file interpretation
5 format(
+ 1X,I2,1X,I2,1X,I2,2X,I2,1X,I2,1X,I2,2X,A7,A1,2X,F5.2,2X,F4.2,
+ 2X,F4.2,2X,F5.2,2X,A2,2X,F4.2,2X,I2,2X,A2,2X,F4.2,2X,I2,2X,A4,
+ 2X,A4,2X,A5,2X,A4,2X,A8,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table4')
write(6,*) '....Loading file: table4'
do i__=1,47
read(1,'(A121)')ar__4
read(ar__4,5)
+ RAh,RAm,RAs,DEd,DEm,DEs,Name_4,m_Name_3,Rmag_3,log_Reff_1,
+ log_a_1,SuBr_3,Diskp_3,b_a_3,e4_3,n_e4_3,b_ae_3,Ampli_3,Bar_3,
+ Disk_3,Spiral_3,Envelope_3,MType_3,Note_3
if(ar__4(32:36) .EQ. '') Rmag_3 = rNULL__
if(ar__4(39:42) .EQ. '') log_Reff_1 = rNULL__
if(ar__4(45:48) .EQ. '') log_a_1 = rNULL__
if(ar__4(51:55) .EQ. '') SuBr_3 = rNULL__
if(ar__4(62:65) .EQ. '') b_a_3 = rNULL__
if(ar__4(68:69) .EQ. '') e4_3 = iNULL__
if(ar__4(76:79) .EQ. '') b_ae_3 = rNULL__
if(ar__4(82:83) .EQ. '') Ampli_3 = iNULL__
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.
c ..............Just test output...........
write(6,5)
+ RAh,RAm,RAs,DEd,DEm,DEs,Name_4,m_Name_3,Rmag_3,log_Reff_1,
+ log_a_1,SuBr_3,Diskp_3,b_a_3,e4_3,n_e4_3,b_ae_3,Ampli_3,Bar_3,
+ Disk_3,Spiral_3,Envelope_3,MType_3,Note_3
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end