Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/636/240 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-Jun-10
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/ApJ/636/240 OGLE II microlensing parameters (Sumi+, 2006)
*================================================================================
*Microlensing optical depth toward the Galactic Bulge using bright sources from
*OGLE-II.
* Sumi T., Wozniak P.R., Udalski A., Szymanski M., Kubiak M., Pietrzynski G.,
* Soszynski I., Zebrun K., Szewczyk O., Wyrzykowski L., Paczynski B.
* <Astrophys. J., 636, 240-260 (2006)>
* =2006ApJ...636..240S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! All 122 candidates with microlensing parameters
fitted with blending
integer*4 nr__
parameter (nr__=122) ! Number of records
character*457 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)
character*11 OGLE (nr__) ! Field number-Identification number
* (SCNN-NNNNNN)
integer*4 Level (nr__) ! Selection criteria level (1)
integer*4 RAh (nr__) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm (nr__) ! (min) Minute of Right Ascension (J2000)
real*4 RAs (nr__) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__) ! Sign of the Declination (J2000)
integer*4 DEd (nr__) ! (deg) Degree of Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000)
real*8 Xpos (nr__) ! (pix) X pixel position
real*8 Ypos (nr__) ! (pix) Y pixel position
real*4 Imag (nr__) ! (mag) The I band magnitude
real*4 v_V_I (nr__) ! (mag) The (V-I) color index
real*4 I0mag (nr__) ! (mag) Extinction corrected I band magnitude
real*4 v_V_I_0 (nr__) ! (mag) Extinction corrected (V-I) color index
real*8 t0 (nr__) ! (d) Time of maximum magnification (HJD-2450000.0)
real*8 tE (nr__) ! (d) Event time-scale
real*8 umin (nr__) ! Minimum impact parameter in units of the
* Einstein radius R_E_
real*8 Amax (nr__) ! Maximum magnification
real*8 Fs (nr__) ! Baseline source star flux fraction in ADUs
real*8 Fb (nr__) ! Baseline blended background star flux in ADUs
real*4 Is (nr__) ! (mag) Best fit source I band magnitude
real*8 chi2 (nr__) ! {Chi}^2^ of the the model fit
integer*4 Ndata (nr__) ! Number of data points in light curve
real*8 redchi2 (nr__) ! Reduced {Chi}^2^ of the the model fit
real*8 sigmax (nr__) ! Maximum Significance of the data points.
real*8 sumsig (nr__) ! Integrated significance of the peak
real*8 e_t0 (nr__) ! (d) 1{sigma} symmetric error of t0
real*8 e1t0 (nr__) ! (d) 1{sigma} asymmetric lower error of t0
real*8 e2t0 (nr__) ! (d) 1{sigma} asymmetric upper error of t0
real*8 e_tE (nr__) ! (d) 1{sigma} symmetric error of tE
real*8 e1tE (nr__) ! (d) 1{sigma} asymmetric lower error of tE
real*8 e2tE (nr__) ! (d) 1{sigma} asymmetric upper error of tE
real*8 e_umin (nr__) ! 1{sigma} symmetric error of umin
real*8 e1umin (nr__) ! 1{sigma} asymmetric lower error of umin
real*8 e2umin (nr__) ! 1{sigma} asymmetric upper error of umin
real*8 e_Fs (nr__) ! 1{sigma} symmetric error of Fs in ADUs
real*8 e1Fs (nr__) ! 1{sigma} asymmetric lower error of Fs
real*8 e2Fs (nr__) ! 1{sigma} asymmetric upper error of Fs
real*8 e_Fb (nr__) ! 1{sigma} symmetric error of Fb in ADUs
real*8 e1Fb (nr__) ! 1{sigma} asymmetric lower error of Fb
real*8 e2Fb (nr__) ! 1{sigma} asymmetric upper error of Fb
real*4 e_Is (nr__) ! (mag) 1{sigma} asymmetric lower error of Is
real*4 E_Is_1 (nr__) ! (mag) 1{sigma} asymmetric upper error of Is
integer*4 Ndatap (nr__) ! Number of data points in the peak within 5tE
real*8 redchi2p (nr__) ! Reduced {Chi}^2^ of the model fit in the peak
real*4 Is0mag (nr__) ! (mag) Extinction corrected I band magnitude of
* source
real*4 e_Is0mag (nr__) ! (mag) 1{sigma} asymmetric lower error of Is0
real*4 E_Is0mag_1 (nr__) ! (mag) 1{sigma} asymmetric upper error of Is0
real*8 fS_1 (nr__) ! Blending fraction Fs/(Fs+Fb)
real*8 e1fS_1 (nr__) ! 1{sigma} asymmetric lower error of fS
real*4 e2fS_1 (nr__) ! 1{sigma} asymmetric upper error of fS
*Note (1): Level from 0 (minimum number of data points) to 6 (serious
* microlensing event), see Table 2 of paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! All 122 candidates with microlensing parameters
fitted without blending
integer*4 nr__1
parameter (nr__1=122) ! Number of records
character*457 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg_1 (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*11 OGLE_1 (nr__1) ! Field number-Identification number
* (SCNN-NNNNNN)
integer*4 Level_1 (nr__1) ! Selection criteria level (1)
integer*4 RAh_1 (nr__1) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm_1 (nr__1) ! (min) Minute of Right Ascension (J2000)
real*4 RAs_1 (nr__1) ! (s) Second of Right Ascension (J2000)
character*1 DE__1 (nr__1) ! Sign of the Declination (J2000)
integer*4 DEd_1 (nr__1) ! (deg) Degree of Declination (J2000)
integer*4 DEm_1 (nr__1) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs_1 (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
real*8 Xpos_1 (nr__1) ! (pix) X pixel position
real*8 Ypos_1 (nr__1) ! (pix) Y pixel position
real*4 Imag_1 (nr__1) ! (mag) The I band magnitude
real*4 v_V_I_1 (nr__1) ! (mag) The (V-I) color index
real*4 I0mag_1 (nr__1) ! (mag) Extinction corrected I band magnitude
real*4 v_V_I_0_1 (nr__1) ! (mag) Extinction corrected (V-I) color index
real*8 t0_1 (nr__1) ! (d) Time of maximum magnification (HJD-2450000.0)
real*8 tE_1 (nr__1) ! (d) Event time-scale
real*8 umin_1 (nr__1) ! Minimum impact parameter in units of the
* Einstein radius R_E_
real*8 Amax_1 (nr__1) ! Maximum magnification
real*8 Fs_2 (nr__1) ! Baseline source star flux fraction in ADUs
real*8 Fb_1 (nr__1) ! Baseline blended background star flux in ADUs
real*4 Is_1 (nr__1) ! (mag) Best fit source I band magnitude
real*8 chi2_1 (nr__1) ! {Chi}^2^ of the the model fit
integer*4 Ndata_1 (nr__1) ! Number of data points in light curve
real*8 redchi2_1 (nr__1) ! Reduced {Chi}^2^ of the the model fit
real*8 sigmax_1 (nr__1) ! Maximum Significance of the data points.
real*8 sumsig_1 (nr__1) ! Integrated significance of the peak
real*8 e_t0_1 (nr__1) ! (d) 1{sigma} symmetric error of t0
real*8 e1t0_1 (nr__1) ! (d) 1{sigma} asymmetric lower error of t0
real*8 e2t0_1 (nr__1) ! (d) 1{sigma} asymmetric upper error of t0
real*8 e_tE_1 (nr__1) ! (d) 1{sigma} symmetric error of tE
real*8 e1tE_1 (nr__1) ! (d) 1{sigma} asymmetric lower error of tE
real*8 e2tE_1 (nr__1) ! (d) 1{sigma} asymmetric upper error of tE
real*8 e_umin_1 (nr__1) ! 1{sigma} symmetric error of umin
real*8 e1umin_1 (nr__1) ! 1{sigma} asymmetric lower error of umin
real*8 e2umin_1 (nr__1) ! 1{sigma} asymmetric upper error of umin
real*8 e_Fs_1 (nr__1) ! 1{sigma} symmetric error of Fs in ADUs
real*8 e1Fs_2 (nr__1) ! 1{sigma} asymmetric lower error of Fs
real*8 e2Fs_2 (nr__1) ! 1{sigma} asymmetric upper error of Fs
real*8 e_Fb_1 (nr__1) ! 1{sigma} symmetric error of Fb in ADUs
real*8 e1Fb_1 (nr__1) ! 1{sigma} asymmetric lower error of Fb
real*8 e2Fb_1 (nr__1) ! 1{sigma} asymmetric upper error of Fb
real*4 e_Is_2 (nr__1) ! (mag) 1{sigma} asymmetric lower error of Is
real*4 E_Is_3 (nr__1) ! (mag) 1{sigma} asymmetric upper error of Is
integer*4 Ndatap_1 (nr__1) ! Number of data points in the peak within 5tE
real*8 redchi2p_1 (nr__1) ! Reduced {Chi}^2^ of the model fit in the peak
real*4 Is0mag_1 (nr__1) ! (mag) Extinction corrected I band magnitude of
* source
real*4 e_Is0mag_2 (nr__1) ! (mag) 1{sigma} asymmetric lower error of Is0
real*4 E_Is0mag_3 (nr__1) ! (mag) 1{sigma} asymmetric upper error of Is0
real*8 fS_3 (nr__1) ! Blending fraction Fs/(Fs+Fb)
real*8 e1fS_3 (nr__1) ! 1{sigma} asymmetric lower error of fS
real*4 e2fS_3 (nr__1) ! 1{sigma} asymmetric upper error of fS
*Note (1): Level from 0 (minimum number of data points) to 6 (serious
* microlensing event), see Table 2 of paper.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! Cross-identification of candidates with other
sources
integer*4 nr__2
parameter (nr__2=122) ! Number of records
character*123 ar__2 ! Full-size record
character*11 OGLE_2 (nr__2) ! Field number-Identification number (SCNN-NNNNNN)
character*111 OName (nr__2) ! Other source identification(s) (1)
*Note (1): Source ID (with separation in arcsec) where:
* MACHO = Thomas et al. 2004, astro-ph/0410341;
* MOA = Sumi et al., 2003ApJ...591..204S;
* SMITH = Smith, Mao, & Wozniak, 2002MNRAS.332..962S;
* EROS = Afonso, C. et al., 2003A&A...404..145A;
* EROS2 = http://dphs10.saclay.cea.fr/Spp/Experiences/EROS/alertes.html;
* JARO = Jaroszynski, M., 2002AcA....52...39J.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea4/*' ! *The light curve data for all 122 candidates
integer*4 nr__3
parameter (nr__3=122) ! Number of records
character*50 ar__3 ! Full-size record
real*8 JD_2450000 (nr__3) ! (d) Julian Date (minus 2450000)
real*8 Flux (nr__3) ! Observed flux (ADU)
real*8 e_Flux (nr__3) ! Error in Flux (ADU)
real*8 Imag_2 (nr__3) ! (mag) Observed I band magnitude
real*8 e_Imag (nr__3) ! (mag) Error in Imag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea5.dat' ! Detection efficiencies for each field with blending
integer*4 nr__4
parameter (nr__4=540) ! Number of records
character*48 ar__4 ! Full-size record
character*4 Field (nr__4) ! Name of file (SCNN)
integer*4 tEi (nr__4) ! Integer value of 10*log of the event time-scale
real*4 log_tE (nr__4) ! ([d]) Log of the event time-scale (1)
real*8 tE_2 (nr__4) ! (d) Event time-scale
integer*4 Nout (nr__4) ! Number of output artificial events
integer*4 Nin (nr__4) ! Number of input artificial events
real*8 Deff (nr__4) ! Detection efficiency
*Note (1): Indicates the mid point of the bin, i.e., the first bin is
* log(tE)=0-0.1.
C=============================================================================
C Loading file 'tablea1.dat' ! All 122 candidates with microlensing parameters
* fitted with blending
C Format for file interpretation
1 format(
+ A11,1X,I1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F7.2,
+ 1X,F7.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F11.6,1X,F11.5,1X,
+ F10.7,1X,F14.3,1X,F14.6,1X,F15.6,1X,F6.3,1X,F10.4,1X,I3,1X,
+ F7.4,1X,F7.3,1X,F8.3,1X,F9.6,1X,F10.6,1X,F10.6,1X,F12.5,1X,
+ F12.5,1X,F9.5,1X,F9.7,1X,F11.7,1X,F9.7,1X,F12.6,1X,F12.6,1X,
+ F11.6,1X,F12.6,1X,F12.6,1X,F11.6,1X,F6.3,1X,F6.3,1X,I3,1X,
+ F8.4,1X,F6.3,1X,F6.3,1X,F6.3,1X,F10.4,1X,F7.4,1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,122
read(1,'(A457)')ar__
read(ar__,1)
+ OGLE(i__),Level(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Xpos(i__),Ypos(i__),Imag(i__),
+ v_V_I(i__),I0mag(i__),v_V_I_0(i__),t0(i__),tE(i__),umin(i__),
+ Amax(i__),Fs(i__),Fb(i__),Is(i__),chi2(i__),Ndata(i__),
+ redchi2(i__),sigmax(i__),sumsig(i__),e_t0(i__),e1t0(i__),
+ e2t0(i__),e_tE(i__),e1tE(i__),e2tE(i__),e_umin(i__),
+ e1umin(i__),e2umin(i__),e_Fs(i__),e1Fs(i__),e2Fs(i__),
+ e_Fb(i__),e1Fb(i__),e2Fb(i__),e_Is(i__),E_Is_1(i__),
+ Ndatap(i__),redchi2p(i__),Is0mag(i__),e_Is0mag(i__),
+ E_Is0mag_1(i__),fS_1(i__),e1fS_1(i__),e2fS_1(i__)
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)
+ OGLE(i__),Level(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Xpos(i__),Ypos(i__),Imag(i__),
+ v_V_I(i__),I0mag(i__),v_V_I_0(i__),t0(i__),tE(i__),umin(i__),
+ Amax(i__),Fs(i__),Fb(i__),Is(i__),chi2(i__),Ndata(i__),
+ redchi2(i__),sigmax(i__),sumsig(i__),e_t0(i__),e1t0(i__),
+ e2t0(i__),e_tE(i__),e1tE(i__),e2tE(i__),e_umin(i__),
+ e1umin(i__),e2umin(i__),e_Fs(i__),e1Fs(i__),e2Fs(i__),
+ e_Fb(i__),e1Fb(i__),e2Fb(i__),e_Is(i__),E_Is_1(i__),
+ Ndatap(i__),redchi2p(i__),Is0mag(i__),e_Is0mag(i__),
+ E_Is0mag_1(i__),fS_1(i__),e1fS_1(i__),e2fS_1(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 'tablea2.dat' ! All 122 candidates with microlensing parameters
* fitted without blending
C Format for file interpretation
2 format(
+ A11,1X,I1,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F7.2,
+ 1X,F7.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F11.6,1X,F11.5,1X,
+ F10.7,1X,F14.3,1X,F14.6,1X,F15.6,1X,F6.3,1X,F10.4,1X,I3,1X,
+ F7.4,1X,F7.3,1X,F8.3,1X,F9.6,1X,F10.6,1X,F10.6,1X,F12.5,1X,
+ F12.5,1X,F9.5,1X,F9.7,1X,F11.7,1X,F9.7,1X,F12.6,1X,F12.6,1X,
+ F11.6,1X,F12.6,1X,F12.6,1X,F11.6,1X,F6.3,1X,F6.3,1X,I3,1X,
+ F8.4,1X,F6.3,1X,F6.3,1X,F6.3,1X,F10.4,1X,F7.4,1X,F6.4)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,122
read(1,'(A457)')ar__1
read(ar__1,2)
+ OGLE_1(i__),Level_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Xpos_1(i__),
+ Ypos_1(i__),Imag_1(i__),v_V_I_1(i__),I0mag_1(i__),
+ v_V_I_0_1(i__),t0_1(i__),tE_1(i__),umin_1(i__),Amax_1(i__),
+ Fs_2(i__),Fb_1(i__),Is_1(i__),chi2_1(i__),Ndata_1(i__),
+ redchi2_1(i__),sigmax_1(i__),sumsig_1(i__),e_t0_1(i__),
+ e1t0_1(i__),e2t0_1(i__),e_tE_1(i__),e1tE_1(i__),e2tE_1(i__),
+ e_umin_1(i__),e1umin_1(i__),e2umin_1(i__),e_Fs_1(i__),
+ e1Fs_2(i__),e2Fs_2(i__),e_Fb_1(i__),e1Fb_1(i__),e2Fb_1(i__),
+ e_Is_2(i__),E_Is_3(i__),Ndatap_1(i__),redchi2p_1(i__),
+ Is0mag_1(i__),e_Is0mag_2(i__),E_Is0mag_3(i__),fS_3(i__),
+ e1fS_3(i__),e2fS_3(i__)
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,2)
+ OGLE_1(i__),Level_1(i__),RAh_1(i__),RAm_1(i__),RAs_1(i__),
+ DE__1(i__),DEd_1(i__),DEm_1(i__),DEs_1(i__),Xpos_1(i__),
+ Ypos_1(i__),Imag_1(i__),v_V_I_1(i__),I0mag_1(i__),
+ v_V_I_0_1(i__),t0_1(i__),tE_1(i__),umin_1(i__),Amax_1(i__),
+ Fs_2(i__),Fb_1(i__),Is_1(i__),chi2_1(i__),Ndata_1(i__),
+ redchi2_1(i__),sigmax_1(i__),sumsig_1(i__),e_t0_1(i__),
+ e1t0_1(i__),e2t0_1(i__),e_tE_1(i__),e1tE_1(i__),e2tE_1(i__),
+ e_umin_1(i__),e1umin_1(i__),e2umin_1(i__),e_Fs_1(i__),
+ e1Fs_2(i__),e2Fs_2(i__),e_Fb_1(i__),e1Fb_1(i__),e2Fb_1(i__),
+ e_Is_2(i__),E_Is_3(i__),Ndatap_1(i__),redchi2p_1(i__),
+ Is0mag_1(i__),e_Is0mag_2(i__),E_Is0mag_3(i__),fS_3(i__),
+ e1fS_3(i__),e2fS_3(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 'tablea3.dat' ! Cross-identification of candidates with other
* sources
C Format for file interpretation
3 format(A11,1X,A111)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,122
read(1,'(A123)')ar__2
read(ar__2,3)OGLE_2(i__),OName(i__)
c ..............Just test output...........
write(6,3)OGLE_2(i__),OName(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea4/*' ! *The light curve data for all 122 candidates
C Format for file interpretation
4 format(F10.5,1X,F12.6,1X,F10.6,1X,F7.3,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'tablea4/*')
write(6,*) '....Loading file: tablea4/*'
do i__=1,122
read(1,'(A50)')ar__3
read(ar__3,4)
+ JD_2450000(i__),Flux(i__),e_Flux(i__),Imag_2(i__),e_Imag(i__)
c ..............Just test output...........
write(6,4)
+ JD_2450000(i__),Flux(i__),e_Flux(i__),Imag_2(i__),e_Imag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea5.dat' ! Detection efficiencies for each field with blending
C Format for file interpretation
5 format(A4,1X,I2,2X,F4.2,2X,F8.4,1X,I7,1X,I6,1X,F9.6)
C Effective file loading
open(unit=1,status='old',file=
+'tablea5.dat')
write(6,*) '....Loading file: tablea5.dat'
do i__=1,540
read(1,'(A48)')ar__4
read(ar__4,5)
+ Field(i__),tEi(i__),log_tE(i__),tE_2(i__),Nout(i__),Nin(i__),
+ Deff(i__)
c ..............Just test output...........
write(6,5)
+ Field(i__),tEi(i__),log_tE(i__),tE_2(i__),Nout(i__),Nin(i__),
+ Deff(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end