Conversion of standardized ReadMe file for
file /./ftp/cats/J/AN/336/590 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/AN/336/590 Bochum Galactic Disk Survey: II (Hackstein+, 2015)
*================================================================================
*The Bochum Survey of the Southern Galactic Disk:
*II. Follow-up measurements and multi-filter photometry for 1323 square degrees
*monitored in 2010-2015.
* Hackstein M., Fein C., Haas M., Ramolla M., Pozo Nunez F.,
* Barr Dominguez A., Kaderhandt L., Thomsch K., Niedworok N., Westhues C.,
* Chini R.
* <Astron. Nachrichten 336, 590 (2015)>
* =2015AN....336..590H (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'varlc.dat' ! Light curves of all variable sources in r and i
integer*4 nr__
parameter (nr__=6945239) ! Number of records
character*128 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*55 Field (nr__) ! GDS Field
real*8 RAdeg (nr__) ! (deg) Right Ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Declination (J2000)
character*19 ID (nr__) ! GDS Variable ID, GDS_JHHMMSSs+DDMMSS
real*8 MJD (nr__) ! (d) Modified Julian Date
character*3 Filter (nr__) ! [r_s i_s] Filter (G1)
real*8 mag (nr__) ! (mag) Magnitude in filter
real*8 e_mag (nr__) ! (mag) Magnitude Error
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'varsum.dat' ! Summary file of variable properties
integer*4 nr__1
parameter (nr__1=64151) ! Number of records
character*297 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*55 Field_1 (nr__1) ! GDS field name
character*19 ID_1 (nr__1) ! GDS Variable ID, GDS_JHHMMSSs+DDMMSS
real*8 RAdeg_1 (nr__1) ! (deg) Mean right ascension (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Mean declination (J2000)
real*8 RArdeg (nr__1) ! (deg) ? r-band right ascension (J2000)
real*8 DErdeg (nr__1) ! (deg) ? r-band declination (J2000)
real*8 RAideg (nr__1) ! (deg) ? i-band right ascension (J2000)
real*8 DEideg (nr__1) ! (deg) ? i-band declination (J2000)
character*1 Grade (nr__1) ! [A-C] GDS grade (1)
integer*4 o_rmag (nr__1) ! Number of measurements in r
integer*4 o_imag (nr__1) ! Number of measurements in i
real*4 rmag (nr__1) ! (mag) ? Median light curve magnitude in r
real*4 s_rmag (nr__1) ! (mag) ? Non-variable equivalent standard deviation
* in r
real*4 e_rmag (nr__1) ! (mag) ? Non-variable equivalent median absolute
* deviation in r
real*4 imag (nr__1) ! (mag) ? Median light curve magnitude in i
real*4 s_imag (nr__1) ! (mag) ? Non-variable equivalent standard deviation
* in i
real*4 e_imag (nr__1) ! (mag) ? Non-variable equivalent median absolute
* deviation in i
integer*4 qflagr (nr__1) ! ? SourceExtractor quality flag in r
integer*4 qflagi (nr__1) ! ? SourceExtractor quality flag in i
real*4 AmpMax (nr__1) ! (mag) Maximum amplitude (2)
integer*4 f_Ampr (nr__1) ! [0/1]? Amplitude flag in r (3)
integer*4 f_Ampi (nr__1) ! [0/1]? Amplitude flag in i (3)
real*4 Ampr (nr__1) ! (mag) ? Light curve amplitude in r
real*4 Ampi (nr__1) ! (mag) ? Light curve amplitude in i
real*4 athrr (nr__1) ! (mag) ? Detection threshold for amplitude method in r
real*4 athri (nr__1) ! (mag) ? Detection threshold for amplitude method in i
real*4 JStetmax (nr__1) ! Maximum Stetson-J (2)
integer*4 f_JStetr (nr__1) ! [0/1]? Stetson-J flag in r (3)
integer*4 f_JSteti (nr__1) ! [0/1]? Stetson-J flag in i (3)
real*4 JStetr (nr__1) ! ? Light curve Stetson-J in r
real*4 JSteti (nr__1) ! ? Light curve Stetson-J in i
real*4 Jthrr (nr__1) ! ? Detection threshold for Stetson-J method in r
real*4 Jthri (nr__1) ! ? Detection threshold for Stetson-J method in i
real*4 sdMax (nr__1) ! (mag) Maximum standard deviation (2)
integer*4 f_sdr (nr__1) ! [0/1]? Standard deviation flag in r (3)
integer*4 f_sdi (nr__1) ! [0/1]? Standard deviation flag in i (3)
real*4 sdr (nr__1) ! (mag) ? Light curve standard deviation in r
real*4 sdi (nr__1) ! (mag) ? Light curve standard deviation in i
real*4 sdthrr (nr__1) ! (mag) ? Detection threshold for standard deviation
* method in r
real*4 sdthri (nr__1) ! (mag) ? Detection threshold for standard deviation
* method in i
integer*4 VSXMatch (nr__1) ! [0/2] Source found in VSX (4)
integer*4 VSXOID (nr__1) ! ? VSX internal identifier
*Note (1): GRADE can be A, B, or C with decreasing variability detection quality.
*Note (2): Maximum value of both filters
*Note (3): Standard deviation flag as follows:
* 0 = a variability detection with method was negative
* 1 = positive
*Note (4): Source found in VSX flag as follows:
* 0 = no match has been found in the AAVSO International Variable Star Intex
* 1 = match in the AAVSO International Variable Star Intex
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'allavg.dat' ! Median magnitudes of all GDS sources in UBVriz
integer*4 nr__2
parameter (nr__2=95330976) ! Number of records
character*49 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination (J2000)
character*3 Filter_1 (nr__2) ! Filter (G1)
real*4 mag_1 (nr__2) ! (mag) ? Magnitude in filter
real*4 e_mag_1 (nr__2) ! (mag) ? Magnitude Error
real*4 Amp (nr__2) ! (mag) ? Maximum amplitude of measurements in filter
integer*4 qflag (nr__2) ! [0/18]? SExtractor quality flag in filter
integer*4 N (nr__2) ! ? Number of measurements in filter
C=============================================================================
C Loading file 'varlc.dat' ! Light curves of all variable sources in r and i
C Format for file interpretation
1 format(
+ A55,1X,F9.5,1X,F9.5,1X,A19,1X,F11.5,1X,A3,1X,F8.5,1X,F7.5)
C Effective file loading
open(unit=1,status='old',file=
+'varlc.dat')
write(6,*) '....Loading file: varlc.dat'
do i__=1,6945239
read(1,'(A128)')ar__
read(ar__,1)
+ Field(i__),RAdeg(i__),DEdeg(i__),ID(i__),MJD(i__),Filter(i__),
+ mag(i__),e_mag(i__)
c ..............Just test output...........
write(6,1)
+ Field(i__),RAdeg(i__),DEdeg(i__),ID(i__),MJD(i__),Filter(i__),
+ mag(i__),e_mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'varsum.dat' ! Summary file of variable properties
C Format for file interpretation
2 format(
+ A55,1X,A19,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,
+ A1,1X,I3,1X,I3,1X,F5.2,1X,F5.3,1X,F5.3,1X,F5.2,1X,F5.3,1X,
+ F5.3,1X,I1,1X,I1,1X,F4.2,1X,I1,1X,I1,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F6.2,1X,I1,1X,I1,1X,F6.2,1X,F6.2,1X,F4.2,1X,F4.2,
+ 1X,F6.4,1X,I1,1X,I1,1X,F6.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,I1,1X,
+ I6)
C Effective file loading
open(unit=1,status='old',file=
+'varsum.dat')
write(6,*) '....Loading file: varsum.dat'
do i__=1,64151
read(1,'(A297)')ar__1
read(ar__1,2)
+ Field_1(i__),ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),RArdeg(i__),
+ DErdeg(i__),RAideg(i__),DEideg(i__),Grade(i__),o_rmag(i__),
+ o_imag(i__),rmag(i__),s_rmag(i__),e_rmag(i__),imag(i__),
+ s_imag(i__),e_imag(i__),qflagr(i__),qflagi(i__),AmpMax(i__),
+ f_Ampr(i__),f_Ampi(i__),Ampr(i__),Ampi(i__),athrr(i__),
+ athri(i__),JStetmax(i__),f_JStetr(i__),f_JSteti(i__),
+ JStetr(i__),JSteti(i__),Jthrr(i__),Jthri(i__),sdMax(i__),
+ f_sdr(i__),f_sdi(i__),sdr(i__),sdi(i__),sdthrr(i__),
+ sdthri(i__),VSXMatch(i__),VSXOID(i__)
if(ar__1(97:105) .EQ. '') RArdeg(i__) = rNULL__
if(ar__1(107:115) .EQ. '') DErdeg(i__) = rNULL__
if(ar__1(117:125) .EQ. '') RAideg(i__) = rNULL__
if(ar__1(127:135) .EQ. '') DEideg(i__) = rNULL__
if(ar__1(147:151) .EQ. '') rmag(i__) = rNULL__
if(ar__1(153:157) .EQ. '') s_rmag(i__) = rNULL__
if(ar__1(159:163) .EQ. '') e_rmag(i__) = rNULL__
if(ar__1(165:169) .EQ. '') imag(i__) = rNULL__
if(ar__1(171:175) .EQ. '') s_imag(i__) = rNULL__
if(ar__1(177:181) .EQ. '') e_imag(i__) = rNULL__
if(ar__1(183:183) .EQ. '') qflagr(i__) = iNULL__
if(ar__1(185:185) .EQ. '') qflagi(i__) = iNULL__
if(ar__1(192:192) .EQ. '') f_Ampr(i__) = iNULL__
if(ar__1(194:194) .EQ. '') f_Ampi(i__) = iNULL__
if(ar__1(196:199) .EQ. '') Ampr(i__) = rNULL__
if(ar__1(201:204) .EQ. '') Ampi(i__) = rNULL__
if(ar__1(206:209) .EQ. '') athrr(i__) = rNULL__
if(ar__1(211:214) .EQ. '') athri(i__) = rNULL__
if(ar__1(223:223) .EQ. '') f_JStetr(i__) = iNULL__
if(ar__1(225:225) .EQ. '') f_JSteti(i__) = iNULL__
if(ar__1(227:232) .EQ. '') JStetr(i__) = rNULL__
if(ar__1(234:239) .EQ. '') JSteti(i__) = rNULL__
if(ar__1(241:244) .EQ. '') Jthrr(i__) = rNULL__
if(ar__1(246:249) .EQ. '') Jthri(i__) = rNULL__
if(ar__1(258:258) .EQ. '') f_sdr(i__) = iNULL__
if(ar__1(260:260) .EQ. '') f_sdi(i__) = iNULL__
if(ar__1(262:267) .EQ. '') sdr(i__) = rNULL__
if(ar__1(269:274) .EQ. '') sdi(i__) = rNULL__
if(ar__1(276:281) .EQ. '') sdthrr(i__) = rNULL__
if(ar__1(283:288) .EQ. '') sdthri(i__) = rNULL__
if(ar__1(292:297) .EQ. '') VSXOID(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ Field_1(i__),ID_1(i__),RAdeg_1(i__),DEdeg_1(i__),RArdeg(i__),
+ DErdeg(i__),RAideg(i__),DEideg(i__),Grade(i__),o_rmag(i__),
+ o_imag(i__),rmag(i__),s_rmag(i__),e_rmag(i__),imag(i__),
+ s_imag(i__),e_imag(i__),qflagr(i__),qflagi(i__),AmpMax(i__),
+ f_Ampr(i__),f_Ampi(i__),Ampr(i__),Ampi(i__),athrr(i__),
+ athri(i__),JStetmax(i__),f_JStetr(i__),f_JSteti(i__),
+ JStetr(i__),JSteti(i__),Jthrr(i__),Jthri(i__),sdMax(i__),
+ f_sdr(i__),f_sdi(i__),sdr(i__),sdi(i__),sdthrr(i__),
+ sdthri(i__),VSXMatch(i__),VSXOID(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'allavg.dat' ! Median magnitudes of all GDS sources in UBVriz
C Format for file interpretation
3 format(F9.5,1X,F9.5,1X,A3,1X,F6.3,1X,F5.3,1X,F5.3,1X,I2,1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'allavg.dat')
write(6,*) '....Loading file: allavg.dat'
do i__=1,95330976
read(1,'(A49)')ar__2
read(ar__2,3)
+ RAdeg_2(i__),DEdeg_2(i__),Filter_1(i__),mag_1(i__),
+ e_mag_1(i__),Amp(i__),qflag(i__),N(i__)
if(ar__2(25:30) .EQ. '') mag_1(i__) = rNULL__
if(ar__2(32:36) .EQ. '') e_mag_1(i__) = rNULL__
if(ar__2(38:42) .EQ. '') Amp(i__) = rNULL__
if(ar__2(44:45) .EQ. '') qflag(i__) = iNULL__
if(ar__2(47:49) .EQ. '') N(i__) = iNULL__
c ..............Just test output...........
write(6,3)
+ RAdeg_2(i__),DEdeg_2(i__),Filter_1(i__),mag_1(i__),
+ e_mag_1(i__),Amp(i__),qflag(i__),N(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end