Conversion of standardized ReadMe file for
file /./ftp/cats/VII/20 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-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__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. VII/20 Catalogue of HII Regions (Sharpless 1959) *================================================================================ *Catalogue of HII Regions * Sharpless S. * <Astrophys. J. Suppl. Ser. 4, 257 (1959)> * =1959ApJS....4..257S C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'catalog.dat' ! The Sharpless (Sh 2) Catalogue integer*4 nr__ parameter (nr__=313) ! Number of records character*57 ar__ ! Full-size record C 1900.0 position composed of: RAh RAm RAds DE- DEd DEm DEs real*8 RAdeg (nr__) ! (deg) Right Ascension 1900.0 real*8 DEdeg (nr__) ! (deg) Declination 1900.0 C ---------------------------------- ! (position vector(s) in degrees) integer*4 Sh2 (nr__) ! [1/313]+ Sharpless HII catalog number integer*4 GlLund (nr__) ! (0.1deg) [0/3600[ Galactic longitude based on * Lund pole integer*4 GbLund (nr__) ! (0.1deg) [-900/+900] Galactic latitude based on * Lund pole integer*4 GLon (nr__) ! (0.1deg) [0/3600[ Galactic longitude (1) integer*4 GLat (nr__) ! (0.1deg) [-900/+900] Galactic latitude (1) integer*4 RAh (nr__) ! (h) Hours RA, 1900.0 integer*4 RAm (nr__) ! (min) Minutes RA, 1900.0 integer*4 RAds (nr__) ! (0.1s) Deci-seconds RA, 1900.0 character*1 DE_ (nr__) ! Sign Dec, 1900.0 integer*4 DEd (nr__) ! (deg) Degrees Dec, 1900.0 integer*4 DEm (nr__) ! (arcmin) Minutes Dec, 1900.0 integer*4 DEs (nr__) ! (arcsec) Seconds Dec, 1900.0 integer*4 RA1950h (nr__) ! (h) Right Ascension (hours) 1950 (1) integer*4 RA1950m (nr__) ! (min) Right Ascension (minutes) 1950 (1) integer*4 RA1950ds (nr__) ! (0.1s) Right Ascension (deci-seconds) 1950 (1) character*1 DE1950_ (nr__) ! Declination (sign) 1950 (1) integer*4 DE1950d (nr__) ! (deg) Declination (degrees) 1950 (1) integer*4 DE1950m (nr__) ! (arcmin) Declination (minutes) 1950 (1) integer*4 DE1950s (nr__) ! (arcsec) Declination (seconds) 1950 (1) integer*4 Diam (nr__) ! (arcmin) Maximun angular diameter of H II region integer*4 Form (nr__) ! [1/3] Classification as to form: * 1=circular; 2=elliptical; 3=irregular integer*4 Struct (nr__) ! [1/3] Classification as to structure, * from 1=amorphous to 3=filamentary integer*4 Bright (nr__) ! [1/3] Classification as to brightness, * from 1=faintest to 3=brightest integer*4 Stars (nr__) ! Number of associated stars with the * H II region *Note (1): Not in the original publication C============================================================================= C Loading file 'catalog.dat' ! The Sharpless (Sh 2) Catalogue C Format for file interpretation 1 format( + I4,I4,I4,I4,I4,I2,I2,I3,A1,I2,I2,I2,I2,I2,I3,A1,I2,I2,I2,I4, + I1,I1,I1,I2) C Effective file loading open(unit=1,status='old',file= +'catalog.dat') write(6,*) '....Loading file: catalog.dat' do i__=1,313 read(1,'(A57)')ar__ read(ar__,1) + Sh2(i__),GlLund(i__),GbLund(i__),GLon(i__),GLat(i__),RAh(i__), + RAm(i__),RAds(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__), + RA1950h(i__),RA1950m(i__),RA1950ds(i__),DE1950_(i__), + DE1950d(i__),DE1950m(i__),DE1950s(i__),Diam(i__),Form(i__), + Struct(i__),Bright(i__),Stars(i__) RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAds(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAds(i__)/2400. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,1) + Sh2(i__),GlLund(i__),GbLund(i__),GLon(i__),GLat(i__),RAh(i__), + RAm(i__),RAds(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__), + RA1950h(i__),RA1950m(i__),RA1950ds(i__),DE1950_(i__), + DE1950d(i__),DE1950m(i__),DE1950s(i__),Diam(i__),Form(i__), + Struct(i__),Bright(i__),Stars(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end