Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/153/96 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/AJ/153/96 Standard Galactic field RR Lyrae. I. Photometry (Monson+, 2017)
*================================================================================
*Standard Galactic field RR Lyrae. I. Optical to mid-infrared phased photometry.
* Monson A.J., Beaton R.L., Scowcroft V., Freedman W.L., Madore B.F.,
* Rich J.A., Seibert M., Kollmeier J.A., Clementini G.
* <Astron. J., 153, 96-96 (2017)>
* =2017AJ....153...96M (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! RR Lyrae (RRL) Galactic calibrators and ephemerides
integer*4 nr__
parameter (nr__=55) ! Number of records
character*78 ar__ ! Full-size record
character*9 Name (nr__) ! Star name
real*8 FPer (nr__) ! (d) [0.25/0.74] Final period
real*8 HJD (nr__) ! (d) Three-hundred MilliMeter Telescope (TMMT) HJD-max
real*4 zeta (nr__) ! (d/yr) [-0.0015/0.0019]? Quadratic O-C shape term
* ({zeta}), if required
character*4 RRL (nr__) ! RR Lyrae class (either RRab, RRc, or RRd)
real*8 BLPer (nr__) ! (d) [7.2/1748.9]? Blazhko effect period
real*4 v_Fe_H_ (nr__) ! ([Sun]) [-2.6/-0.07] Metallicity (1)
character*1 r__Fe_H_ (nr__) ! Reference flag on [Fe/H] (2)
character*3 HIP (nr__) ! [HIP] Parallax from Hipparcos (HIP) (3)
character*3 BW (nr__) ! Parallax from Baade-Wesselink (3)
character*3 HST (nr__) ! [HST] Parallax from HST (HST) (3)
*Note (1): Unless otherwise noted, values are taken from Feast et al. 2008
* (Cat. J/MNRAS/386/2115), but the measurements were first compiled by
* Fernley et al. 1998 (Cat. J/A+A/330/515) and references therein and are on a
* metallicity scale defined by Fernley & Barnes 1997 (Cat. J/A+AS/125/313) and
* references therein.
*Note (2): Reference codesa are defined as follows:
* c = Fernley et al. (1998MNRAS.293L..61F);
* d = Clementini et al. (2000AJ....120.2054C).
*Note (3): The codes indicating where the parallax is derived from are defined as
* follows:
* HIP = Hipparcos (see Cat. I/239);
* HST = Hubble Space Telescope (see Benedict et al. 2011, Cat. J/AJ/142/187);
* 1 = Liu & Janes 1990ApJ...354..273L;
* 2 = Jones et al. 1992ApJ...386..646J;
* 3 = Skillen et al. 1993MNRAS.265..301S;
* 4 = Fernley et al. 1989MNRAS.236..447F;
* 5 = Jones et al. 1988ApJ...332..206J;
* 6 = Cacciari et al. 1989A&A...209..154C;
* 7 = Fernley et al. 1990MNRAS.242..685F;
* 8 = Skillen et al. 1989MNRAS.241..281S.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Three-hundred MilliMeter Telescope (TMMT)
photometry and phased archival data
integer*4 nr__1
parameter (nr__1=90243) ! Number of records
character*54 ar__1 ! Full-size record
character*9 Name_1 (nr__1) ! Star name
character*2 Flt (nr__1) ! Filter (G1)
real*4 mag (nr__1) ! (mag) [6.1/20.3] Vega magnitude in Flt
real*4 e_mag (nr__1) ! (mag) [0.001/3.9]?=99.005 Uncertainty in mag
* ({sigma}_phot_)
real*4 err (nr__1) ! (mag) [0.003/0.05] Systematic error in the zero-point
* determination ({sigma}_sys_)
real*8 MJD (nr__1) ! (d) Heliocentric Modified Julian Day (HJD-2400000.5)
real*4 Phase (nr__1) ! (deg) [0/1] Final derived phase ({Phi})
integer*4 Ref (nr__1) ! [0/999]? References (1)
*Note (1): Reference codes are defined as follows:
* 0 = Three-hundred MilliMeter Telescope (TMMT), this work;
* 1 = Spitzer, this work;
* 4 = Skillen et al. (1993MNRAS.265..301S);
* 5 = Barnes et al. (1992PASP..104..514B);
* 7 = Liu & Janes (1989ApJS...69..593L);
* 8 = Liu & Janes (1990ApJ...354..273L);
* 9 = Barcza & Benko (2014MNRAS.442.1863B);
* 10 = Paczynski (1965AcA....15..115P);
* 11 = 2MASS (Cutri et al. 2003, Cat. II/246);
* 17 = ASAS Pojmanski (1997AcA....47..467P);
* 15 = Jones et al. (1992ApJ...386..646J);
* 19 = IBVS Broglia & Conconi (1992IBVS.3748....1B);
* 31 = Fernley et al. (1990MNRAS.242..685F);
* 41 = Fernley et al. (1989MNRAS.236..447F);
* 98 = Clementini et al. (1990A&AS...85..865C);
* 99 = Clementini et al. (2000AJ....120.2054C);
* 999 = Three-hundred MilliMeter Telescope (TMMT), modified for Blazhko
* effect.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Gaussian-windowed LOcal regrESSion (GLOESS) light
curves
integer*4 nr__2
parameter (nr__2=94208) ! Number of records
character*27 ar__2 ! Full-size record
character*9 Name_2 (nr__2) ! Star name
character*2 Flt_1 (nr__2) ! Filter (G1)
real*8 Phase_1 (nr__2) ! (deg) [0/1] Phase
real*4 mag_1 (nr__2) ! (mag) [6.3/14.1] Gaussian-windowed LOcal regrESSion
* (GLOESS) Vega magnitude (1)
*Note (1): Mean magnitudes were determined by computing the mean intensity of the
* evenly sampled GLOESS fit points to the light curve, then converting back
* to a magnitude. The light curve must be sampled with enough data points to
* capture all the nuances of the shape for an accurate mean; in our case we
* sampled with 256 data points spaced every 1/256 in phase. Generally, GLOESS
* fits were determined only for those stars and bands that had more than 20
* individual data points over a reasonable portion of phase (i.e., 20 data
* points only spanning {Phi}~0.1 would not have a GLOESS fit nor a mean
* magnitude).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Intensity mean magnitudes from Gaussian-windowed
LOcal regrESSion (GLOESS) light curves
integer*4 nr__3
parameter (nr__3=55) ! Number of records
character*139 ar__3 ! Full-size record
character*9 Name_3 (nr__3) ! Star name
real*4 v_Umag_ (nr__3) ! (mag) [8/13]? Intensity mean Johnson U band magnitude
real*4 e__Umag_ (nr__3) ! (mag) [0.01/0.021]? Uncertainty in Umag (4)
real*4 v_Bmag_ (nr__3) ! (mag) [8/14]? Intensity mean Johnson B band magnitude
real*4 e__Bmag_ (nr__3) ! (mag) [0.003/0.021]? Uncertainty in Bmag (4)
real*4 v_Vmag_ (nr__3) ! (mag) [7.7/13.4]? Intensity mean Johnson V band
* magnitude
real*4 e__Vmag_ (nr__3) ! (mag) [0.003/0.021]? Uncertainty in Vmag (4)
real*4 v_Rmag_ (nr__3) ! (mag) [9.2/11.9]? Intensity mean Kron-Cousins R_C_
* band magnitude
real*4 e__Rmag_ (nr__3) ! (mag) [0.02/0.02]? Uncertainty in Rmag (4)
real*4 v_Imag_ (nr__3) ! (mag) [7.2/12.88]? Intensity mean Kron-Cousins I_C_
* band magnitude
real*4 e__Imag_ (nr__3) ! (mag) [0.003/0.028]? Uncertainty in Imag (4)
real*4 v_Jmag_ (nr__3) ! (mag) [6.74/11.21]? Intensity mean 2MASS J band
* magnitude
real*4 e__Jmag_ (nr__3) ! (mag) [0.009/0.02]? Uncertainty in Jmag (4)
real*4 v_Hmag_ (nr__3) ! (mag) [6.5/11]? Intensity mean 2MASS H band magnitude
real*4 e__Hmag_ (nr__3) ! (mag) [0.013/0.013]? Uncertainty in Hmag (4)
real*4 v_Ksmag_ (nr__3) ! (mag) [6.49/10.89]? Intensity mean 2MASS K_s_ band
* magnitude
real*4 e__Ksmag_ (nr__3) ! (mag) [0.009/0.009]? Uncertainty in Ksmag (4)
real*4 v__3_6__ (nr__3) ! (mag) [6.47/12.27]? Intensity mean Spitzer/IRAC
* 3.6{mu}m band magnitude
real*4 e___3_6__ (nr__3) ! (mag) [0.008/0.011]? Uncertainty in [3.6] (4)
real*4 v__4_5__ (nr__3) ! (mag) [6.46/12.25]? Intensity mean Spitzer/IRAC
* 4.5{mu}m band magnitude
real*4 e___4_5__ (nr__3) ! (mag) [0.008/0.011]? Uncertainty in [4.5] (4)
*Note (4):
* The random uncertainty of the Gaussian-windowed LOcal regrESSion
* (GLOESS)-derived mean magnitude is simply the error on the mean of data
* points going into GLOESS fitting. Thus, stars with more data points will
* have a smaller uncertainty in GLOESS mean magnitude. The systematic
* uncertainty is determined by the photometric transformations, either in
* transforming our Three-hundred MilliMeter Telescope (TMMT) photometry onto
* an absolute system (see Figure 2) or as reported in the literature and in
* transforming from other filter systems (as described in Section 4). The
* final reported error is:
* {sigma}^2^=1/{sum}(1/{sigma}^2^_phot_)+1/{sum}(1/{sigma}^2^_sys_),
* where the sum over {sigma}_sys_ includes only the unique entries from each
* reference, i.e., it is not counted for every measurement.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Amplitudes, minimum magnitude and rise times from
Gaussian-windowed LOcal regrESSion (GLOESS) light
curves
integer*4 nr__4
parameter (nr__4=55) ! Number of records
character*199 ar__4 ! Full-size record
character*9 Name_4 (nr__4) ! Star name
real*4 Uamp (nr__4) ! (mag) [0.46/1.74]? Amplitude in Johnson U band (5)
real*4 Umin (nr__4) ! (mag) [8.5/13.4]? Minimum U band magnitude at HJDmax
real*4 RTU (nr__4) ! [0.1/0.44]? U band rise time (5)
real*4 Bamp (nr__4) ! (mag) [0.46/1.7]? Amplitude in Johnson B band (5)
real*4 Bmin (nr__4) ! (mag) [8.5/14.1]? Minimum B band magnitude at HJDmax
real*4 RTB (nr__4) ! [0.1/0.51]? B band rise time (5)
real*4 Vamp (nr__4) ! (mag) [0.35/1.31]? Amplitude in Johnson V band (5)
real*4 Vmin (nr__4) ! (mag) [8/13.6]? Minimum V band magnitude at HJDmax
real*4 RTV (nr__4) ! [0.1/0.5]? V band rise time (5)
real*4 Ramp (nr__4) ! (mag) [0.3/1.1]? Amplitude in Kron-Cousins R_C_ band (5)
real*4 Rmin (nr__4) ! (mag) [9.5/12.3]? Minimum R_C_ band magnitude at HJDmax
real*4 RTR (nr__4) ! [0.13/0.45]? R_C_ band rise time (5)
real*4 Iamp (nr__4) ! (mag) [0.1/0.9]? Amplitude in Kron-Cousins I_C_ band (5)
real*4 Imin (nr__4) ! (mag) [7.5/13.1]? Minimum I_C_ band magnitude at HJDmax
real*4 RTI (nr__4) ! [0.1/0.48]? I_C_ band rise time (5)
real*4 Jamp (nr__4) ! (mag) [0.2/0.6]? Amplitude in 2MASS J band (5)
real*4 Jmin (nr__4) ! (mag) [6.9/11.5]? Minimum J band magnitude at HJDmax
real*4 RTJ (nr__4) ! [0.12/0.5]? J band rise time (5)
real*4 Hamp (nr__4) ! (mag) [0.23/0.38]? Amplitude in 2MASS H band (5)
real*4 Hmin (nr__4) ! (mag) [6.7/11.3]? Minimum H band magnitude at HJDmax
real*4 RTH (nr__4) ! [0.18/0.52]? H band rise time (5)
real*4 Ksamp (nr__4) ! (mag) [0.09/0.34]? Amplitude in 2MASS Ks band (5)
real*4 Ksmin (nr__4) ! (mag) [6.6/11.2]? Minimum Ks band magnitude at HJDmax
real*4 RTKs (nr__4) ! [0.19/0.56]? Ks band rise time (5)
real*4 v3_6amp (nr__4) ! (mag) [0.05/0.35]? Amplitude in Spitzer/IRAC 3.6{mu}m
* band (5)
real*4 v3_6min (nr__4) ! (mag) [6.6/12.4]? Minimum IRAC 3.6{mu}m magnitude at
* HJDmax
real*4 RT3_6 (nr__4) ! [0.17/0.6]? IRAC 3.6{mu}m rise time (5)
real*4 v4_5amp (nr__4) ! (mag) [0.06/0.36]? Amplitude in IRAC 4.5{mu}m band (5)
real*4 v4_5min (nr__4) ! (mag) [6.6/12.4]? Minimum IRAC 4.5{mu}m band magnitude
* at HJDmax
real*4 RT4_5 (nr__4) ! [0.16/0.93]? IRAC 4.5{mu}m band rise time
*Note (5): Amplitude (a_{lambda}_) and rise time (RT_{lambda}_) are defined as
* the difference in magnitude (mag) and phase (dimensionless), respectively,
* between the minimum and maximum of the Gaussian-windowed LOcal regrESSion
* (GLOESS) light curve.
C=============================================================================
C Loading file 'table1.dat' ! RR Lyrae (RRL) Galactic calibrators and ephemerides
C Format for file interpretation
1 format(
+ A9,1X,F11.9,1X,F12.4,1X,E10.3,1X,A4,1X,F7.2,1X,F5.2,1X,A1,1X,
+ A3,1X,A3,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,55
read(1,'(A78)')ar__
read(ar__,1)
+ Name(i__),FPer(i__),HJD(i__),zeta(i__),RRL(i__),BLPer(i__),
+ v_Fe_H_(i__),r__Fe_H_(i__),HIP(i__),BW(i__),HST(i__)
if(ar__(36:45) .EQ. '') zeta(i__) = rNULL__
if(ar__(52:58) .EQ. '') BLPer(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),FPer(i__),HJD(i__),zeta(i__),RRL(i__),BLPer(i__),
+ v_Fe_H_(i__),r__Fe_H_(i__),HIP(i__),BW(i__),HST(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Three-hundred MilliMeter Telescope (TMMT)
* photometry and phased archival data
C Format for file interpretation
2 format(A9,1X,A2,1X,F6.3,1X,F6.3,1X,F5.3,1X,F10.4,1X,F5.3,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,90243
read(1,'(A54)')ar__1
read(ar__1,2)
+ Name_1(i__),Flt(i__),mag(i__),e_mag(i__),err(i__),MJD(i__),
+ Phase(i__),Ref(i__)
if(ar__1(51:54) .EQ. '') Ref(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Flt(i__),mag(i__),e_mag(i__),err(i__),MJD(i__),
+ Phase(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Gaussian-windowed LOcal regrESSion (GLOESS) light
* curves
C Format for file interpretation
3 format(A9,1X,A2,1X,F7.5,1X,F6.3)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,94208
read(1,'(A27)')ar__2
read(ar__2,3)Name_2(i__),Flt_1(i__),Phase_1(i__),mag_1(i__)
c ..............Just test output...........
write(6,3)Name_2(i__),Flt_1(i__),Phase_1(i__),mag_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Intensity mean magnitudes from Gaussian-windowed
* LOcal regrESSion (GLOESS) light curves
C Format for file interpretation
4 format(
+ A9,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,55
read(1,'(A139)')ar__3
read(ar__3,4)
+ Name_3(i__),v_Umag_(i__),e__Umag_(i__),v_Bmag_(i__),
+ e__Bmag_(i__),v_Vmag_(i__),e__Vmag_(i__),v_Rmag_(i__),
+ e__Rmag_(i__),v_Imag_(i__),e__Imag_(i__),v_Jmag_(i__),
+ e__Jmag_(i__),v_Hmag_(i__),e__Hmag_(i__),v_Ksmag_(i__),
+ e__Ksmag_(i__),v__3_6__(i__),e___3_6__(i__),v__4_5__(i__),
+ e___4_5__(i__)
if(ar__3(11:16) .EQ. '') v_Umag_(i__) = rNULL__
if(ar__3(18:22) .EQ. '') e__Umag_(i__) = rNULL__
if(ar__3(24:29) .EQ. '') v_Bmag_(i__) = rNULL__
if(ar__3(31:35) .EQ. '') e__Bmag_(i__) = rNULL__
if(ar__3(37:42) .EQ. '') v_Vmag_(i__) = rNULL__
if(ar__3(44:48) .EQ. '') e__Vmag_(i__) = rNULL__
if(ar__3(50:55) .EQ. '') v_Rmag_(i__) = rNULL__
if(ar__3(57:61) .EQ. '') e__Rmag_(i__) = rNULL__
if(ar__3(63:68) .EQ. '') v_Imag_(i__) = rNULL__
if(ar__3(70:74) .EQ. '') e__Imag_(i__) = rNULL__
if(ar__3(76:81) .EQ. '') v_Jmag_(i__) = rNULL__
if(ar__3(83:87) .EQ. '') e__Jmag_(i__) = rNULL__
if(ar__3(89:94) .EQ. '') v_Hmag_(i__) = rNULL__
if(ar__3(96:100) .EQ. '') e__Hmag_(i__) = rNULL__
if(ar__3(102:107) .EQ. '') v_Ksmag_(i__) = rNULL__
if(ar__3(109:113) .EQ. '') e__Ksmag_(i__) = rNULL__
if(ar__3(115:120) .EQ. '') v__3_6__(i__) = rNULL__
if(ar__3(122:126) .EQ. '') e___3_6__(i__) = rNULL__
if(ar__3(128:133) .EQ. '') v__4_5__(i__) = rNULL__
if(ar__3(135:139) .EQ. '') e___4_5__(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3(i__),v_Umag_(i__),e__Umag_(i__),v_Bmag_(i__),
+ e__Bmag_(i__),v_Vmag_(i__),e__Vmag_(i__),v_Rmag_(i__),
+ e__Rmag_(i__),v_Imag_(i__),e__Imag_(i__),v_Jmag_(i__),
+ e__Jmag_(i__),v_Hmag_(i__),e__Hmag_(i__),v_Ksmag_(i__),
+ e__Ksmag_(i__),v__3_6__(i__),e___3_6__(i__),v__4_5__(i__),
+ e___4_5__(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Amplitudes, minimum magnitude and rise times from
* Gaussian-windowed LOcal regrESSion (GLOESS) light
* curves
C Format for file interpretation
5 format(
+ A9,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,
+ 1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,55
read(1,'(A199)')ar__4
read(ar__4,5)
+ Name_4(i__),Uamp(i__),Umin(i__),RTU(i__),Bamp(i__),Bmin(i__),
+ RTB(i__),Vamp(i__),Vmin(i__),RTV(i__),Ramp(i__),Rmin(i__),
+ RTR(i__),Iamp(i__),Imin(i__),RTI(i__),Jamp(i__),Jmin(i__),
+ RTJ(i__),Hamp(i__),Hmin(i__),RTH(i__),Ksamp(i__),Ksmin(i__),
+ RTKs(i__),v3_6amp(i__),v3_6min(i__),RT3_6(i__),v4_5amp(i__),
+ v4_5min(i__),RT4_5(i__)
if(ar__4(11:15) .EQ. '') Uamp(i__) = rNULL__
if(ar__4(17:22) .EQ. '') Umin(i__) = rNULL__
if(ar__4(24:28) .EQ. '') RTU(i__) = rNULL__
if(ar__4(30:34) .EQ. '') Bamp(i__) = rNULL__
if(ar__4(36:41) .EQ. '') Bmin(i__) = rNULL__
if(ar__4(43:47) .EQ. '') RTB(i__) = rNULL__
if(ar__4(49:53) .EQ. '') Vamp(i__) = rNULL__
if(ar__4(55:60) .EQ. '') Vmin(i__) = rNULL__
if(ar__4(62:66) .EQ. '') RTV(i__) = rNULL__
if(ar__4(68:72) .EQ. '') Ramp(i__) = rNULL__
if(ar__4(74:79) .EQ. '') Rmin(i__) = rNULL__
if(ar__4(81:85) .EQ. '') RTR(i__) = rNULL__
if(ar__4(87:91) .EQ. '') Iamp(i__) = rNULL__
if(ar__4(93:98) .EQ. '') Imin(i__) = rNULL__
if(ar__4(100:104) .EQ. '') RTI(i__) = rNULL__
if(ar__4(106:110) .EQ. '') Jamp(i__) = rNULL__
if(ar__4(112:117) .EQ. '') Jmin(i__) = rNULL__
if(ar__4(119:123) .EQ. '') RTJ(i__) = rNULL__
if(ar__4(125:129) .EQ. '') Hamp(i__) = rNULL__
if(ar__4(131:136) .EQ. '') Hmin(i__) = rNULL__
if(ar__4(138:142) .EQ. '') RTH(i__) = rNULL__
if(ar__4(144:148) .EQ. '') Ksamp(i__) = rNULL__
if(ar__4(150:155) .EQ. '') Ksmin(i__) = rNULL__
if(ar__4(157:161) .EQ. '') RTKs(i__) = rNULL__
if(ar__4(163:167) .EQ. '') v3_6amp(i__) = rNULL__
if(ar__4(169:174) .EQ. '') v3_6min(i__) = rNULL__
if(ar__4(176:180) .EQ. '') RT3_6(i__) = rNULL__
if(ar__4(182:186) .EQ. '') v4_5amp(i__) = rNULL__
if(ar__4(188:193) .EQ. '') v4_5min(i__) = rNULL__
if(ar__4(195:199) .EQ. '') RT4_5(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4(i__),Uamp(i__),Umin(i__),RTU(i__),Bamp(i__),Bmin(i__),
+ RTB(i__),Vamp(i__),Vmin(i__),RTV(i__),Ramp(i__),Rmin(i__),
+ RTR(i__),Iamp(i__),Imin(i__),RTI(i__),Jamp(i__),Jmin(i__),
+ RTJ(i__),Hamp(i__),Hmin(i__),RTH(i__),Ksamp(i__),Ksmin(i__),
+ RTKs(i__),v3_6amp(i__),v3_6min(i__),RT3_6(i__),v4_5amp(i__),
+ v4_5min(i__),RT4_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end