Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/618/A44 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-Jun-13
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/A+A/618/A44 Predicted microlensing events from Gaia DR2 (Bramich, 2018)
*================================================================================
*Predicted microlensing events from analysis of Gaia Data Release 2.
* Bramich D.M.
* <Astron. Astrophys. 618, A44 (2018)>
* =2018A&A...618A..44B (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea7.dat' ! Characteristics of the microlensing events ME34-ME76
integer*4 nr__
parameter (nr__=43) ! Number of records
character*460 ar__ ! Full-size record
C Position composed of: RAdeg DEdeg (Epoch=2015.5)
character*4 Event ! Event name
character*5 LensSpType ! Lens spectral type
character*19 LensID ! Lens GDR2 ID
real*8 RALdeg ! (deg) Lens right ascension (ICRS, at Ep=2015.5)
real*8 e_RALdeg ! (mas) Error on lens right ascension
real*8 DELdeg ! (deg) Lens declination (ICRS, at Ep=2015.5)
real*8 e_DELdeg ! (mas) Error on lens declination
real*8 LenspmRA ! (mas/yr) Lens proper motion in right ascension
real*4 e_LenspmRA ! (mas/yr) Error on lens proper motion in
* right ascension
real*8 LenspmDE ! (mas/yr) Lens proper motion in declination
real*4 e_LenspmDE ! (mas/yr) Error on lens proper motion in declination
real*8 Lensplx ! (mas) Lens parallax
real*4 e_Lensplx ! (mas) Error on lens parallax
real*8 LensGmag ! (mag) Lens G mean magnitude
real*4 e_LensGmag ! (mag) Error on lens G mean magnitude
real*8 LensBPmag ! (mag) Lens G_BP mean magnitude
real*4 e_LensBPmag ! (mag) Error on lens G_BP mean magnitude
real*8 LensRPmag ! (mag) Lens G_RP mean magnitude
real*4 e_LensRPmag ! (mag) Error on lens G_RP mean magnitude
real*4 LensMass ! (Msun) Lens mass
integer*8 Source ! Source GDR2 ID
real*8 RAdeg ! (deg) Source right ascension (ICRS, at Ep=2015.5)
real*4 e_RAdeg ! (mas) Error on source right ascension
real*8 DEdeg ! (deg) Source declination (ICRS, at Ep=2015.5)
real*4 e_DEdeg ! (mas) Error on source declination
real*8 pmRA ! (mas/yr) ? Source proper motion in right ascension
real*4 e_pmRA ! (mas/yr) ? Error on source proper motion in
* right ascension
real*8 pmDE ! (mas/yr) ? Source proper motion in declination
real*4 e_pmDE ! (mas/yr) ? Error on source proper motion in
* declination
real*4 plx ! (mas) ? Source parallax
real*4 e_plx ! (mas) ? Error on source parallax
real*8 Gmag ! (mag) ? Source G mean magnitude
real*4 e_Gmag ! (mag) ? Error on source G mean magnitude
real*8 BPmag ! (mag) ? Source G_BP mean magnitude
real*4 e_BPmag ! (mag) ? Error on source G_BP mean magnitude
real*8 RPmag ! (mag) ? Source G_RP mean magnitude
real*4 e_RPmag ! (mag) ? Error on source G_RP mean magnitude
real*4 thetaE ! (mas) Einstein radius theta_E
real*4 e_thetaE ! (mas) Error on Einstein radius
real*4 u0t ! Minimum normalised source-lens separation
* u_0 in theta_E units
real*4 b_u0t ! 15.9 percentile of the minimum normalised
* source-lens separation u_0
* in theta_E units
real*4 B_u0t_1 ! 84.1 percentile of the minimum normalised
* source-lens separation u_0
* in theta_E units
real*8 u0 ! (mas) Minimum source-lens separation u_0
real*8 b_u0 ! (mas) 15.9 percentile of the minimum
* source-lens separation u_0
real*8 B_u0_1 ! (mas) 84.1 percentile of the minimum
* source-lens separation u_0
real*8 t0 ! (yr) Epoch of the event peak t_0 in Julian year
character*3 n_t0 ! [(ab) ] Note on t0 (1)
real*8 b_t0 ! (yr) 15.9 percentile of the epoch of the event
* peak t_0 in Julian year
real*8 B_t0_1 ! (yr) 84.1 percentile of the epoch of the event
* peak t_0 in Julian year
real*4 dtheta2 ! (mas) Difference between the minimum and maximum
* astrometric shifts of an event over the
* 12-year baseline of an extended
* Gaia mission
real*4 b_dtheta2 ! (mas) 15.9 percentile of the difference between
* the minimum and maximum astrometric shifts
* of an event over the 12-year baseline of
* an extended Gaia mission
real*4 B_dtheta2_1 ! (mas) 84.1 percentile of the difference between
* the minimum and maximum astrometric shifts
* of an event over the 12-year baseline of
* an extended Gaia mission
real*8 Tdtheta2 ! (d) Amount of time that an event spends with
* its astrometric shift above halfway
* between the minimum and maximum shift
real*8 b_Tdtheta2 ! (d) 15.9 percentile of the amount of time that
* an event spends with its astrometric
* shift above halfway between the minimum
* and maximum shift
real*8 B_Tdtheta2_1 ! (d) 84.1 percentile of the amount of time that
* an event spends with its astrometric
* shift above halfway between the minimum
* and maximum shift
*Note (1): Note as follows:
* (a) = Event peaks before this date and the event properties are only
* computed using the time period considered in this paper.
* (b) = Event peaks after this date. For event properties computed from a
* more appropriate time period, see Bramich & Nielsen (2018).
C=============================================================================
C Loading file 'tablea7.dat' ! Characteristics of the microlensing events ME34-ME76
C Format for file interpretation
1 format(
+ A4,1X,A5,1X,A19,1X,F14.10,2X,F7.4,1X,F14.10,2X,F7.4,1X,F9.3,
+ 1X,F5.3,1X,F8.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F8.4,1X,F6.4,1X,
+ F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F4.2,1X,I19,1X,F14.10,1X,F6.3,
+ 1X,F14.10,1X,F6.3,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.3,
+ 1X,F5.3,1X,F6.2,1X,F6.2,1X,F6.2,1X,F7.2,1X,F7.2,1X,F7.2,1X,
+ F10.5,A3,1X,F10.5,1X,F10.5,1X,F5.3,1X,F5.3,1X,F5.3,1X,F7.2,1X,
+ F7.2,1X,F7.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea7.dat')
write(6,*) '....Loading file: tablea7.dat'
do i__=1,43
read(1,'(A460)')ar__
read(ar__,1)
+ Event,LensSpType,LensID,RALdeg,e_RALdeg,DELdeg,e_DELdeg,
+ LenspmRA,e_LenspmRA,LenspmDE,e_LenspmDE,Lensplx,e_Lensplx,
+ LensGmag,e_LensGmag,LensBPmag,e_LensBPmag,LensRPmag,
+ e_LensRPmag,LensMass,Source,RAdeg,e_RAdeg,DEdeg,e_DEdeg,pmRA,
+ e_pmRA,pmDE,e_pmDE,plx,e_plx,Gmag,e_Gmag,BPmag,e_BPmag,RPmag,
+ e_RPmag,thetaE,e_thetaE,u0t,b_u0t,B_u0t_1,u0,b_u0,B_u0_1,t0,
+ n_t0,b_t0,B_t0_1,dtheta2,b_dtheta2,B_dtheta2_1,Tdtheta2,
+ b_Tdtheta2,B_Tdtheta2_1
if(ar__(240:246) .EQ. '') pmRA = rNULL__
if(ar__(248:252) .EQ. '') e_pmRA = rNULL__
if(ar__(254:260) .EQ. '') pmDE = rNULL__
if(ar__(262:266) .EQ. '') e_pmDE = rNULL__
if(ar__(268:273) .EQ. '') plx = rNULL__
if(ar__(275:279) .EQ. '') e_plx = rNULL__
if(ar__(281:287) .EQ. '') Gmag = rNULL__
if(ar__(289:294) .EQ. '') e_Gmag = rNULL__
if(ar__(296:302) .EQ. '') BPmag = rNULL__
if(ar__(304:309) .EQ. '') e_BPmag = rNULL__
if(ar__(311:317) .EQ. '') RPmag = rNULL__
if(ar__(319:324) .EQ. '') e_RPmag = rNULL__
c ..............Just test output...........
write(6,1)
+ Event,LensSpType,LensID,RALdeg,e_RALdeg,DELdeg,e_DELdeg,
+ LenspmRA,e_LenspmRA,LenspmDE,e_LenspmDE,Lensplx,e_Lensplx,
+ LensGmag,e_LensGmag,LensBPmag,e_LensBPmag,LensRPmag,
+ e_LensRPmag,LensMass,Source,RAdeg,e_RAdeg,DEdeg,e_DEdeg,pmRA,
+ e_pmRA,pmDE,e_pmDE,plx,e_plx,Gmag,e_Gmag,BPmag,e_BPmag,RPmag,
+ e_RPmag,thetaE,e_thetaE,u0t,b_u0t,B_u0t_1,u0,b_u0,B_u0_1,t0,
+ n_t0,b_t0,B_t0_1,dtheta2,b_dtheta2,B_dtheta2_1,Tdtheta2,
+ b_Tdtheta2,B_Tdtheta2_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end