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 2024-Mar-29 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