Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJS/244/29 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 2026-Aug-16
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/ApJS/244/29 Microlensing events toward the Galactic bulge (Mroz+, 2019)
*================================================================================
*Microlensing optical depth and event rate toward the Galactic bulge from 8yr of
*OGLE-IV observations.
* Mroz P., Udalski A., Skowron J., Szymanski M.K., Soszynski I.,
* Wyrzykowski L., Pietrukowicz P., Kozlowski S., Poleski R., Ulaczyk K.,
* Rybicki K., Iwanek P.
* <Astrophys. J. Suppl. Ser., 244, 29 (2019)>
* =2019ApJS..244...29M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Best-fitting parameters of the analyzed
microlensing events in low-cadence OGLE fields
integer*4 nr__
parameter (nr__=5790) ! Number of records
character*294 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*16 Name ! Star name
character*6 Field ! OGLE-IV field name
integer*4 sField ! OGLE-IV subfield name within Field
integer*4 Star ! [50/208567] Star ID in the OGLE-IV database
integer*4 RAh ! (h) [17/18] Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! [-] Sign of the Declination (J2000)
integer*4 DEd ! (deg) [19/38] Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*8 RAdeg ! (deg) [258.28/284.31] Right Ascension (J2000)
real*8 DEdeg ! (deg) [-38.4/-19.2] Declination (J2000)
real*8 GLON ! (deg) [-8.6/9.9] Galactic longitude
real*8 GLAT ! (deg) [-14.9/6.3] Galactic latitude
real*8 t0_best ! (d) Best-fit t0, time of closest approach,
* Heliocentric Julian Date
real*8 tE_best ! (d) [1.58/303] Best-fit tE, Einstein timescale
real*4 u0_best ! [0/1] Best-fit u0, projected separation
real*4 Is_best ! (mag) [13.9/21] Best-fit I-band source magnitude
real*4 fs_best ! [0.01/3.8] Best-fit blending parameter (1)
real*8 t0_med ! (d) Median t0, time of closest approach,
* Heliocentric Julian Date
real*8 e_t0_med ! (d) [-37.2/0] Negative error bar, t0-med
real*4 E_t0_med_1 ! (d) [0/43] Positive error bar, t0-med
real*8 tE_med ! (d) [1.7/8812] Median tE, Einstein timescale
real*8 e_tE_med ! (d) [-4821/-0.04] Negative error bar, tE-med
real*8 E_tE_med_1 ! (d) [0.06/19953] Positive error bar, tE-med
real*4 u0_med ! [0/1.3] Median u0, projected separation
real*4 e_u0_med ! [-0.8/0] Negative error bar, u0-med
real*4 E_u0_med_1 ! [0/1.5] Positive error bar, u0-med
real*4 Is_med ! (mag) [13.9/26.5] Median I-band source magnitude
real*4 e_Is_med ! (mag) [-6.8/-0.005] Negative error bar, Is-med
real*4 E_Is_med_1 ! (mag) [0.006/5.2] Positive error bar, Is-med
real*4 fs_med ! [0/3.5] Median blending parameter (1)
real*4 e_fs_med ! [-0.6/0] Negative error bar, fs-med
real*4 E_fs_med_1 ! [0/0.9] Positive error bar, fs-med
real*4 Weight ! [2/129.2] Inverse of the detection efficiency
character*37 EWS ! OGLE EWS ID(s)
*Note (1): The blending parameter is defined as the ratio between the source
* flux and the total flux of the detected object (i.e., the sum of the
* fluxes of the source and unrelated blends).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! Surface density of stars in OGLE-IV subfields
calculated using image-level simulations
integer*4 nr__1
parameter (nr__1=3842) ! Number of records
character*66 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 Field_1 ! OGLE-IV field name
integer*4 sField_1 ! OGLE-IV subfield name within Field
real*8 RAdeg_1 ! (deg) [258.28/286.67] Right ascension (J2000)
real*8 DEdeg_1 ! (deg) [-38.4/-19.2] Declination (J2000)
real*4 GLON_1 ! (deg) [-8.6/10.6] Galactic longitude
real*8 GLAT_1 ! (deg) [-16.5/6.4] Galactic latitude
real*4 sigma18 ! (arcmin-2) [3.7/369.7] Star surface density in
* the range 14<I<18
real*4 sigma21 ! (arcmin-2) [42.8/4609.2] Star surface density in
* the range 14<I<21
integer*4 N18 ! [603/60940] Number of stars in
* the range 14<I<18
integer*4 N21 ! [7183/756286] Number of stars in
* the range 14<I<21
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! Basic information about analyzed fields
integer*4 nr__2
parameter (nr__2=121) ! Number of records
character*51 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*6 Field_2 ! OGLE-IV field name
real*8 RAdeg_2 ! (deg) [259/286] Right ascension (J2000)
real*8 DEdeg_2 ! (deg) [-38/-19.7] Declination (J2000)
real*8 GLON_2 ! (deg) [-8/10] Galactic longitude
real*8 GLAT_2 ! (deg) [-16/6] Galactic latitude
real*4 Ns ! (10+6) [0.5/7.1] Number of monitored stars
integer*4 Nep ! [36/12099] Number of epochs
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Microlensing optical depth and event rates in the
OGLE-IV fields (averaged over sources brighter
than I=21).
integer*4 nr__3
parameter (nr__3=121) ! Number of records
character*73 ar__3 ! Full-size record
C Position composed of: GLON GLAT
character*6 Field_3 ! OGLE-IV field name
real*8 GLON_3 ! (deg) [-8/10] Galactic longitude
real*8 GLAT_3 ! (deg) [-16/6] Galactic latitude
real*4 tau ! (10-6) [0/3.7] Optical depth {tau}_300_
real*4 e_tau ! (10-6) [0/1.3] Uncertainty of tau
real*4 Gam ! (10-6/yr) [0/26.5] Event rate {Gamma} per star
* per 10^6^ year
real*4 e_Gam ! (10-6/yr) [0/7.1] Uncertainty of gam
real*4 Gam_d ! (deg-2/yr) [0/265.3] Event rate {Gamma}_deg2_
* per unit area (in square degree)
* per year
real*4 e_Gam_d ! (deg-2/yr) [0/16.3] Uncertainty of gam-deg2
real*4 tE_mean ! (d) [0/79.6] Mean Einstein timescale <t_E_>
real*4 e_tE_mean ! (d) [0/79.6] Uncertainty of tE-mean
integer*4 Nev ! [0/441] Number of detected events
real*4 N21_1 ! (10+6) [0.4/19.1] Number of sources brighter
* than I=21
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'phot/*' ! Calibrated I-band photometry of individual events,
in ASCII format
integer*4 nr__4
parameter (nr__4=5790) ! Number of records
character*26 ar__4 ! Full-size record
real*8 HJD ! (d) Heliocentric Julian Date
real*4 Imag ! (mag) [12.2/23.8] I-band magnitude
real*4 e_Imag ! (mag) [0.003/4] I-band magnitude uncertainty
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'eff/*' ! Event detection efficiencies in individual fields
integer*4 nr__5
parameter (nr__5=121) ! Number of records
character*22 ar__5 ! Full-size record
real*4 logtE_min ! ([d]) [-1/2.4] Minimal log(t_E_)
real*4 logtE_max ! ([d]) [-0.9/2.5] Maximal log(t_E_)
real*8 Eff ! [0.0005/0.41] Detection efficiency
C=============================================================================
C Loading file 'table3.dat' ! Best-fitting parameters of the analyzed
* microlensing events in low-cadence OGLE fields
C Format for file interpretation
1 format(
+ A16,1X,A6,1X,I2,1X,I6,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
+ F4.1,1X,F9.5,1X,F9.5,1X,F8.5,1X,F9.5,1X,F11.3,1X,F7.3,1X,F5.3,
+ 1X,F6.3,1X,F5.3,1X,F11.3,1X,F7.3,1X,F6.3,1X,F8.3,1X,F9.3,1X,
+ F9.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F5.3,
+ 1X,F6.3,1X,F5.3,1X,F6.2,1X,A37)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,5790
read(1,'(A294)')ar__
read(ar__,1)
+ Name,Field,sField,Star,RAh,RAm,RAs,DE_,DEd,DEm,DEs,RAdeg,
+ DEdeg,GLON,GLAT,t0_best,tE_best,u0_best,Is_best,fs_best,
+ t0_med,e_t0_med,E_t0_med_1,tE_med,e_tE_med,E_tE_med_1,u0_med,
+ e_u0_med,E_u0_med_1,Is_med,e_Is_med,E_Is_med_1,fs_med,
+ e_fs_med,E_fs_med_1,Weight,EWS
c ..............Just test output...........
write(6,1)
+ Name,Field,sField,Star,RAh,RAm,RAs,DE_,DEd,DEm,DEs,RAdeg,
+ DEdeg,GLON,GLAT,t0_best,tE_best,u0_best,Is_best,fs_best,
+ t0_med,e_t0_med,E_t0_med_1,tE_med,e_tE_med,E_tE_med_1,u0_med,
+ e_u0_med,E_u0_med_1,Is_med,e_Is_med,E_Is_med_1,fs_med,
+ e_fs_med,E_fs_med_1,Weight,EWS
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! Surface density of stars in OGLE-IV subfields
* calculated using image-level simulations
C Format for file interpretation
2 format(
+ A6,1X,I2,1X,F7.3,1X,F7.3,1X,F6.3,1X,F7.3,1X,F5.1,1X,F6.1,1X,
+ I5,1X,I6)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,3842
read(1,'(A66)')ar__1
read(ar__1,2)
+ Field_1,sField_1,RAdeg_1,DEdeg_1,GLON_1,GLAT_1,sigma18,
+ sigma21,N18,N21
c ..............Just test output...........
write(6,2)
+ Field_1,sField_1,RAdeg_1,DEdeg_1,GLON_1,GLAT_1,sigma18,
+ sigma21,N18,N21
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! Basic information about analyzed fields
C Format for file interpretation
3 format(A6,1X,F8.4,1X,F8.4,1X,F7.4,1X,F8.4,1X,F3.1,1X,I5)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,121
read(1,'(A51)')ar__2
read(ar__2,3)Field_2,RAdeg_2,DEdeg_2,GLON_2,GLAT_2,Ns,Nep
c ..............Just test output...........
write(6,3)Field_2,RAdeg_2,DEdeg_2,GLON_2,GLAT_2,Ns,Nep
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Microlensing optical depth and event rates in the
* OGLE-IV fields (averaged over sources brighter
* than I=21).
C Format for file interpretation
4 format(
+ A6,1X,F7.4,1X,F8.4,1X,F4.2,1X,F4.2,1X,F4.1,1X,F3.1,1X,F5.1,1X,
+ F4.1,1X,F4.1,1X,F4.1,1X,I3,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,121
read(1,'(A73)')ar__3
read(ar__3,4)
+ Field_3,GLON_3,GLAT_3,tau,e_tau,Gam,e_Gam,Gam_d,e_Gam_d,
+ tE_mean,e_tE_mean,Nev,N21_1
c ..............Just test output...........
write(6,4)
+ Field_3,GLON_3,GLAT_3,tau,e_tau,Gam,e_Gam,Gam_d,e_Gam_d,
+ tE_mean,e_tE_mean,Nev,N21_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'phot/*' ! Calibrated I-band photometry of individual events,
* in ASCII format
C Format for file interpretation
5 format(F13.5,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'phot/*')
write(6,*) '....Loading file: phot/*'
do i__=1,5790
read(1,'(A26)')ar__4
read(ar__4,5)HJD,Imag,e_Imag
c ..............Just test output...........
write(6,5)HJD,Imag,e_Imag
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'eff/*' ! Event detection efficiencies in individual fields
C Format for file interpretation
6 format(F6.3,1X,F6.3,1X,F8.6)
C Effective file loading
open(unit=1,status='old',file=
+'eff/*')
write(6,*) '....Loading file: eff/*'
do i__=1,121
read(1,'(A22)')ar__5
read(ar__5,6)logtE_min,logtE_max,Eff
c ..............Just test output...........
write(6,6)logtE_min,logtE_max,Eff
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end