Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/107/1022 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-20
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/AJ/107/1022 UBVRI photometry of SN 1993J (Richmond+, 1994)
*================================================================================
*UBVRI photometry of SN 1993J in M81: the first 120 days
* RICHMOND M.W., TREFFERS R.R., FILIPPENKO A.V., PAIK Y., LEIBUNDGUT B.,
* SCHULMAN E., COX C.V.
* <Astron. J. 107, 1022 (1994)>
* =1994AJ....107.1022R (SIMBAD/NED Reference)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2' ! Leuschner 76-cm observations of SN 1993J
integer*4 nr__
parameter (nr__=86) ! Number of records
character*62 ar__ ! Full-size record
character*9 UT (nr__) ! UT date (Month day)
real*8 JD (nr__) ! (d) Julian Date
real*4 Umag (nr__) ! (mag) []? U band magnitude (1)
character*1 u_Umag (nr__) ! [:+] Uncertainty flag for U (2)
real*4 Bmag (nr__) ! (mag) []? B band magnitude
character*1 u_Bmag (nr__) ! [:+] Uncertainty flag for B (2)
real*4 Vmag (nr__) ! (mag) []? V band magnitude
character*1 u_Vmag (nr__) ! [:+] Uncertainty flag for V (2)
real*4 Rmag (nr__) ! (mag) []? R band magnitude
character*1 u_Rmag (nr__) ! [:+] Uncertainty flag for R (2)
real*4 Imag (nr__) ! (mag) []? I band magnitude
character*1 u_Imag (nr__) ! [:+] Uncertainty flag for I (2)
character*4 Notes (nr__) ! [ CP?] C = cloudy, P = poor seeing
character*1 Note7 (nr__) ! [IP] I = Isaac Newton Telescope;
* P = poor conditions (3)
*Note (1): The 76 cm telescope's U-band values do not transform well onto the
* standard color system, but may at least give a qualitative idea of
* the behavior of SN 1993J at near-UV wavelengths.
*Note (2): ":" denotes an uncertain value, and
* "+" one which has been derived using color terms from the previous
* or following night.
*Note (3): In table7, all data are from the 1-m Jacobus Kapteyn Telescope unless
* marked I.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3' ! Leuschner 50-cm observations of SN 1993J
integer*4 nr__1
parameter (nr__1=47) ! Number of records
character*62 ar__1 ! Full-size record
character*9 UT_1 (nr__1) ! UT date (Month day)
real*8 JD_1 (nr__1) ! (d) Julian Date
real*4 Umag_1 (nr__1) ! (mag) []? U band magnitude (1)
character*1 u_Umag_1 (nr__1) ! [:+] Uncertainty flag for U (2)
real*4 Bmag_1 (nr__1) ! (mag) []? B band magnitude
character*1 u_Bmag_1 (nr__1) ! [:+] Uncertainty flag for B (2)
real*4 Vmag_1 (nr__1) ! (mag) []? V band magnitude
character*1 u_Vmag_1 (nr__1) ! [:+] Uncertainty flag for V (2)
real*4 Rmag_1 (nr__1) ! (mag) []? R band magnitude
character*1 u_Rmag_1 (nr__1) ! [:+] Uncertainty flag for R (2)
real*4 Imag_1 (nr__1) ! (mag) []? I band magnitude
character*1 u_Imag_1 (nr__1) ! [:+] Uncertainty flag for I (2)
character*4 Notes_1 (nr__1) ! [ CP?] C = cloudy, P = poor seeing
character*1 Note7_1 (nr__1) ! [IP] I = Isaac Newton Telescope;
* P = poor conditions (3)
*Note (1): The 76 cm telescope's U-band values do not transform well onto the
* standard color system, but may at least give a qualitative idea of
* the behavior of SN 1993J at near-UV wavelengths.
*Note (2): ":" denotes an uncertain value, and
* "+" one which has been derived using color terms from the previous
* or following night.
*Note (3): In table7, all data are from the 1-m Jacobus Kapteyn Telescope unless
* marked I.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table14' ! Summed optical fluxes of SN 1993J
integer*4 nr__2
parameter (nr__2=95) ! Number of records
character*54 ar__2 ! Full-size record
real*8 JD_2 (nr__2) ! (d) Julian Date
real*4 Flux (nr__2) ! (mW/m2) Fluxes from measured light curves
real*4 LowExt (nr__2) ! (mW/m2) Fluxes for extinction [E(B-V)=0.08 mag]
real*4 ModExt (nr__2) ! (mW/m2) Fluxes for extinction [E(B-V)=0.32 mag]
character*2 Site (nr__2) ! Observing site (1)
character*5 Bands (nr__2) ! Bands summed to determine fluxes
*Note (1): Site: LP = La Palma; KP = Kitt Peak; L = Leuschner
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table15' ! Blackbody temp. and luminosity of SN 1993J
integer*4 nr__3
parameter (nr__3=95) ! Number of records
character*45 ar__3 ! Full-size record
real*8 JD_3 (nr__3) ! (d) Julian Date
integer*4 LowTemp (nr__3) ! (K) Temperature from low extinction
* [E(B-V)=0.08 mag]
real*4 LowLum (nr__3) ! (10+35W) Bolometric luminosity from LowTemp
integer*4 ModTemp (nr__3) ! (K) Temperature from moderate extinction
* [E(B-V)=0.32 mag]
real*4 ModLum (nr__3) ! (10+35W) Bolometric luminosity from ModTemp
character*2 Site_1 (nr__3) ! Observing site (1)
character*5 Bands_1 (nr__3) ! Bands used to determine luminosities
*Note (1): Site: LP = La Palma; KP = Kitt Peak; L = Leuschner
C=============================================================================
C Loading file 'table2' ! Leuschner 76-cm observations of SN 1993J
C Format for file interpretation
1 format(
+ A9,1X,F10.2,1X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A1,1X,
+ F5.2,A1,1X,A4,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table2')
write(6,*) '....Loading file: table2'
do i__=1,86
read(1,'(A62)')ar__
read(ar__,1)
+ UT(i__),JD(i__),Umag(i__),u_Umag(i__),Bmag(i__),u_Bmag(i__),
+ Vmag(i__),u_Vmag(i__),Rmag(i__),u_Rmag(i__),Imag(i__),
+ u_Imag(i__),Notes(i__),Note7(i__)
if(ar__(22:26) .EQ. '') Umag(i__) = rNULL__
if(ar__(29:33) .EQ. '') Bmag(i__) = rNULL__
if(ar__(36:40) .EQ. '') Vmag(i__) = rNULL__
if(ar__(43:47) .EQ. '') Rmag(i__) = rNULL__
if(ar__(50:54) .EQ. '') Imag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ UT(i__),JD(i__),Umag(i__),u_Umag(i__),Bmag(i__),u_Bmag(i__),
+ Vmag(i__),u_Vmag(i__),Rmag(i__),u_Rmag(i__),Imag(i__),
+ u_Imag(i__),Notes(i__),Note7(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3' ! Leuschner 50-cm observations of SN 1993J
C Format for file interpretation
2 format(
+ A9,1X,F10.2,1X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A1,1X,F5.2,A1,1X,
+ F5.2,A1,1X,A4,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3')
write(6,*) '....Loading file: table3'
do i__=1,47
read(1,'(A62)')ar__1
read(ar__1,2)
+ UT_1(i__),JD_1(i__),Umag_1(i__),u_Umag_1(i__),Bmag_1(i__),
+ u_Bmag_1(i__),Vmag_1(i__),u_Vmag_1(i__),Rmag_1(i__),
+ u_Rmag_1(i__),Imag_1(i__),u_Imag_1(i__),Notes_1(i__),
+ Note7_1(i__)
if(ar__1(22:26) .EQ. '') Umag_1(i__) = rNULL__
if(ar__1(29:33) .EQ. '') Bmag_1(i__) = rNULL__
if(ar__1(36:40) .EQ. '') Vmag_1(i__) = rNULL__
if(ar__1(43:47) .EQ. '') Rmag_1(i__) = rNULL__
if(ar__1(50:54) .EQ. '') Imag_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ UT_1(i__),JD_1(i__),Umag_1(i__),u_Umag_1(i__),Bmag_1(i__),
+ u_Bmag_1(i__),Vmag_1(i__),u_Vmag_1(i__),Rmag_1(i__),
+ u_Rmag_1(i__),Imag_1(i__),u_Imag_1(i__),Notes_1(i__),
+ Note7_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table14' ! Summed optical fluxes of SN 1993J
C Format for file interpretation
3 format(F10.2,2X,E9.3,2X,E9.3,2X,E9.3,2X,A2,2X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table14')
write(6,*) '....Loading file: table14'
do i__=1,95
read(1,'(A54)')ar__2
read(ar__2,3)
+ JD_2(i__),Flux(i__),LowExt(i__),ModExt(i__),Site(i__),
+ Bands(i__)
c ..............Just test output...........
write(6,3)
+ JD_2(i__),Flux(i__),LowExt(i__),ModExt(i__),Site(i__),
+ Bands(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table15' ! Blackbody temp. and luminosity of SN 1993J
C Format for file interpretation
4 format(F10.2,1X,I5,2X,F4.2,1X,I5,1X,F5.2,2X,A2,2X,A5)
C Effective file loading
open(unit=1,status='old',file=
+'table15')
write(6,*) '....Loading file: table15'
do i__=1,95
read(1,'(A45)')ar__3
read(ar__3,4)
+ JD_3(i__),LowTemp(i__),LowLum(i__),ModTemp(i__),ModLum(i__),
+ Site_1(i__),Bands_1(i__)
c ..............Just test output...........
write(6,4)
+ JD_3(i__),LowTemp(i__),LowLum(i__),ModTemp(i__),ModLum(i__),
+ Site_1(i__),Bands_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end