Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/590/A116 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 2026-Apr-23
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/590/A116 Nainital-Cape Survey. IV. (Joshi+, 2016)
*================================================================================
*The Nainital-Cape Survey.
*IV. A Search for pulsational variability in 108 chemically peculiar stars.
* Joshi S., Martinez P., Chowdhury S., Chakradhari N.K., Joshi Y.C.,
* van Heerden P., Medupe T., Kumar Y. B., Kuhn R.B.
* <Astron. Astrophys. 590, A116 (2016)>
* =2016A&A...590A.116J
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The known roAp stars and their physical parameters
integer*4 nr__
parameter (nr__=61) ! Number of records
character*109 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 Seq (nr__) ! Serial Number of a star in the table
integer*4 HD (nr__) ! ? HD no. of the star
character*12 OName (nr__) ! Other (popular) name of the star
integer*4 RAh (nr__) ! (h) Right ascension (J2000)
integer*4 RAm (nr__) ! (min) Right ascension (J2000)
integer*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)
character*13 SpType (nr__) ! MK spectral type
real*4 b_y (nr__) ! (mag) ?=- b-y colour index
real*4 m1 (nr__) ! (mag) ?=- m1 colour index
real*4 c1 (nr__) ! (mag) ?=- c1 colour index
real*4 Dm1 (nr__) ! (mag) ?=- Delta m1 spectral index
real*4 Dc1 (nr__) ! (mag) ?=- Delta c1 spectral index
real*4 beta (nr__) ! (mag) ?=- H{beta} index
integer*4 Teff (nr__) ! (K) Effective temperature of the star
character*8 Ref (nr__) ! Reference detailed in refs.dat file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Additional parameters for the known roAp stars
integer*4 nr__1
parameter (nr__1=61) ! Number of records
character*98 ar__1 ! Full-size record
integer*4 Seq_1 (nr__1) ! Serial Number of a star in the table
character*12 Name (nr__1) ! Name of the star
real*4 Vmag (nr__1) ! (mag) Visual magnitude of the star
real*4 plx (nr__1) ! (mas) ?=- Parallax of the star
real*4 Dist (nr__1) ! (pc) ?=- Distance of the star
real*4 VMAG_1 (nr__1) ! (mag) ?=- Absolute magnitude of the star
real*4 logLstar (nr__1) ! ([Lsun]) ?=- Luminosity of the star relative to the Sun
real*4 Ppul (nr__1) ! (min) Pulsation period of the star
integer*4 Deltamu (nr__1) ! (muHz) ?=- Frequency separation in micro-hertz units
real*4 Amax (nr__1) ! (mmag) ?=- Maximum photometric amplitude variation
real*4 RVmax (nr__1) ! (km/s) ?=- Maximum radial velocity variation
character*1 l_Prot (nr__1) ! Limit flag on Prot
real*4 Prot (nr__1) ! (d) ?=- Rotation Period of the star
real*4 logg (nr__1) ! ([cm/s2]) ?=- Surface gravity of the star
real*4 Mstar (nr__1) ! (Msun) ?=- Mass of the star relative to Sun
real*4 Rstar (nr__1) ! (Rsun) ?=- Radius of the star relative to Sun
real*4 B (nr__1) ! (kgauss) ?=- Mean longitudinal magnetic field of
* the star
real*4 vsini (nr__1) ! (km/s) ?=- Projected rotational velocity of the star
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Stars classified as null results during the
Nainital-Cape survey (108 stars)
integer*4 nr__2
parameter (nr__2=108) ! Number of records
character*251 ar__2 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__2) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__2) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
integer*4 Seq_2 (nr__2) ! Serial Number of a star in the table
integer*4 HD_1 (nr__2) ! HD no. of the star
integer*4 RAh_1 (nr__2) ! (h) Right ascension (J2000)
integer*4 RAm_1 (nr__2) ! (min) Right ascension (J2000)
integer*4 RAs_1 (nr__2) ! (s) Right ascension (J2000)
character*1 DE__1 (nr__2) ! Declination sign (J2000)
integer*4 DEd_1 (nr__2) ! (deg) Declination (J2000)
integer*4 DEm_1 (nr__2) ! (arcmin) Declination (J2000)
integer*4 DEs_1 (nr__2) ! (arcsec) Declination (J2000)
real*4 Vmag_2 (nr__2) ! (mag) Visual magnitude of the star
character*7 SpType_1 (nr__2) ! MK spectral type
real*4 plx_1 (nr__2) ! (mas) ?=- Parallax of the star
real*4 e_plx (nr__2) ! (mas) ?=- Error in parallax
real*4 b_y_1 (nr__2) ! (mag) ?=- b-y colour index
real*4 m1_1 (nr__2) ! (mag) ?=- m1 colour index
real*4 c1_1 (nr__2) ! (mag) ?=- c1 colour index
real*4 beta_1 (nr__2) ! (mag) ?=- H Beta index
real*4 Dm1_1 (nr__2) ! (mag) ?=- Delta m1 spectral index
real*4 Dc1_1 (nr__2) ! (mag) ?=- Delta c1 spectral index
integer*4 Teff_1 (nr__2) ! (K) ?=- Effective temperature of the star
real*4 logLstar_1 (nr__2) ! ([Lsun]) ?=- Luminosity of the star relative to the Sun
real*4 Deltat1 (nr__2) ! (h) Duration of observation 1
integer*4 HJD1 (nr__2) ! (d) Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year1 (nr__2) ! Year when the star was observed
real*4 Deltat2 (nr__2) ! (h) ? Duration of observation 2
integer*4 HJD2 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year2 (nr__2) ! ? Year when the star was observed
real*4 Deltat3 (nr__2) ! (h) ? Duration of observation 3
integer*4 HJD3 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year3 (nr__2) ! ? Year when the star was observed
real*4 Deltat4 (nr__2) ! (h) ? Duration of observation 4
integer*4 HJD4 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year4 (nr__2) ! ? Year when the star was observed
real*4 Deltat5 (nr__2) ! (h) ? Duration of observation 5
integer*4 HJD5 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year5 (nr__2) ! ? Year when the star was observed
real*4 Deltat6 (nr__2) ! (h) ? Duration of observation 6
integer*4 HJD6 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year6 (nr__2) ! ? Year when the star was observed
real*4 Deltat7 (nr__2) ! (h) ? Duration of observation 7
integer*4 HJD7 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year7 (nr__2) ! ? Year when the star was observed
real*4 Deltat8 (nr__2) ! (h) ? Duration of observation 8
integer*4 HJD8 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year8 (nr__2) ! ? Year when the star was observed
real*4 Deltat9 (nr__2) ! (h) ? Duration of observation 9
integer*4 HJD9 (nr__2) ! (d) ? Julian date (2450000+) when the star was
* observed
integer*4 Obs_Year9 (nr__2) ! ? Year when the star was observed
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References used in table1.dat
integer*4 nr__3
parameter (nr__3=54) ! Number of records
character*80 ar__3 ! Full-size record
integer*4 Ref_1 (nr__3) ! Reference number
character*19 BibCode (nr__3) ! BibCode
character*23 Aut (nr__3) ! Author's name
character*32 Com (nr__3) ! Comments
C=============================================================================
C Loading file 'table1.dat' ! The known roAp stars and their physical parameters
C Format for file interpretation
1 format(
+ I2,3X,I6,1X,A12,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,A13,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,I4,1X,A8)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,61
read(1,'(A109)')ar__
read(ar__,1)
+ Seq(i__),HD(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),SpType(i__),b_y(i__),
+ m1(i__),c1(i__),Dm1(i__),Dc1(i__),beta(i__),Teff(i__),Ref(i__)
if(ar__(6:11) .EQ. '') HD(i__) = iNULL__
if (idig(ar__(59:63)).EQ.0) b_y(i__) = rNULL__
if (idig(ar__(65:69)).EQ.0) m1(i__) = rNULL__
if (idig(ar__(71:75)).EQ.0) c1(i__) = rNULL__
if (idig(ar__(77:82)).EQ.0) Dm1(i__) = rNULL__
if (idig(ar__(84:89)).EQ.0) Dc1(i__) = rNULL__
if (idig(ar__(91:95)).EQ.0) beta(i__) = rNULL__
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)
+ Seq(i__),HD(i__),OName(i__),RAh(i__),RAm(i__),RAs(i__),
+ DE_(i__),DEd(i__),DEm(i__),DEs(i__),SpType(i__),b_y(i__),
+ m1(i__),c1(i__),Dm1(i__),Dc1(i__),beta(i__),Teff(i__),Ref(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 'table2.dat' ! Additional parameters for the known roAp stars
C Format for file interpretation
2 format(
+ I2,1X,A12,1X,F5.2,1X,F5.2,1X,F6.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ I2,1X,F5.2,1X,F4.2,A1,F6.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2,1X,
+ F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,61
read(1,'(A98)')ar__1
read(ar__1,2)
+ Seq_1(i__),Name(i__),Vmag(i__),plx(i__),Dist(i__),VMAG_1(i__),
+ logLstar(i__),Ppul(i__),Deltamu(i__),Amax(i__),RVmax(i__),
+ l_Prot(i__),Prot(i__),logg(i__),Mstar(i__),Rstar(i__),B(i__),
+ vsini(i__)
if (idig(ar__1(23:27)).EQ.0) plx(i__) = rNULL__
if (idig(ar__1(29:34)).EQ.0) Dist(i__) = rNULL__
if (idig(ar__1(36:39)).EQ.0) VMAG_1(i__) = rNULL__
if (idig(ar__1(41:44)).EQ.0) logLstar(i__) = rNULL__
if (Deltamu(i__) .EQ. 45) Deltamu(i__) = iNULL__
if (idig(ar__1(55:59)).EQ.0) Amax(i__) = rNULL__
if (idig(ar__1(61:64)).EQ.0) RVmax(i__) = rNULL__
if (idig(ar__1(66:71)).EQ.0) Prot(i__) = rNULL__
if (idig(ar__1(73:76)).EQ.0) logg(i__) = rNULL__
if (idig(ar__1(78:81)).EQ.0) Mstar(i__) = rNULL__
if (idig(ar__1(83:86)).EQ.0) Rstar(i__) = rNULL__
if (idig(ar__1(88:92)).EQ.0) B(i__) = rNULL__
if (idig(ar__1(94:98)).EQ.0) vsini(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Seq_1(i__),Name(i__),Vmag(i__),plx(i__),Dist(i__),VMAG_1(i__),
+ logLstar(i__),Ppul(i__),Deltamu(i__),Amax(i__),RVmax(i__),
+ l_Prot(i__),Prot(i__),logg(i__),Mstar(i__),Rstar(i__),B(i__),
+ vsini(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Stars classified as null results during the
* Nainital-Cape survey (108 stars)
C Format for file interpretation
3 format(
+ I3,1X,I6,1X,I2,1X,I2,1X,I2,1X,A1,I2,1X,I2,1X,I2,1X,F5.2,2X,A7,
+ F5.2,1X,F4.2,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F6.3,
+ 1X,I6,1X,F5.2,2X,F4.2,1X,I4,1X,I4,2X,F4.2,1X,I4,1X,I4,2X,F4.2,
+ 1X,I4,1X,I4,2X,F4.2,1X,I4,1X,I4,2X,F4.2,1X,I4,1X,I4,2X,F4.2,
+ 1X,I4,1X,I4,2X,F4.2,1X,I4,1X,I4,2X,F4.2,1X,I4,1X,I4,2X,F4.2,
+ 1X,I4,1X,I4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,108
read(1,'(A251)')ar__2
read(ar__2,3)
+ Seq_2(i__),HD_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Vmag_2(i__),
+ SpType_1(i__),plx_1(i__),e_plx(i__),b_y_1(i__),m1_1(i__),
+ c1_1(i__),beta_1(i__),Dm1_1(i__),Dc1_1(i__),Teff_1(i__),
+ logLstar_1(i__),Deltat1(i__),HJD1(i__),Obs_Year1(i__),
+ Deltat2(i__),HJD2(i__),Obs_Year2(i__),Deltat3(i__),HJD3(i__),
+ Obs_Year3(i__),Deltat4(i__),HJD4(i__),Obs_Year4(i__),
+ Deltat5(i__),HJD5(i__),Obs_Year5(i__),Deltat6(i__),HJD6(i__),
+ Obs_Year6(i__),Deltat7(i__),HJD7(i__),Obs_Year7(i__),
+ Deltat8(i__),HJD8(i__),Obs_Year8(i__),Deltat9(i__),HJD9(i__),
+ Obs_Year9(i__)
if (idig(ar__2(45:49)).EQ.0) plx_1(i__) = rNULL__
if (idig(ar__2(51:54)).EQ.0) e_plx(i__) = rNULL__
if (idig(ar__2(56:61)).EQ.0) b_y_1(i__) = rNULL__
if (idig(ar__2(63:68)).EQ.0) m1_1(i__) = rNULL__
if (idig(ar__2(70:74)).EQ.0) c1_1(i__) = rNULL__
if (idig(ar__2(76:80)).EQ.0) beta_1(i__) = rNULL__
if (idig(ar__2(82:87)).EQ.0) Dm1_1(i__) = rNULL__
if (idig(ar__2(89:94)).EQ.0) Dc1_1(i__) = rNULL__
if (Teff_1(i__) .EQ. 45) Teff_1(i__) = iNULL__
if (idig(ar__2(103:107)).EQ.0) logLstar_1(i__) = rNULL__
if(ar__2(126:129) .EQ. '') Deltat2(i__) = rNULL__
if(ar__2(131:134) .EQ. '') HJD2(i__) = iNULL__
if(ar__2(136:139) .EQ. '') Obs_Year2(i__) = iNULL__
if(ar__2(142:145) .EQ. '') Deltat3(i__) = rNULL__
if(ar__2(147:150) .EQ. '') HJD3(i__) = iNULL__
if(ar__2(152:155) .EQ. '') Obs_Year3(i__) = iNULL__
if(ar__2(158:161) .EQ. '') Deltat4(i__) = rNULL__
if(ar__2(163:166) .EQ. '') HJD4(i__) = iNULL__
if(ar__2(168:171) .EQ. '') Obs_Year4(i__) = iNULL__
if(ar__2(174:177) .EQ. '') Deltat5(i__) = rNULL__
if(ar__2(179:182) .EQ. '') HJD5(i__) = iNULL__
if(ar__2(184:187) .EQ. '') Obs_Year5(i__) = iNULL__
if(ar__2(190:193) .EQ. '') Deltat6(i__) = rNULL__
if(ar__2(195:198) .EQ. '') HJD6(i__) = iNULL__
if(ar__2(200:203) .EQ. '') Obs_Year6(i__) = iNULL__
if(ar__2(206:209) .EQ. '') Deltat7(i__) = rNULL__
if(ar__2(211:214) .EQ. '') HJD7(i__) = iNULL__
if(ar__2(216:219) .EQ. '') Obs_Year7(i__) = iNULL__
if(ar__2(222:225) .EQ. '') Deltat8(i__) = rNULL__
if(ar__2(227:230) .EQ. '') HJD8(i__) = iNULL__
if(ar__2(232:235) .EQ. '') Obs_Year8(i__) = iNULL__
if(ar__2(238:241) .EQ. '') Deltat9(i__) = rNULL__
if(ar__2(243:246) .EQ. '') HJD9(i__) = iNULL__
if(ar__2(248:251) .EQ. '') Obs_Year9(i__) = iNULL__
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,3)
+ Seq_2(i__),HD_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Vmag_2(i__),
+ SpType_1(i__),plx_1(i__),e_plx(i__),b_y_1(i__),m1_1(i__),
+ c1_1(i__),beta_1(i__),Dm1_1(i__),Dc1_1(i__),Teff_1(i__),
+ logLstar_1(i__),Deltat1(i__),HJD1(i__),Obs_Year1(i__),
+ Deltat2(i__),HJD2(i__),Obs_Year2(i__),Deltat3(i__),HJD3(i__),
+ Obs_Year3(i__),Deltat4(i__),HJD4(i__),Obs_Year4(i__),
+ Deltat5(i__),HJD5(i__),Obs_Year5(i__),Deltat6(i__),HJD6(i__),
+ Obs_Year6(i__),Deltat7(i__),HJD7(i__),Obs_Year7(i__),
+ Deltat8(i__),HJD8(i__),Obs_Year8(i__),Deltat9(i__),HJD9(i__),
+ Obs_Year9(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 'refs.dat' ! References used in table1.dat
C Format for file interpretation
4 format(I2,1X,A19,1X,A23,2X,A32)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,54
read(1,'(A80)')ar__3
read(ar__3,4)Ref_1(i__),BibCode(i__),Aut(i__),Com(i__)
c ..............Just test output...........
write(6,4)Ref_1(i__),BibCode(i__),Aut(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end