Conversion of standardized ReadMe file for
file /./ftp/cats/J/AcA/63/159 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-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. J/AcA/63/159 LMC OGLE-III Shallow Survey variable stars (Ulaczyk+, 2013) *================================================================================ *Variable stars from the OGLE-III Shallow Survey in the Large Magellanic Cloud. * Ulaczyk K., Szymanski M.K., Udalski A., Kubiak M., Pietrzynski G., * Soszynski I., Wyrzykowski L., Poleski R., Gieren W., Walker A.R., * Garcia-varela A. * <Acta Astron., 63, 159-179 (2013)> * =2013AcA....63..159U (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'cep.dat' ! Basic parameters of Cepheids derived with OGLE Shallow Survey photometry (table 1) integer*4 nr__ parameter (nr__=28) ! Number of records character*189 ar__ ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg ! (deg) Right Ascension J2000.0 real*8 DEdeg ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*18 Name ! Name (OGLE-LMC-CEP-NNNN or BRIGHT-LMC-CEP-N * for new objects) or BRIGHT-LMC-RRLYR-N integer*4 RAh ! (h) Right ascension (J2000.0) integer*4 RAm ! (min) Right ascension (J2000.0) real*4 RAs ! (s) Right ascension (J2000.0) character*1 DE_ ! Declination sign (J2000.0) integer*4 DEd ! (deg) Declination (J2000.0) integer*4 DEm ! (arcmin) Declination (J2000.0) integer*4 DEs ! (arcsec) Declination (J2000.0) character*8 Subfield ! Subfield designation (LMCNNN.N) real*8 Per ! (d) Period real*4 Imag ! (mag) I-band mean magnitude real*4 Vmag ! (mag) V-band mean magnitude real*4 V_I ! (mag) V-I colour index real*8 TmaxI ! (d) I-band Time of maximum (HJD-2450000) real*4 AI ! (mag) I-band Amplitude real*4 R21I ! I-band Fourier R21 coefficient real*4 phi21I ! I-band Fourier phi21 coefficient real*4 R31I ! I-band Fourier R31 coefficient real*4 phi31I ! I-band Fourier phi31 coefficient real*8 TmaxV ! (d) V-band Time of maximum (HJD-2450000) real*4 AV ! (mag) V-band Amplitude real*4 R21V ! V-band Fourier R21 coefficient real*4 phi21V ! V-band Fourier phi21 coefficient real*4 R31V ! V-band Fourier R31 coefficient real*4 phi31V ! V-band Fourier phi31 coefficient integer*4 o_Imag ! Number of observations in I-band integer*4 o_Vmag ! Number of observations in V-band character*19 OName ! Other name c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'rrlyr.dat' ! Basic parameters of RR Lyrae integer*4 nr__1 parameter (nr__1=3) ! Number of records character*189 ar__1 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_1 ! (deg) Right Ascension J2000.0 real*8 DEdeg_1 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*18 Name_1 ! Name (OGLE-LMC-CEP-NNNN or BRIGHT-LMC-CEP-N * for new objects) or BRIGHT-LMC-RRLYR-N integer*4 RAh_1 ! (h) Right ascension (J2000.0) integer*4 RAm_1 ! (min) Right ascension (J2000.0) real*4 RAs_1 ! (s) Right ascension (J2000.0) character*1 DE__1 ! Declination sign (J2000.0) integer*4 DEd_1 ! (deg) Declination (J2000.0) integer*4 DEm_1 ! (arcmin) Declination (J2000.0) integer*4 DEs_1 ! (arcsec) Declination (J2000.0) character*8 Subfield_1 ! Subfield designation (LMCNNN.N) real*8 Per_1 ! (d) Period real*4 Imag_1 ! (mag) I-band mean magnitude real*4 Vmag_1 ! (mag) V-band mean magnitude real*4 V_I_1 ! (mag) V-I colour index real*8 TmaxI_1 ! (d) I-band Time of maximum (HJD-2450000) real*4 AI_1 ! (mag) I-band Amplitude real*4 R21I_1 ! I-band Fourier R21 coefficient real*4 phi21I_1 ! I-band Fourier phi21 coefficient real*4 R31I_1 ! I-band Fourier R31 coefficient real*4 phi31I_1 ! I-band Fourier phi31 coefficient real*8 TmaxV_1 ! (d) V-band Time of maximum (HJD-2450000) real*4 AV_1 ! (mag) V-band Amplitude real*4 R21V_1 ! V-band Fourier R21 coefficient real*4 phi21V_1 ! V-band Fourier phi21 coefficient real*4 R31V_1 ! V-band Fourier R31 coefficient real*4 phi31V_1 ! V-band Fourier phi31 coefficient integer*4 o_Imag_1 ! Number of observations in I-band integer*4 o_Vmag_1 ! Number of observations in V-band character*19 OName_1 ! Other name c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'ecl.dat' ! Eclipsing variable stars catalog (table 2) integer*4 nr__2 parameter (nr__2=82) ! Number of records character*134 ar__2 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_2 ! (deg) Right Ascension J2000.0 real*8 DEdeg_2 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*17 Name_2 ! Name (BRIGHT-LMC-ECL-NN) integer*4 RAh_2 ! (h) Right ascension (J2000.0) integer*4 RAm_2 ! (min) Right ascension (J2000.0) real*4 RAs_2 ! (s) Right ascension (J2000.0) character*1 DE__2 ! Declination sign (J2000.0) integer*4 DEd_2 ! (deg) Declination (J2000.0) integer*4 DEm_2 ! (arcmin) Declination (J2000.0) integer*4 DEs_2 ! (arcsec) Declination (J2000.0) character*8 Subfield_2 ! Subfield designation (LMCNNN.N) real*8 Per_2 ! (d) Period real*4 Imag_2 ! (mag) I-band mean magnitude real*4 Vmag_2 ! (mag) V-band mean magnitude real*4 V_I_2 ! (mag) V-I colour index integer*4 o_Imag_2 ! Number of observations in I-band integer*4 o_Vmag_2 ! Number of observations in V-band real*4 e_Imag ! (mag) Dispersion of magnitudes for the I-band real*4 e_Vmag ! (mag) Dispersion of magnitudes for the V-band character*5 SubType ! Subtype (EA, EB, EW, Ell) character*30 OName_2 ! Other name c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'lpv.dat' ! Long-period variable stars catalog (table 3) integer*4 nr__3 parameter (nr__3=110) ! Number of records character*160 ar__3 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_3 ! (deg) Right Ascension J2000.0 real*8 DEdeg_3 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*18 Name_3 ! Name (BRIGHT-LMC-LPV-NNN) integer*4 RAh_3 ! (h) Right ascension (J2000.0) integer*4 RAm_3 ! (min) Right ascension (J2000.0) real*4 RAs_3 ! (s) Right ascension (J2000.0) character*1 DE__3 ! Declination sign (J2000.0) integer*4 DEd_3 ! (deg) Declination (J2000.0) integer*4 DEm_3 ! (arcmin) Declination (J2000.0) integer*4 DEs_3 ! (arcsec) Declination (J2000.0) character*8 Subfield_3 ! Subfield designation (LMCNNN.N) real*8 Per_3 ! (d) Dominant period real*4 Imag_3 ! (mag) I-band mean magnitude real*4 Vmag_3 ! (mag) V-band mean magnitude real*4 V_I_3 ! (mag) V-I colour index integer*4 o_Imag_3 ! Number of observations in I-band integer*4 o_Vmag_3 ! Number of observations in V-band real*4 e_Imag_1 ! (mag) Dispersion uncertainty on Imag real*4 e_Vmag_1 ! (mag) Dispersion uncertainty on Vmag character*4 SubType_1 ! Subtype (MIRA, SRV) real*4 d ! (arcsec) ? Separation distance real*4 Jmag ! (mag) ?=- 2MASS J magnitude real*4 Hmag ! (mag) ?=- 2MASS H magnitude real*4 Kmag ! (mag) ?=- 2MASS K magnitude character*31 OName_3 ! Other name c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'misc.dat' ! Miscellaneous stars integer*4 nr__4 parameter (nr__4=558) ! Number of records character*144 ar__4 ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg_4 ! (deg) Right Ascension J2000.0 real*8 DEdeg_4 ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*19 Name_4 ! Name (BRIGHT-LMC-MISC-NNN) integer*4 RAh_4 ! (h) Right ascension (J2000.0) integer*4 RAm_4 ! (min) Right ascension (J2000.0) real*4 RAs_4 ! (s) Right ascension (J2000.0) character*1 DE__4 ! Declination sign (J2000.0) integer*4 DEd_4 ! (deg) Declination (J2000.0) integer*4 DEm_4 ! (arcmin) Declination (J2000.0) integer*4 DEs_4 ! (arcsec) Declination (J2000.0) character*8 Subfield_4 ! Subfield designation (LMCNNN.N) real*4 Imag_4 ! (mag) I-band mean magnitude real*4 Vmag_4 ! (mag) V-band mean magnitude real*4 V_I_4 ! (mag) V-I colour index integer*4 o_Imag_4 ! Number of observations in I-band integer*4 o_Vmag_4 ! Number of observations in V-band real*4 e_Imag_2 ! (mag) Dispersion uncertainty on Imag real*4 e_Vmag_2 ! (mag) Dispersion uncertainty on Vmag real*4 d_1 ! (arcsec) ? Separation distance real*4 Jmag_1 ! (mag) ?=- 2MASS J magnitude real*4 Hmag_1 ! (mag) ?=- 2MASS H magnitude real*4 Kmag_1 ! (mag) ?=- 2MASS K magnitude character*31 OName_4 ! Other name C============================================================================= C Loading file 'cep.dat' ! Basic parameters of Cepheids derived with OGLE * Shallow Survey photometry (table 1) C Format for file interpretation 1 format( + A18,1X,I1,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A8,2X,F10.6, + 1X,F6.3,1X,F6.3,1X,F6.3,2X,F10.5,1X,F4.2,1X,F5.3,1X,F5.3,1X, + F5.3,1X,F5.3,3X,F10.5,1X,F4.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3, + 1X,I2,1X,I2,1X,A19) C Effective file loading open(unit=1,status='old',file= +'cep.dat') write(6,*) '....Loading file: cep.dat' do i__=1,28 read(1,'(A189)')ar__ read(ar__,1) + Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Subfield,Per,Imag,Vmag,V_I, + TmaxI,AI,R21I,phi21I,R31I,phi31I,TmaxV,AV,R21V,phi21V,R31V, + phi31V,o_Imag,o_Vmag,OName 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. if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg c ..............Just test output........... write(6,1) + Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Subfield,Per,Imag,Vmag,V_I, + TmaxI,AI,R21I,phi21I,R31I,phi31I,TmaxV,AV,R21V,phi21V,R31V, + phi31V,o_Imag,o_Vmag,OName write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'rrlyr.dat' ! Basic parameters of RR Lyrae C Format for file interpretation 2 format( + A18,1X,I1,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A8,2X,F10.6, + 1X,F6.3,1X,F6.3,1X,F6.3,2X,F10.5,1X,F4.2,1X,F5.3,1X,F5.3,1X, + F5.3,1X,F5.3,3X,F10.5,1X,F4.2,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3, + 1X,I2,1X,I2,1X,A19) C Effective file loading open(unit=1,status='old',file= +'rrlyr.dat') write(6,*) '....Loading file: rrlyr.dat' do i__=1,3 read(1,'(A189)')ar__1 read(ar__1,2) + Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Subfield_1, + Per_1,Imag_1,Vmag_1,V_I_1,TmaxI_1,AI_1,R21I_1,phi21I_1,R31I_1, + phi31I_1,TmaxV_1,AV_1,R21V_1,phi21V_1,R31V_1,phi31V_1, + o_Imag_1,o_Vmag_1,OName_1 RAdeg_1 = rNULL__ DEdeg_1 = 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 .GT. -180) RAdeg_1=RAh_1*15. if(RAm_1 .GT. -180) RAdeg_1=RAdeg_1+RAm_1/4. if(RAs_1 .GT. -180) RAdeg_1=RAdeg_1+RAs_1/240. if(DEd_1 .GE. 0) DEdeg_1=DEd_1 if(DEm_1 .GE. 0) DEdeg_1=DEdeg_1+DEm_1/60. if(DEs_1 .GE. 0) DEdeg_1=DEdeg_1+DEs_1/3600. if(DE__1.EQ.'-'.AND.DEdeg_1.GE.0) DEdeg_1=-DEdeg_1 c ..............Just test output........... write(6,2) + Name_1,RAh_1,RAm_1,RAs_1,DE__1,DEd_1,DEm_1,DEs_1,Subfield_1, + Per_1,Imag_1,Vmag_1,V_I_1,TmaxI_1,AI_1,R21I_1,phi21I_1,R31I_1, + phi31I_1,TmaxV_1,AV_1,R21V_1,phi21V_1,R31V_1,phi31V_1, + o_Imag_1,o_Vmag_1,OName_1 write(6,'(6H Pos: 2F8.4)') RAdeg_1,DEdeg_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'ecl.dat' ! Eclipsing variable stars catalog (table 2) C Format for file interpretation 3 format( + A17,1X,I1,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A8,2X,F10.6, + 1X,F6.3,1X,F6.3,1X,F6.3,1X,I2,1X,I2,1X,F5.3,1X,F5.3,1X,A5,1X, + A30) C Effective file loading open(unit=1,status='old',file= +'ecl.dat') write(6,*) '....Loading file: ecl.dat' do i__=1,82 read(1,'(A134)')ar__2 read(ar__2,3) + Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,Subfield_2, + Per_2,Imag_2,Vmag_2,V_I_2,o_Imag_2,o_Vmag_2,e_Imag,e_Vmag, + SubType,OName_2 RAdeg_2 = rNULL__ DEdeg_2 = rNULL__ c Derive coordinates RAdeg_2 and DEdeg_2 from input data c (RAdeg_2 and DEdeg_2 are set to rNULL__ when unknown) if(RAh_2 .GT. -180) RAdeg_2=RAh_2*15. if(RAm_2 .GT. -180) RAdeg_2=RAdeg_2+RAm_2/4. if(RAs_2 .GT. -180) RAdeg_2=RAdeg_2+RAs_2/240. if(DEd_2 .GE. 0) DEdeg_2=DEd_2 if(DEm_2 .GE. 0) DEdeg_2=DEdeg_2+DEm_2/60. if(DEs_2 .GE. 0) DEdeg_2=DEdeg_2+DEs_2/3600. if(DE__2.EQ.'-'.AND.DEdeg_2.GE.0) DEdeg_2=-DEdeg_2 c ..............Just test output........... write(6,3) + Name_2,RAh_2,RAm_2,RAs_2,DE__2,DEd_2,DEm_2,DEs_2,Subfield_2, + Per_2,Imag_2,Vmag_2,V_I_2,o_Imag_2,o_Vmag_2,e_Imag,e_Vmag, + SubType,OName_2 write(6,'(6H Pos: 2F8.4)') RAdeg_2,DEdeg_2 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'lpv.dat' ! Long-period variable stars catalog (table 3) C Format for file interpretation 4 format( + A18,1X,I1,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A8,2X,F10.6, + 1X,F6.3,1X,F6.3,1X,F6.3,1X,I2,1X,I2,1X,F5.3,1X,F5.3,1X,A4,1X, + F4.2,1X,F6.3,1X,F6.3,1X,F6.3,A31) C Effective file loading open(unit=1,status='old',file= +'lpv.dat') write(6,*) '....Loading file: lpv.dat' do i__=1,110 read(1,'(A160)')ar__3 read(ar__3,4) + Name_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3,Subfield_3, + Per_3,Imag_3,Vmag_3,V_I_3,o_Imag_3,o_Vmag_3,e_Imag_1,e_Vmag_1, + SubType_1,d,Jmag,Hmag,Kmag,OName_3 if(ar__3(105:108) .EQ. '') d = rNULL__ if (idig(ar__3(110:115)).EQ.0) Jmag = rNULL__ if (idig(ar__3(117:122)).EQ.0) Hmag = rNULL__ if (idig(ar__3(124:129)).EQ.0) Kmag = rNULL__ RAdeg_3 = rNULL__ DEdeg_3 = rNULL__ c Derive coordinates RAdeg_3 and DEdeg_3 from input data c (RAdeg_3 and DEdeg_3 are set to rNULL__ when unknown) if(RAh_3 .GT. -180) RAdeg_3=RAh_3*15. if(RAm_3 .GT. -180) RAdeg_3=RAdeg_3+RAm_3/4. if(RAs_3 .GT. -180) RAdeg_3=RAdeg_3+RAs_3/240. if(DEd_3 .GE. 0) DEdeg_3=DEd_3 if(DEm_3 .GE. 0) DEdeg_3=DEdeg_3+DEm_3/60. if(DEs_3 .GE. 0) DEdeg_3=DEdeg_3+DEs_3/3600. if(DE__3.EQ.'-'.AND.DEdeg_3.GE.0) DEdeg_3=-DEdeg_3 c ..............Just test output........... write(6,4) + Name_3,RAh_3,RAm_3,RAs_3,DE__3,DEd_3,DEm_3,DEs_3,Subfield_3, + Per_3,Imag_3,Vmag_3,V_I_3,o_Imag_3,o_Vmag_3,e_Imag_1,e_Vmag_1, + SubType_1,d,Jmag,Hmag,Kmag,OName_3 write(6,'(6H Pos: 2F8.4)') RAdeg_3,DEdeg_3 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'misc.dat' ! Miscellaneous stars C Format for file interpretation 5 format( + A19,1X,I1,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A8,1X,F6.3,1X, + F6.3,1X,F6.3,1X,I2,1X,I2,1X,F5.3,1X,F5.3,1X,F4.2,1X,F6.3,1X, + F6.3,1X,F6.3,A31) C Effective file loading open(unit=1,status='old',file= +'misc.dat') write(6,*) '....Loading file: misc.dat' do i__=1,558 read(1,'(A144)')ar__4 read(ar__4,5) + Name_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4,Subfield_4, + Imag_4,Vmag_4,V_I_4,o_Imag_4,o_Vmag_4,e_Imag_2,e_Vmag_2,d_1, + Jmag_1,Hmag_1,Kmag_1,OName_4 if(ar__4(89:92) .EQ. '') d_1 = rNULL__ if (idig(ar__4(94:99)).EQ.0) Jmag_1 = rNULL__ if (idig(ar__4(101:106)).EQ.0) Hmag_1 = rNULL__ if (idig(ar__4(108:113)).EQ.0) Kmag_1 = rNULL__ RAdeg_4 = rNULL__ DEdeg_4 = rNULL__ c Derive coordinates RAdeg_4 and DEdeg_4 from input data c (RAdeg_4 and DEdeg_4 are set to rNULL__ when unknown) if(RAh_4 .GT. -180) RAdeg_4=RAh_4*15. if(RAm_4 .GT. -180) RAdeg_4=RAdeg_4+RAm_4/4. if(RAs_4 .GT. -180) RAdeg_4=RAdeg_4+RAs_4/240. if(DEd_4 .GE. 0) DEdeg_4=DEd_4 if(DEm_4 .GE. 0) DEdeg_4=DEdeg_4+DEm_4/60. if(DEs_4 .GE. 0) DEdeg_4=DEdeg_4+DEs_4/3600. if(DE__4.EQ.'-'.AND.DEdeg_4.GE.0) DEdeg_4=-DEdeg_4 c ..............Just test output........... write(6,5) + Name_4,RAh_4,RAm_4,RAs_4,DE__4,DEd_4,DEm_4,DEs_4,Subfield_4, + Imag_4,Vmag_4,V_I_4,o_Imag_4,o_Vmag_4,e_Imag_2,e_Vmag_2,d_1, + Jmag_1,Hmag_1,Kmag_1,OName_4 write(6,'(6H Pos: 2F8.4)') RAdeg_4,DEdeg_4 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