Conversion of standardized ReadMe file for
file /./ftp/cats/I/284 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.5, on 2013-Jun-19
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. I/284 The USNO-B1.0 Catalog (Monet+ 2003)
*================================================================================
*The USNO-B Catalog
* Monet D.G., Levine S.E., Casian B., et al.
* <Astron. J., 125, 984 (2003)>
* =2003AJ....125..984M
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'out.sam' ! Sample of data (around 0+0)
integer*4 nr__
parameter (nr__=100) ! Number of records
character*250 ar__ ! Full-size record
character*12 USNO_B1_0 (nr__) ! Designation of the object (1)
character*12 Tycho_2 (nr__) ! Designation in the Tycho-2 Catalog <I/259>
real*8 RAdeg (nr__) ! (deg) Right Ascension at Eq=J2000, Ep=J2000 (2)
real*8 DEdeg (nr__) ! (deg) Declination at Eq=J2000, Ep=J2000 (2)
integer*4 e_RAdeg (nr__) ! (mas) Mean error on RAdeg*cos(DEdeg) at Epoch
integer*4 e_DEdeg (nr__) ! (mas) Mean error on DEdeg at Epoch
real*4 Epoch (nr__) ! (yr) Mean epoch of observation (2)
integer*4 pmRA (nr__) ! (mas/yr) Proper motion in RA (relative to YS4.0)
integer*4 pmDE (nr__) ! (mas/yr) Proper motion in DE (relative to YS4.0)
integer*4 muPr (nr__) ! (0.1) ? Total Proper Motion probability (7)
integer*4 e_pmRA (nr__) ! (mas/yr) Mean error on pmRA
integer*4 e_pmDE (nr__) ! (mas/yr) Mean error on pmDE
integer*4 fit_RA (nr__) ! (100mas) Mean error on RA fit
integer*4 fit_DE (nr__) ! (100mas) Mean error on DE fit
integer*4 Ndet (nr__) ! [0,5] Number of detections (7)
character*3 Flags (nr__) ! [MsY.] Flags on object (3)
real*4 B1mag (nr__) ! (mag) ? First blue magnitude
integer*4 B1C (nr__) ! ? source of photometric calibration (4)
integer*4 B1S (nr__) ! ? Survey number (see "Surveys" section below)
integer*4 B1f (nr__) ! ? Field number in survey
integer*4 B1s_g (nr__) ! ? Star-galaxy separation (6)
real*4 B1xi (nr__) ! (arcsec) ? Residual in X direction (5)
real*4 B1eta (nr__) ! (arcsec) ? Residual in Y direction (5)
real*4 R1mag (nr__) ! (mag) ? First red magnitude
integer*4 R1C (nr__) ! ? source of photometric calibration (4)
integer*4 R1S (nr__) ! ? Survey number (see "Surveys" section below)
integer*4 R1f (nr__) ! ? Field number in survey
integer*4 R1s_g (nr__) ! ? Star-galaxy separation (6)
real*4 R1xi (nr__) ! (arcsec) ? Residual in X direction (5)
real*4 R1eta (nr__) ! (arcsec) ? Residual in Y direction (5)
real*4 B2mag (nr__) ! (mag) ? Second blue magnitude
integer*4 B2C (nr__) ! ? source of photometric calibration (4)
integer*4 B2S (nr__) ! ? Survey number (see "Surveys" section below)
integer*4 B2f (nr__) ! ? Field number in survey
integer*4 B2s_g (nr__) ! ? Star-galaxy separation (6)
real*4 B2xi (nr__) ! (arcsec) ? Residual in X direction (5)
real*4 B2eta (nr__) ! (arcsec) ? Residual in Y direction (5)
real*4 R2mag (nr__) ! (mag) ? Second red magnitude
integer*4 R2C (nr__) ! ? source of photometric calibration (4)
integer*4 R2S (nr__) ! ? Survey number (see "Surveys" section below)
integer*4 R2f (nr__) ! ? Field number in survey
integer*4 R2s_g (nr__) ! ? Star-galaxy separation (6)
real*4 R2xi (nr__) ! (arcsec) ? Residual in X direction (5)
real*4 R2eta (nr__) ! (arcsec) ? Residual in Y direction (5)
real*4 Imag (nr__) ! (mag) ? Infrared (N) magnitude
integer*4 IC (nr__) ! ? source of photometric calibration (4)
integer*4 IS (nr__) ! ? Survey number (see "Surveys" section below)
integer*4 If (nr__) ! ? Field number in survey
integer*4 Is_g (nr__) ! ? Star-galaxy separation (6)
real*4 Ixi (nr__) ! (arcsec) ? Residual in X direction (5)
real*4 Ieta (nr__) ! (arcsec) ? Residual in Y direction (5)
*Note (1): The USNO-B1.0 is arranged in zones of 0.1deg in Declination,
* with objects ordered by Right Ascension in each zone.
* The USNO-B1.0 is made of the zone number (from 0000 in South Pole
* to 1799 in North Pole), followed by a sequential number.
*Note (2): The proper motion was applied to compute the RAdeg and DEdeg
* values. Please note that the large uncertainties in the proper motions
* mean that the RAdeg and DEdeg are less accurate than the mean errors
* e_RAdeg and e_DEdeg which apply to the position at the mean Epoch.
*Note (3): . denotes the absence of any flag
* M = Existence in a proper motion catalog,
* s = object on a diffraction spike
* Y = Correlation with YS4.0 catalog (Monet, in prep.)
*Note (4): the photometric calibration is represented by a number:
* 0 = bright photometric standards on the plate
* 1 = faint photometric standard on the plate
* 2 = faint photometric standard one plate away (on overlap plate)
* 3 = faint photometric standard two plate away (on overlap of overlap)
* etc
*Note (5): Distance, along the x- and y- direction, of the object position
* compared to the mean epoch.
*Note (6): The star/galaxy separation is a measure of the similarity
* of the point-spread function to a stellar profile:
* 0 means quite dissimilar -- i.e. a non-stellar object
* 11 means quite similar -- i.e. a stellar object
*Note (7): For Tycho-2 stars, the Total Proper Motion probability is
* not given, and the number of detections Ndet is set to zero.
* For other stars, Ndet is 2 or more.
C=============================================================================
C Loading file 'out.sam' ! Sample of data (around 0+0)
C Format for file interpretation
1 format(
+ A12,1X,A12,1X,F10.6,F10.6,1X,I3,1X,I3,1X,F6.1,1X,I6,1X,I6,1X,
+ I1,1X,I3,1X,I3,1X,I1,1X,I1,1X,I1,1X,A3,2X,F5.2,1X,I1,1X,I1,1X,
+ I3,1X,I2,1X,F6.2,F6.2,2X,F5.2,1X,I1,1X,I1,1X,I3,1X,I2,1X,F6.2,
+ F6.2,2X,F5.2,1X,I1,1X,I1,1X,I3,1X,I2,1X,F6.2,F6.2,2X,F5.2,1X,
+ I1,1X,I1,1X,I3,1X,I2,1X,F6.2,F6.2,2X,F5.2,1X,I1,1X,I1,1X,I3,
+ 1X,I2,1X,F6.2,F6.2)
C Effective file loading
open(unit=1,file='out.sam', status='old')
write(6,*) '....Loading file: out.sam'
do i__=1,100
read(1,'(A250)')ar__
read(ar__,1)
+ USNO_B1_0(i__),Tycho_2(i__),RAdeg(i__),DEdeg(i__),
+ e_RAdeg(i__),e_DEdeg(i__),Epoch(i__),pmRA(i__),pmDE(i__),
+ muPr(i__),e_pmRA(i__),e_pmDE(i__),fit_RA(i__),fit_DE(i__),
+ Ndet(i__),Flags(i__),B1mag(i__),B1C(i__),B1S(i__),B1f(i__),
+ B1s_g(i__),B1xi(i__),B1eta(i__),R1mag(i__),R1C(i__),R1S(i__),
+ R1f(i__),R1s_g(i__),R1xi(i__),R1eta(i__),B2mag(i__),B2C(i__),
+ B2S(i__),B2f(i__),B2s_g(i__),B2xi(i__),B2eta(i__),R2mag(i__),
+ R2C(i__),R2S(i__),R2f(i__),R2s_g(i__),R2xi(i__),R2eta(i__),
+ Imag(i__),IC(i__),IS(i__),If(i__),Is_g(i__),Ixi(i__),Ieta(i__)
if(ar__(77:77) .EQ. '') muPr(i__) = iNULL__
if(ar__(98:102) .EQ. '') B1mag(i__) = rNULL__
if(ar__(104:104) .EQ. '') B1C(i__) = iNULL__
if(ar__(106:106) .EQ. '') B1S(i__) = iNULL__
if(ar__(108:110) .EQ. '') B1f(i__) = iNULL__
if(ar__(112:113) .EQ. '') B1s_g(i__) = iNULL__
if(ar__(115:120) .EQ. '') B1xi(i__) = rNULL__
if(ar__(121:126) .EQ. '') B1eta(i__) = rNULL__
if(ar__(129:133) .EQ. '') R1mag(i__) = rNULL__
if(ar__(135:135) .EQ. '') R1C(i__) = iNULL__
if(ar__(137:137) .EQ. '') R1S(i__) = iNULL__
if(ar__(139:141) .EQ. '') R1f(i__) = iNULL__
if(ar__(143:144) .EQ. '') R1s_g(i__) = iNULL__
if(ar__(146:151) .EQ. '') R1xi(i__) = rNULL__
if(ar__(152:157) .EQ. '') R1eta(i__) = rNULL__
if(ar__(160:164) .EQ. '') B2mag(i__) = rNULL__
if(ar__(166:166) .EQ. '') B2C(i__) = iNULL__
if(ar__(168:168) .EQ. '') B2S(i__) = iNULL__
if(ar__(170:172) .EQ. '') B2f(i__) = iNULL__
if(ar__(174:175) .EQ. '') B2s_g(i__) = iNULL__
if(ar__(177:182) .EQ. '') B2xi(i__) = rNULL__
if(ar__(183:188) .EQ. '') B2eta(i__) = rNULL__
if(ar__(191:195) .EQ. '') R2mag(i__) = rNULL__
if(ar__(197:197) .EQ. '') R2C(i__) = iNULL__
if(ar__(199:199) .EQ. '') R2S(i__) = iNULL__
if(ar__(201:203) .EQ. '') R2f(i__) = iNULL__
if(ar__(205:206) .EQ. '') R2s_g(i__) = iNULL__
if(ar__(208:213) .EQ. '') R2xi(i__) = rNULL__
if(ar__(214:219) .EQ. '') R2eta(i__) = rNULL__
if(ar__(222:226) .EQ. '') Imag(i__) = rNULL__
if(ar__(228:228) .EQ. '') IC(i__) = iNULL__
if(ar__(230:230) .EQ. '') IS(i__) = iNULL__
if(ar__(232:234) .EQ. '') If(i__) = iNULL__
if(ar__(236:237) .EQ. '') Is_g(i__) = iNULL__
if(ar__(239:244) .EQ. '') Ixi(i__) = rNULL__
if(ar__(245:250) .EQ. '') Ieta(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ USNO_B1_0(i__),Tycho_2(i__),RAdeg(i__),DEdeg(i__),
+ e_RAdeg(i__),e_DEdeg(i__),Epoch(i__),pmRA(i__),pmDE(i__),
+ muPr(i__),e_pmRA(i__),e_pmDE(i__),fit_RA(i__),fit_DE(i__),
+ Ndet(i__),Flags(i__),B1mag(i__),B1C(i__),B1S(i__),B1f(i__),
+ B1s_g(i__),B1xi(i__),B1eta(i__),R1mag(i__),R1C(i__),R1S(i__),
+ R1f(i__),R1s_g(i__),R1xi(i__),R1eta(i__),B2mag(i__),B2C(i__),
+ B2S(i__),B2f(i__),B2s_g(i__),B2xi(i__),B2eta(i__),R2mag(i__),
+ R2C(i__),R2S(i__),R2f(i__),R2s_g(i__),R2xi(i__),R2eta(i__),
+ Imag(i__),IC(i__),IS(i__),If(i__),Is_g(i__),Ixi(i__),Ieta(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end