Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-21
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/AcA/51/175 OGLE-II DIA microlensing events (Wozniak+, 2001)
*================================================================================
*Difference Image Analysis of the OGLE-II bulge data. II. Microlensing events.
* Wozniak P.R., Udalski A., Szymanski M., Kubiak M., Pietrzynski G.,
* Soszynski I., Zebrun K.
* <Acta Astron., 51, 175 (2001)>
* =2001AcA....51..175W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'crossid.dat' ! DIA-OGLE corresponding numbers
integer*4 nr__
parameter (nr__=520) ! Number of records
character*30 ar__ ! Full-size record
character*4 Field (nr__) ! Field number (SCNN)
integer*4 DIA (nr__) ! DIA star number in the field
integer*4 OGLE (nr__) ! ? OGLE number
real*4 Dist (nr__) ! (pix) ? Distance between identified positions (1pix=0.4")
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! OGLE-II bulge fields
integer*4 nr__1
parameter (nr__1=49) ! Number of records
character*47 ar__1 ! Full-size record
real*8 RA_ (nr__1) ! (deg) Right Ascension J2000
real*8 Dec (nr__1) ! (deg) Declination J2000
character*4 Field_1 (nr__1) ! Field number (scNN)
integer*4 RAh (nr__1) ! (h) Right ascension (J2000)
integer*4 RAm (nr__1) ! (min) Right ascension (J2000)
real*4 RAs (nr__1) ! (s) Right ascension (J2000)
character*1 DE_ (nr__1) ! Declination sign (J2000)
integer*4 DEd (nr__1) ! (deg) Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Declination (J2000)
integer*4 DEs (nr__1) ! (arcsec) Declination (J2000)
real*4 GLON (nr__1) ! (deg) [-12/11] Galactic longitude
real*4 GLAT (nr__1) ! (deg) Galactic latitude
integer*4 Nobs (nr__1) ! Number of observations
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablesa.dat' ! All microlensing events, positions
integer*4 nr__2
parameter (nr__2=520) ! Number of records
character*70 ar__2 ! Full-size record
real*8 RA__1 (nr__2) ! (deg) Right Ascension J2000
real*8 Dec_1 (nr__2) ! (deg) Declination J2000
character*4 Field_2 (nr__2) ! Field designation (scNN)
integer*4 DIA_1 (nr__2) ! DIA star number in the field (1)
real*8 Xpos (nr__2) ! (pix) X position (2)
real*8 Ypos (nr__2) ! (pix) Y position (2)
integer*4 RAh_1 (nr__2) ! (h) Right ascension (J2000)
integer*4 RAm_1 (nr__2) ! (min) Right ascension (J2000)
real*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)
real*4 DEs_1 (nr__2) ! (arcsec) Declination (J2000)
real*4 Dist1 (nr__2) ! Centroid distance between the variable and
* the nearest object on the reference frame
* detected with DoPhot
real*8 Flux1 (nr__2) ! Flux of the nearest star on the reference frame
real*4 m_OGLE (nr__2) ! (mag) Zero point offset between DIA fluxes and OGLE
* database (varies with X, Y)
* m(OGLE)=m(Zero)-2.5*log(flux)
character*1 Det (nr__2) ! [234] Flag for detection method (3)
*Note (1): Number in the database of candidate variable objects from DIA
*Note (2): 0.0 for the center of the pixel in lower left corner
*Note (3): Method of detection as follows:
* 2 = Microlensing events found in the algorithmic search
* 3 = Microlensing events recovered in the visual search
* 4 = Events recovered by cross identification with the OGLE catalog
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablesb.dat' ! All microlensing events, parameters
integer*4 nr__3
parameter (nr__3=520) ! Number of records
character*68 ar__3 ! Full-size record
character*4 Field_3 (nr__3) ! Field number (scNN)
integer*4 DIA_2 (nr__3) ! DIA star number in the field
real*8 F0 (nr__3) ! Baseline difference flux
real*8 F1 (nr__3) ! Lensed flux
real*8 u0 (nr__3) ! Impact parameter (units of Einstein radius)
real*8 t0 (nr__3) ! (s) Einstein ring crossing time
real*8 tmax (nr__3) ! (d) Moment of maximum light (HJD-2450000)
character*1 l_I0 (nr__3) ! Limit flag on I0
real*4 I0 (nr__3) ! (mag) Baseline I magnitude
real*8 f (nr__3) ! Ratio of the lensed flux to the total baseline flux
C=============================================================================
C Loading file 'crossid.dat' ! DIA-OGLE corresponding numbers
C Format for file interpretation
1 format(A4,4X,I4,4X,I6,3X,F5.3)
C Effective file loading
open(unit=1,file='crossid.dat', status='old')
write(6,*) '....Loading file: crossid.dat'
do i__=1,520
read(1,'(A30)')ar__
read(ar__,1)Field(i__),DIA(i__),OGLE(i__),Dist(i__)
if(ar__(17:22) .EQ. '') OGLE(i__) = iNULL__
if(ar__(26:30) .EQ. '') Dist(i__) = rNULL__
c ..............Just test output...........
write(6,1)Field(i__),DIA(i__),OGLE(i__),Dist(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table1.dat' ! OGLE-II bulge fields
C Format for file interpretation
2 format(
+ 4X,A4,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,F6.2,1X,
+ F5.2,2X,I3)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,49
read(1,'(A47)')ar__1
read(ar__1,2)
+ Field_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Nobs(i__)
c Derive coordinates RA_ and Dec from input data
c (RA_ and Dec are set to rNULL_ when unknown)
RA_(i__) = RAh(i__)
if(RAh(i__) .GE. 0) RA_(i__)=RAh(i__)*15.
if(RAm(i__) .GE. 0) RA_(i__)=RA_(i__)+RAm(i__)/4.
if(RAs(i__) .GE. 0) RA_(i__)=RA_(i__)+RAs(i__)/240.
Dec(i__) = DEd(i__)
if(DEm(i__) .GE. 0) Dec(i__)=Dec(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) Dec(i__)=Dec(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.Dec(i__).GE.0) Dec(i__)=-Dec(i__)
c ..............Just test output...........
write(6,2)
+ Field_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__),
+ DEm(i__),DEs(i__),GLON(i__),GLAT(i__),Nobs(i__)
write(6,'(6H Pos: 2F8.4)') RA_(i__),Dec(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablesa.dat' ! All microlensing events, positions
C Format for file interpretation
3 format(
+ A4,1X,I4,1X,F7.2,1X,F7.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,
+ 1X,F4.1,1X,F4.2,1X,F7.2,1X,F5.2,1X,A1)
C Effective file loading
open(unit=1,file='tablesa.dat', status='old')
write(6,*) '....Loading file: tablesa.dat'
do i__=1,520
read(1,'(A70)')ar__2
read(ar__2,3)
+ Field_2(i__),DIA_1(i__),Xpos(i__),Ypos(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),Dist1(i__),Flux1(i__),m_OGLE(i__),Det(i__)
c Derive coordinates RA__1 and Dec_1 from input data
c (RA__1 and Dec_1 are set to rNULL_ when unknown)
RA__1(i__) = RAh_1(i__)
if(RAh_1(i__) .GE. 0) RA__1(i__)=RAh_1(i__)*15.
if(RAm_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAm_1(i__)/4.
if(RAs_1(i__) .GE. 0) RA__1(i__)=RA__1(i__)+RAs_1(i__)/240.
Dec_1(i__) = DEd_1(i__)
if(DEm_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEm_1(i__)/60.
if(DEs_1(i__) .GE. 0) Dec_1(i__)=Dec_1(i__)+DEs_1(i__)/3600.
if(DE__1(i__).EQ.'-'.AND.Dec_1(i__).GE.0) Dec_1(i__)=-Dec_1(i__)
c ..............Just test output...........
write(6,3)
+ Field_2(i__),DIA_1(i__),Xpos(i__),Ypos(i__),RAh_1(i__),
+ RAm_1(i__),RAs_1(i__),DE__1(i__),DEd_1(i__),DEm_1(i__),
+ DEs_1(i__),Dist1(i__),Flux1(i__),m_OGLE(i__),Det(i__)
write(6,'(6H Pos: 2F8.4)') RA__1(i__),Dec_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablesb.dat' ! All microlensing events, parameters
C Format for file interpretation
4 format(
+ A4,1X,I4,1X,F8.2,1X,F8.2,1X,F7.5,1X,F7.2,1X,F7.2,1X,A1,F6.3,
+ 1X,F8.3)
C Effective file loading
open(unit=1,file='tablesb.dat', status='old')
write(6,*) '....Loading file: tablesb.dat'
do i__=1,520
read(1,'(A68)')ar__3
read(ar__3,4)
+ Field_3(i__),DIA_2(i__),F0(i__),F1(i__),u0(i__),t0(i__),
+ tmax(i__),l_I0(i__),I0(i__),f(i__)
c ..............Just test output...........
write(6,4)
+ Field_3(i__),DIA_2(i__),F0(i__),F1(i__),u0(i__),t0(i__),
+ tmax(i__),l_I0(i__),I0(i__),f(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end