Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/708/560 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/ApJ/708/560 Spectroscopy of UMa II and Coma Ber (Frebel+, 2010) *================================================================================ *High-resolution spectroscopy of extremely metal-poor stars in the least evolved *galaxies: Ursa Major II and Coma Berenices. * Frebel A., Simon J.D., Geha M., Willman B. * <Astrophys. J., 708, 560-583 (2010)> * =2010ApJ...708..560F C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'stars.dat' ! Photometry and other stellar parameters (table 1 and 5 of the paper) integer*4 nr__ parameter (nr__=8) ! Number of records character*88 ar__ ! Full-size record character*14 Gal (nr__) ! Dwarf galaxy name (Coma Berenices or * Ursa Major II) character*14 SDSS (nr__) ! Star SDSS designation (JHHMMSS+DDMMSS) character*9 Name (nr__) ! Star name character*1 f_Name (nr__) ! [f] Non-member star (1) real*4 Vmag (nr__) ! (mag) ? V-band magnitude (2) real*4 rmag (nr__) ! (mag) ? SDSS-DR7 r-band magnitude (2) real*4 g_r (nr__) ! (mag) ? SDSS-DR7 g-r color index (2) real*4 rMag_1 (nr__) ! (mag) ? Absolute r-band magnitude (2) integer*4 Teff (nr__) ! (K) Effective temperature (3) real*4 logg (nr__) ! ([cm/s2]) Surface gravity real*4 v_Fe_H_ (nr__) ! ([-]) Fe/H abundance real*4 Vt (nr__) ! (km/s) Microturbulent velocity {xi} *Note (1): * f = UMa II-NM was observed but was later determined not to be a member of * UMa II. *Note (2): The r magnitudes and g-r colors were obtained from the updated SDSS * DR7 (Abazajian et al. 2009, Cat. II/294), while the V magnitudes were * determined from g and r using the conversions given by Smith et al. * (2002AJ....123.2121S). *Note (3): Temperature are rounded to the nearest 10K. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Observing details integer*4 nr__1 parameter (nr__1=8) ! Number of records character*92 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__1) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*9 Name_1 (nr__1) ! Star name integer*4 RAh (nr__1) ! (h) Hour of right ascension (J2000) integer*4 RAm (nr__1) ! (min) Minute of right ascension (J2000) real*4 RAs (nr__1) ! (s) Second of right ascension (J2000) character*1 DE_ (nr__1) ! Sign of declination (J2000) integer*4 DEd (nr__1) ! (deg) Degree of declination (J2000) integer*4 DEm (nr__1) ! (arcmin) Arcminute of declination (J2000) real*4 DEs (nr__1) ! (arcsec) Arcsecond of declination (J2000) real*8 JD (nr__1) ! (d) Julian Date of observation real*4 Exp (nr__1) ! (h) Exposure time integer*4 S_N5 (nr__1) ! 5000{AA} signal to noise ratio (4) integer*4 S_N6 (nr__1) ! 6000{AA} signal to noise ratio (4) integer*4 S_N65 (nr__1) ! 6500{AA} signal to noise ratio (4) real*4 RVH (nr__1) ! (km/s) HIRES radial velocity real*4 e_RVH (nr__1) ! (km/s) RVH uncertainty real*4 RVD (nr__1) ! (km/s) ? DEIMOS radial velocity real*4 e_RVD (nr__1) ! (km/s) ? RVD uncertainty character*8 Comm (nr__1) ! Comment *Note (4): The S/N measurements are for ~22m{AA} ("blue" CCD), ~26m{AA} * ("green" CCD), and ~28m{AA} ("green" CCD) pixel sizes, respectively. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'ew.dat' ! Equivalent width measurements (table 4) integer*4 nr__2 parameter (nr__2=2191) ! Number of records character*52 ar__2 ! Full-size record character*9 Name_2 (nr__2) ! Star name character*4 Ion (nr__2) ! Species identification real*8 lambda (nr__2) ! (0.1nm) Wavelength in Angstroms real*4 ExPot (nr__2) ! (eV) ? Excitation potential real*8 log_gf (nr__2) ! ([-]) ? Log of the oscillator strength character*1 l_EW (nr__2) ! Limit flag on EW real*4 EW (nr__2) ! (0.1pm) ? Equivalent width in milli-Angstroms character*1 f_EW (nr__2) ! [s] s: indicates a spectrum synthesis measurement character*1 l_log_e (nr__2) ! Limit flag on log(e) real*4 log_e (nr__2) ! ([-]) ? log{epsilon} abundance (H=12) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table6.dat' ! HIRES abundances of the UMa II stars integer*4 nr__3 parameter (nr__3=84) ! Number of records character*47 ar__3 ! Full-size record character*9 Name_3 (nr__3) ! Star name character*5 Ion_1 (nr__3) ! Species identification character*1 l_log_e_1 (nr__3) ! Limit flag on loge real*4 log_e_1 (nr__3) ! ([-]) log{epsilon} abundance (X) character*1 f_log_e (nr__3) ! [a] distorted lines (1) character*1 l__X_H_ (nr__3) ! Limit flag on [X/H] real*4 v_X_H_ (nr__3) ! ([Sun]) X/H abundance character*1 l__X_Fe_ (nr__3) ! Limit flag on [X/Fe] real*4 v_X_Fe_ (nr__3) ! ([Sun]) ? X/Fe abundance (G1) integer*4 N (nr__3) ! Number of spectra used real*4 Sig (nr__3) ! ([-]) ? Abundance uncertainty *Note (1): * a = The Zn lines are somewhat distorted; this may lead to an overestimated * Zn abundance. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table7.dat' ! HIRES abundances of the Coma Ber stars integer*4 nr__4 parameter (nr__4=84) ! Number of records character*47 ar__4 ! Full-size record character*9 Name_4 (nr__4) ! Star name character*5 Ion_2 (nr__4) ! Species identification character*1 l_log_e_1 (nr__4) ! Limit flag on loge real*4 log_e_1 (nr__4) ! ([-]) log{epsilon} abundance (X) character*1 f_log_e_1 (nr__4) ! [a] distorted lines (1) character*1 l__X_H__1 (nr__4) ! Limit flag on [X/H] real*4 v_X_H__1 (nr__4) ! ([Sun]) X/H abundance character*1 l__X_Fe__1 (nr__4) ! Limit flag on [X/Fe] real*4 v_X_Fe__1 (nr__4) ! ([Sun]) ? X/Fe abundance (G1) integer*4 N_1 (nr__4) ! Number of spectra used real*4 Sig_1 (nr__4) ! ([-]) ? Abundance uncertainty *Note (1): * a = The Zn lines are somewhat distorted; this may lead to an overestimated * Zn abundance. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table9.dat' ! HIRES abundances of standard star HD 122563 integer*4 nr__5 parameter (nr__5=54) ! Number of records character*54 ar__5 ! Full-size record character*9 Name_5 (nr__5) ! Star name character*5 Ion_3 (nr__5) ! Species identification real*4 log_e_X0 (nr__5) ! ([-]) ? Log {epsilon} abundance (X)_0_ character*1 l_log_e_1 (nr__5) ! Limit flag on loge real*4 log_e_1 (nr__5) ! ([-]) ? Log {epsilon} abundance (X) character*1 f_log_e_1 (nr__5) ! [b] neutron capture element (1) character*1 l__X_H__2 (nr__5) ! Limit flag on [X/H] real*4 v_X_H__2 (nr__5) ! ([Sun]) ? X/H abundance character*1 l__X_Fe__2 (nr__5) ! Limit flag on [X/Fe] real*4 v_X_Fe__2 (nr__5) ! ([Sun]) ? X/Fe abundance (G1) integer*4 N_2 (nr__5) ! ? Number of spectra used real*4 Sig_2 (nr__5) ! ([-]) ? Abundance uncertainty character*1 ref (nr__5) ! [cd] source of the data: c=this paper, * d=Aoki et al. (2007ApJ...660..747A) *Note (1): b = Neutron-capture abundances are taken from Honda et al. * (2006ApJ...643.1180H). C============================================================================= C Loading file 'stars.dat' ! Photometry and other stellar parameters (table 1 * and 5 of the paper) C Format for file interpretation 1 format( + A14,6X,A14,1X,A9,1X,A1,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,I4, + 1X,F3.1,1X,F5.2,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'stars.dat') write(6,*) '....Loading file: stars.dat' do i__=1,8 read(1,'(A88)')ar__ read(ar__,1) + Gal(i__),SDSS(i__),Name(i__),f_Name(i__),Vmag(i__),rmag(i__), + g_r(i__),rMag_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),Vt(i__) if(ar__(48:52) .EQ. '') Vmag(i__) = rNULL__ if(ar__(54:58) .EQ. '') rmag(i__) = rNULL__ if(ar__(60:63) .EQ. '') g_r(i__) = rNULL__ if(ar__(65:69) .EQ. '') rMag_1(i__) = rNULL__ c ..............Just test output........... write(6,1) + Gal(i__),SDSS(i__),Name(i__),f_Name(i__),Vmag(i__),rmag(i__), + g_r(i__),rMag_1(i__),Teff(i__),logg(i__),v_Fe_H_(i__),Vt(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Observing details C Format for file interpretation 2 format( + A9,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F5.2,1X,F9.1,1X,F4.2, + 1X,I3,1X,I3,1X,I3,1X,F6.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,A8) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,8 read(1,'(A92)')ar__1 read(ar__1,2) + Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),JD(i__),Exp(i__),S_N5(i__),S_N6(i__), + S_N65(i__),RVH(i__),e_RVH(i__),RVD(i__),e_RVD(i__),Comm(i__) if(ar__1(74:79) .EQ. '') RVD(i__) = rNULL__ if(ar__1(81:83) .EQ. '') e_RVD(i__) = rNULL__ 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(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. 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,2) + Name_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),DEd(i__), + DEm(i__),DEs(i__),JD(i__),Exp(i__),S_N5(i__),S_N6(i__), + S_N65(i__),RVH(i__),e_RVH(i__),RVD(i__),e_RVD(i__),Comm(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'ew.dat' ! Equivalent width measurements (table 4) C Format for file interpretation 3 format( + A9,1X,A4,1X,F7.2,1X,F4.2,1X,F7.3,1X,A1,F5.1,1X,A1,1X,A1,F6.2) C Effective file loading open(unit=1,status='old',file= +'ew.dat') write(6,*) '....Loading file: ew.dat' do i__=1,2191 read(1,'(A52)')ar__2 read(ar__2,3) + Name_2(i__),Ion(i__),lambda(i__),ExPot(i__),log_gf(i__), + l_EW(i__),EW(i__),f_EW(i__),l_log_e(i__),log_e(i__) if(ar__2(24:27) .EQ. '') ExPot(i__) = rNULL__ if(ar__2(29:35) .EQ. '') log_gf(i__) = rNULL__ if(ar__2(38:42) .EQ. '') EW(i__) = rNULL__ if(ar__2(47:52) .EQ. '') log_e(i__) = rNULL__ c ..............Just test output........... write(6,3) + Name_2(i__),Ion(i__),lambda(i__),ExPot(i__),log_gf(i__), + l_EW(i__),EW(i__),f_EW(i__),l_log_e(i__),log_e(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table6.dat' ! HIRES abundances of the UMa II stars C Format for file interpretation 4 format( + A9,1X,A5,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,F5.2,1X,I3,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table6.dat') write(6,*) '....Loading file: table6.dat' do i__=1,84 read(1,'(A47)')ar__3 read(ar__3,4) + Name_3(i__),Ion_1(i__),l_log_e_1(i__),log_e_1(i__), + f_log_e(i__),l__X_H_(i__),v_X_H_(i__),l__X_Fe_(i__), + v_X_Fe_(i__),N(i__),Sig(i__) if(ar__3(34:38) .EQ. '') v_X_Fe_(i__) = rNULL__ if(ar__3(44:47) .EQ. '') Sig(i__) = rNULL__ c ..............Just test output........... write(6,4) + Name_3(i__),Ion_1(i__),l_log_e_1(i__),log_e_1(i__), + f_log_e(i__),l__X_H_(i__),v_X_H_(i__),l__X_Fe_(i__), + v_X_Fe_(i__),N(i__),Sig(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table7.dat' ! HIRES abundances of the Coma Ber stars C Format for file interpretation 5 format( + A9,1X,A5,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,F5.2,1X,I3,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table7.dat') write(6,*) '....Loading file: table7.dat' do i__=1,84 read(1,'(A47)')ar__4 read(ar__4,5) + Name_4(i__),Ion_2(i__),l_log_e_1(i__),log_e_1(i__), + f_log_e_1(i__),l__X_H__1(i__),v_X_H__1(i__),l__X_Fe__1(i__), + v_X_Fe__1(i__),N_1(i__),Sig_1(i__) if(ar__4(34:38) .EQ. '') v_X_Fe__1(i__) = rNULL__ if(ar__4(44:47) .EQ. '') Sig_1(i__) = rNULL__ c ..............Just test output........... write(6,5) + Name_4(i__),Ion_2(i__),l_log_e_1(i__),log_e_1(i__), + f_log_e_1(i__),l__X_H__1(i__),v_X_H__1(i__),l__X_Fe__1(i__), + v_X_Fe__1(i__),N_1(i__),Sig_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table9.dat' ! HIRES abundances of standard star HD 122563 C Format for file interpretation 6 format( + A9,1X,A5,1X,F4.2,1X,A1,F5.2,1X,A1,1X,A1,F5.2,1X,A1,F5.2,1X,I3, + 1X,F4.2,1X,A1) C Effective file loading open(unit=1,status='old',file= +'table9.dat') write(6,*) '....Loading file: table9.dat' do i__=1,54 read(1,'(A54)')ar__5 read(ar__5,6) + Name_5(i__),Ion_3(i__),log_e_X0(i__),l_log_e_1(i__), + log_e_1(i__),f_log_e_1(i__),l__X_H__2(i__),v_X_H__2(i__), + l__X_Fe__2(i__),v_X_Fe__2(i__),N_2(i__),Sig_2(i__),ref(i__) if(ar__5(17:20) .EQ. '') log_e_X0(i__) = rNULL__ if(ar__5(23:27) .EQ. '') log_e_1(i__) = rNULL__ if(ar__5(32:36) .EQ. '') v_X_H__2(i__) = rNULL__ if(ar__5(39:43) .EQ. '') v_X_Fe__2(i__) = rNULL__ if(ar__5(45:47) .EQ. '') N_2(i__) = iNULL__ if(ar__5(49:52) .EQ. '') Sig_2(i__) = rNULL__ c ..............Just test output........... write(6,6) + Name_5(i__),Ion_3(i__),log_e_X0(i__),l_log_e_1(i__), + log_e_1(i__),f_log_e_1(i__),l__X_H__2(i__),v_X_H__2(i__), + l__X_Fe__2(i__),v_X_Fe__2(i__),N_2(i__),Sig_2(i__),ref(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end