Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/493/1171 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 2024-Apr-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/A+A/493/1171 Mutual phenomena of Galilean satellites PHEMU03 (Arlot+, 2009)
*================================================================================
*The PHEMU03 catalogue of observations of the mutual phenomena of the
*Galilean satellites of Jupiter.
* Arlot J.E., Thuillot W., Ruatti C., Ahmad A., Amosse A., Anbazhagan P.,
* Andreyev M., Antov A., Appakutty M., Asher D., Aubry S., Baron N.,
* Bassiere N., Berthe M., Bogdanovski R., Bosq F., Bredner E., Buettner D.,
* Buromsky M., Cammarata S., Casas R., Chis G.D., Christou A.A.,
* Coquerel J.-P., Corlan R., Cremaschini C., Crussaire D., Cuypers J.,
* Dennefeld M., Descamps P., Devyatkin A., Dimitrov D., Dorokhova T.N.,
* Dorokhov N.I., Dourneau G., Duenas M., Dumitrescu A., Emelianov N.,
* Ferrara D., Fiel D., Fienga A., Flatres T., Foglia S., Garlitz J.,
* Gerbos J., Gilbert R., Goncalves R.M.D., Gonzales D., , Gorda S.Yu.,
* Gorshanov D.L., Hansen M.W., Harrington M., Irsmambetova T.R., Ito Y.,
* Ivanova V., Izmailov I.S., Khovritchev M.Yu., Khrutskaya E.V., Kieken J.,
* Kiseleva T.P., Kuppuswamy K., Lainey V., Lavayssiere M., Lazzarotti P.,
* Le Campion J.-F., Lellouch E., Li Z.L., Lo Savio E., Lou M., Magny E.,
* Manek J., Marinello W., Marino G., McAuliffe J.P., Michelli M.,
* Moldovan D., Montagnac S., Moorthy V., Nickel O., Nier J.M., Noel T.,
* Noyelles B., Oksanen A., Parrat D., Pauwels T., Peng Q.Y., Pizzetti G.,
* Priban V., Ramachandran B., Rambaux N., Rapaport M., Rapavy P., Rau G.,
* Sacre J.-J., Sada P.V., Salvaggio F., Sarlin P., Sciuto C., Selvakumar G.,
* Sergeyev A., Sidorov M., Sorescu S., Spampinato S.A., Stellmacher I.,
* Trunkovsky E., Tejfel V., Tudose V., Turcu V., Ugarte I., Vantyghem P.,
* Vasundhara R., Vaubaillon J., Velu C., Venkataramana A.K., Vidal-Sainz J.,
* Vienne A., Vilar J., Vingerhoets P., Vollman W.
* <Astron. Astrophys. 493, 1171 (2009)>
* =2009A&A...493.1171A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'events.dat' ! List of mutual events
integer*4 nr__
parameter (nr__=119) ! Number of records
character*66 ar__ ! Full-size record
character*10 Date (nr__) ! ("Y/M/D") UT Observation date
character*1 n_Date (nr__) ! [*] Phenomenon observed at several locations
character*3 Phen (nr__) ! Phenomenon designation (G1)
character*1 n_Phen (nr__) ! [APT] Annular, Partial or Total phenomenon
character*8 Obs_Beg (nr__) ! ("h:m:s") Predicted time of the beginning of the event
character*8 Obs_Max (nr__) ! ("h:m:s") Predicted time of the maximum of the event
character*8 Obs_End (nr__) ! ("h:m:s") Predicted time of the end of the event
real*4 dmag (nr__) ! (mag) Magnitude drop
real*4 Ph (nr__) ! (s) Phase angle
real*4 Dist (nr__) ! Distance from the satellites to the center
* of Jupiter in Jovian radii
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'phemu03.dat' ! PHEMU03 observations (table4 of paper)
integer*4 nr__1
parameter (nr__1=377) ! Number of records
character*118 ar__1 ! Full-size record
character*10 Date_1 (nr__1) ! ("Y/M/D") Observation date (1)
character*1 n_Date_1 (nr__1) ! [*] Phenomenon observed at several locations
character*3 Phen_1 (nr__1) ! Phenomenon designation (G1)
character*1 n_Phen_1 (nr__1) ! [APT] Annular, Partial or Total phenomenon
character*14 Site (nr__1) ! Name of the site of observation
character*8 Obs_Max_1 (nr__1) ! ("h:m:s") Time of the observed maximum of the event
integer*4 e_Obs_Max (nr__1) ! (s) Error on Obs.Max
real*4 dmag_1 (nr__1) ! (mag) Magnitude drop
real*4 e_dmag (nr__1) ! (mag) rms uncertainty on magdrop
integer*4 C_O (nr__1) ! (s) Calculated minus observed time
character*1 Inst (nr__1) ! [TL] Code for the instrument used (2)
integer*4 Ap (nr__1) ! (cm) Aperture of the instrument in centimeters
character*5 Recept (nr__1) ! Code for the receptor
integer*4 EJup (nr__1) ! (deg) Elevation of Jupiter upon the horizon
integer*4 ESun (nr__1) ! (deg) Elevation of the Sun upon the horizon
integer*4 Qobs (nr__1) ! [0/3]?=- Conditions of observation code (3)
character*1 Filt (nr__1) ! [BCGRV-] Filter code
real*4 Tint (nr__1) ! (s) ?=- Integration time
real*4 Tint2 (nr__1) ! (s) ? Second value of integration time when
* interval
integer*4 Size (nr__1) ! (arcsec) ?=- Size of the diaphragm or window
character*3 Sat (nr__1) ! Order number of satellites in diaphragm,
* i.e. the satellites, the global magnitude
* drop of which was observed
character*13 File (nr__1) ! Name of observation file (4)
*Note (1): General informations of observation are in table times.dat,
* for the same Date, n_Date and Phen.
*Note (2): Code for the instrument used as follows:
* T = reflector
* L = refractor (i.e. Lunette)
*Note (3): Conditions of observation code as follows:
* 1 = Very good conditions
* 2 = Acceptable conditions
* 3 = Very difficult conditions
*Note (4): In "lc" subdirectory, with .txt extension for ascii files,
* in "fig" subdirectory, with .png extension for image files.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! Summary of photometric light curves
integer*4 nr__2
parameter (nr__2=372) ! Number of records
character*147 ar__2 ! Full-size record
character*13 File_1 (nr__2) ! File name (1)
character*13 Phen_2 (nr__2) ! Phenomenon
character*10 Date_2 (nr__2) ! ("D/M/Y") Observation date
character*28 Site_1 (nr__2) ! Site of observation
integer*4 Npts (nr__2) ! Number of points
real*8 alpha (nr__2) ! (h) Right ascension of Jupiter at 0h
* (in decimal hours)
real*8 delta (nr__2) ! (rad) Declination of Jupiter at 0h
* (in decimal radians)
real*8 JD (nr__2) ! (d) Julian date at 0h
real*8 Lon (nr__2) ! (h) Longitude of the observatory (in decimal hours)
* counted positively towards west
real*8 Lat (nr__2) ! (rad) Latitude of the observatory (in decimal radians)
* counted positively towards north
real*4 KA (nr__2) ! ? KA parameter (no used)
real*4 KX (nr__2) ! ? KX parameter (no used)
real*4 Z (nr__2) ! ? Zenithal distance
*Note (1): The light curves are provided in the files dXXXYYY.ZZZ where
* XXX is the coded name of the site of observation on 3 characters,
* YYY the type of phenomena, ZZZ the date.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'lc/*' ! Individual light curves in ascii
integer*4 nr__3
parameter (nr__3=372) ! Number of records
character*44 ar__3 ! Full-size record
character*2 Com (nr__3) ! Comments
integer*4 Obs_h (nr__3) ! (h) Hour of observation
integer*4 Obs_min (nr__3) ! (min) Minute of observation
real*4 Obs_s (nr__3) ! (s) Second of observation
real*4 mag (nr__3) ! (mag) Magnitude drop
real*8 e_mag (nr__3) ! (mag) rms uncertainty on mag
C=============================================================================
C Loading file 'events.dat' ! List of mutual events
C Format for file interpretation
1 format(
+ A10,A1,1X,A3,1X,A1,3X,A8,3X,A8,2X,A8,3X,F4.3,1X,F4.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'events.dat')
write(6,*) '....Loading file: events.dat'
do i__=1,119
read(1,'(A66)')ar__
read(ar__,1)
+ Date(i__),n_Date(i__),Phen(i__),n_Phen(i__),Obs_Beg(i__),
+ Obs_Max(i__),Obs_End(i__),dmag(i__),Ph(i__),Dist(i__)
c ..............Just test output...........
write(6,1)
+ Date(i__),n_Date(i__),Phen(i__),n_Phen(i__),Obs_Beg(i__),
+ Obs_Max(i__),Obs_End(i__),dmag(i__),Ph(i__),Dist(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'phemu03.dat' ! PHEMU03 observations (table4 of paper)
C Format for file interpretation
2 format(
+ A10,A1,1X,A3,1X,A1,2X,A14,1X,A8,1X,I3,1X,F5.3,2X,F4.3,1X,I4,
+ 1X,A1,1X,I3,1X,A5,1X,I3,2X,I3,1X,I1,1X,A1,2X,F4.2,1X,F3.1,1X,
+ I1,1X,A3,1X,A13)
C Effective file loading
open(unit=1,status='old',file=
+'phemu03.dat')
write(6,*) '....Loading file: phemu03.dat'
do i__=1,377
read(1,'(A118)')ar__1
read(ar__1,2)
+ Date_1(i__),n_Date_1(i__),Phen_1(i__),n_Phen_1(i__),Site(i__),
+ Obs_Max_1(i__),e_Obs_Max(i__),dmag_1(i__),e_dmag(i__),
+ C_O(i__),Inst(i__),Ap(i__),Recept(i__),EJup(i__),ESun(i__),
+ Qobs(i__),Filt(i__),Tint(i__),Tint2(i__),Size(i__),Sat(i__),
+ File(i__)
if (Qobs(i__) .EQ. 45) Qobs(i__) = iNULL__
if (idig(ar__1(91:94)).EQ.0) Tint(i__) = rNULL__
if(ar__1(96:98) .EQ. '') Tint2(i__) = rNULL__
if (Size(i__) .EQ. 45) Size(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Date_1(i__),n_Date_1(i__),Phen_1(i__),n_Phen_1(i__),Site(i__),
+ Obs_Max_1(i__),e_Obs_Max(i__),dmag_1(i__),e_dmag(i__),
+ C_O(i__),Inst(i__),Ap(i__),Recept(i__),EJup(i__),ESun(i__),
+ Qobs(i__),Filt(i__),Tint(i__),Tint2(i__),Size(i__),Sat(i__),
+ File(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'list.dat' ! Summary of photometric light curves
C Format for file interpretation
3 format(
+ A13,1X,A13,1X,A10,1X,A28,1X,I5,1X,F7.4,1X,F7.5,F9.1,2X,F9.6,
+ 1X,F9.7,1X,F3.2,1X,F6.2,1X,E16.10)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,372
read(1,'(A147)')ar__2
read(ar__2,3)
+ File_1(i__),Phen_2(i__),Date_2(i__),Site_1(i__),Npts(i__),
+ alpha(i__),delta(i__),JD(i__),Lon(i__),Lat(i__),KA(i__),
+ KX(i__),Z(i__)
if(ar__2(121:123) .EQ. '') KA(i__) = rNULL__
if(ar__2(125:130) .EQ. '') KX(i__) = rNULL__
if(ar__2(132:147) .EQ. '') Z(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ File_1(i__),Phen_2(i__),Date_2(i__),Site_1(i__),Npts(i__),
+ alpha(i__),delta(i__),JD(i__),Lon(i__),Lat(i__),KA(i__),
+ KX(i__),Z(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'lc/*' ! Individual light curves in ascii
C Format for file interpretation
4 format(3X,A2,1X,I2,1X,I2,1X,F5.2,2X,E16.11,2X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'lc/*')
write(6,*) '....Loading file: lc/*'
do i__=1,372
read(1,'(A44)')ar__3
read(ar__3,4)
+ Com(i__),Obs_h(i__),Obs_min(i__),Obs_s(i__),mag(i__),
+ e_mag(i__)
c ..............Just test output...........
write(6,4)
+ Com(i__),Obs_h(i__),Obs_min(i__),Obs_s(i__),mag(i__),
+ e_mag(i__)
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