Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/503/913 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.5, on 2013-May-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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/A+A/503/913 Synthetic spectrophotometry for C-rich giants (Aringer+, 2009)
*================================================================================
*Synthetic photometry for carbon rich giants. I. Hydrostatic dust-free models.
* Aringer B., Girardi L., Nowotny W., Marigo P., Lederer M.T.
* <Astron. Astrophys., 503, 913-928 (2009)>
* =2009A&A...503..913A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'photom.dat' ! List of photometric system files
integer*4 nr__
parameter (nr__=36) ! Number of records
character*141 ar__ ! Full-size record
character*15 FileName (nr__) ! Name of the file in subdirectory phot
character*125 Filt (nr__) ! Filters for which bolometric corrections are
* presented in the file
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'spectra.dat' ! List of spectra
integer*4 nr__1
parameter (nr__1=746) ! Number of records
character*89 ar__1 ! Full-size record
integer*4 Mod (nr__1) ! [1/746] Model number
integer*4 Ver (nr__1) ! [3/4] Version of COMARCS model
integer*4 Teff (nr__1) ! (K) [2400/4000] Effective temperature
real*4 logg (nr__1) ! ([cm/s2]) [-1.00/+0.00] Surface gravity
real*4 Mass (nr__1) ! (solMass) [1/99] Mass
real*4 Z (nr__1) ! (Sun) (0.10, 0.33, 1.00) Metallicity
* (in solar metallicity)
real*4 C_O (nr__1) ! [1.05/5.00] C/O abundance ratio
real*4 xi (nr__1) ! (km/s) [2.5] Microturbulent velocity
character*49 SpFileName (nr__1) ! Name of the spectrum in subdirectory spec
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'phot/*' ! Individual photometric files
with bolometric corrections
integer*4 nr__2
parameter (nr__2=36) ! Number of records
character*368 ar__2 ! Full-size record
integer*4 Mod_1 (nr__2) ! [1/746] Model number
integer*4 Teff_1 (nr__2) ! (K) Effective temperature
real*4 logg_1 (nr__2) ! ([cm/s2]) Surface gravity
real*4 Mass_1 (nr__2) ! (solMass) Mass
real*4 Z_1 (nr__2) ! ([Sun]) Metallicity (in solar metallicity)
real*4 C_O_1 (nr__2) ! C/O abundance ratio
real*4 logL (nr__2) ! ([solLum]) Luminosity
real*4 xi_1 (nr__2) ! (km/s) microturbulent velocity
integer*4 Ver_1 (nr__2) ! [3/4] Version of COMARCS model
character*8 Band1 (nr__2) ! Band name
real*8 BC1 (nr__2) ! (mag) ?=- Bolometric correction in Band1 filter (1)
character*8 Band2 (nr__2) ! Band name
real*8 BC2 (nr__2) ! (mag) ?=- Bolometric correction in Band2 filter (1)
character*8 Band3 (nr__2) ! Band name
real*8 BC3 (nr__2) ! (mag) ?=- Bolometric correction in Band3 filter (1)
character*8 Band4 (nr__2) ! Band name
real*8 BC4 (nr__2) ! (mag) ?=- Bolometric correction in Band4 filter (1)
character*8 Band5 (nr__2) ! Band name
real*8 BC5 (nr__2) ! (mag) ?=- Bolometric correction in Band5 filter (1)
character*8 Band6 (nr__2) ! Band name
real*8 BC6 (nr__2) ! (mag) ?=- Bolometric correction in Band6 filter (1)
character*8 Band7 (nr__2) ! Band name
real*8 BC7 (nr__2) ! (mag) ?=- Bolometric correction in Band7 filter (1)
character*8 Band8 (nr__2) ! Band name
real*8 BC8 (nr__2) ! (mag) ?=- Bolometric correction in Band8 filter (1)
character*8 Band9 (nr__2) ! Band name
real*8 BC9 (nr__2) ! (mag) ?=- Bolometric correction in Band9 filter (1)
character*8 Band10 (nr__2) ! Band name
real*8 BC10 (nr__2) ! (mag) ?=- Bolometric correction in Band10 filter (1)
character*9 Band11 (nr__2) ! Band name
real*8 BC11 (nr__2) ! (mag) ?=- Bolometric correction in Band11 filter (1)
character*9 Band12 (nr__2) ! Band name
real*8 BC12 (nr__2) ! (mag) ?=- Bolometric correction in Band12 filter (1)
character*9 Band13 (nr__2) ! Band name
real*8 BC13 (nr__2) ! (mag) ?=- Bolometric correction in Band13 filter (1)
character*8 Band14 (nr__2) ! Band name
real*8 BC14 (nr__2) ! (mag) ?=- Bolometric correction in Band14 filter (1)
character*8 Band15 (nr__2) ! Band name
real*8 BC15 (nr__2) ! (mag) ?=- Bolometric correction in Band15 filter (1)
character*8 Band16 (nr__2) ! Band name
real*8 BC16 (nr__2) ! (mag) ?=- Bolometric correction in Band16 filter (1)
character*8 Band17 (nr__2) ! Band name
real*8 BC17 (nr__2) ! (mag) ?=- Bolometric correction in Band17 filter (1)
*Note (1): If a certain filter curve is not covered by the synthetic spectra
* (below 4444 and above 250000{AA}) a value of --- is listed for the
* bolometric correction.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'spec/*' ! Individual synthetic spectra
integer*4 nr__3
parameter (nr__3=746) ! Number of records
character*31 ar__3 ! Full-size record
real*8 lambda (nr__3) ! (0.1nm) Wavelength
real*4 nFlux (nr__3) ! [0/2] Continuum normalized flux (1)
real*4 nuL_nu (nr__3) ! (10-7W) Frequency times specific luminosity
* {nu}.L_{nu}_
*Note (1): Normalized to calculation without atomic and molecular line opacities
C=============================================================================
C Loading file 'photom.dat' ! List of photometric system files
C Format for file interpretation
1 format(A15,1X,A125)
C Effective file loading
open(unit=1,file='photom.dat', status='old')
write(6,*) '....Loading file: photom.dat'
do i__=1,36
read(1,'(A141)')ar__
read(ar__,1)FileName(i__),Filt(i__)
c ..............Just test output...........
write(6,1)FileName(i__),Filt(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'spectra.dat' ! List of spectra
C Format for file interpretation
2 format(
+ I4,2X,I1,1X,I4,1X,F5.2,1X,F5.2,1X,F4.2,2X,F4.2,1X,F3.1,1X,A49)
C Effective file loading
open(unit=1,file='spectra.dat', status='old')
write(6,*) '....Loading file: spectra.dat'
do i__=1,746
read(1,'(A89)')ar__1
read(ar__1,2)
+ Mod(i__),Ver(i__),Teff(i__),logg(i__),Mass(i__),Z(i__),
+ C_O(i__),xi(i__),SpFileName(i__)
c ..............Just test output...........
write(6,2)
+ Mod(i__),Ver(i__),Teff(i__),logg(i__),Mass(i__),Z(i__),
+ C_O(i__),xi(i__),SpFileName(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'phot/*' ! Individual photometric files
* with bolometric corrections
C Format for file interpretation
3 format(
+ I4,5X,I4,1X,F5.2,1X,F4.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.1,2X,
+ I1,A8,1X,F7.3,2X,A8,2X,F7.3,2X,A8,2X,F7.3,2X,A8,2X,F7.3,2X,A8,
+ 2X,F7.3,2X,A8,2X,F7.3,2X,A8,1X,F7.3,2X,A8,2X,F7.3,2X,A8,2X,
+ F7.3,2X,A8,2X,F7.3,1X,A9,2X,F7.3,1X,A9,2X,F7.3,1X,A9,2X,F7.3,
+ 2X,A8,2X,F7.3,2X,A8,2X,F7.3,2X,A8,2X,F7.3,2X,A8,1X,F7.3)
C Effective file loading
open(unit=1,file='phot/*', status='old')
write(6,*) '....Loading file: phot/*'
do i__=1,36
read(1,'(A368)')ar__2
read(ar__2,3)
+ Mod_1(i__),Teff_1(i__),logg_1(i__),Mass_1(i__),Z_1(i__),
+ C_O_1(i__),logL(i__),xi_1(i__),Ver_1(i__),Band1(i__),BC1(i__),
+ Band2(i__),BC2(i__),Band3(i__),BC3(i__),Band4(i__),BC4(i__),
+ Band5(i__),BC5(i__),Band6(i__),BC6(i__),Band7(i__),BC7(i__),
+ Band8(i__),BC8(i__),Band9(i__),BC9(i__),Band10(i__),BC10(i__),
+ Band11(i__),BC11(i__),Band12(i__),BC12(i__),Band13(i__),
+ BC13(i__),Band14(i__),BC14(i__),Band15(i__),BC15(i__),
+ Band16(i__),BC16(i__),Band17(i__),BC17(i__)
if (BC1(i__) .EQ. 2.2233e-322) BC1(i__) = rNULL__
if (BC2(i__) .EQ. 2.2233e-322) BC2(i__) = rNULL__
if (BC3(i__) .EQ. 2.2233e-322) BC3(i__) = rNULL__
if (BC4(i__) .EQ. 2.2233e-322) BC4(i__) = rNULL__
if (BC5(i__) .EQ. 2.2233e-322) BC5(i__) = rNULL__
if (BC6(i__) .EQ. 2.2233e-322) BC6(i__) = rNULL__
if (BC7(i__) .EQ. 2.2233e-322) BC7(i__) = rNULL__
if (BC8(i__) .EQ. 2.2233e-322) BC8(i__) = rNULL__
if (BC9(i__) .EQ. 2.2233e-322) BC9(i__) = rNULL__
if (BC10(i__) .EQ. 2.2233e-322) BC10(i__) = rNULL__
if (BC11(i__) .EQ. 2.2233e-322) BC11(i__) = rNULL__
if (BC12(i__) .EQ. 2.2233e-322) BC12(i__) = rNULL__
if (BC13(i__) .EQ. 2.2233e-322) BC13(i__) = rNULL__
if (BC14(i__) .EQ. 2.2233e-322) BC14(i__) = rNULL__
if (BC15(i__) .EQ. 2.2233e-322) BC15(i__) = rNULL__
if (BC16(i__) .EQ. 2.2233e-322) BC16(i__) = rNULL__
if (BC17(i__) .EQ. 2.2233e-322) BC17(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Mod_1(i__),Teff_1(i__),logg_1(i__),Mass_1(i__),Z_1(i__),
+ C_O_1(i__),logL(i__),xi_1(i__),Ver_1(i__),Band1(i__),BC1(i__),
+ Band2(i__),BC2(i__),Band3(i__),BC3(i__),Band4(i__),BC4(i__),
+ Band5(i__),BC5(i__),Band6(i__),BC6(i__),Band7(i__),BC7(i__),
+ Band8(i__),BC8(i__),Band9(i__),BC9(i__),Band10(i__),BC10(i__),
+ Band11(i__),BC11(i__),Band12(i__),BC12(i__),Band13(i__),
+ BC13(i__),Band14(i__),BC14(i__),Band15(i__),BC15(i__),
+ Band16(i__),BC16(i__),Band17(i__),BC17(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'spec/*' ! Individual synthetic spectra
C Format for file interpretation
4 format(1X,F10.3,3X,F6.4,1X,E10.5)
C Effective file loading
open(unit=1,file='spec/*', status='old')
write(6,*) '....Loading file: spec/*'
do i__=1,746
read(1,'(A31)')ar__3
read(ar__3,4)lambda(i__),nFlux(i__),nuL_nu(i__)
c ..............Just test output...........
write(6,4)lambda(i__),nFlux(i__),nuL_nu(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end