Conversion of standardized ReadMe file for
file /./ftp/cats/J/AN/336/590 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 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/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 ! GDS Field
real*8 RAdeg ! (deg) Right Ascension (J2000)
real*8 DEdeg ! (deg) Declination (J2000)
character*19 ID ! GDS Variable ID, GDS_JHHMMSSs+DDMMSS
real*8 MJD ! (d) Modified Julian Date
character*3 Filter ! [r_s i_s] Filter (G1)
real*8 mag ! (mag) Magnitude in filter
real*8 e_mag ! (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 ! GDS field name
character*19 ID_1 ! GDS Variable ID, GDS_JHHMMSSs+DDMMSS
real*8 RAdeg_1 ! (deg) Mean right ascension (J2000)
real*8 DEdeg_1 ! (deg) Mean declination (J2000)
real*8 RArdeg ! (deg) ? r-band right ascension (J2000)
real*8 DErdeg ! (deg) ? r-band declination (J2000)
real*8 RAideg ! (deg) ? i-band right ascension (J2000)
real*8 DEideg ! (deg) ? i-band declination (J2000)
character*1 Grade ! [A-C] GDS grade (1)
integer*4 o_rmag ! Number of measurements in r
integer*4 o_imag ! Number of measurements in i
real*4 rmag ! (mag) ? Median light curve magnitude in r
real*4 s_rmag ! (mag) ? Non-variable equivalent standard deviation
* in r
real*4 e_rmag ! (mag) ? Non-variable equivalent median absolute
* deviation in r
real*4 imag ! (mag) ? Median light curve magnitude in i
real*4 s_imag ! (mag) ? Non-variable equivalent standard deviation
* in i
real*4 e_imag ! (mag) ? Non-variable equivalent median absolute
* deviation in i
integer*4 qflagr ! ? SourceExtractor quality flag in r
integer*4 qflagi ! ? SourceExtractor quality flag in i
real*4 AmpMax ! (mag) Maximum amplitude (2)
integer*4 f_Ampr ! [0/1]? Amplitude flag in r (3)
integer*4 f_Ampi ! [0/1]? Amplitude flag in i (3)
real*4 Ampr ! (mag) ? Light curve amplitude in r
real*4 Ampi ! (mag) ? Light curve amplitude in i
real*4 athrr ! (mag) ? Detection threshold for amplitude method in r
real*4 athri ! (mag) ? Detection threshold for amplitude method in i
real*4 JStetmax ! Maximum Stetson-J (2)
integer*4 f_JStetr ! [0/1]? Stetson-J flag in r (3)
integer*4 f_JSteti ! [0/1]? Stetson-J flag in i (3)
real*4 JStetr ! ? Light curve Stetson-J in r
real*4 JSteti ! ? Light curve Stetson-J in i
real*4 Jthrr ! ? Detection threshold for Stetson-J method in r
real*4 Jthri ! ? Detection threshold for Stetson-J method in i
real*4 sdMax ! (mag) Maximum standard deviation (2)
integer*4 f_sdr ! [0/1]? Standard deviation flag in r (3)
integer*4 f_sdi ! [0/1]? Standard deviation flag in i (3)
real*4 sdr ! (mag) ? Light curve standard deviation in r
real*4 sdi ! (mag) ? Light curve standard deviation in i
real*4 sdthrr ! (mag) ? Detection threshold for standard deviation
* method in r
real*4 sdthri ! (mag) ? Detection threshold for standard deviation
* method in i
integer*4 VSXMatch ! [0/2] Source found in VSX (4)
integer*4 VSXOID ! ? 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 ! (deg) Right Ascension (J2000)
real*8 DEdeg_2 ! (deg) Declination (J2000)
character*3 Filter_1 ! Filter (G1)
real*4 mag_1 ! (mag) ? Magnitude in filter
real*4 e_mag_1 ! (mag) ? Magnitude Error
real*4 Amp ! (mag) ? Maximum amplitude of measurements in filter
integer*4 qflag ! [0/18]? SExtractor quality flag in filter
integer*4 N ! ? 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,RAdeg,DEdeg,ID,MJD,Filter,mag,e_mag
c ..............Just test output...........
write(6,1)Field,RAdeg,DEdeg,ID,MJD,Filter,mag,e_mag
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,ID_1,RAdeg_1,DEdeg_1,RArdeg,DErdeg,RAideg,DEideg,
+ Grade,o_rmag,o_imag,rmag,s_rmag,e_rmag,imag,s_imag,e_imag,
+ qflagr,qflagi,AmpMax,f_Ampr,f_Ampi,Ampr,Ampi,athrr,athri,
+ JStetmax,f_JStetr,f_JSteti,JStetr,JSteti,Jthrr,Jthri,sdMax,
+ f_sdr,f_sdi,sdr,sdi,sdthrr,sdthri,VSXMatch,VSXOID
if(ar__1(97:105) .EQ. '') RArdeg = rNULL__
if(ar__1(107:115) .EQ. '') DErdeg = rNULL__
if(ar__1(117:125) .EQ. '') RAideg = rNULL__
if(ar__1(127:135) .EQ. '') DEideg = rNULL__
if(ar__1(147:151) .EQ. '') rmag = rNULL__
if(ar__1(153:157) .EQ. '') s_rmag = rNULL__
if(ar__1(159:163) .EQ. '') e_rmag = rNULL__
if(ar__1(165:169) .EQ. '') imag = rNULL__
if(ar__1(171:175) .EQ. '') s_imag = rNULL__
if(ar__1(177:181) .EQ. '') e_imag = rNULL__
if(ar__1(183:183) .EQ. '') qflagr = iNULL__
if(ar__1(185:185) .EQ. '') qflagi = iNULL__
if(ar__1(192:192) .EQ. '') f_Ampr = iNULL__
if(ar__1(194:194) .EQ. '') f_Ampi = iNULL__
if(ar__1(196:199) .EQ. '') Ampr = rNULL__
if(ar__1(201:204) .EQ. '') Ampi = rNULL__
if(ar__1(206:209) .EQ. '') athrr = rNULL__
if(ar__1(211:214) .EQ. '') athri = rNULL__
if(ar__1(223:223) .EQ. '') f_JStetr = iNULL__
if(ar__1(225:225) .EQ. '') f_JSteti = iNULL__
if(ar__1(227:232) .EQ. '') JStetr = rNULL__
if(ar__1(234:239) .EQ. '') JSteti = rNULL__
if(ar__1(241:244) .EQ. '') Jthrr = rNULL__
if(ar__1(246:249) .EQ. '') Jthri = rNULL__
if(ar__1(258:258) .EQ. '') f_sdr = iNULL__
if(ar__1(260:260) .EQ. '') f_sdi = iNULL__
if(ar__1(262:267) .EQ. '') sdr = rNULL__
if(ar__1(269:274) .EQ. '') sdi = rNULL__
if(ar__1(276:281) .EQ. '') sdthrr = rNULL__
if(ar__1(283:288) .EQ. '') sdthri = rNULL__
if(ar__1(292:297) .EQ. '') VSXOID = iNULL__
c ..............Just test output...........
write(6,2)
+ Field_1,ID_1,RAdeg_1,DEdeg_1,RArdeg,DErdeg,RAideg,DEideg,
+ Grade,o_rmag,o_imag,rmag,s_rmag,e_rmag,imag,s_imag,e_imag,
+ qflagr,qflagi,AmpMax,f_Ampr,f_Ampi,Ampr,Ampi,athrr,athri,
+ JStetmax,f_JStetr,f_JSteti,JStetr,JSteti,Jthrr,Jthri,sdMax,
+ f_sdr,f_sdi,sdr,sdi,sdthrr,sdthri,VSXMatch,VSXOID
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,DEdeg_2,Filter_1,mag_1,e_mag_1,Amp,qflag,N
if(ar__2(25:30) .EQ. '') mag_1 = rNULL__
if(ar__2(32:36) .EQ. '') e_mag_1 = rNULL__
if(ar__2(38:42) .EQ. '') Amp = rNULL__
if(ar__2(44:45) .EQ. '') qflag = iNULL__
if(ar__2(47:49) .EQ. '') N = iNULL__
c ..............Just test output...........
write(6,3)
+ RAdeg_2,DEdeg_2,Filter_1,mag_1,e_mag_1,Amp,qflag,N
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end