FORTRAN Generation
(/./ftp/cats/J/A_A/604/L6)

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/604/L6 EBLM J0555-57 photometry and RV (von Boetticher+, 2017) *================================================================================ *The EBLM Project. *III. A Saturn-size low-mass star at the hydrogen-burning limit. * von Boetticher A., Triaud A.H.M.J, Queloz D., Gill S., Lendl M., Delrez L., * Anderson D.R., Collier Cameron A., Faedi F., Gillon M., * Gomez Maqueo Chew Y., Hebb L., Hellier C., Jehin E., Maxted P.F.L., * Martin D.V., Pepe F., Pollacco D., Segransan D., Smalley B., Udry S., * West R. * <Astron. Astrophys. 604, L6 (2017)> * =2017A&A...604L...6V (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'coralie1.dat' ! CORALIE radial velocities I integer*4 nr__ parameter (nr__=20) ! Number of records character*86 ar__ ! Full-size record real*8 BJD ! (d) Barycentric Julian date (BJD-2400000) real*8 RV ! (km/s) Radial velocity real*8 e_RV ! (km/s) 1 sigma photon noise error bars on RV real*8 FWHM ! (km/s) Full-width at half-maximum of * cross-correlation function real*4 Cont ! Contrast of cross-correlation function real*8 BisSpan ! (km/s) Span of the bisector inverse slope real*4 pclean ! ?=- Probability that spectrum is NOT * contaminated by blend integer*4 Rej ! [0/1] Spectrum rejected: * 1: rejected, 0: accepted (if pclean<0.01) character*23 Com ! Comments c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'coralie2.dat' ! CORALIE radial velocities II integer*4 nr__1 parameter (nr__1=11) ! Number of records character*86 ar__1 ! Full-size record real*8 BJD_1 ! (d) Barycentric Julian date (BJD-2400000) real*8 RV_1 ! (km/s) Radial velocity real*8 e_RV_1 ! (km/s) 1 sigma photon noise error bars on RV real*8 FWHM_1 ! (km/s) Full-width at half-maximum of * cross-correlation function real*4 Cont_1 ! Contrast of cross-correlation function real*8 BisSpan_1 ! (km/s) Span of the bisector inverse slope real*4 pclean_1 ! ?=- Probability that spectrum is NOT * contaminated by blend integer*4 Rej_1 ! [0/1] Spectrum rejected: * 1: rejected, 0: accepted (if pclean<0.01) character*23 Com_1 ! Comments c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'eulerl.dat' ! Euler z band photometry, large aperture including A and B components integer*4 nr__2 parameter (nr__2=526) ! Number of records character*64 ar__2 ! Full-size record real*8 HJD ! (d) Heliocentric Julian Date (HJD-2450000) UTC real*8 Flux ! Normalised flux real*8 e_Flux ! 1 sigma error bar on the flux real*4 xshift ! (pix) Centroid shift on detector real*4 yshift ! (pix) Centroid shift on detector real*4 AM ! Airmass real*4 FWHM_2 ! (pix) Full-width at half-maximum of centroid real*4 Sky ! Sky background flux real*4 ExpTime ! (s) Exposure times c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'eulers.dat' ! Euler z band photometry, small aperture on A component integer*4 nr__3 parameter (nr__3=526) ! Number of records character*64 ar__3 ! Full-size record real*8 HJD_1 ! (d) Heliocentric Julian Date (HJD-2450000) UTC real*8 Flux_1 ! Normalised flux real*8 e_Flux_1 ! 1 sigma error bar on the flux real*4 xshift_1 ! (pix) Centroid shift on detector real*4 yshift_1 ! (pix) Centroid shift on detector real*4 AM_1 ! Airmass real*4 FWHM_3 ! (pix) Full-width at half-maximum of centroid real*4 Sky_1 ! Sky background flux real*4 ExpTime_1 ! (s) Exposure times c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'trappist.dat' ! TRAPPIST z band photometry, aperture including A and B components integer*4 nr__4 parameter (nr__4=1188) ! Number of records character*230 ar__4 ! Full-size record real*8 BJD_2 ! (d) Barycentric Julian Date (BJD-2450000) UTC real*8 Flux_2 ! Normalised flux real*4 e_Flux_2 ! 1 sigma error bar on the flux real*4 xshift_2 ! (pix) Centroid shift on detector real*4 yshift_2 ! (pix) Centroid shift on detector real*8 FWHMx ! (pix) Centroid FWHM in x real*8 FWHMy ! (pix) Centroid FWHM in y real*8 FWHM_4 ! (pix) Centroid FWHM real*8 Sky_2 ! Sky background flux real*8 AM_2 ! Airmass real*8 ExpTime_2 ! (s) Exposure time C============================================================================= C Loading file 'coralie1.dat' ! CORALIE radial velocities I C Format for file interpretation 1 format( + F12.6,1X,F8.5,1X,F7.5,1X,F8.5,1X,F6.3,1X,F8.5,1X,F5.3,1X,I1, + 1X,A23) C Effective file loading open(unit=1,status='old',file= +'coralie1.dat') write(6,*) '....Loading file: coralie1.dat' do i__=1,20 read(1,'(A86)')ar__ read(ar__,1)BJD,RV,e_RV,FWHM,Cont,BisSpan,pclean,Rej,Com if (idig(ar__(56:60)).EQ.0) pclean = rNULL__ c ..............Just test output........... write(6,1)BJD,RV,e_RV,FWHM,Cont,BisSpan,pclean,Rej,Com c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'coralie2.dat' ! CORALIE radial velocities II C Format for file interpretation 2 format( + F12.6,1X,F8.5,1X,F7.5,1X,F8.5,1X,F6.3,1X,F8.5,1X,F5.3,1X,I1, + 1X,A23) C Effective file loading open(unit=1,status='old',file= +'coralie2.dat') write(6,*) '....Loading file: coralie2.dat' do i__=1,11 read(1,'(A86)')ar__1 read(ar__1,2) + BJD_1,RV_1,e_RV_1,FWHM_1,Cont_1,BisSpan_1,pclean_1,Rej_1,Com_1 if (idig(ar__1(56:60)).EQ.0) pclean_1 = rNULL__ c ..............Just test output........... write(6,2) + BJD_1,RV_1,e_RV_1,FWHM_1,Cont_1,BisSpan_1,pclean_1,Rej_1,Com_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'eulerl.dat' ! Euler z band photometry, large aperture * including A and B components C Format for file interpretation 3 format( + F11.6,1X,F8.6,1X,F8.6,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2, + 1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'eulerl.dat') write(6,*) '....Loading file: eulerl.dat' do i__=1,526 read(1,'(A64)')ar__2 read(ar__2,3) + HJD,Flux,e_Flux,xshift,yshift,AM,FWHM_2,Sky,ExpTime c ..............Just test output........... write(6,3) + HJD,Flux,e_Flux,xshift,yshift,AM,FWHM_2,Sky,ExpTime c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'eulers.dat' ! Euler z band photometry, small aperture * on A component C Format for file interpretation 4 format( + F11.6,1X,F8.6,1X,F8.6,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F5.2, + 1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'eulers.dat') write(6,*) '....Loading file: eulers.dat' do i__=1,526 read(1,'(A64)')ar__3 read(ar__3,4) + HJD_1,Flux_1,e_Flux_1,xshift_1,yshift_1,AM_1,FWHM_3,Sky_1, + ExpTime_1 c ..............Just test output........... write(6,4) + HJD_1,Flux_1,e_Flux_1,xshift_1,yshift_1,AM_1,FWHM_3,Sky_1, + ExpTime_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'trappist.dat' ! TRAPPIST z band photometry, aperture including * A and B components C Format for file interpretation 5 format( + F18.13,1X,F19.17,1X,E23.18,1X,E24.17,1X,E26.19,1X,F18.16,1X, + F19.16,1X,F18.16,1X,F19.15,1X,F18.16,1X,F18.16) C Effective file loading open(unit=1,status='old',file= +'trappist.dat') write(6,*) '....Loading file: trappist.dat' do i__=1,1188 read(1,'(A230)')ar__4 read(ar__4,5) + BJD_2,Flux_2,e_Flux_2,xshift_2,yshift_2,FWHMx,FWHMy,FWHM_4, + Sky_2,AM_2,ExpTime_2 c ..............Just test output........... write(6,5) + BJD_2,Flux_2,e_Flux_2,xshift_2,yshift_2,FWHMx,FWHMy,FWHM_4, + Sky_2,AM_2,ExpTime_2 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