Conversion of standardized ReadMe file for
file /./ftp/cats/J/AZh/79/867 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-18 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/AZh/79/867 Abundances of HII regions in blue galaxies (Holovatyy+, 2002) *================================================================================ *The chemical composition of HII regions in blue compact dwarf galaxies. * Holovatyy V.V., Melekh B.Ya. * <Astron. Zh. 79, 867 (2002)> * =2002AZh....79..867H * =2002ARep...46..779H C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'position.dat' ! HII region positions from Simbad integer*4 nr__ parameter (nr__=41) ! Number of records character*63 ar__ ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg ! (deg) Right Ascension J2000.0 real*8 DEdeg ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*14 HII ! Name of the HII region integer*4 RAh ! (h) Right ascension (J2000.0) integer*4 RAm ! (min) Right ascension (J2000.0) real*4 RAs ! (s) Right ascension (J2000.0) character*1 DE_ ! Declination sign (J2000.0) integer*4 DEd ! (deg) Declination (J2000.0) integer*4 DEm ! (arcmin) Declination (J2000.0) integer*4 DEs ! (arcsec) ? Declination (J2000.0) character*27 Com ! Comments c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Empirical expression for determining the compositions of nebular gas and values of (A/H)calc/(A/H)mod for HII regions in BCD galaxies integer*4 nr__1 parameter (nr__1=17) ! Number of records character*102 ar__1 ! Full-size record character*3 TN ! Analytical expression designation (T1-T17) character*8 r_TN ! Reference numbers, in file "refs.dat" character*83 Expression ! Expression real*4 D_A_H ! Calculated to model (A/H) abundance ratio c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'refs.dat' ! References from table1 integer*4 nr__2 parameter (nr__2=12) ! Number of records character*64 ar__2 ! Full-size record integer*4 RefNo ! Reference number character*60 Text ! Author and reference (bibcode) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Composition of HII regions in BCD galaxies, found using ICFs in table1.dat integer*4 nr__3 parameter (nr__3=41) ! Number of records character*110 ar__3 ! Full-size record integer*4 N ! Running number character*14 HII_1 ! Name of HII region real*4 v_He_H_T1_8 ! (10-2) Abundance [He/H]*100 using T1 and * recombination coefficient from * Brocklehurst, 1972MNRAS.157..211B real*4 v_He_H_T2_8 ! (10-2) ? Abundance [He/H]*100 using T2 and * recombination coefficient from * Brocklehurst, 1972MNRAS.157..211B real*4 v_He_H_T1_9 ! (10-2) Abundance [He/H]*100 using T1 and * recombination coefficient from * Smits, 1996MNRAS.278..683S real*4 v_He_H_T2_9 ! (10-2) ? Abundance [He/H]*100 using T2 and * recombination coefficient from * Smits, 1996MNRAS.278..683S real*4 v_N_H_T6 ! (10-6) ? Abundance [N/H]*10^6^ using T6 real*4 v_O_H_T7 ! (10-5) Abundance [O/H]*10^5^ using T7 real*4 v_O_H_T8 ! (10-5) ? Abundance [O/H]*10^5^ using T8 real*4 v_O_H_T9 ! (10-5) ? Abundance [O/H]*10^5^ using T9 real*4 v_O_H_T10 ! (10-5) Abundance [O/H]*10^5^ using T10 real*4 v_Ne_H_T11 ! (10-6) ? Abundance [Ne/H]*10^6^ using T11 real*4 v_S_H_T13 ! (10-6) ? Abundance [S/H]*10^6^ using T13 real*4 v_Ar_H_T14 ! (10-7) ? Abundance [Ar/H]*10^7^ using T14 real*4 v_Ar_H_T16 ! (10-7) ? Abundance [Ar/H]*10^7^ using T16 real*4 v_Ar_H_T17 ! (10-7) ? Abundance [Ar/H]*10^7^ using T17 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Composition of HII regions in BCD galaxies, found using ICFs from grids of photoionization models from Stasinska, 1990A&AS...83..501S integer*4 nr__4 parameter (nr__4=41) ! Number of records character*104 ar__4 ! Full-size record integer*4 N_1 ! Running number character*14 HII_2 ! Name of HII region real*4 v_He_H_8 ! (10-2) ? Abundance [He/H]*100 using T1 and * recombination coefficient from * Brocklehurst, 1972MNRAS.157..211B real*4 e__He_H_8 ! (10-2) ? Uncertainty in [He/H]8 real*4 v_He_H_9 ! (10-2) ? Abundance [He/H]*100 using recombination * coefficient from Smits, 1996MNRAS.278..683S real*4 e__He_H_9 ! (10-2) ? Uncertainty in [He/H]9 real*4 v_O_H_ ! (10-5) Abundance [O/H]*10^5^, using G3, G4, G26 (1) real*4 e__O_H_ ! (10-5) Uncertainty in [O/H] real*4 v_N_H_ ! (10-6) ? Abundance [N/H]*10^6^, using G6 (1) real*4 e__N_H_ ! (10-6) ? Uncertainty in [N/H] real*4 v_Ne_H_ ! (10-6) ? Abundance [Ne/H]*10^6^, using G29 (1) real*4 e__Ne_H_ ! (10-6) ? Uncertainty in [Ne/H] real*4 v_S_H_ ! (10-6) ? Abundance [S/H]*10^6^, using G30 (1) real*4 e__S_H_ ! (10-6) ? Uncertainty in [S/H] real*4 v_Ar_H_ ! (10-7) ? Abundance [Ar/H]*10^7^, * using G10, G11, G31, G32 (1) real*4 e__Ar_H_ ! (10-7) ? Uncertainty in [Ar/H] *Note (1): G1-G32 are ionization correction factors (ICF) expression of table3 * in Holovatyy and Melekh, Kinematica Fiz. Nebesnykh Tel (2002) (in press). C============================================================================= C Loading file 'position.dat' ! HII region positions from Simbad C Format for file interpretation 1 format(A14,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,A27) C Effective file loading open(unit=1,status='old',file= +'position.dat') write(6,*) '....Loading file: position.dat' do i__=1,41 read(1,'(A63)')ar__ read(ar__,1)HII,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Com if(ar__(34:35) .EQ. '') DEs = iNULL__ RAdeg = rNULL__ DEdeg = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh .GT. -180) RAdeg=RAh*15. if(RAm .GT. -180) RAdeg=RAdeg+RAm/4. if(RAs .GT. -180) RAdeg=RAdeg+RAs/240. if(DEd .GE. 0) DEdeg=DEd if(DEm .GE. 0) DEdeg=DEdeg+DEm/60. if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600. if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg c ..............Just test output........... write(6,1)HII,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Com write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table1.dat' ! Empirical expression for determining the * compositions of nebular gas and values of * (A/H)calc/(A/H)mod for HII regions in BCD * galaxies C Format for file interpretation 2 format(A3,1X,A8,1X,A83,2X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,17 read(1,'(A102)')ar__1 read(ar__1,2)TN,r_TN,Expression,D_A_H c ..............Just test output........... write(6,2)TN,r_TN,Expression,D_A_H c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'refs.dat' ! References from table1 C Format for file interpretation 3 format(I2,2X,A60) C Effective file loading open(unit=1,status='old',file= +'refs.dat') write(6,*) '....Loading file: refs.dat' do i__=1,12 read(1,'(A64)')ar__2 read(ar__2,3)RefNo,Text c ..............Just test output........... write(6,3)RefNo,Text c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Composition of HII regions in BCD galaxies, * found using ICFs in table1.dat C Format for file interpretation 4 format( + I2,1X,A14,2X,F5.2,2X,F5.2,2X,F5.2,2X,F5.2,2X,F4.2,2X,F5.2,2X, + F5.2,2X,F5.2,2X,F5.2,2X,F5.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,41 read(1,'(A110)')ar__3 read(ar__3,4) + N,HII_1,v_He_H_T1_8,v_He_H_T2_8,v_He_H_T1_9,v_He_H_T2_9, + v_N_H_T6,v_O_H_T7,v_O_H_T8,v_O_H_T9,v_O_H_T10,v_Ne_H_T11, + v_S_H_T13,v_Ar_H_T14,v_Ar_H_T16,v_Ar_H_T17 if(ar__3(27:31) .EQ. '') v_He_H_T2_8 = rNULL__ if(ar__3(41:45) .EQ. '') v_He_H_T2_9 = rNULL__ if(ar__3(48:51) .EQ. '') v_N_H_T6 = rNULL__ if(ar__3(61:65) .EQ. '') v_O_H_T8 = rNULL__ if(ar__3(68:72) .EQ. '') v_O_H_T9 = rNULL__ if(ar__3(82:86) .EQ. '') v_Ne_H_T11 = rNULL__ if(ar__3(89:92) .EQ. '') v_S_H_T13 = rNULL__ if(ar__3(95:98) .EQ. '') v_Ar_H_T14 = rNULL__ if(ar__3(101:104) .EQ. '') v_Ar_H_T16 = rNULL__ if(ar__3(107:110) .EQ. '') v_Ar_H_T17 = rNULL__ c ..............Just test output........... write(6,4) + N,HII_1,v_He_H_T1_8,v_He_H_T2_8,v_He_H_T1_9,v_He_H_T2_9, + v_N_H_T6,v_O_H_T7,v_O_H_T8,v_O_H_T9,v_O_H_T10,v_Ne_H_T11, + v_S_H_T13,v_Ar_H_T14,v_Ar_H_T16,v_Ar_H_T17 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Composition of HII regions in BCD galaxies, * found using ICFs from grids of photoionization * models from Stasinska, 1990A&AS...83..501S C Format for file interpretation 5 format( + I2,1X,A14,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F5.2,2X,F4.2,2X, + F4.2,2X,F4.2,2X,F5.2,2X,F4.2,2X,F4.2,2X,F4.2,2X,F4.2,3X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,41 read(1,'(A104)')ar__4 read(ar__4,5) + N_1,HII_2,v_He_H_8,e__He_H_8,v_He_H_9,e__He_H_9,v_O_H_, + e__O_H_,v_N_H_,e__N_H_,v_Ne_H_,e__Ne_H_,v_S_H_,e__S_H_, + v_Ar_H_,e__Ar_H_ if(ar__4(20:23) .EQ. '') v_He_H_8 = rNULL__ if(ar__4(26:29) .EQ. '') e__He_H_8 = rNULL__ if(ar__4(32:35) .EQ. '') v_He_H_9 = rNULL__ if(ar__4(38:41) .EQ. '') e__He_H_9 = rNULL__ if(ar__4(57:60) .EQ. '') v_N_H_ = rNULL__ if(ar__4(63:66) .EQ. '') e__N_H_ = rNULL__ if(ar__4(69:73) .EQ. '') v_Ne_H_ = rNULL__ if(ar__4(76:79) .EQ. '') e__Ne_H_ = rNULL__ if(ar__4(82:85) .EQ. '') v_S_H_ = rNULL__ if(ar__4(88:91) .EQ. '') e__S_H_ = rNULL__ if(ar__4(94:97) .EQ. '') v_Ar_H_ = rNULL__ if(ar__4(101:104) .EQ. '') e__Ar_H_ = rNULL__ c ..............Just test output........... write(6,5) + N_1,HII_2,v_He_H_8,e__He_H_8,v_He_H_9,e__He_H_9,v_O_H_, + e__O_H_,v_N_H_,e__N_H_,v_Ne_H_,e__Ne_H_,v_S_H_,e__S_H_, + v_Ar_H_,e__Ar_H_ c .......End.of.Just test output........... end do close(1) C============================================================================= stop end