Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/229/362 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 2026-Oct-11
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/229/362 Small-scale structures in HI hole of M33 (Deul+, 1990)
*================================================================================
*Small-scale structures in the HI distribution of M 33.
* Deul E.R., den Hartog R.H.
* <Astron. Astrophys. 229, 362 (1990)>
* =1990A&A...229..362D
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Properties of the H I holes in the catalogue
integer*4 nr__
parameter (nr__=148) ! Number of records
character*106 ar__ ! Full-size record
C B1950 position composed of: RAdeg DEdeg
integer*4 v_DD90_ (nr__) ! Sequential number of the hole
real*4 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (B1950)
real*4 DEdeg (nr__) ! (deg) Declination in decimal degrees (B1950)
real*4 R (nr__) ! (kpc) Distance to the HI hole with respect to the
* M33 center (1)
real*4 theta (nr__) ! (deg) Orientation angle of the HI hole in the
* galaxy with respect to the north
* equatorial pole, with positive sense in
* the anti-clockwise direction
integer*4 HRV (nr__) ! (km/s) Radial heliocentric velocity of maximum
* contrast
integer*4 DV (nr__) ! (km/s) Expansion velocity (to an accuracy of 2km/s)
real*4 Cont (nr__) ! Hole contrast defined as
* Cont=I(surr)-I(50%)/I(surr)
integer*4 T (nr__) ! [1/3] Type of Hole (in Fig.1 of paper):
* 1=gap, 2=shift, 3=depression
integer*4 Q (nr__) ! Quality of detection (-1 met for the hole)
character*1 C (nr__) ! [cp] complete/partial classification code (2)
integer*4 Min (nr__) ! (pc) Extent of the minor axis of the ellipse
* that can be fitted to the 50% intensity
* contour of a hole
integer*4 Maj (nr__) ! (pc) Extent of the major axis of the ellipse
* that can be fitted to the 50% intensity
* contour of a hole
integer*4 Diam (nr__) ! (pc) Diameter of the hole, Diam=(Min+Maj)^1/2^
real*4 eps (nr__) ! Eccentricity of the hole
integer*4 PA (nr__) ! (deg) [-90/90] Position angle of the major axis
* of the ellipse fitted
real*4 Age (nr__) ! (Gyr) Kinematic age assuming a constant
* expansion velocity
real*4 nHI (nr__) ! (cm-3) Ambiant density derived from the observed
* total column densities for the HI material
* surrounding the hole
real*4 Mass (nr__) ! (10+4Msun) Mass assuming that nHi equals the pre-hole
* density and a spherical hole.
real*4 Energy (nr__) ! (10+43J) Energy needed to create the hole (x10^50^ergs)
*Note (1): We have used the central position from table1 of Deul and
* van der Hulst (1987A&AS...67..509D)
*Note (2): Classification code as follows:
* c = holes that show a clear elliptical shape over more than 90% of
* the outline
* p = partial
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Overlap properties of the H I holes with respect
to the H II regions and OB associations
integer*4 nr__1
parameter (nr__1=202) ! Number of records
character*28 ar__1 ! Full-size record
integer*4 v_DD90__1 (nr__1) ! Sequential number of the hole
character*6 HII (nr__1) ! Name of the HII region, CPSDP, Courtes et al.,
* 1987, Cat. <J/A+AS/131/317>
integer*4 EM (nr__1) ! [1-4]? Emission measure code (1)
integer*4 HI (nr__1) ! [1-3]? HI-HII code (2)
integer*4 OB (nr__1) ! ? OB association number, [HS80], Humphreys &
* Sandage, 1980ApJS...44..319H
integer*4 HI_OB (nr__1) ! [1-3]? HI-OB code (2)
integer*4 HII_OB (nr__1) ! [1-3]? HII-OB code (2)
*Note (1): The EM code is taken from the H II region catalogue of Boulesteix
* et al. (1974A&A....37...33B). 1 (faint) to 4 (extremely strong) on a
* logarithmic scale
*Note (2): The HI-HII, HI-OB, and HII-OB overlap codes are defined as follows :
* 1 = 100% overlap (e.g. generally the HII regions lies within the HI hole)
* 2 = between 100% and 50% overlap
* 3 = between 50% and no overlap
C=============================================================================
C Loading file 'table1.dat' ! Properties of the H I holes in the catalogue
C Format for file interpretation
1 format(
+ I3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.1,1X,I4,1X,I3,1X,F4.2,1X,I2,
+ 1X,I2,2X,A1,2X,I3,1X,I4,1X,I4,2X,F5.3,1X,I3,2X,F4.1,1X,F5.2,
+ 1X,F6.1,2X,F6.1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,148
read(1,'(A106)')ar__
read(ar__,1)
+ v_DD90_(i__),RAdeg(i__),DEdeg(i__),R(i__),theta(i__),HRV(i__),
+ DV(i__),Cont(i__),T(i__),Q(i__),C(i__),Min(i__),Maj(i__),
+ Diam(i__),eps(i__),PA(i__),Age(i__),nHI(i__),Mass(i__),
+ Energy(i__)
c ..............Just test output...........
write(6,1)
+ v_DD90_(i__),RAdeg(i__),DEdeg(i__),R(i__),theta(i__),HRV(i__),
+ DV(i__),Cont(i__),T(i__),Q(i__),C(i__),Min(i__),Maj(i__),
+ Diam(i__),eps(i__),PA(i__),Age(i__),nHI(i__),Mass(i__),
+ Energy(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Overlap properties of the H I holes with respect
* to the H II regions and OB associations
C Format for file interpretation
2 format(I3,2X,A6,2X,I1,2X,I1,2X,I3,2X,I1,2X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,202
read(1,'(A28)')ar__1
read(ar__1,2)
+ v_DD90__1(i__),HII(i__),EM(i__),HI(i__),OB(i__),HI_OB(i__),
+ HII_OB(i__)
if(ar__1(14:14) .EQ. '') EM(i__) = iNULL__
if(ar__1(17:17) .EQ. '') HI(i__) = iNULL__
if(ar__1(20:22) .EQ. '') OB(i__) = iNULL__
if(ar__1(25:25) .EQ. '') HI_OB(i__) = iNULL__
if(ar__1(28:28) .EQ. '') HII_OB(i__) = iNULL__
c ..............Just test output...........
write(6,2)
+ v_DD90__1(i__),HII(i__),EM(i__),HI(i__),OB(i__),HI_OB(i__),
+ HII_OB(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end