Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/474/1081 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 2024-Apr-16
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+A/474/1081 Stellar population gradients in bulges. I (Gorgas+, 2007)
*================================================================================
*Stellar population gradients in bulges along the Hubble sequence. I. The data.
* Gorgas J., Jablonka P., Goudfrooij P.
* <Astron. Astrophys., 474, 1081-1092 (2007)>
* =2007A&A...474.1081G
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Galaxy sample
integer*4 nr__
parameter (nr__=32) ! Number of records
character*87 ar__ ! Full-size record
character*11 Galaxy ! Galaxy name
character*5 MType ! Morphological type
real*4 TType ! Hubble T type (1)
real*4 Btmag ! (mag) Total apparent B magnitude (1)
integer*4 i ! (deg) Galaxy inclination (1)
character*12 Bulge ! Bulge type, from Lutticke et al. (2000,
* Cat. J/A+AS/145/405)
character*11 Env ! Environment
character*3 Run ! Run number(s) in which each galaxy was observed
integer*4 PA ! (deg) Position angle along the slit
integer*4 ExpTime ! (s) Total exposure time
real*4 reff ! (arcsec) Effective radius of the bulge, computed as
* explained in Sect. 4
real*4 e_reff ! (arcsec) rms uncertainty on reff
real*4 Vmax ! (km/s) ? Maximum rotational velocity along the disk
* major axis, corrected from inclination
* when necessary (1)
*Note (1): Data have been taken from HyperLeda (http://leda.univ-lyon1.fr/)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'data1.dat' ! Line-strength indices in flux-calibrated spectra,
not transformed to the Lick/IDS system
integer*4 nr__1
parameter (nr__1=983) ! Number of records
character*424 ar__1 ! Full-size record
character*11 Galaxy_1 ! Galaxy Name
integer*4 Sigma0 ! (km/s) Velocity dispersion of the spectra (1)
real*8 Rad ! (arcsec) [-116/62] Distance to the galaxy center
real*8 HdA ! (0.1nm) ? HdeltaA index value
real*8 e_HdA ! (0.1nm) ? Uncertainty of HdeltaA index value
real*8 HdF ! (0.1nm) ? HdeltaF index value
real*8 e_HdF ! (0.1nm) ? Uncertainty of HdeltaF index value
real*8 CN1 ! (mag) ? CN1 index value
real*8 e_CN1 ! (mag) ? Uncertainty of CN1 index value
real*8 CN2 ! (mag) ? CN2 index value
real*8 e_CN2 ! (mag) ? Uncertainty of CN2 index value
real*8 Ca4227 ! (0.1nm) ? Ca4227 index value
real*8 e_Ca4227 ! (0.1nm) ? Uncertainty of Ca4227 index value
real*8 G4300 ! (0.1nm) ? G4300 index value
real*8 e_G4300 ! (0.1nm) ? Uncertainty of G4300 index value
real*8 HgA ! (0.1nm) ? HGammaA index value
real*8 e_HgA ! (0.1nm) ? Uncertainty of HGammaA index value
real*8 HgF ! (0.1nm) ? HgammaF index value
real*8 e_HgF ! (0.1nm) ? Uncertainty of HgammaF index value
real*8 Fe4383 ! (0.1nm) ? Fe4383 index value
real*8 e_Fe4383 ! (0.1nm) ? Uncertainty of Fe4383 index value
real*8 Ca4455 ! (0.1nm) ? Ca4455 index value
real*8 e_Ca4455 ! (0.1nm) ? Uncertainty of Ca4455 index value
real*8 Fe4531 ! (0.1nm) ? Fe4531 index value
real*8 e_Fe4531 ! (0.1nm) ? Uncertainty of Fe4531 index value
real*8 C4668 ! (0.1nm) ? C4668 index value
real*8 e_C4668 ! (0.1nm) ? Uncertainty of C4668 index value
real*8 HB4861 ! (0.1nm) ? HB4861 index value
real*8 e_HB4861 ! (0.1nm) ? Uncertainty of HB4861 index value
real*8 Fe5015 ! (0.1nm) ? Fe5015 index value
real*8 e_Fe5015 ! (0.1nm) ? Uncertainty of Fe5015 index value
real*8 Mg1 ! (mag) ? Mg1 index value
real*8 e_Mg1 ! (mag) ? Uncertainty of Mg1 index value
real*8 Mg2 ! (mag) ? Mg2 index value
real*8 e_Mg2 ! (mag) ? Uncertainty of Mg2 index value
real*8 Mgb ! (0.1nm) ? Mgb index value
real*8 e_Mgb ! (0.1nm) ? Uncertainty of Mgb index value
real*8 Fe5270 ! (0.1nm) ? Fe5270 index value
real*8 e_Fe5270 ! (0.1nm) ? Uncertainty of Fe5270 index value
real*8 Fe5335 ! (0.1nm) ? Fe5335 index value
real*8 e_Fe5335 ! (0.1nm) ? Uncertainty of Fe5335 index value
real*8 Fe5406 ! (0.1nm) ? Fe5406 index value
real*8 e_Fe5406 ! (0.1nm) ? Uncertainty of Fe5406 index value
real*8 Fe5709 ! (0.1nm) ? Fe5709 index value
real*8 e_Fe5709 ! (0.1nm) ? Uncertainty of Fe5709 index value
real*8 Fe5782 ! (0.1nm) ? Fe5782 index value
real*8 e_Fe5782 ! (0.1nm) ? Uncertainty of Fe5782 index value
real*8 Na5895 ! (0.1nm) ? Na5895 index value
real*8 e_Na5895 ! (0.1nm) ? Uncertainty of Na5895 index value
real*8 TiO1 ! (mag) ? TiO1 index value
real*8 e_TiO1 ! (mag) ? Uncertainty of TiO1 index value
real*8 TiO2 ! (mag) ? TiO2 index value
real*8 e_TiO2 ! (mag) ? Uncertainty of TiO2 index value
*Note (1): This indicates the velocity dispersion at which the indices have been
* measured
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'data2.dat' ! Dynamical parameters and indices transformed to
the Lick/IDS system
integer*4 nr__2
parameter (nr__2=960) ! Number of records
character*480 ar__2 ! Full-size record
character*11 Galaxy_2 ! Galaxy Name
real*4 Rad_1 ! (arcsec) [-62/62] Distance to the galaxy center
real*8 Gamma ! ? Relative line-strength parameter
real*8 e_Gamma ! ? Uncertainty on relative line-strength
real*8 RV ! (km/s) ? Radial velocity
real*8 e_RV ! (km/s) ? Uncertainty of radial velocity
real*8 Sigma ! (km/s) ? Velocity dispersion
real*8 e_Sigma ! (km/s) ? Uncertainty of velocity dispersion
real*8 HdA_1 ! (0.1nm) ? HdeltaA index value
real*8 e_HdA_1 ! (0.1nm) ? Uncertainty of HdeltaA index value
real*8 HdF_1 ! (0.1nm) ? HdeltaF index value
real*8 e_HdF_1 ! (0.1nm) ? Uncertainty of HdeltaF index value
real*8 CN1_1 ! (mag) ? CN1 index value
real*8 e_CN1_1 ! (mag) ? Uncertainty of CN1 index value
real*8 CN2_1 ! (mag) ? CN2 index value
real*8 e_CN2_1 ! (mag) ? Uncertainty of CN2 index value
real*8 Ca4227_1 ! (0.1nm) ? Ca4227 index value
real*8 e_Ca4227_1 ! (0.1nm) ? Uncertainty of Ca4227 index value
real*8 G4300_1 ! (0.1nm) ? G4300 index value
real*8 e_G4300_1 ! (0.1nm) ? Uncertainty of G4300 index value
real*8 HgA_1 ! (0.1nm) ? HGammaA index value
real*8 e_HgA_1 ! (0.1nm) ? Uncertainty of HGammaA index value
real*8 HgF_1 ! (0.1nm) ? HgammaF index value
real*8 e_HgF_1 ! (0.1nm) ? Uncertainty of HgammaF index value
real*8 Fe4383_1 ! (0.1nm) ? Fe4383 index value
real*8 e_Fe4383_1 ! (0.1nm) ? Uncertainty of Fe4383 index value
real*8 Ca4455_1 ! (0.1nm) ? Ca4455 index value
real*8 e_Ca4455_1 ! (0.1nm) ? Uncertainty of Ca4455 index value
real*8 Fe4531_1 ! (0.1nm) ? Fe4531 index value
real*8 e_Fe4531_1 ! (0.1nm) ? Uncertainty of Fe4531 index value
real*8 C4668_1 ! (0.1nm) ? C4668 index value
real*8 e_C4668_1 ! (0.1nm) ? Uncertainty of C4668 index value
real*8 HB4861_1 ! (0.1nm) ? HB4861 index value
real*8 e_HB4861_1 ! (0.1nm) ? Uncertainty of HB4861 index value
real*8 Fe5015_1 ! (0.1nm) ? Fe5015 index value
real*8 e_Fe5015_1 ! (0.1nm) ? Uncertainty of Fe5015 index value
real*8 Mg1_1 ! (mag) ? Mg1 index value
real*8 e_Mg1_1 ! (mag) ? Uncertainty of Mg1 index value
real*8 Mg2_1 ! (mag) ? Mg2 index value
real*8 e_Mg2_1 ! (mag) ? Uncertainty of Mg2 index value
real*8 Mgb_1 ! (0.1nm) ? Mgb index value
real*8 e_Mgb_1 ! (0.1nm) ? Uncertainty of Mgb index value
real*8 Fe5270_1 ! (0.1nm) ? Fe5270 index value
real*8 e_Fe5270_1 ! (0.1nm) ? Uncertainty of Fe5270 index value
real*8 Fe5335_1 ! (0.1nm) ? Fe5335 index value
real*8 e_Fe5335_1 ! (0.1nm) ? Uncertainty of Fe5335 index value
real*8 Fe5406_1 ! (0.1nm) ? Fe5406 index value
real*8 e_Fe5406_1 ! (0.1nm) ? Uncertainty of Fe5406 index value
real*8 Fe5709_1 ! (0.1nm) ? Fe5709 index value
real*8 e_Fe5709_1 ! (0.1nm) ? Uncertainty of Fe5709 index value
real*8 Fe5782_1 ! (0.1nm) ? Fe5782 index value
real*8 e_Fe5782_1 ! (0.1nm) ? Uncertainty of Fe5782 index value
real*8 Na5895_1 ! (0.1nm) ? Na5895 index value
real*8 e_Na5895_1 ! (0.1nm) ? Uncertainty of Na5895 index value
real*8 TiO1_1 ! (mag) ? TiO1 index value
real*8 e_TiO1_1 ! (mag) ? Uncertainty of TiO1 index value
real*8 TiO2_1 ! (mag) ? TiO2 index value
real*8 e_TiO2_1 ! (mag) ? Uncertainty of TiO2 index value
real*8 EW_OIII_ ! (0.1nm) ? OIII 5007 equivalent width
real*8 e_EW_OIII_ ! (0.1nm) ? Uncertainty of [OIII] 5007 equivalent width
C=============================================================================
C Loading file 'table1.dat' ! Galaxy sample
C Format for file interpretation
1 format(
+ A11,2X,A5,1X,F4.1,1X,F5.2,1X,I2,1X,A12,1X,A11,1X,A3,1X,I3,1X,
+ I5,1X,F4.1,1X,F4.1,1X,F5.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,32
read(1,'(A87)')ar__
read(ar__,1)
+ Galaxy,MType,TType,Btmag,i,Bulge,Env,Run,PA,ExpTime,reff,
+ e_reff,Vmax
if(ar__(83:87) .EQ. '') Vmax = rNULL__
c ..............Just test output...........
write(6,1)
+ Galaxy,MType,TType,Btmag,i,Bulge,Env,Run,PA,ExpTime,reff,
+ e_reff,Vmax
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'data1.dat' ! Line-strength indices in flux-calibrated spectra,
* not transformed to the Lick/IDS system
C Format for file interpretation
2 format(
+ A11,I6,F7.1,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'data1.dat')
write(6,*) '....Loading file: data1.dat'
do i__=1,983
read(1,'(A424)')ar__1
read(ar__1,2)
+ Galaxy_1,Sigma0,Rad,HdA,e_HdA,HdF,e_HdF,CN1,e_CN1,CN2,e_CN2,
+ Ca4227,e_Ca4227,G4300,e_G4300,HgA,e_HgA,HgF,e_HgF,Fe4383,
+ e_Fe4383,Ca4455,e_Ca4455,Fe4531,e_Fe4531,C4668,e_C4668,HB4861,
+ e_HB4861,Fe5015,e_Fe5015,Mg1,e_Mg1,Mg2,e_Mg2,Mgb,e_Mgb,Fe5270,
+ e_Fe5270,Fe5335,e_Fe5335,Fe5406,e_Fe5406,Fe5709,e_Fe5709,
+ Fe5782,e_Fe5782,Na5895,e_Na5895,TiO1,e_TiO1,TiO2,e_TiO2
if(ar__1(25:32) .EQ. '') HdA = rNULL__
if(ar__1(33:40) .EQ. '') e_HdA = rNULL__
if(ar__1(41:48) .EQ. '') HdF = rNULL__
if(ar__1(49:56) .EQ. '') e_HdF = rNULL__
if(ar__1(57:64) .EQ. '') CN1 = rNULL__
if(ar__1(65:72) .EQ. '') e_CN1 = rNULL__
if(ar__1(73:80) .EQ. '') CN2 = rNULL__
if(ar__1(81:88) .EQ. '') e_CN2 = rNULL__
if(ar__1(89:96) .EQ. '') Ca4227 = rNULL__
if(ar__1(97:104) .EQ. '') e_Ca4227 = rNULL__
if(ar__1(105:112) .EQ. '') G4300 = rNULL__
if(ar__1(113:120) .EQ. '') e_G4300 = rNULL__
if(ar__1(121:128) .EQ. '') HgA = rNULL__
if(ar__1(129:136) .EQ. '') e_HgA = rNULL__
if(ar__1(137:144) .EQ. '') HgF = rNULL__
if(ar__1(145:152) .EQ. '') e_HgF = rNULL__
if(ar__1(153:160) .EQ. '') Fe4383 = rNULL__
if(ar__1(161:168) .EQ. '') e_Fe4383 = rNULL__
if(ar__1(169:176) .EQ. '') Ca4455 = rNULL__
if(ar__1(177:184) .EQ. '') e_Ca4455 = rNULL__
if(ar__1(185:192) .EQ. '') Fe4531 = rNULL__
if(ar__1(193:200) .EQ. '') e_Fe4531 = rNULL__
if(ar__1(201:208) .EQ. '') C4668 = rNULL__
if(ar__1(209:216) .EQ. '') e_C4668 = rNULL__
if(ar__1(217:224) .EQ. '') HB4861 = rNULL__
if(ar__1(225:232) .EQ. '') e_HB4861 = rNULL__
if(ar__1(233:240) .EQ. '') Fe5015 = rNULL__
if(ar__1(241:248) .EQ. '') e_Fe5015 = rNULL__
if(ar__1(249:256) .EQ. '') Mg1 = rNULL__
if(ar__1(257:264) .EQ. '') e_Mg1 = rNULL__
if(ar__1(265:272) .EQ. '') Mg2 = rNULL__
if(ar__1(273:280) .EQ. '') e_Mg2 = rNULL__
if(ar__1(281:288) .EQ. '') Mgb = rNULL__
if(ar__1(289:296) .EQ. '') e_Mgb = rNULL__
if(ar__1(297:304) .EQ. '') Fe5270 = rNULL__
if(ar__1(305:312) .EQ. '') e_Fe5270 = rNULL__
if(ar__1(313:320) .EQ. '') Fe5335 = rNULL__
if(ar__1(321:328) .EQ. '') e_Fe5335 = rNULL__
if(ar__1(329:336) .EQ. '') Fe5406 = rNULL__
if(ar__1(337:344) .EQ. '') e_Fe5406 = rNULL__
if(ar__1(345:352) .EQ. '') Fe5709 = rNULL__
if(ar__1(353:360) .EQ. '') e_Fe5709 = rNULL__
if(ar__1(361:368) .EQ. '') Fe5782 = rNULL__
if(ar__1(369:376) .EQ. '') e_Fe5782 = rNULL__
if(ar__1(377:384) .EQ. '') Na5895 = rNULL__
if(ar__1(385:392) .EQ. '') e_Na5895 = rNULL__
if(ar__1(393:400) .EQ. '') TiO1 = rNULL__
if(ar__1(401:408) .EQ. '') e_TiO1 = rNULL__
if(ar__1(409:416) .EQ. '') TiO2 = rNULL__
if(ar__1(417:424) .EQ. '') e_TiO2 = rNULL__
c ..............Just test output...........
write(6,2)
+ Galaxy_1,Sigma0,Rad,HdA,e_HdA,HdF,e_HdF,CN1,e_CN1,CN2,e_CN2,
+ Ca4227,e_Ca4227,G4300,e_G4300,HgA,e_HgA,HgF,e_HgF,Fe4383,
+ e_Fe4383,Ca4455,e_Ca4455,Fe4531,e_Fe4531,C4668,e_C4668,HB4861,
+ e_HB4861,Fe5015,e_Fe5015,Mg1,e_Mg1,Mg2,e_Mg2,Mgb,e_Mgb,Fe5270,
+ e_Fe5270,Fe5335,e_Fe5335,Fe5406,e_Fe5406,Fe5709,e_Fe5709,
+ Fe5782,e_Fe5782,Na5895,e_Na5895,TiO1,e_TiO1,TiO2,e_TiO2
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'data2.dat' ! Dynamical parameters and indices transformed to
* the Lick/IDS system
C Format for file interpretation
3 format(
+ A11,F6.1,F8.3,F8.3,F9.1,F7.1,F8.1,F7.1,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,
+ F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3,F8.3)
C Effective file loading
open(unit=1,status='old',file=
+'data2.dat')
write(6,*) '....Loading file: data2.dat'
do i__=1,960
read(1,'(A480)')ar__2
read(ar__2,3)
+ Galaxy_2,Rad_1,Gamma,e_Gamma,RV,e_RV,Sigma,e_Sigma,HdA_1,
+ e_HdA_1,HdF_1,e_HdF_1,CN1_1,e_CN1_1,CN2_1,e_CN2_1,Ca4227_1,
+ e_Ca4227_1,G4300_1,e_G4300_1,HgA_1,e_HgA_1,HgF_1,e_HgF_1,
+ Fe4383_1,e_Fe4383_1,Ca4455_1,e_Ca4455_1,Fe4531_1,e_Fe4531_1,
+ C4668_1,e_C4668_1,HB4861_1,e_HB4861_1,Fe5015_1,e_Fe5015_1,
+ Mg1_1,e_Mg1_1,Mg2_1,e_Mg2_1,Mgb_1,e_Mgb_1,Fe5270_1,e_Fe5270_1,
+ Fe5335_1,e_Fe5335_1,Fe5406_1,e_Fe5406_1,Fe5709_1,e_Fe5709_1,
+ Fe5782_1,e_Fe5782_1,Na5895_1,e_Na5895_1,TiO1_1,e_TiO1_1,
+ TiO2_1,e_TiO2_1,EW_OIII_,e_EW_OIII_
if(ar__2(18:25) .EQ. '') Gamma = rNULL__
if(ar__2(26:33) .EQ. '') e_Gamma = rNULL__
if(ar__2(34:42) .EQ. '') RV = rNULL__
if(ar__2(43:49) .EQ. '') e_RV = rNULL__
if(ar__2(50:57) .EQ. '') Sigma = rNULL__
if(ar__2(58:64) .EQ. '') e_Sigma = rNULL__
if(ar__2(65:72) .EQ. '') HdA_1 = rNULL__
if(ar__2(73:80) .EQ. '') e_HdA_1 = rNULL__
if(ar__2(81:88) .EQ. '') HdF_1 = rNULL__
if(ar__2(89:96) .EQ. '') e_HdF_1 = rNULL__
if(ar__2(97:104) .EQ. '') CN1_1 = rNULL__
if(ar__2(105:112) .EQ. '') e_CN1_1 = rNULL__
if(ar__2(113:120) .EQ. '') CN2_1 = rNULL__
if(ar__2(121:128) .EQ. '') e_CN2_1 = rNULL__
if(ar__2(129:136) .EQ. '') Ca4227_1 = rNULL__
if(ar__2(137:144) .EQ. '') e_Ca4227_1 = rNULL__
if(ar__2(145:152) .EQ. '') G4300_1 = rNULL__
if(ar__2(153:160) .EQ. '') e_G4300_1 = rNULL__
if(ar__2(161:168) .EQ. '') HgA_1 = rNULL__
if(ar__2(169:176) .EQ. '') e_HgA_1 = rNULL__
if(ar__2(177:184) .EQ. '') HgF_1 = rNULL__
if(ar__2(185:192) .EQ. '') e_HgF_1 = rNULL__
if(ar__2(193:200) .EQ. '') Fe4383_1 = rNULL__
if(ar__2(201:208) .EQ. '') e_Fe4383_1 = rNULL__
if(ar__2(209:216) .EQ. '') Ca4455_1 = rNULL__
if(ar__2(217:224) .EQ. '') e_Ca4455_1 = rNULL__
if(ar__2(225:232) .EQ. '') Fe4531_1 = rNULL__
if(ar__2(233:240) .EQ. '') e_Fe4531_1 = rNULL__
if(ar__2(241:248) .EQ. '') C4668_1 = rNULL__
if(ar__2(249:256) .EQ. '') e_C4668_1 = rNULL__
if(ar__2(257:264) .EQ. '') HB4861_1 = rNULL__
if(ar__2(265:272) .EQ. '') e_HB4861_1 = rNULL__
if(ar__2(273:280) .EQ. '') Fe5015_1 = rNULL__
if(ar__2(281:288) .EQ. '') e_Fe5015_1 = rNULL__
if(ar__2(289:296) .EQ. '') Mg1_1 = rNULL__
if(ar__2(297:304) .EQ. '') e_Mg1_1 = rNULL__
if(ar__2(305:312) .EQ. '') Mg2_1 = rNULL__
if(ar__2(313:320) .EQ. '') e_Mg2_1 = rNULL__
if(ar__2(321:328) .EQ. '') Mgb_1 = rNULL__
if(ar__2(329:336) .EQ. '') e_Mgb_1 = rNULL__
if(ar__2(337:344) .EQ. '') Fe5270_1 = rNULL__
if(ar__2(345:352) .EQ. '') e_Fe5270_1 = rNULL__
if(ar__2(353:360) .EQ. '') Fe5335_1 = rNULL__
if(ar__2(361:368) .EQ. '') e_Fe5335_1 = rNULL__
if(ar__2(369:376) .EQ. '') Fe5406_1 = rNULL__
if(ar__2(377:384) .EQ. '') e_Fe5406_1 = rNULL__
if(ar__2(385:392) .EQ. '') Fe5709_1 = rNULL__
if(ar__2(393:400) .EQ. '') e_Fe5709_1 = rNULL__
if(ar__2(401:408) .EQ. '') Fe5782_1 = rNULL__
if(ar__2(409:416) .EQ. '') e_Fe5782_1 = rNULL__
if(ar__2(417:424) .EQ. '') Na5895_1 = rNULL__
if(ar__2(425:432) .EQ. '') e_Na5895_1 = rNULL__
if(ar__2(433:440) .EQ. '') TiO1_1 = rNULL__
if(ar__2(441:448) .EQ. '') e_TiO1_1 = rNULL__
if(ar__2(449:456) .EQ. '') TiO2_1 = rNULL__
if(ar__2(457:464) .EQ. '') e_TiO2_1 = rNULL__
if(ar__2(465:472) .EQ. '') EW_OIII_ = rNULL__
if(ar__2(473:480) .EQ. '') e_EW_OIII_ = rNULL__
c ..............Just test output...........
write(6,3)
+ Galaxy_2,Rad_1,Gamma,e_Gamma,RV,e_RV,Sigma,e_Sigma,HdA_1,
+ e_HdA_1,HdF_1,e_HdF_1,CN1_1,e_CN1_1,CN2_1,e_CN2_1,Ca4227_1,
+ e_Ca4227_1,G4300_1,e_G4300_1,HgA_1,e_HgA_1,HgF_1,e_HgF_1,
+ Fe4383_1,e_Fe4383_1,Ca4455_1,e_Ca4455_1,Fe4531_1,e_Fe4531_1,
+ C4668_1,e_C4668_1,HB4861_1,e_HB4861_1,Fe5015_1,e_Fe5015_1,
+ Mg1_1,e_Mg1_1,Mg2_1,e_Mg2_1,Mgb_1,e_Mgb_1,Fe5270_1,e_Fe5270_1,
+ Fe5335_1,e_Fe5335_1,Fe5406_1,e_Fe5406_1,Fe5709_1,e_Fe5709_1,
+ Fe5782_1,e_Fe5782_1,Na5895_1,e_Na5895_1,TiO1_1,e_TiO1_1,
+ TiO2_1,e_TiO2_1,EW_OIII_,e_EW_OIII_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end